Search for R-Parity Violating Decays of Supersymmetric Particles in e+e- Collisions at Centre-of-Mass Energies from 189GeV to 202GeV

Searches for the production of supersymmetric particles under the assumption that R-parity is violated via a single dominant LLE, LQD or UDD coupling were performed. These use the data collected by the ALEPH detector at LEP at centre-of-mass energies from 188.6 to 201.6 GeV. The numbers of candidate events observed in the data are consistent with Standard Model expectations. Upper limits on the production cross sections and lower limits on the masses of charginos, sleptons, squarks and sneutrinos are derived.


Introduction
Minimal supersymmetric extensions of the Standard Model (MSSM) [1] usually make the assumption that R-parity, R p = −1 3B+L+2S , is conserved [2], where B denotes the baryon number, L the lepton number and S the spin of a field.The conservation of R-parity is not required theoretically and models in which R-parity is violated can be constructed which are compatible with existing experimental constraints.
The R-parity violating terms of the superpotential considered here are [3] where D, Ū ( Ē) are the down-like and up-like quark (lepton) singlet superfields, and Q (L) is the quark (lepton) doublet superfield respectively; λ, λ ′ and λ ′′ are Yukawa couplings and i, j, k = 1, 2, 3 are generation indices.The presence of such R-parity violating terms imply that the lightest supersymmetric particle (LSP) is no longer stable and that sparticles can be produced singly.The sparticle decays which proceed directly to standard model particles are called direct decays.Decays in which the sparticle first decays, conserving R-parity, to the lightest neutralino are referred to as indirect decays.Both cases are illustrated in Fig. 1.
Other cascade decays are possible but not considered in the following.
In this paper a new search for the resonant production of single sneutrinos decaying indirectly is presented.In addition, previously reported searches for both direct and indirect decays of pair produced sparticles at 183 GeV [4] are extended and applied to new data at higher energies.In particular, new selections for indirect decays of sleptons and squarks via the LQ D operator are developed.Table 1 summarises the possible decays and indicates those addressed in this paper.Other collaborations at LEP have published similar searches at lower energies [5,6,7,8,9].
The following assumptions are made throughout: • All three terms in Equation ( 1) are addressed, however only one term for a specific set of indices (i, j and k) is considered non zero.Unless otherwise stated the derived limits correspond to the choice of indices for the coupling giving the worst limit.
• The lifetime of the sparticles can be neglected, i.e. the mean flight path is less than 1 cm.
• Results are interpreted within the framework of the MSSM.Gaugino mass unification at the electroweak scale is assumed, giving the condition M 1 = 5 3 M 2 tan 2 θ W . • For the case of the charginos and neutralinos, only large values of the universal scalar mass m 0 are considered; this implies the direct decays of the lightest chargino and the next-to-lightest neutralino are suppressed.It also implies three-body decay kinematics for the lightest neutralino.
1 Table 1: For each sparticle the table lists whether the decay mode is searched for (•), possible but not considered (⋄), or not possible (×).Those processes marked with † were not considered in the 183 GeV results [4].
The search results reported here use data collected by the ALEPH detector in 1998 and 1999 from e + e − collisions at centre-of-mass energies between 188.6 GeV and 201.6 GeV.The total data sample used corresponds to an integrated recorded luminosity of 173.6 pb −1 at 188.6 GeV, 29.0 pb −1 at 191.6 GeV, 80.1 pb −1 at 195.5 GeV, 85.9 pb −1 at 199.5 GeV and 41.9 pb −1 at 201.6 GeV.This paper is organised as follows: after a brief description of the ALEPH detector in Section 2, the Monte Carlo samples used for signal and background generation are detailed in Section 3. Sections 4, 5 and 6 give the results and interpretations for each of the R-parity violating couplings, and finally Section 7 gives a summary of the results.

The ALEPH Detector
The ALEPH detector is described in detail in Ref. [10].An account of the performance of the detector and a description of the standard analysis algorithms can be found in Ref. [11].Here, only a brief description of the detector components and the algorithms relevant for this analysis is given.
The trajectories of charged particles are measured with a silicon vertex detector, a cylindrical drift chamber, and a large time projection chamber (TPC).The central detectors are immersed in a 1.5 T axial magnetic field provided by a superconducting solenoidal coil.The electromagnetic calorimeter (ECAL), placed between the TPC and the coil, is a highly segmented sampling calorimeter which is used to identify electrons and photons and to measure their energies.The luminosity monitors extend the calorimetric coverage down to 34 mrad from the beam axis.The hadron calorimeter (HCAL) consists of the iron return yoke of the magnet instrumented with streamer tubes.It provides a measurement of hadronic energy and, together with the external muon chambers, muon identification.The calorimetric and tracking information are combined in an energy flow algorithm which gives a measure of the total energy, and therefore the missing energy, with an uncertainty of (0.6 √ E + 0.6) GeV.
Electron identification is primarily based upon the matching between the measured momentum of the charged track and the energy deposited in the ECAL.Additional information from the shower profile in the ECAL and the measured rate of specific ionisation energy loss in the TPC are also used.Muons are separated from hadrons by their characteristic pattern in HCAL and the presence of associated hits in the muon chambers.

Monte Carlo Samples and Efficiencies
The signal topologies were simulated using the SUSYGEN Monte Carlo program [12] modified as described in Ref. [4].The events were subsequently passed through either a full simulation or a faster simplified simulation of the ALEPH detector.Where the fast simulation was used a subselection of these were also passed through the full simulation to verify the accuracy of the fast simulation.
Samples of all major backgrounds were generated and passed through the full simulation, corresponding to at least 10 times the collected luminosity in the data.The PYTHIA generator [13] was used to produce qq events and four-fermion final states from Weν, ZZ and Zee, with a vector-boson invariant mass cut of 0.2 GeV/c 2 for ZZ and Weν, and 2 GeV/c 2 for Zee.Pairs of W bosons were generated with KORALW [14].The KORALW cross sections were adjusted to agree with the most recent theoretical calculations [15].Pair production of leptons was simulated with UNIBAB [16] (electrons) and KORALZ [17] (muons and taus).The γγ → f f processes were generated with PHOT02 [18].
The selections were optimised to give the minimum expected 95% C.L. excluded cross section in the absence of a signal for masses close to the high end of the expected sensitivity.Selection efficiencies were determined as a function of the SUSY particle masses and the generation structure of the R-parity violating couplings λ ijk , λ ′ ijk and λ ′′ ijk .The cross section limits were evaluated at the highest centre-of-mass energy.Where data taken at a range of centre-of-mass energies contributed to the exclusions the data were weighted with the expected evolution of the cross-section with √ s.
The systematic uncertainties on the selection efficiencies are of order of 4-5% and are dominated by the statistical uncertainty of the Monte Carlo signal samples, with small additional contributions from lepton identification and energy flow reconstruction.They were taken into account by reducing the selection efficiencies by one standard deviation of the statistical error.
When setting the limits, background subtraction was performed for two-and fourfermion final states according to the prescription given in Ref. [19].To take into account the uncertainties on the background estimates, the amount of background subtracted is reduced.For two-fermion processes it was reduced by its statistical error.The contribution from WW and ZZ processes were reduced by the statistical error added in quadrature with 1% of its estimate.The components from Weν and Zee processes were reduced by 20% of their estimate.No background is subtracted for the γγ → f f processes.

4
Decays via a dominant LL Ē coupling Under the assumption of a dominant LL Ē coupling, the decay topologies can consist of as little as two acoplanar leptons in the simplest case (direct slepton decay or single resonant sneutrino production), or they may consist of as many as six leptons plus four neutrinos in the most complicated case (indirect chargino decay).In addition to the purely leptonic topologies, the MSSM cascade decays of charginos into lighter neutralinos may produce multi-jet and multi-lepton final states.No direct decays are possible for the squarks.
The absolute lower limit on the mass of the lightest neutralino of 23 GeV/c 2 obtained in Ref. [20], which is valid for any choice of µ, M 2 , m 0 and generational indices (i, j and k), is used to restrict the range of neutralino mass considered for the indirect decays.
The various selections addressing the above topologies, the expected backgrounds, and the numbers of candidates selected in the data at √ s = 188.6-201.6GeV are summarised in Table 2. Details of the "6 Leptons + E", the "4 Leptons + E" and the "4 Leptons" analyses are given in Ref. [20].The "Acoplanar Leptons" and "Leptons and Hadrons" selections are described in Ref. [4], although the "Leptons and Hadrons" has been updated for the increased centre-of-mass energy as described in section 4.2.Wherever the "Acoplanar Leptons" selection is used the expected combination of final state flavours (e, µ and τ ) for each process is used to set the exclusion.

Single resonant sneutrino production
Single resonant production of sneutrinos [21] can occur for the specific couplings λ 121 and λ 131 .The sneutrino may decay indirectly through the diagram shown in Fig. 2 or directly to e + e − .Since the production cross-section is a function of |λ 1j1 | 2 , limits can be set on the magnitude of λ 1j1 as a function of the sneutrino mass.The best sensitivity is obtained for the case where the sneutrino is produced exactly on shell, √ s = M ν .For centre-of-mass energies above M ν initial state radiation from the e + e − system allows a radiative return to the sneutrino resonance.There is also some sensitivity for M ν > √ s via the production of a virtual sneutrino.
Since the final state consists of two leptons and two neutrinos the "Acoplanar Leptons" selection is used to select these events.Figure 3 shows the excluded values of λ 121 and λ 131 as a function of the mass of the sneutrinos; all data taken in the range √ s = 130 to 189 GeV are used.The data and background numbers for the lower energies are given in Ref. [4,20].Also shown are the results for the direct decays from electroweak fits [22] and the exclusion from low energy measurements [23].

Charginos and neutralinos decaying via LL Ē
Depending on the masses of the gauginos and on the lepton flavour composition in the decay, the indirect decays of charginos to neutralinos and of heavier neutralinos to lighter neutralinos populate different regions in track multiplicity, visible mass and leptonic energy.For this reason three different subselections were developed [4], covering topologies with large leptonic energies and at least two jets (Subselection I), topologies with small multiplicities and large leptonic energy fractions (Subselection II), and topologies with a moderate leptonic energy fraction (Subselection III).The combination of the three subselections is defined as the "Leptons and Hadrons" selection.The complete set of cuts, updated for √ s > 184 GeV is shown in Table 3.
Interpreting the results, shown in Table 2, in the framework of the MSSM, 95% C.L. exclusion limits are derived in the (µ,M 2 ) plane and shown in Fig. 4(a) for large scalar masses m 0 = 500 GeV/c 2 .The corresponding lower limit on the mass of the lightest chargino is essentially at the kinematic limit for pair production.
The searches for the lightest and second lightest neutralino do not extend the excluded region in the (µ,M 2 ) plane beyond that achieved with the chargino search alone.
Table 3: The list of cuts for the "Leptons and Hadrons" selection, which is used for charginos and squarks decaying indirectly via the LL Ē operator.The event variables are defined in Ref. [4].

Squarks decaying via LL Ē
Although squarks cannot decay directly with an LL Ē coupling, they may decay indirectly to the lightest neutralino.This topology is searched for by means of the "Leptons and Hadrons" selection.The 95% C.L. squark mass limits are presented as functions of M χ in Fig. 5 for the case of t1 and b1 squarks.The following limits upon the right-handed squarks can be derived: M ũR > 90 GeV/c 2 and M dR > 89 GeV/c 2 for any λ ijk .

Sleptons decaying via LL Ē
A right-handed slepton can decay directly via the LL Ē coupling to a lepton and antineutrino, hence the acoplanar lepton selection is used.For a given choice of generation indices the decay will produce two final states equally; for the coupling λ ijk these decays are lk R → l i νl j or νl i l j .Excluded cross sections are shown in Fig. 6(a) for the different mixtures of acoplanar lepton states.The MSSM production cross section for right-handed smuon pairs and selectron pairs at µ = −200 GeV/c 2 and tan β = 2 are superimposed.The cross section limit translates into a lower bound on the smuon (or stau) mass of M μR ,τ R > 81 GeV/c 2 and M ẽR > 92 GeV/c 2 (µ = −200 GeV/c 2 , tan β = 2) for the direct decays and the worst case coupling.

Sneutrinos decaying via LL Ē
In pair production each sneutrino can decay directly into pairs of charged leptons giving the final states eeee, eeµµ, eeτ τ , µµµµ, µµτ τ and τ τ τ τ .The different final states correspond to different choices of generation indices.The "Four Lepton" selection was used to derive exclusion limits on the sneutrino pair production cross section shown in Fig. 7(a).These limits translate into a lower bound on the electron sneutrino mass of M νe > 98 GeV/c 2 (µ = −200 GeV/c 2 , tan β = 2) and the muon sneutrino mass of M νµ > 86 GeV/c 2 for direct decays and the worst case coupling.
Indirect decays of sneutrinos are selected using the "Four Leptons + E" selection.The limits in the (M χ , M ν ) plane corresponding to this case are shown in Fig. 7(b) and (c).Using the bound M χ > 23 GeV/c 2 this limit can be interpreted as M νµ,τ > 83 GeV/c 2 and M νe > 94 GeV/c 2 for the worst case coupling, where the cross section for the electron sneutrino is evaluated at µ = −200 GeV/c 2 and tan β = 2.

Decays via a dominant LQ D Coupling
For a dominant LQ D operator the event topologies are mainly characterised by large hadronic activity, possibly with some leptons and/or missing energy.In the simplest case the topology consists of four-jet final states, and in more complicated scenarios of multi-jet and multi-lepton and/or multi-neutrino states.A summary of the results of the various selections is given in Table 4.
The acoplanar jet selection (AJ-H) and the four jets and missing energy selection (4JH) developed for R-parity conserving SUSY searches are used [24].The "MultiJets + Leptons" and the "2J+2τ " selections are updated from Ref. [4] and reoptimised for the higher centreof-mass energy as described below.The "Jets-HM" selection is a retuned version of the "MultiJets + Leptons" subselection I for a high visible mass system.The "4 Jets + 2τ " is unchanged from Ref. [25].Two new selections for five jets and one isolated lepton and four jets and two isolated leptons, "5 Jets + 1 Iso.l" and "4 Jets + 2 Iso.l", are described below.

Charginos and neutralinos decaying via LQ D
Three subselections were developed to select the chargino indirect topologies [4]; some cuts have been reoptimised for the higher centre-of-mass energy.Subselection I is designed to select final states based on hadronic activity, e.g.χ + χ − → qqqqχχ; subselection II is designed for decays such as χ + χ − → lνqqχχ where the leptonic energy is larger, and subselection III is designed to select the decays χ + χ − → lνlνχχ.The combination of the three subselections is defined as the "Multi-jets plus Leptons" selection.The complete set of cuts is shown in Table 5. 5 Jets + 4 4.2 1 0.7 1 1.9 0 1.9 0 0.9 6 9.6 1 Iso.l 4 Jets + 0 3.1 0 0.5 1 1.3 3 1.5 1 0.7 5 7.1 2 Iso.l Table 5: The list of cuts for the "Multi-jets plus Leptons" selection used to select indirect chargino decays via the LQ D operator.The "Jets-HM" selection is also listed; this is used to select intermediate ∆M indirect squark decays with high visible mass.Primed event variables are calculated from physical quantities excluding identified leptons.The event variables are defined in Ref. [4].

Subselection I
Subselection II Subselection III Jets-HM Interpreting these results in the framework of the MSSM, 95% C.L. exclusion limits are derived in the (µ,M 2 ) plane and shown in Fig. 4(b) for large scalar masses (m 0 = 500 GeV/c 2 ).The corresponding lower limit on the mass of the lightest chargino is essentially at the kinematic limit for pair production.
The searches for the lightest and second lightest neutralino do not extend the excluded region in the (µ,M 2 ) plane beyond that achieved with the chargino search alone.

Squarks decaying via LQ D
A squark can decay directly to a quark and either a lepton or a neutrino leading to topologies with acoplanar jets and up to two leptons.Couplings with electrons or muons in the final state are not considered as existing limits from the Tevatron [26] exclude the possibility of seeing such a signal at LEP.To select q → qτ and q → qν, the "2J+2τ " and the "AJ-H" selections are used.The "2J+2τ " selection is unchanged from Ref. [4] except that the sliding mass window is now 10 GeV/c 2 wide centred on the squark mass.Examples of limits for squark production are shown in Fig. 8.In particular, for a dominant λ ′ 33k coupling, which implies Br( tL → qτ ) = 100%, a lower limit of 93 GeV/c 2 is obtained for M tL .
Indirect decays of squarks via the LQ D operator will give six jets and up to two charged leptons.These are selected by the "4 Jets + 2τ " selection if M χ ≤ 20 GeV/c 2 , either the "5 Jets + 1 Iso.l" or "Multi-jets plus Leptons" if (M q − M χ ) ≤ 15 GeV/c 2 , otherwise the Table 6: The list of cuts for the "5 Jets + 1 Isolated Leptons" and "4 Jets + 2 Isolated Leptons" selections, as used for the indirect squark and slepton searches decaying via an LQ D operator.M ′′ vis is the visible mass after the two leading leptons are removed.Other variables are defined in Ref. [4].
5 Jets + 1 Iso.l 4 Jets + 2 Iso.l N ch > 9 "Jets-HM" selection is used.The cuts used for the "5 Jets + 1 Iso.l" selection are listed in Table 6.Limits for the optimistic case of left-handed squarks are shown in Fig. 9.The following limits for tL and bL are derived: M tL > 84 GeV/c 2 and M bL > 74 GeV/c 2 .

Sleptons and sneutrinos decaying via LQ D
Direct decays of sleptons and sneutrinos via the LQ D operator lead to four jet final states.The "Four-Jets" selection from Ref. [4] is applied.The distributions of the di-jet masses for data and Monte Carlo are shown in Fig. 10(a).Fewer events are observed around the W mass peak region than expected.Limits are derived by sliding a mass window of 5 GeV/c 2 across the di-jet mass distribution.The results are shown in Fig. 10(b) and imply M νµ > 77 GeV/c 2 and M μL > 79 GeV/c 2 .
Indirect decays of the sleptons via the LQ D operator will give two, three or four leptons and four jets in the final state; two leptons will be of the same flavour as the initial sleptons.The indirect decays of sneutrinos will give a final state with four jets, up to two leptons and missing energy.For selectrons and smuons the "4 Jets + 2 Iso.l" selection is used except for the special case of λ ′ 3jk = 0 and (M lR − M χ ) < 10 GeV/c 2 where the "4 Jets + 2τ " selection is used.The cuts used in the "4 Jets + 2 Iso.l" are listed in Table 6.Indirect stau decays are selected with the "5 Jets + 1 Iso.l" selection if M χ > 20 GeV/c 2 and either λ ′ 2jk = 0 or λ ′ 1jk = 0, otherwise the inclusive combination of the "5 Jets + 1 Iso.l" and the "Leptons and Hadrons" selections is used.The sneutrinos are selected with the "4JH" selection for M χ > 20 GeV/c 2 and "AJ-H" (acoplanar jets) otherwise.Limits for these decays are shown in Fig. 11 and Fig. 12.The limits are M ẽR > 89 GeV/c 2 , Table 7: The list of cuts for the "Four Jets + E" and "Many Jets + E" selections, as used for the indirect sneutrino searches via a Ū D D operator.

Decays via a dominant Ū D D coupling
For a dominant Ū D D operator the final states are characterised by topologies having many hadronic jets, possibly associated with leptons and missing energy.A number of selections are used: "Four Jets", "Four Jets-Broad", "Many Jets", "Many Jets + Leptons", "Four Jets + 2 Leptons", "Many Jets + 2 Leptons", "Four Jets + E" and "Many Jets + E", all introduced in Ref. [4].The latter two selections, which require missing energy, have been slightly modified for the higher energy and are summarised in Table 7.These selections rely mainly on two characteristics of the events: mass reconstruction of the pair produced sparticles and/or the presence of many jets in the event.Table 8 gives a list of all the selections and the numbers of observed and expected events.

Charginos and neutralinos decaying via Ū D D
The decay of the lightest neutralino and direct decays of the chargino both lead to six hadronic jets in the final state.The indirect chargino decays give rise to a variety of final states depending on the W * decay, they range from ten hadronic jets to six jets associated with leptons and missing energy.The "Many Jets", "Four Jets" and "Many Jets + Lepton" selections are used to cover these topologies.
Interpreting these results in the framework of the MSSM, Fig. 4(c) shows the 95% C.L. exclusion in the (µ, M 2 ) plane.As for the LL Ē and LQ D couplings the lower limit on the The searches for the lightest and second lightest neutralino do not extend the excluded region in the (µ,M 2 ) plane beyond that achieved with the chargino search alone.

Squarks decaying via Ū D D
The direct decay of pair produced squarks leads to four-quark final states.The "Four Jet" selection is therefore used to extract the mass limits.As shown in Fig. 10(b) the mass limits are 82 GeV/c 2 for up-type squarks and 68 GeV/c 2 for down-type squarks.
For the indirect squark decays, which lead to eight-jet topologies, the "Four Jets-Broad" selection is used.The resulting 95% C.L. exclusion in the (M χ , M q) plane for left-handed stop and sbottom are shown in Fig. 13.The corresponding mass limits are M tL > 71.5 GeV/c 2 and M bL > 71.5 GeV/c 2 .

Sleptons decaying via Ū D D
No direct slepton decays are possible via the Ū D D coupling.For the indirect decays of selectron and smuon pairs, which lead to six-jets plus two-lepton final states, the "Four Jets + 2 Leptons" selection is used for large mass differences between the slepton and neutralino, and the "Many Jets + 2 Leptons" for the low mass difference region.In addition, for the very low mass difference region the leptons are very soft and the "Four Jets" selection is used.
Figure 14 shows the 95% C.L. exclusion in the (M χ , M l) plane for selectrons and smuons.The selectron cross section is evaluated at µ = −200 GeV/c 2 and tan β = 2.The shape of the limits at M χ ≈ 20 GeV/c 2 is due to the switch between selections.For M l − M χ > 10 GeV/c 2 this yields M ẽR > 88.5 GeV/c 2 and M μR > 82.5 GeV/c 2 .

Sneutrinos decaying via Ū D D
No direct sneutrino decays are possible via the Ū D D coupling.Sneutrinos decaying indirectly lead to six-jet final states containing two neutrinos.Large mass differences between the sneutrino and neutralino lead to event topologies with significant missing energy, and the "Four Jets + E" selection is used.For small mass differences the "Many Jets + E" selection is used.

Summary
A number of searches were developed to select R-parity violating decay topologies for single and pair production of SUSY particles.It has been assumed that the LSP has a negligible lifetime, and that only one coupling λ ijk , λ ′ ijk or λ ′′ ijk is nonzero.These searches found no evidence for R-parity violating supersymmetry in the data collected at √ s = 188.6-201.6GeV, and various limits were set within the framework of the MSSM.
From searches for singly produced sneutrinos, upper limits on the values of the λ 121 and λ 131 couplings were set as a function of M νµ and M ντ , respectively.
For the indirect decays of sfermions, the limits listed in Table 9  Assuming large m 0 the chargino mass limit is given by the kinematic limit M χ + > 100 GeV/c 2 , irrespective of the R-parity violating operator.
2 and M νµ > 75 GeV/c 2 for the worst case couplings; the selectron and electron sneutrino cross sections are evaluated at µ = −200 GeV/c 2 and tan β = 2.

Table 2 :
The observed numbers of events in the data and the corresponding Standard Model background expectations for the LL Ē selections.

Table 4 :
The observed numbers of events in the data and the Standard Model background expectations for the LQ D selections.

Table 8 :
The observed numbers of events in the data and the corresponding Standard Model background expectations for the Ū D D selections, quoted with any sliding mass cuts removed.

Table 9 :
The 95% confidence level lower mass limits for sparticles decaying indirectly for each of the three R-parity violating couplings.