Search for Supersymmetric Particles with R-Parity Violating Decays in e+e- Collisions at sqrt{s} up to 209GeV

Searches for the pair production of supersymmetric particles under the assumption that R-parity is violated via a single dominant LLEbar, LQDbar or UbarDbarDbar coupling are performed using the data collected by the ALEPH detector at LEP at centre-of-mass energies from 189 to 209GeV. The numbers of observed candidate events in the data are in agreement with the Standard Model expectation, and limits on the production cross sections and on the masses of charginos, sleptons, squarks and sneutrinos are derived.


Introduction
The search for supersymmetric particles addressed in this paper is performed in the framework of the minimal supersymmetric extension of the Standard Model (MSSM) [1] with R-parity violation.Conservation of R-parity [2] is usually assumed in order to prevent experimentally forbidden low energy processes (e.g.fast proton decay).Nevertheless this symmetry is not required theoretically, and models with R-parity violation can be constructed, which are compatible with experimental constraints.
A generic model can be built from the R-parity violating terms of the superpotential [3] where D, Ū and Ē are down-like quark, up-like quark and lepton singlet superfields, Q and L are the quark and lepton doublet superfields; λ, λ ′ and λ ′′ are Yukawa couplings and i, j, k = 1, 2, 3 are generation indices.The presence of such R-parity violating terms implies that the lightest supersymmetric particle (LSP) is unstable and that supersymmetric particles can decay directly to Standard Model particles.
The sparticle decays which proceed directly to Standard Model particles are called direct decays (Fig. 1).Decays in which the sparticle first decays, conserving R-parity, to the lightest neutralino are referred to as indirect decays (Fig. 2).Other cascade decays are possible but not considered in the following.
The following assumptions are made throughout: • All three terms in Equation (1) are addressed, however only one term at a time is considered to be nonzero for a specific set of indices (i, j and k).Unless otherwise stated the derived limits correspond to the choice of indices for the coupling giving the least stringent limit.
• The lifetime of the sparticles can be neglected, i.e. the mean flight path is less than 1 cm.This assumption restricts the sensitivity of the search to R-parity violating couplings greater than 10 −4 for gauginos and 10 −7 for direct decays of sfermions, and constrains the mass of the neutralino to be above 10 GeV/c 2 [4].
• Results are interpreted within the framework of the MSSM with R-parity violation.
In addition to R-parity violating couplings, the parameters are the gaugino mass terms (M i ), the sfermion masses (m f ), the ratio of the Higgs doublet vacuum expectation values (tanβ), the higgsino mass term (µ) and the trilinear couplings (A i ).
Gaugino mass term unification at the GUT scale is assumed, giving the condition M 1 = 5 3 M 2 tan 2 θ W at the electroweak scale.• For chargino and neutralino decays, only large values of sfermion masses are considered, with the consequence that the direct decays of the lightest chargino and the next-to-lightest neutralino are suppressed.It also implies three-body kinematics for the lightest neutralino decay.
The searches presented in this paper cover all types of sparticle pair production; the case of single sneutrino production is addressed in [5].The results reported are based on all data collected by the ALEPH detector in the years 1998, 1999 and 2000.The new data collected in year 2000 at centre-of-mass energies up to 209 GeV is grouped in two samples of 81.6 pb −1 and 133.7 pb −1 luminosity at √ s = 204.9GeV and 206.5 GeV, respectively.Those two samples will be called 205 GeV and 207 GeV in the following.This paper is organized as follows: after a brief description of the ALEPH detector, the Monte Carlo generators used for signal and background are listed in Section 2. The selections used for all topologies and the results obtained when those selections are applied on data and Monte Carlo events are presented in Section 3. Section 4 gives the interpretations, within the MSSM framework, of the absence of any supersymmetric signal in the data.A summary of the results is presented in Section 5.

The ALEPH Detector and Monte Carlo generators
The ALEPH detector is described in detail in Ref. [6].An account of the performance of the detector and a description of the standard analysis algorithms can be found in Ref. [7].
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 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 luminosity monitors extend the calorimetric coverage down to 34 mrad from the beam axis.The calorimetric and tracking information are combined in an energy flow algorithm which gives a measure of the total energy with an uncertainty of (0.6 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.
The signal topologies were simulated using the SUSYGEN Monte Carlo program [8], modified as described in Ref. [9].The events were subsequently passed either through a full simulation, or through a faster simplified simulation of the ALEPH detector for interpolation purposes.
Samples of all major backgrounds were generated and passed through the full simulation.The PYTHIA generator [10] was used to produce qq events and four-fermion final states from Weν, ZZ and Zee.Pairs of W bosons were generated with KORALW [11].Pair production of leptons was simulated with BHWIDE [12] (electrons) and KORALZ [13] (muons and taus).The γγ → f f processes were generated with PHOT02 [14].

Selections and results
The selections were optimized 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 .Two new selections have been added with respect to previous publications, the others are unchanged from those used in Refs.[4,9,15,16,17], except for centre-of-mass energy rescaling.The new selections : "Many jets + Taus" and "Four jets + Taus" address topologies with many jets and taus and are used in the search for indirect stau decays.The set of cuts is shown in Table 1.Tau production is tagged by missing energy and a low multiplicity jet.The corresponding event variable (N min jet−ch ) is the number of charged tracks in the lowest multiplicity jet when forcing the event into four jets.The decay topologies consist either of purely leptonic final states -as few as two acoplanar leptons in the simplest case (direct slepton decay) or as many as six leptons plus four neutrinos in the most complicated case (indirect chargino decay) -or of multi-jet and multi-lepton final states.
The various selections addressing the above topologies are summarized in Table 2 together with the references of the papers in which detailed descriptions of the selection cuts can be found.The numbers of selected data candidates and the expected backgrounds are also given in this table.

Selections for a dominant LQ D Coupling
For a dominant LQ D operator the event topologies are mainly characterized by large hadronic activity, possibly with some leptons and some missing energy.In the simplest case the topology consists of four jets, and in the more complicated scenarios, of several jets or leptons, with or without missing mass.The various selections are listed in Table 3 together with the corresponding numbers of observed data candidates and expected background events.

Selections for a dominant Ū D D coupling
For a dominant Ū D D operator the final states are characterized by topologies having many hadronic jets, possibly associated with leptons or taus and missing energy.
These selections rely mainly on two characteristics of the signal events: the reconstructed mass of pair produced sparticles and the presence of many jets in the final state.In Table 4 the list of all the selections is given, together with the numbers of data candidates and expected background events.

Table 3:
The observed numbers of events in the year 2000 data sample and the corresponding background expectations for the LQ D selections.The selections are given together with the references to the papers in which they are described.

Selection
Ref.

Interpretation within the MSSM framework
For all selections, the number of candidate events observed in the data is in agreement with the Standard Model background expectations.The results of the selections have been used to set limits on the MSSM parameter space.
The cross-section limits were evaluated at 208 GeV.Data taken at lower centre-of-mass energies also contribute to the limits with a reduced weight.The weight was calculated from the expected evolution of the cross section with √ s.
The systematic uncertainties on the selection efficiencies are of the order of 4-5% and are dominated by the statistics of the simulated signal samples, with small additional contributions from lepton identification and energy flow simulation.They were taken into account by reducing the selection efficiencies by the estimated systematic uncertainty.
When setting the limits, background subtraction was performed for two-and fourfermion final states according to the prescription described in Ref. [18].The systematic uncertainty on the expected Standard Model background has been evaluated by detailed comparison of the simulation with the data on control samples obtained with relaxed cuts.The subtracted background has been reduced by the systematic uncertainty derived from these comparisons (typically a few percent depending on the analysis).For the Weν and Zee processes the subtracted background has been further reduced by 20% due to the poor knowledge of the production cross section in the kinematic region selected by this analysis.No background is subtracted for the γγ → f f process.
The absolute lower limit on the mass of the lightest neutralino of 23 GeV/c 2 obtained in Ref. [4] for a dominant LL Ē coupling, which is valid for any choice of µ, M 2 , m 0 (the unified sfermion mass term at the GUT scale) and generational indices (i, j and k), is used to restrict the range of neutralino mass considered for the indirect decays of the LL Ē searches.
The limits on sfermion masses derived from searches and indirect constraints obtained at LEP1 are discussed in Ref. [9].They range from 40 to 45 GeV/c 2 and are indicated on the exclusion plots.

Charginos and neutralinos decaying via LL Ē
The results are interpreted assuming large scalar masses (m 0 = 500 GeV/c 2 ).Depending on the masses of the gauginos and on the lepton flavour composition in the decay, the indirect decays of charginos to neutralinos populate different regions in track multiplicity, visible mass and leptonic energy.The "Leptons and Hadrons" selection, optimized for each possible topology, is used.
In the framework of the MSSM, 95% C.L. exclusion limits are derived in the (µ,M 2 ) plane as shown in Fig. 3a for tan β = 1.41.The corresponding lower limit on the mass of the lightest chargino is 103 GeV/c 2 .
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 LL Ē
Squarks cannot decay directly with an LL Ē coupling but they may decay indirectly via the lightest neutralino.Because the resulting topology is close to that arising from the indirect chargino decay, the "Leptons and Hadrons" selection is used.The 95% C.L. squark mass limits are presented as functions of m χ in Fig. 4 for the case of t1 and b1 squarks.The results are displayed for left-handed squarks and for the values of the mixing angle for which the coupling to the Z vanishes.In the case of purely left-handed squarks, the following limits can be derived: m tL > 91 GeV/c 2 and m bL > 90 GeV/c 2 for any λ ijk .

Sleptons decaying via LL Ē
A slepton can decay directly via the LL Ē coupling to a lepton and anti-neutrino, hence the "Acoplanar Leptons" selection is used.For a given choice of LL Ē coupling, the decay of a right-handed slepton produces two different final states, lk R → ℓ i νℓ j or νℓ i ℓ j , with equal branching ratios.Excluded cross sections are shown in Fig. 5a  Indirect decays of sleptons are selected using the "Six Leptons + E" selection.Limits corresponding to this case are shown in Figs.5b, 5c and 5d.Using the bound of m χ > 23 GeV/c 2 these limits can be interpreted as the mass limits m ẽR > 96 GeV/c 2 (µ = −200 GeV/c 2 , tan β = 2), m μR > 96 GeV/c 2 and m τR > 95 GeV/c 2 .

Sneutrinos decaying via LL Ē
Sneutrinos can decay directly via LL Ē into pairs of charged leptons.For pair produced sneutrinos, the final states, depending on the generation indices, are eeee, eeµµ, eeτ τ , µµµµ, µµτ τ and τ τ τ τ , and can be selected with the "Four Lepton" selection.The exclusion limits on the sneutrino pair production cross section are shown in Fig. 6a.These limits translate into a lower bound on the electron sneutrino mass of m νe > 100 GeV/c 2 (µ = −200 GeV/c 2 , tan β = 2) and the muon sneutrino mass of m νµ > 90 GeV/c 2 .
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 Figs.6b and 6c.Using the bound m χ > 23 GeV/c 2 this limit can be interpreted as m νµ,τ > 89 GeV/c 2 and m νe > 98 GeV/c 2 , where the cross section for the electron sneutrino is evaluated at µ = −200 GeV/c 2 and tan β = 2.

Charginos and neutralinos decaying via LQ D
The results are interpreted assuming large scalar masses (m 0 = 500 GeV/c 2 ).For the various topologies produced in the indirect decays of the chargino pairs via an LQ D coupling, the "MultiJets + Leptons" selection is used.
In the framework of the MSSM, 95% C.L. exclusion limits are derived in the (µ,M 2 ) plane as shown in Fig. 3b.The corresponding lower limit on the mass of the lightest chargino is 103 GeV/c 2 .
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 via LQ D 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 [19] 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.Examples of limits for squark production are shown in Fig. 7.In particular, for a dominant λ ′ 33k coupling, which implies Br( tL → qτ ) = 100%, a lower limit of m tL > 97 GeV/c 2 is obtained.
Limits for left-handed squarks are shown in Fig. 8.The following limits for tL and bL are derived: m tL > 85 GeV/c 2 and m bL > 80 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 is applied.The distributions of the di-jet masses for data and Monte Carlo are shown in Figs.9a and 9b.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. 9c and imply m νµ > 79 GeV/c 2 and m μL > 81 GeV/c 2 .
Indirect decays of the sleptons via the LQ D operator yield two, three or four leptons and four jets in the final state; two leptons are of the same flavour as the initial sleptons.The indirect decays of sneutrinos produce a final state with four jets, up to two leptons and missing energy.For selectrons and smuons the "4 Jets + 2 Iso.ℓ" 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.Indirect stau decays are selected with the "5 Jets + 1 Iso.ℓ" selection if m χ > 20 GeV/c 2 and either λ ′ 2jk = 0 or λ ′ 1jk = 0.The combination of the "5 Jets + 1 Iso.ℓ"

Charginos and neutralinos decaying via Ū D D
The decay of the lightest neutralino leads to six hadronic jets in the final state.The indirect decays of a chargino or the second lightest neutralino give rise to a variety of final states which 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.
In the framework of the MSSM, 95% C.L. exclusion limits in the (µ, M 2 ) plane are obtained as shown in Fig. 3c.The lower limit on the lightest chargino mass is 103 GeV/c 2 .
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. 9c the mass limits are 82.5 GeV/c 2 for up-type squarks and 77 GeV/c 2 for down-type squarks.
For indirect squark decays, which lead to eight-jet topologies, the "Four Jets Broad" and "Four Jets" selections are used.Figure 12 shows the 95% C.L. exclusion limits in the (m χ , m q) plane for left-handed stop and sbottom.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 pair produced selectrons and smuons, which lead to six-jet 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.For indirect stau decays, which lead to six-jet plus two tau final states, the "Four Jets + Taus" and "Many Jets + Taus" selections are used for large and low mass differences between the stau and the neutralino, respectively.

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 plus two neutrinos.For small mass differences between the sneutrino and neutralino, the six jets are well separated and the "Many Jets + E" selection is used.For large mass differences the event is characterized by a significant missing energy, and the "Four Jets + E" selection is used.

Summary
Pair production of supersymmetric particles, followed by direct or indirect decays involving R-parity violating couplings, has been searched for in the data collected with the ALEPH detector at LEP at centre-of-mass energies between 189 and 209 GeV.It has been assumed that the LSP has a negligible lifetime, and that only one λ ijk , λ ′ ijk or λ ′′ ijk coupling is nonzero.Several selections covering all the possible final states have been applied.No evidence for a signal has been found and various limits have been set within the framework of the MSSM with R-parity violating couplings.These results improve on those previously published by ALEPH [15] and by the other LEP collaborations [20].
The limits obtained for direct decays of sfermions are • for an LQ D coupling: • for a Ū D D coupling: -m ũL > 82.5 GeV/c 2 , -m dL > 77.0 GeV/c 2 .
Table 5: The 95% confidence level lower mass limits for indirect sparticle decays for each of the three R-parity violating couplings, assuming m l,ν − m χ > 10 GeV/c 2 for Ū D D and µ = −200 GeV/c 2 and tan β = 2 for ẽ and νe .
Lower mass limit (GeV/c For the indirect decays of sfermions, mass limits are listed in Table 5.For large sfermion masses, an absolute limit of 103 GeV/c 2 has been set on the chargino mass, irrespective of the R-parity violating operator. for the different mixtures of acoplanar lepton states.The MSSM production cross sections for right-handed smuon pairs, and for selectron pairs at µ = −200 GeV/c 2 and tan β = 2, are superimposed.The cross section limits translate into lower bounds of m μR ,τ R > 87 GeV/c 2 and m ẽR > 96 GeV/c 2 (µ = −200 GeV/c 2 , tan β = 2).

Figure 1 :Figure 2 :
Figure 1: Direct R-parity violating decays of supersymmetric particles via the λ, λ ′ and λ ′′ couplings.The points mark the R-parity violating vertex in the decay.

Figure 9 :
Figure 9: The distributions of the reconstructed jet-pair invariant masses after forcing each event into four jets.The points are the data taken in year 2000, for (a) the 205 GeV sample and (b) the 207 GeV sample.The solid histogram is the predicted Standard Model background.In (c), the 95% C.L. cross section upper limit for sleptons (via LQ D), sneutrinos (via LQ D) and squarks (via Ū D D) decaying directly to four jets is shown.The MSSM cross sections for pair production of muon sneutrinos, left-handed smuons and right-handed squarks are superimposed.

Table 1 :
The list of cuts for the "Four Jets + Taus" and "Many Jets + Taus" selections, as used

Table 2 :
The observed numbers of events in the year 2000 data sample and the corresponding background expectations for the LL Ē selections.The selections are given together with the references to the papers in which they are described.

Table 4 :
The observed numbers of events in the year 2000 data sample and the corresponding background expectations for the Ū D D selections.The selections are given together with the references to the papers in which they are described.