Search for R-parity violating decays of a top squark in proton-proton collisions at sqrt(s) = 8 TeV

The results of a search for a supersymmetric partner of the top quark (top squark), pair-produced in proton-proton collisions at sqrt(s) = 8 TeV, are presented. The search, which focuses on R-parity violating, chargino-mediated decays of the top squark, is performed in final states with low missing transverse momentum, two oppositely charged electrons or muons, and at least five jets. The analysis uses a data sample corresponding to an integrated luminosity of 19.7 inverse femtobarns collected with the CMS detector at the LHC in 2012. The data are found to be in agreement with the standard model expectation, and upper limits are placed on the top squark pair production cross section at 95% confidence level. Assuming a 100% branching fraction for the top squark decay chain, top squark to t + chargino, chargino to lepton + jj, top squark masses less than 890 (1000) GeV for the electron (muon) channel are excluded for the first time in models with a single nonzero R-parity violating coupling lambda'[ijk] (i,j,k<= 2), where i,j,k correspond to the three generations.


Introduction
Supersymmetry (SUSY) [1,2] is an extension of the standard model (SM) that may provide a solution to the hierarchy problem [3,4]. In the SUSY framework, quadratically divergent radiative corrections to the Higgs boson mass, dominated by loops involving the top quark, are canceled by loops with a supersymmetric partner of the top quark (top squark). The mass of the top squark is expected to be within a few hundred GeV of the top quark mass, and the supersymmetric Higgs boson partners are also expected to have masses less than 1 TeV [5,6].
Searches for SUSY are performed in many decay channels and are classified into R-parity conserving (RPC) and R-parity violating (RPV) scenarios. The quantum number, R-parity, P R = (−1) 3B+L+2s has a value +1 for SM particles and −1 for superpartners, where B, L, and s are baryon number, lepton number, and spin, respectively [7]. In RPC models the top squark is expected to decay into the lightest SUSY particle, which escapes detection. This results in an event signature with substantial missing transverse momentum. Recent searches performed at the LHC at CERN in events with high missing transverse momentum have reduced the parameter space available for a low mass top squark [8][9][10][11][12][13]. However, R-parity may not be conserved, in which case searches for SUSY particle decaying to SM particles without substantial missing transverse momentum are important.
The superpotential terms that result in R-parity violation are given by: where λ ijk , λ ijk , and λ ijk are three trilinear Yukawa couplings; i, j, k = 1, 2, 3 are generation indices; L and Q are the SU(2) L doublet superfields of the lepton and quark; H u is the Higgs field that gives mass to the up-type quarks; µ i are the bilinear terms that mix lepton and Higgs superfields, and E, D, and U are the SU(2) L singlet superfields of the charged lepton, downtype quark, and up-type quark. The third term violates the conservation of baryon number, while the first two violate the conservation of lepton number. If baryon number and lepton number were both violated, proton decay would proceed at a rate excluded by experimental observations [14,15]. To avoid these experimental constraints and to simplify the interpretation of results, it is commonly assumed that only one of the λ ijk , λ ijk , or λ ijk couplings is different from zero. In this analysis only λ ijk couplings with (i, j, k) ≤ 2 are considered.
In RPV SUSY models with the chargino χ ± 1 lighter than the top squark and nonzero λ ijk , the top squark t can decay via t → b χ ± 1 , with subsequent decay of the chargino to a lepton and two jets via an off-shell sneutrino ( χ ± 1 → ± + jj) [16], as depicted in Fig. 1. The branching fraction of decay χ ± 1 → ν + jj via an off-shell slepton will be negligible unless the slepton and sneutrino masses are comparable. The decay χ ± 1 → W ± χ 0 1 is suppressed for models with χ ± 1 and χ 0 1 almost degenerate in mass. We perform a search for top squark decays, as depicted in Fig. 1, using proton-proton (pp) collisions at a center-of-mass energy of 8 TeV, corresponding to an integrated luminosity of 19.7 fb −1 , collected with the CMS detector at the LHC in 2012. As top squarks are expected to be dominantly pair-produced at the LHC [17], the search is performed using events with exactly two oppositely charged electrons (e ± e ∓ ) or muons (µ ± µ ∓ ), at least five jets of which one or more jet is identified as arising from hadronization of a bottom quark (b-tagged jet), and high S T , where S T is defined as the scalar sum of the transverse momenta of leptons and jets. As a consequence of the assumption that only one of the λ ijk couplings is nonzero, the two leptons must have opposite charge and the same flavor. Details of the event selection are described in Figure 1: Diagram for the R-parity violating, chargino-mediated decay of a top squark. The chargino decays to a lepton and two jets via an off-shell sneutrino with nonzero λ ijk coupling.
The sensitivity of the e ± e ∓ (µ ± µ ∓ ) search does not depend on which of the four RPV couplings associated with the second operator LQD (L i Q j D k ) in Eq. (1): λ 111 , λ 112 , λ 121 , and λ 122 (λ 211 , λ 212 , λ 221 , and λ 222 ) is non-zero, because the final states and kinematic distributions are the same in each case. We expect that the searches have some sensitivity to models with thirdgeneration couplings λ 311 , λ 312 , λ 321 , and λ 322 , via leptonic τ decays; however, we do not include this possible extra contribution in this paper. The difference ∆M t, χ ± 1 between top squark mass M t , and chargino mass M χ ± 1 , is chosen to be 100 GeV, since this value is representative of the bulk of the M t -M χ ± 1 parameter space where the signal reconstruction efficiency is slowly varying. This analysis does not attempt to quantify the decrease in efficiency (and signal sensitivity) in the regions of parameter space where either ∆M t, χ ± 1 or M χ ± 1 is very small (<100 GeV).
Several searches for R-parity violating top squark decays via LQD couplings have been performed by the CMS [18-20] and ATLAS [21] collaborations. These searches have focused on top squark pairs decaying via λ i32 couplings into final states of two leptons (e ± or µ ± ) and two jets or two leptons (e ± or µ ± ) and six jets, four of which are b-tagged jets [20, 21]; via λ 3jk couplings into a final state including two tau leptons and two b-tagged jets [19]; and via the λ 233 coupling into a final state including three leptons and additional jets [18]. The analysis described in this paper is the first search for R-parity violating top squark decays via purely first-or second-generation LQD couplings; in this case, the final states are two leptons (e ± or µ ± ) and six jets, two of which are b-tagged jets.

The CMS detector
A detailed description of the CMS detector, together with a definition of the coordinate system used, can be found elsewhere [22]. A notable feature of the CMS detector is its 6 m internal diameter superconducting solenoid magnet that provides a field of 3.8 T. Within the field volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter, and a brass and scintillator hadron calorimeter. Muon detectors based on gas ionization chambers are embedded in a steel flux-return yoke located outside the solenoid. Events are collected by a two-layer trigger system based on a hardware level-1 trigger, followed by a software-based high-level trigger.
The pseudorapidity range covered by the tracking system is |η| < 2.5, the muon detector extends up to |η| < 2.4, and the calorimeters cover a region with |η| < 3.0. The region of 3 < |η| < 5 is instrumented with steel and quartz fiber forward calorimeters. The hermeticity of the detector up to large values of |η| permits accurate measurement of the momentum balance transverse to the beam direction.

Trigger and event selection
Events are selected using a trigger that requires at least one electron (muon) with a transverse momentum (p T ) threshold of 27 (24) GeV, and |η| < 2.5 (2.1). All objects are reconstructed using a particle-flow (PF) algorithm [23, 24], which uses information from all subsystems to reconstruct photons, electrons, muons, charged hadrons, and neutral hadrons.
To reduce the background from jets containing leptons, we impose isolation constraints on the transverse energy E T,cone from charged-particle tracks or deposits in the calorimeter within a cone ∆R = √ (∆η) 2 + (∆φ) 2 = 0.3 (0.4) around the trajectory of the electron (muon), where φ is the azimuthal angle. The energy from the reconstructed lepton and the average transverse energy density from pileup are subtracted from E T,cone , where pileup is defined as additional inelastic pp collisions within the same or the adjacent LHC bunch crossing. Tracking information together with calorimeter information is used to identify and subtract hadronic energy depositions from charged particles originating from pileup. The contributions to the neutral hadron and photon energy components due to pileup are also computed and subtracted. In the electron channel, the contributions to the neutral hadron and photon energy components due to pileup interactions are subtracted from E T,cone using the jet area technique [25], which computes the transverse energy density of neutral particles from the median of the neutral energy distribution in jets with p T > 3 GeV on an event-by-event basis. In the muon channel, the method assumes the pileup energy density from neutral particles to be half of that from charged hadrons, based on measurements performed in jets [24].
Electrons are reconstructed by matching an energy cluster in the ECAL with a track reconstructed using a Gaussian sum filter [26]. Electrons are required to have p T > 50 GeV and |η| < 2.5. The transition region between the ECAL barrel and endcap is excluded (1.444 < |η| < 1.566) because the calorimeter is not well modeled in this region. Electrons are identified using a multivariate identification algorithm [26], whose input variables are sensitive to bremsstrahlung along the electron path, matching between tracks and ECAL energy deposits, and shower-shape variables. The algorithm is trained with a sample of simulated Drell-Yan (DY) events that contains true electrons and a data sample enriched in misidentified electrons. In addition, the transverse impact parameter of the electron track is required to be less than 2 mm. To reduce backgrounds that arise from photon conversions in the inner pixel detector, at least one pixel hit in the innermost pixel layer is required and the electron must be inconsistent with the hypothesis that it resulted from photon pair creation. We ensure that the electron is isolated from other activity in the event by requiring that E T,cone be less than 10% of the electron p T .
Muon tracks are reconstructed using the information from the muon chambers and the silicon tracker and are required to be consistent with the reconstructed primary vertex. The tracks are required to have at least one hit in both the pixel tracker and muon detector, and at least six hits in the silicon strip tracker. Muons are required to have p T > 50 GeV and |η| < 2.1. Most cosmic ray muons are rejected by requiring that the transverse (longitudinal) impact parameter be less than 2 (5) mm relative to the primary vertex, defined as the vertex with the largest sum of the p 2 T from all tracks associated with it. Only muons with at least ten hits in the silicon strip tracker and at least one hit in the pixel detector are considered, which ensures a precise momentum measurement. Isolation is imposed by the requirement that E T,cone be less than 12% The differences in lepton reconstruction and trigger efficiencies between data and simulation are corrected in simulation in bins of p T and η, using a tag-and-probe method [28].
Jets are reconstructed from PF objects [29] using the anti-k T clustering algorithm [30] with a distance parameter of 0.5. The tracker and ECAL granularity are exploited to precisely measure the charged particles, and hence to determine jet directions at the production vertex. To remove jets arising from instrumental and non-collision backgrounds, additional criteria on charged and neutral hadron energy are applied.
The energy and momentum of each jet are corrected as a function of the jet p T and η to account for the combined response function of the calorimeters. The average energy from pileup is subtracted from the jet [31]. Only jets within |η| < 2.4 are considered. The corrected jet p T must be at least 100 GeV for the leading jet, 50 GeV for the second-leading jet, and 30 GeV for the remaining jets. At least five jets are required in the event.
Events with at least one b-tagged jet are selected. The combined secondary-vertex algorithm [32] uses information from the track impact parameter and vertex information to discriminate between jets that originate from b quarks and jets from light-flavor quarks and gluons. The algorithm correctly identifies jets produced by the hadronization of a b quark (b jets) with an efficiency of approximately 70% and misidentifies jets from light-flavor quarks or gluons (charm quarks) at a rate of approximately 1% (20%) [32]. The b-tagging efficiency in the simulation is scaled to match the measured efficiency in data as a function of p T , η, and the flavor of the jet.
The missing transverse momentum p miss T in the event is defined as the projection of the negative vector sum of the momenta of all reconstructed PF candidates on the plane perpendicular to the beams. The magnitude of p miss T in the event is referred to as E miss T . To suppress leptonic tt decays that often have significant E miss T because of the presence of neutrinos in the final state, E miss T is required to be less than 100 GeV. The dilepton mass M , computed from the two lepton four-momenta, is required to be greater than 130 GeV, based on an optimization to reduce the contribution from low-mass resonances and Z boson decays.
To enhance the statistical significance, for each lepton flavor the sample is divided into three exclusive categories of jet multiplicity: N jets = 5, 6, or ≥ 7. To improve the sensitivity to signal decays, we compute an S T threshold S min T optimized for each top squark mass hypothesis and for each N jets bin. The S min T is determined by maximizing the value of S/ √ S + B, where S and B are the number of expected signal and background events above S min T , respectively.

Simulation of background and signal events
Monte Carlo (MC) simulations of background and signal events are used to optimize the selection criteria for maximum signal sensitivity and to estimate backgrounds. The simulation of the hard-scattering event is performed using the leading-order (LO) matrix element event generator MADGRAPH 5 [33], unless noted otherwise. The CTEQ6L1 [34] set of parton distribution functions (PDF) is used to describe the proton structure. The simulation of the hard-scattering event is then passed to PYTHIA 6.426 [35] with the Z2* tune [36] to model the parton shower, hadronization, and the underlying event. A full simulation of the response of the CMS detector is performed using GEANT4 [37]. Additional simulated minimum bias events are overlaid to reproduce the effects of pileup.
The main SM backgrounds for this search are DY and tt pair production. Additional SM backgrounds, which include diboson (WW, WZ, and ZZ) and single top quark production, are small. The tt sample is generated with up to three additional partons, the DY events are produced with up to four additional partons, and the diboson samples are generated with up to two additional partons. Single top quark production (t-, s-, and tW-channels) is simulated with POWHEG v1.0 [38][39][40][41][42]. Simulated samples of tt and DY are normalized using cross sections computed at next-to-next-to-leading-order (NNLO) [43,44]. Cross sections computed at next-to-leading-order (NLO) [45,46] are used to normalize the single top quark and diboson samples.
The signal samples are generated using MADGRAPH 5, PYTHIA 6.426, and the CTEQ6L1 PDF set. The top squark pair production cross section is computed at NLO as a function of M t , including soft gluon resummation at next-to-leading logarithm (NLL) [47][48][49][50]. The uncertainty in the cross section includes uncertainties associated with the renormalization and factorization scale, and the PDF set [51].

Background estimation
Corrections to the normalization of tt and DY simulations are estimated by examining background enriched samples in data. A summary of the selection criteria for the signal search region and the control regions, including selections on the dilepton mass, is presented in Table 1. Diboson and single top quark production yield small contributions to the background and are estimated from simulation. In simulated tt sample, events are reweighted so that the p T of the top quark matches the data in a dedicated control sample [52]. Table 1: Summary of the selection criteria for the signal region and the control regions. Data in the control regions described as tt, DY normalization, and DY shape are used to estimate SM backgrounds in the signal region.
We use this control sample to compute correction factors for the tt simulation for different jet multiplicities in the signal region. The e ± µ ∓ control sample is well modeled by the simulation, thus correction factors are statistically consistent with unity.
The Drell-Yan production constitutes approximately 8% of the SM background in the signal region, and is reduced by requiring at least one b-tagged jet. The contribution from this source is estimated using a control sample of two oppositely charged same-flavor leptons, which have an invariant mass M in the range 50-130 GeV. We perform a fit to the M distribution to estimate the number of DY events. The DY shape is obtained from background-subtracted data using a DY-enriched sample with no b-tagged jets. The background from diboson decays including leptonic Z boson decays is estimated from simulation and is constrained in the fit. The M shape for the remaining backgrounds does not exhibit a Z boson mass peak, and is described by a linear function. The fit determines the number N DY of DY events and the number of all other background events. To check that the procedure is insensitive to a potential signal contamination, we performed a fit with signal events included, and observed that the obtained N DY is independent of the presence of the potential signal in the control sample. The ratio of N DY from the fit to the simulated number of DY events is calculated for each value of N jets and is used to correct the simulation. This correction factor ranges from 1.2 ± 0.1 to 2.1 ± 0.6 and increases with jet multiplicity.
We checked that the corrections to the DY normalization are valid in the signal region with M > 130 GeV. We compared the numbers of events in different mass ranges using a DYenriched sample with at least five jets and no b-tagged jets. The ratio of the number of events with M in the Z-peak (normalization region) to the number with M in the high-mass tail (signal region) is predicted from simulation to be 11.8 ± 0.4 and observed to be 14.0 ± 3.5 in data, in reasonable agreement.

Systematic uncertainties
We evaluate systematic uncertainties related to each background and to the signal reconstruction efficiency; these are summarized in Table 2. Since the tt correction factor for the simulated sample is estimated from a control sample of e ± µ ∓ events in data, the systematic uncertainty in this background is given by the statistical uncertainty in the control sample. This uncertainty ranges from 10 to 50%, depending on the value of N jets and of S min T . The uncertainties related to lepton trigger, identification, and isolation are negligible. For the small DY background, we take 50 (100)% of correction factor as the systematic uncertainty on the correction in 5 (≥6) jet bin(s). We assign a 30% uncertainty to the diboson and single top quark background contributions to account for the difference between the NLO theoretical calculation and the CMS measurements of the WW and ZZ cross sections [53] and the single top cross sections [54]. The statistical uncertainty due to the finite size of the simulated background samples is 10-30%, depending on the N jets bin and S min T value. The following systematic uncertainties in the signal efficiency are included: jet energy scale (5%) [31], jet b-tagging efficiencies (3%), integrated luminosity (2.6%) [55], lepton identification and reconstruction efficiency (3%), electron energy scale (2%), muon momentum scale (0.9%), and trigger efficiency (1%). Noting the effect of the b-tagging uncertainty on the signal prediction is evaluated by varying the efficiency and misidentification rates by their uncertainties [32,56], and the effect on the signal prediction. The uncertainty related to the lepton isolation requirement for signal events with many jets is estimated using a tt control sample selected as shown in Table 1, but with ≥ 7 jets, and is determined to be 5%. The uncertainty due to the limited size of the simulated signal sample varies from 2 to 7%. The impact of uncertainties related to the PDF set choice, modeling of the top quark p T spectrum, and pileup modeling is determined to be negligible. Figure 2 shows the observed distributions of jet multiplicity, the estimated background distributions, and the expected distributions for signals with a mass M t of either 300 GeV or 900 GeV. In Tables 3 and 4  We use these results to determine 95% confidence level (CL) limits, as a function of M t , on the product of the top-squark pair-production cross section and the square of the branching 7 Results fraction B for the decay t → b ± qq. We use the modified frequentist CL s method [57] with profiling of nuisance parameters. For each M t hypothesis, the Poisson likelihoods of the three N jets bins are combined. Systematic uncertainties are incorporated into the test statistic as nuisance parameters. The nuisance parameter probability density function (pdf) for the tt background normalization, which is estimated from background control regions containing limited numbers of events in high N jets bins, is described by a gamma function. All other uncertainties are treated with log-normal pdfs. With the exception of uncertainties related to the finite size of a control sample, we assume the systematic uncertainties are fully correlated across different N jets bins.

Results
The observed and expected limits on the product of the cross section and the branching fraction squared are shown in Fig. 3. The green (yellow) band corresponds to a variation of one (two) standard deviation(s) on the expected limit. The dotted curve shows the signal cross section, with the width of the associated band showing the sensitivity to uncertainties in the renormalization and factorization scales and the PDF uncertainties [51]. Comparing the observed cross section limits to the signal cross section, we exclude top squarks with masses less than 890 (1000) GeV for the electron (muon) channel. The expected mass exclusion is 950 (970) GeV for the electron (muon) channel.
These cross section limits strictly apply to models with mass difference ∆M t, χ ± 1 = 100 GeV; however, the sensitivities for models with ∆M t, χ ± 1 > 50 GeV are similar. The mass exclusions assume B = 100%. As described earlier, the limits for the electron channel apply equally to models with nonzero λ 111 , λ 112 , λ 121 , or λ 122 and the limits for the muon channel apply equally to models with nonzero λ 211 , λ 212 , λ 221 , or λ 222 . Because the coupling strength does not affect the production cross section and the branching fraction is assumed to be 100%, the value of λ ijk is not important as long as it is sufficiently large to ensure that the sneutrino decays promptly. For coupling values smaller than 10 −5 , the decay lengths are of order 1 mm or greater, resulting in a decreased signal reconstruction efficiency and sensitivity. These are the first limits on chargino-mediated top squark decays via a single LQD coupling λ ijk with (i, j, k ≤ 2).     Figure 3: Observed and expected 95% CL upper limits on the product of the cross section and the branching fraction (B) squared, for e ± e ∓ (left) and µ ± µ ∓ (right). The green (inner) and yellow (outer) bands show the 1 s.d. and 2 s.d. uncertainty ranges in the expected limits, respectively. The dotted curve shows the expected top squark cross section computed at NLO+NLL. The difference M t -M χ ± 1 is assumed to be 100 GeV for the signal model.

Summary
A search for new phenomena using events with two oppositely charged electrons or muons, at least five jets, with at least one b-tagged jet, and low missing transverse momentum has been performed. No excess over the estimated background is observed. The results are interpreted in the framework of chargino-mediated, R-parity violating top squark decays, assuming a 100% branching fraction for the top squark decay chain, t → b χ ± 1 , χ ± 1 → ± + jj. In models with a single nonzero λ ijk coupling with (i, j, k ≤ 2), the results exclude top squarks with mass less than 890 (1000) GeV for the electron (muon) channel at 95% confidence level. These limits are the first obtained for this model. Table 3: Observed events, estimated background, and expected signal yields, for N jets = 5, 6, and ≥7, along with the optimized value of S min T , for different M t in the electron channel. The signal and background uncertainties include both statistical and systematic contributions.    [14] Super-Kamiokande Collaboration, "Search for Proton Decay via p → e + π 0 and p → µ + π 0 in a Large Water Cherenkov Detector", Phys. [20] CMS Collaboration, "Search for pair production of first and second generation leptoquarks in proton-proton collisions at √ s = 8 TeV", (2015). arXiv:1509.03744. In press: Phys. Rev. D.

A Supplementary Material
Figure A.1 shows the number of observed and predicted events in each optimized S min T selection for 5th, 6th, and ≥ 7th jets for electron and muon channels. The event counts correspond to Tables 3 and 4 respectively.
[GeV]  for 5th, 6th, and ≥7th jets for e ± e ∓ (left) and µ ± µ ∓ (right). The S min T selections are optimized for M t ranging from 300 to 1100 GeV. The asymmetric error bars indicate the central confidence intervals for Poisson-distributed data.

B The CMS Collaboration
INFN Sezione di Napoli a , Università di Napoli 'Federico II' b , Napoli, Italy, Università della