Search for R -parity-violating supersymmetry in events with four or more leptons in √ s = 7 TeV pp collisions with the ATLAS detector

: A search for new phenomena in ﬁnal states with four or more leptons (electrons or muons) is presented. The analysis is based on 4.7 fb − 1 of √ s = 7 TeV proton-proton collisions delivered by the Large Hadron Collider and recorded with the ATLAS detector. Observations are consistent with Standard Model expectations in two signal regions: one that requires moderate values of missing transverse momentum and another that requires large eﬀective mass. The results are interpreted in a simpliﬁed model of R -parity-violating supersymmetry in which a 95% CL exclusion region is set for charged wino masses up to 540 GeV. In an R -parity-violating MSUGRA/CMSSM model, values of m 1 / 2 up to 820 GeV are excluded for 10 < tan β < 40.


Introduction
Events with four or more leptons are rarely produced in Standard Model (SM) processes while being predicted by a variety of theories for physics beyond the SM.These include supersymmetry (SUSY) [1][2][3][4][5][6][7][8][9], technicolour [10], and models with extra bosons [11] or heavy neutrinos [12].This paper presents a search with the ATLAS detector for anomalous production of events with four or more leptons."Leptons" refers to electrons or muons, including those from τ decays, but does not include τ leptons that decay hadronically.The analysis is based on 4.7 fb −1 of proton-proton collisions delivered by the Large Hadron Collider (LHC) at a centre-of-mass energy √ s = 7 TeV between March and October 2011.The results are interpreted in the context of R-parity-violating (RPV) SUSY.Similar searches have been conducted at the Tevatron Collider [13,14] and by CMS [15].
2 R-parity-violating supersymmetry SUSY postulates the existence of SUSY particles, or "sparticles", each with spin (S) differing by one-half unit from that of its SM partner.Gauge-invariant and renormalisable interactions introduced in SUSY models can violate the conservation of baryon (B) and lepton (L) number and lead to a proton lifetime shorter than current experimental limits [16].This is usually solved by assuming that R-parity, defined by P R = (−1) 2S+3B+L , is conserved [17][18][19][20][21], which makes the lightest supersymmetric particle (LSP) stable.In P R -conserving models where the LSP is neutral and weakly interacting, sparticle production is characterised by large missing transverse momentum (E miss T ) due to LSPs escaping detection.Many SUSY searches at hadron colliders rely on this large E miss T signature.Alternatively, proton decay can be prevented by imposing other symmetries [22] that require the conservation of either lepton or baryon number, while allowing R-parity violation.Such models can accommodate non-zero neutrino masses and neutrino-mixing angles consistent with the observation of neutrino oscillations [23].If R-parity is violated, the LSP decays into SM particles and the signature of large E miss T may be lost.Trilinear lepton-number-violating RPV interactions can generate both charged leptons and neutrinos during the LSP decay, and therefore lead to a characteristic signature with high lepton multiplicity and moderate values of E miss T compared to R-parity-conserving models.

R-parity-violating models
The results of this analysis are interpreted in two Minimal Supersymmetric Standard Model (MSSM) scenarios with an RPV superpotential term given by W RP V = λ ijk L i L j Ēk .The i, j, k indices of the λ ijk Yukawa couplings refer to the lepton generations.The lepton SU (2) doublet superfields are denoted by L i , while the corresponding singlet superfields are given by E k .Single coupling dominance is assumed with λ 121 as the only non-zero coupling.The λ 121 coupling is chosen as a representative model with multiple electrons and muons in the final state.Comparable signal yields are expected with a choice of λ 122 single coupling dominance.
The first scenario is a simplified model [24] where the lightest chargino and neutralino are the only sparticles with masses below the TeV scale.Charginos ( χ± i , i = 1, 2) and neutralinos ( χ0 j , j = 1, 2, 3, 4) are the mass eigenstates formed from the linear superposition of the SUSY partners of the Higgs and electroweak gauge bosons.These are the Higgsinos, winos and bino.Wino-like charginos are pair-produced and each decays into a W boson and a bino-like χ0 1 , which is the LSP.The LSP then undergoes the three-body decay χ0 1 → eµν e or χ0 1 → eeν µ through a virtual slepton or sneutrino, with a branching fraction of 50% each, as shown in figure 1(a-c).The width of the χ0 1 is fixed at a value of 100 MeV resulting in prompt decays.
The second scenario is taken from ref. [25] and is the (m 1/2 , tan β) slice of the minimal SUper GRAvity/Constrained MSSM (MSUGRA/CMSSM) containing the BC1 point of ref. [26].The unification parameters m 0 and A 0 are zero and µ is positive1 .The RPV coupling λ 121 is set to 0.032 at the unification scale [27].Both strong and weak processes contribute to SUSY pair production, where weak processes are dominant above m 1/2 ∼ 600 GeV.The lighter stau, τ1 , is the LSP over most of the parameter-space, and decays with equal probability through the four-body decays τ1 → τ eµν e or τ1 → τ eeν µ , via a virtual neutralino and sneutrino or slepton, as shown in figure 1(d).While the Higgs mass values in this model are lower than those of the recently observed Higgs-like resonance [28,29], the MSUGRA/CMSSM scenario considered nevertheless remains an instructive benchmark model.

Detector description
ATLAS [30] is a multipurpose particle detector with forward-backward symmetric cylindrical geometry.It includes an inner tracker (ID) immersed in a 2 T axial magnetic field providing precision tracking of charged particles for pseudorapidities 2 |η| < 2.5.Sampling calorimeter systems with either liquid argon or scintillator tiles as the active media provide energy measurements over the range |η| < 4.9.The muon detectors are positioned outside the calorimeters and are contained in an air-core toroidal magnetic field produced by superconducting magnets with field integrals varying from 1 T•m to 8 T•m.They provide trigger and high-precision tracking capabilities for |η| < 2.4 and |η| < 2.7, respectively.
of the two Higgs vacuum expectation values in the MSSM, A0 is the trilinear coupling and µ is the Higgs mixing parameter, all defined at the unification scale.
2 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe.The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upward.Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the beam pipe.The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2).

Monte Carlo simulation
Several Monte Carlo (MC) generators are used to simulate SM processes and new physics signals relevant for this analysis.SHERPA [31] is used to simulate the diboson processes W W , W Z and ZZ, where Z also includes virtual photons.These diboson samples correspond to all SM diagrams leading to the ℓνℓ ′ ν ′ , ℓℓℓ ′ ν ′ , and ℓℓℓ ′ ℓ ′ final states, where ℓ, ℓ ′ = e, µ, τ and ν, ν ′ = ν e , ν µ , ν τ .Interference between the diagrams is taken into account.MadGraph [32] is used for the t tW , t tW W , t tZ, W γ and Zγ processes.MC@NLO [33] is chosen for the simulation of single and pair production of top quarks, and ALPGEN [34] is used to simulate W +jets and Z+jets processes.Expected diboson yields are normalised using next-to-leading-order (NLO) QCD predictions obtained with MCFM [35,36].The topquark pair-production contribution is normalised to approximate next-to-next-to-leadingorder calculations (NNLO) [37] and the t tW , t tW W , t tZ contributions are normalised to NLO predictions [38,39].The W γ and Zγ yields are normalised to be consistent with the ATLAS cross-section measurements [40].The QCD NNLO FEWZ [41,42] cross-sections are used for normalisation of the inclusive W +jets and Z+jets processes.The choice of the parton distribution functions (PDFs) depends on the generator.The CTEQ6L1 [43] PDFs are used with MadGraph and ALPGEN, and the CT10 [44] PDFs with MC@NLO and SHERPA.
The simplified model samples are produced with Herwig++ [45] and the MSUGRA/ CMSSM BC1-like samples are produced with HERWIG [46].The yields of the SUSY samples are normalised to the NLO cross-sections obtained from PROSPINO [47] for weak processes, and to next-to-leading-logarithmic accuracy (NLL) for strong processes.
Fragmentation and hadronisation for the ALPGEN and MC@NLO samples are performed with HERWIG, while for MadGraph, PYTHIA [48] is used, and for SHERPA these are performed internally.JIMMY [49] is interfaced to HERWIG for simulation of the underlying event.For all MC samples, the propagation of particles through the ATLAS detector is modelled using GEANT4 [50,51].The effect of multiple proton-proton collisions from the same or different bunch crossings is incorporated into the simulation by overlaying additional minimum-bias events generated by PYTHIA onto hard-scatter events.Simulated events are weighted to match the distribution of the number of interactions per bunch crossing observed in data.Simulated data are reconstructed in the same manner as the data.

Event reconstruction and preselection
The data sample was collected with an inclusive selection of single-lepton and double-lepton triggers.For single-lepton triggers, at least one reconstructed muon (electron) is required to have transverse momentum p µ T (transverse energy E e T ) above 20 GeV (25 GeV).For dilepton triggers, at least two reconstructed leptons are required to have triggered the event, with transverse energy or momentum above threshold.The two muons are each required to have p µ T > 12 GeV for dimuon triggers, and the two electrons to have E e T > 17 GeV for dielectron triggers, while the thresholds for electron-muon triggers are E e T > 15 GeV and p µ T > 10 GeV.These thresholds are chosen such that the overall trigger efficiency is high, typically in excess of 90%, and independent of the transverse momentum of the triggerable objects within uncertainties.
Events recorded during normal running conditions are analysed if the primary vertex has five or more tracks associated to it.The primary vertex of an event is identified as the vertex with the highest Σp 2 T of associated tracks.Electrons must satisfy "medium" identification criteria [52] and fulfil |η| < 2.47 and E T > 10 GeV, where E T and |η| are determined from the calibrated clustered energy deposits in the electromagnetic calorimeter and the matched ID track, respectively.Muons are reconstructed by combining tracks in the ID and tracks in the muon spectrometer [53].Reconstructed muons are considered as candidates if they have transverse momentum p T > 10 GeV and |η| < 2.4.
Jets are reconstructed with the anti-k t algorithm [54] with a radius parameter of R = 0.4 using clustered energy deposits calibrated at the electromagnetic scale 3 .The jet energy is corrected to account for the non-compensating nature of the calorimeter using correction factors parameterised as a function of the jet E T and η [55].The correction factors were obtained from simulation and have been refined and validated using data.Jets considered in this analysis have E T > 20 GeV and |η| < 2.5.The p T -weighted fraction of the tracks in the jet that are associated with the primary vertex is required to be larger than 0.75.
Events containing jets failing the quality criteria described in ref. [55] are rejected to suppress both SM and beam-induced background.Jets are identified as containing b-hadron decays, and thus called "b-tagged", using a multivariate technique based on quantities such as the impact parameters of the tracks associated to a reconstructed secondary vertex.The chosen working point of the b-tagging algorithm [56] correctly identifies b-quark jets in simulated top-quark decays with an efficiency of 60% and misidentifies the jets initiated by light-flavour quarks or gluons with a rate of < 1%, for jets with E T > 20 GeV and |η| < 2.5.
The missing transverse momentum, E miss T , is the magnitude of the vector sum of the transverse momentum or transverse energy of all p T > 10 GeV muons, E T > 20 GeV electrons, E T > 20 GeV jets, and calibrated calorimeter energy clusters with |η| < 4.9 not associated to these objects [57].
In this analysis, "tagged" leptons are leptons separated from each other and from candidate jets as described below.If two candidate electrons are reconstructed with ∆R ≡ (∆φ) 2 + (∆η) 2 < 0.1, the lower energy one is discarded.Candidate jets within ∆R = 0.2 of an electron candidate are rejected.To suppress leptons originating from semi-leptonic decays of c-and b-quarks, all lepton candidates within ∆R = 0.4 of any remaining jet candidates are removed.Muons undergoing bremsstrahlung can be reconstructed with an overlapping electron candidate.To reject these, tagged electrons and muons separated from jets and reconstructed within ∆R = 0.1 of each other are both discarded.Events containing one or more tagged muons that have transverse impact parameter with respect to the primary vertex |d 0 | > 0.2 mm or longitudinal impact parameter with respect to the primary vertex |z 0 | > 1 mm are rejected to suppress cosmic muon background.
"Signal" leptons are tagged leptons for which the scalar sum of the transverse momenta of tracks within a cone of ∆R = 0.2 around the lepton candidate, and excluding the lepton candidate track itself, is less than 10% of the lepton E T for electrons and less than 1.8 GeV for muons.Tracks selected for the electron and muon isolation requirement defined above are those which have p T > 1 GeV and are associated to the primary vertex of the event.Signal electrons must also pass "tight" identification criteria [52].

Signal region selection
Selected events must contain four or more signal leptons.The invariant mass of any sameflavour opposite-sign (SFOS) lepton pair, m SFOS , must be above 20 GeV, otherwise the lepton pair is discarded to suppress background from low-mass resonances.Events which contain a SFOS lepton pair inside the [81.2, 101.2]GeV interval are vetoed to reject Z-boson candidates (Z-veto).
Two signal regions are then defined: a signal region with E miss T > 50 GeV (SR1) and one with effective mass m eff > 300 GeV (SR2).The effective mass is defined by the scalar sum shown in (eq.1), where p µ T (E e T ) is the transverse momentum of the signal muons (electrons) and E j T is the transverse energy of jets with E T > 40 GeV: SR1 is optimised for regions of the parameter space with values of m 1/2 below ∼700 GeV or m χ± 1 below ∼300 GeV.SR2 targets the production of heavier sparticles.In general, SR1 is sensitive to models with E miss T originating from neutrinos and SR2 to scenarios with a large multiplicity of high-p T objects.
Three further regions are defined to validate the expected background against data.The validation regions are described in section 10.Table 1 summarises the signal and validation region definitions.

Standard Model background estimation
Several SM processes, which are classified into irreducible and reducible components below, contribute to the background in the signal regions.The dominant sources are ZZ, W Z, and t t production in both SR1 and SR2.

Irreducible background processes
A background process is considered "irreducible" if it leads to events with four real, isolated leptons, referred to as "real" leptons below.These include ZZ, t tZ and t tW W production, where a gauge boson may be produced off-mass-shell.These contributions are determined using the corresponding MC samples, for which lepton and jet selection efficiencies [58][59][60][61], and trigger efficiencies are corrected to account for differences with respect to data.

Reducible background processes
A "reducible" process has at least one "fake" lepton, that is either a lepton from a semileptonic decay of a b-or c-quark, referred to as heavy-flavour, or an electron from an isolated single-track photon conversion.The contribution from misidentified light-flavour quark or gluon jets is found to be negligible based on studies in simulation.The reducible background includes W Z, t t, t tW , W W , single t-quark, or single Z-boson production, in all cases produced in association with jets or photons.The yield of W bosons with three fake leptons is negligible.The t t and the W Z backgrounds correspond respectively to 59% and 35% of the reducible background in SR1, and to 46% each in SR2.In both SR1 and SR2, fake leptons are predominantly fake electrons (99%), originating either from b-quarks in t t candidate events or from conversions in W Z candidate events.Fake muons from b-quark decays in t t events are suppressed by the object separation scheme described in section 6. Regardless of the origin, the misidentification probability decreases as the lepton p T increases.
The reducible background is estimated using a weighting method applied to events containing signal leptons (ℓ S ) and loose leptons (ℓ L ), which are tagged leptons failing the signal lepton requirements.Since the reducible background is dominated by events with at most two fake leptons, it is estimated as: where the second term corrects for the double counting of reducible-background events with two fake leptons in the first term.The term N data (3ℓ S + ℓ L ) is the total number of events with three signal and one loose lepton, while N MCirr (3ℓ S + ℓ L ) is the irreducible contribution of events obtained from simulation.The definitions of N data (2ℓ and N MCirr (2ℓ S + ℓ L 1 + ℓ L 2 ) are analogous.As a conservative approach, the potential signal contamination in the 3ℓ S + ℓ L and 2ℓ S + ℓ L 1 + ℓ L 2 loose lepton data samples is not taken into account.The average "fake ratio" F depends on the flavour and kinematics of the loose lepton ℓ L and it is defined as: where i is the type of fake (heavy-flavour leptons or conversion electrons) and j is the process category the fake originates from (top quark or W/Z boson).The fake ratios f ij are defined as the ratios of the probabilities that fake tagged leptons are identified as signal leptons to the probabilities that they are identified as loose leptons.The f ij are determined for each relevant fake type and for each reducible-background process, and they are parameterised in muon (electron) p T (E T ) and η.The fake ratios are weighted according to the fractional contribution of the process they originate from through R ij fractions.Both f ij and R ij are determined in simulation.Each correction factor α i is the fake ratio measured in data divided by that in simulation, in control samples described below.The fake ratios F are estimated to vary from 0.8 to 0.05 for muons when the p T increases from 10 to 100 GeV.For electrons, there is little p T dependence, and the F values are ∼0.T .The remaining non-b b background (∼1% level) is subtracted from the control sample in data using MC predictions.In this control sample, the fake ratio from heavy-flavour decays is calculated using the ratio of tagged leptons passing signal lepton requirements to those that fail.The correction factors are found to be 0.97 ± 0.13 and 0.86 ± 0.13 for electrons and muons respectively.
The correction factor for electron candidates originating from photon conversions is determined in a sample of photons radiated from a muon in Z → µµ decays.Events with two opposite-sign muons and one tagged electron are selected and the invariant mass of the µµe triplet is required to lie within 10 GeV of the nominal Z-boson mass value.In this control sample, the fake ratio for conversion electrons is calculated using the ratio of tagged electrons identified as signal electrons to those identified as loose electrons.The correction factor is found to be 1.24 ± 0.13.

Systematic uncertainties
Several sources of systematic uncertainty are considered in the signal, control and validation regions.Correlations of systematic uncertainties between processes and regions are accounted for.
MC-based sources of systematic uncertainty affect the irreducible background, the R ij fractions of the average fake ratios, and the signal yields.The MC-based systematic sources include the acceptance uncertainty due to the PDFs and the theoretical cross-section uncertainties due to the renormalisation/factorisation scale and PDFs.The uncertainty due to the PDF set is determined using the error set of the original PDF and the scale uncertainties are calculated by varying the factorisation and renormalisation scales.Additional systematic uncertainties are those resulting from the jet energy scale [58] and resolution [59], the lepton efficiencies [60,61], energy scales and resolutions, and uncertainties in b-tagging rates [62][63][64][65].The choice of MC generator is also included for the irreducible background.The systematic uncertainty on the luminosity (3.9%) [66,67] affects only the yields of the irreducible background and the signal yield.
In SR1, the total uncertainty on the irreducible background is 70%.This is dominated by the uncertainty on the efficiency times acceptance of the signal region selection for the ZZ MC event generator, determined by comparing the SHERPA and POWHEG [68][69][70][71] generators and found to be 65% of the SHERPA ZZ yield.The next largest uncertainties on the irreducible background are due to the jet energy scale (53%) and that due to the limited number of MC events generated (34%).All the remaining uncertainties in this signal region lie in the range of 0.1-5%.In SR2, the uncertainties are similar, except for smaller uncertainties due to the jet energy scale (3%) and the limited number of generated events (17%).
For the average fake ratio F , the R ij fractions are varied between 0% and 100% to account for the uncertainty on the R ij fractions from the sources listed above.The fake ratios f ij are generally found to be similar across processes and types of fakes in most of the kinematic regions.Also included in the uncertainty on the reducible background is the uncertainty from the dependence of the fake ratio on E miss T (10%) and the uncertainty on the fake-ratio correction factors (15-34%).In SR1, the dominant uncertainty on the reducible background is that due to the R ij fractions (75%).This relatively large uncertainty is due to the limited number of events in the loose lepton data sample, and the fact that those leptons fall in kinematic regions in which the fake ratio has larger variations between processes and types of fakes.The next to leading uncertainties are those from the statistical uncertainty on the data (60%) and the dependence of the fake ratio on E miss T (25%).In SR2, the uncertainty on the reducible background is dominated by that from the limited number of data events (140%), followed by that due to the R ij fractions (6%) and the dependence of the fake ratio on E miss T (4%).
The total uncertainties on the signal yields are 10-20% and are calculated using the method described in ref. [72].Signal cross-sections are calculated to NLO in the strong coupling constant using PROSPINO4 .An envelope of cross-section predictions is defined using the 68% CL intervals of the CTEQ6.6 [77] (including the α S uncertainty) and MSTW [78] PDF sets, together with variations of the factorisation and renormalisation scales by factors of two or one half.The nominal cross-section value is taken to be the midpoint of the envelope and the uncertainty assigned is half the full width of the envelope, closely following the PDF4LHC recommendations [79].
The background predictions have been verified in three validation regions (VR).A region (VR1) selects events with three signal leptons, E miss T > 50 GeV, and vetos events with Zboson candidates described in section 7.In VR1 the reducible background is dominated by t t production.The contribution from ZZ production is validated in a region (VR2) defined by events with four leptons containing a Z-boson candidate.A region (VR3) containing events with E miss T < 50 GeV, m eff < 300 GeV, four signal leptons, and no Z-boson candidate is used to validate the sum of the reducible and irreducible backgrounds.The data and predictions are in agreement within the quoted statistical and systematic uncertainties, as shown in table 2. Reported are the probabilities that the background fluctuates to the observed number of events or higher (p 0 -value) and the corresponding number of standard deviations (σ).No significant excess of events is found in the signal regions.Upper limits on the visible cross-section of new physics processes are calculated, defined by the product of production cross-section, acceptance and efficiency, and placed at 95% CL with the modified frequentist CL s prescription [80].All systematic uncertainties and their correlations are taken into account via nuisance parameters in a profile likelihood fit [81].Observed 95% CL limits on the visible cross-section of new physics processes are placed at 1.3 fb in SR1 and 1.1 fb in SR2.The corresponding expected limits are 0.8 fb and 0.7 fb, respectively.
The results of the analysis are interpreted in two RPV SUSY scenarios and shown in fig.3, where the limits are calculated choosing the signal region with the best expected limit for each of the model points.The uncertainties on the signal cross-section are not included in the limit calculation but their impact on the observed limit is shown by the red dotted lines.
The main features of the exclusion limits can be explained in broad terms as follows.In the simplified model the sensitivity is governed by the production cross-section, which decreases at large m χ± 1 values, and by the efficiency of the signal region selection cuts.The requirements on the minimum value of the SFOS dilepton mass, and on the minimum lepton-lepton separation reduce the acceptance and selection efficiency for leptons from light χ0 1 , hence only values of 10 GeV and above are considered.As the mass of the χ0 1 approaches zero, the phase space for leptonic decays is greatly reduced, while values of m χ0 1 = 0 GeV are not possible without the LSP becoming stable.In the MSUGRA/CMSSM model, the production cross-section decreases as m 1/2 increases, leading to smaller sensi-tivities in the high-mass region.Due to lower values of the four-body decay branching ratio and an increased τ1 lifetime, the sensitivity drops for tan β values above 40.In the region of the parameter space with m 1/2 ∼ 800 GeV, weak gaugino production contributes 50% to the total SUSY production cross-section, while ẽ/μ/ν (τ 1 ) production contributes 10-20% (10-30%), and strong production 20%.Similar fractional contributions are also seen in the two signal regions.

Summary
Results from a search for new phenomena in the final state with four or more leptons (electrons or muons) and either moderate values of missing transverse momentum or large effective mass are reported.The analysis is based on 4.7 fb −1 of proton-proton collision data delivered by the LHC at √ s = 7 TeV.No significant excess of events is found in data.
Observed 95% CL limits on the visible cross-section are placed at 1.3 fb and 1.1 fb in the two signal regions, respectively.The null result is interpreted in a simplified model of chargino pair-production in which each chargino cascades to the lightest neutralino that decays into two charged leptons (ee or eµ) and a neutrino via an RPV coupling.In the simplified model of RPV supersymmetry, chargino masses up to 540 GeV are excluded for LSP masses above 300 GeV.Limits are also set in an RPV MSUGRA model with a τ1 LSP that promptly decays into a τ lepton, two charged leptons and a neutrino, where values of m 1/2 below 820 GeV are excluded when 10 < tan β < 40.
[GeV]    MSUGRA/CMSSM.The expected and observed limits are calculated without signal cross-section uncertainty taken into account.The yellow band is the ±1σ experimental uncertainty on the expected limit (black dashed line).The red dotted lines are the ±1σ signal theory uncertainty on the observed limit (red solid line).Linear interpolation is used to account for the discreteness of the signal grids.The exclusion contours are optimised by applying in each signal grid point the CL values from the more sensitive signal region (lowest expected CL).It should be noted that the y-axis does not start at zero for either plot.The LEP limit corresponds to the limit on the χ± 1 mass in ref. [82], while the D0 limit corresponds to the limit in the χ± 1 , χ0 1 mass plane in ref. [13].The hard m τ1 < 40 GeV limit is from the LEP Z width measurement [83].

Figure 1 .
Figure 1.Illustration of the χ0 1 decay in the chosen simplified model (a-c), and the four-body τ1 decay in the MSUGRA/CMSSM model (d).The charge conjugate decay of the χ0 1 decay is implied.

3 .
The correction factor for heavy-flavour fakes is measured in a b b-dominated control sample.This is defined by selecting events with only one b-tagged jet (containing a muon) and a tagged lepton.The non-b b contributions from the single and pair production of top quarks and W bosons produced in association with b-quarks are suppressed with E miss T < 40 GeV and transverse mass m T = 2 • E miss T • p ℓ T • (1 − cos ∆φ ℓ,E miss T ) < 40 GeV requirements, where ∆φ ℓ,E miss T is the azimuthal angle between the tagged lepton ℓ and the E miss

Figure 2 . 1 =
Figure 2. Distributions of (a) E miss T and (b) m eff for events with at least four leptons and no Z-boson candidates (before either the E miss T or m eff requirements).The uncertainty band includes both statistical and systematic uncertainty.The yields of two benchmark SUSY models are shown for illustration purposes."SUSY ref. point 1" is a simplified model point defined by m χ± 1 = 500 GeV and m χ0 1 = 300 GeV, while "SUSY ref. point 2" is a MSUGRA model point defined by m 1/2 = 860 GeV and tan β = 37.The signal distributions are not stacked on top of the expected background.

Figure 3 .
Figure 3. Observed and expected 95% CL limit contours for (a) simplified model and (b) MSUGRA/CMSSM.The expected and observed limits are calculated without signal cross-section uncertainty taken into account.The yellow band is the ±1σ experimental uncertainty on the expected limit (black dashed line).The red dotted lines are the ±1σ signal theory uncertainty on the observed limit (red solid line).Linear interpolation is used to account for the discreteness of the signal grids.The exclusion contours are optimised by applying in each signal grid point the CL values from the more sensitive signal region (lowest expected CL).It should be noted that the y-axis does not start at zero for either plot.The LEP limit corresponds to the limit on the χ± 1 mass in ref.[82], while the D0 limit corresponds to the limit in the χ± 1 , χ0 1 mass plane in ref.[13].The hard m τ1 < 40 GeV limit is from the LEP Z width measurement[83].

Table 1 .
The selection requirements for the two signal (SR) and three validation regions (VR).The Z-candidate veto (requirement) rejects (selects) events that have a same-flavour opposite-sign lepton pair with mass inside the [81.2, 101.2]GeV interval.

Table 2 .
Expected number of events from SUSY signals, SM backgrounds, and observed number of events in data in validation regions VR1, VR2 and VR3 (4.7 fb −1 )."SUSY ref. point 1" refers to a simplified model with m χ± 1 = 500 GeV, m χ0 1 = 300 GeV, while "SUSY ref. point 2" refers to a MSUGRA model with m 1/2 = 860 GeV and tan β = 37.The uncertainties quoted include statistical and systematic effects.The processes labelled with ( †) are classified as irreducible background in the three-lepton validation regions.In the four-lepton validation regions they are included in the reducible background.The subtraction of the double-counted contributions can lead to negative central values for the reducible background.

Table 3 .
Expected number of events from SUSY signals, SM backgrounds, and observed number of events in data in signal regions SR1 and SR2 (4.7 fb −1 )."SUSY ref. point 1" refers to a simplified model with m χ± 1 = 500 GeV, m χ0 1 = 300 GeV, while "SUSY ref. point 2" refers to a MSUGRA model with m 1/2 = 860 GeV and tan β = 37.The uncertainties quoted include statistical and systematic effects.Upper limits on the observed (expected) visible production cross-section of new physics processes at 95% CL are also shown.