Search for SUSY with R-parity violating LL\bar{E} couplings at sqrt(s) = 189 GeV

Searches for pair production of supersymmetric particles under the assumption that R-parity is not conserved are presented, based on data recorded by the DELPHI detector in 1998 from e+e- collisions at a centre-of-mass energy of 189 GeV. Only one R-parity violating LL\bar{E} term (i.e. one lambda coupling), which couples scalar leptons to leptons, is considered to be dominant at a time. Moreover, it is assumed that the strength of the R-parity violating couplings is such that the lifetimes can be neglected. The search for pair production of neutralinos, charginos and sleptons has been performed for both direct R-parity violating decays and indirect cascade decays. The results are in agreement with Standard Model expectations, and are used to update the constraints on the MSSM parameter values and the mass limits previously derived at sqrt(s) = 183 GeV. The present 95% C.L. limits on supersymmetric particle masses are: m_{\tilde{\chi}^0}>30 GeV/c^2 and m_{\tilde{\chi}^\pm}>94 GeV/c^2; m_{\tilde{\nu}}>76.5 GeV/c^2 (direct and indirect decays); m_{\tilde{\ell}_R}>83 GeV/c^2 (indirect decay only).


Motivations
The R-parity symmetry plays an essential role in the construction of supersymmetric theories of interactions, such as the Minimal Supersymmetric extension of the Standard Model (MSSM) [1]. The conservation of R-parity is closely related to the conservation of lepton (L) and baryon (B) numbers and the multiplicative quantum number associated to the R-parity symmetry is defined by R p = (−1) 3B+L+2S for a particle with spin S [2]. Standard particles have even R-parity, and the corresponding superpartners have odd R-parity. The conservation of R-parity guarantees that the new spin-0 sfermions cannot be directly exchanged between standard fermions. It implies that the new sparticles (R p = −1) can only be pair-produced, and that the decay of a sparticle should lead to another one, or an odd number of them. Then, it ensures the stability of the Lightest Supersymmetric Particle (LSP). The MSSM is designed to conserve R-parity: it is phenomenologically justified by proton decay constraints, and by the hope that a neutral LSP will provide a good dark matter candidate.
One of the major consequences of the R-parity violation is obviously that the LSP is no longer stable since it is allowed to decay to standard fermions. This fact modifies the signatures of the supersymmetric particle production compared to the expected signatures in case of R-parity conservation. In any case, whether it turns out to be absolutely conserved or not, R-parity plays an essential role in the study of the phenomenological implications of supersymmetric theories.
In complementarity with the searches for supersymmetric particles in the hypothesis of R-parity conservation, direct searches for R-parity violation ( R p ) signatures in sparticle production have been performed by the LEP2 experiments [3,4]. No evidence for supersymmetric particle production has been observed so far, independently of the hypothesis on R-parity. In 1998, the LEP centre-of-mass energy reached 189 GeV, and an integrated luminosity of 158 pb −1 was collected by the DELPHI experiment. The results of the searches for pair production of supersymmetric particles under the hypothesis of R-parity violating couplings between sleptons and leptons, performed with the data collected by DELPHI in 1997 at a centre-of-mass energy of 183 GeV [3], are updated by the analyses of the data recorded in 1998 presented in this paper.

R-parity violation in the MSSM
The R p superpotential [5] contains three trilinear terms, two violating L conservation, and one violating B conservation. We consider here only the λ ijk L i L jĒk term (i, j, k are generation indices, L (Ē) denote the lepton doublet (singlet) superfields) which couples the sleptons to the leptons; since λ ijk = −λ jik , there are nine independent λ ijk couplings. Upper limits on the λ ijk couplings can be derived from indirect searches of R-parity violating effects [6]- [8], assuming that only one λ ijk is dominant at a time. For example, charged current universality allows a limit on λ 122 to be derived: λ 122 < 0.049 × mẽ R 100 GeV/c 2 and the upper limits on the neutrino mass are used to derive a limit on λ 133 : λ 133 < 0.006 × mτ R 100 GeV/c 2 [9]. Taking into account recent data on neutrino masses and mixings, smaller values of the upper limits on several λ ijk have been derived, all being over 0.0007 (for ml = 100 GeV/c 2 ) [10]. In the analyses described here, only one λ ijk was assumed to be dominant and its upper bound has been taken into account.
The relevant MSSM parameters for these R p searches are: M 1 , M 2 , the U(1) and SU (2) gaugino mass at the electroweak scale, m 0 , the scalar common mass at the GUT scale, µ, the mixing mass term of the Higgs doublets at the electroweak scale and tanβ, the ratio of the vacuum expectation values of the two Higgs doublets. The unification of the gaugino masses at the GUT scale, which implies M 1 = 5 3 tan 2 θ W M 2 ≃ 1 2 M 2 at the electroweak scale, is assumed in the study of production and/or decay processes involving neutralinos and charginos.
We assume that the running of the R p couplings from the GUT scale to the electroweak does not have a significant effect on the evolution of the gaugino and fermion masses. This is an assumption that will be reconsidered once detailed theoretical calculations on this subject become available.

R-parity violating decays
This paper presents the searches for pair produced gauginos and sfermions. In case of pair production R p is conserved at the production vertex; the cross-sections do not depend on the R p couplings. The R-parity violation affects only the decay of sparticles. Two types of supersymmetric particle decays are considered. First, the direct decay, corresponding to the sfermion R p direct decay into two standard fermions, or to the neutralino (chargino) decay into a fermion and a virtual sfermion which then decays into two standard fermions. Second, the indirect decay corresponding to the supersymmetric particle cascade decay through R-parity conserving vertices to on-shell supersymmetric particles down to a lighter supersymmetric particle decaying via one LLĒ coupling.
The direct decay of a neutralino or a chargino via a dominant λ ijk coupling leads to purely leptonic decay products, with or without neutrinos (ℓℓ ′ ν, ℓℓ ′ ℓ ′′ , ℓνν). The indirect decay of a heavier neutralino or a chargino adds jets and/or leptons to the leptons produced in the LSP decay.
The sneutrino direct decay gives two charged leptons: via λ ijk only theν i andν j are allowed to decay directly to ℓ ± j ℓ ∓ k and ℓ ± i ℓ ∓ k respectively. The charged slepton direct decay gives one neutrino and one charged lepton (the lepton flavour may be different from the slepton one): the supersymmetric partner of the right-handed leptonl kR decays directly into ν iL ℓ jL or ℓ iL ν jL , and the supersymmetric partner of the left-handed leptoñ ℓ i(j)L decays intoν j(i)L ℓ kR .
The indirect decay of a sneutrino (resp. charged slepton) into the lightest neutralino and a neutrino (resp. charged lepton) leads to a purely leptonic final state: two charged leptons and two neutrinos (resp. three charged leptons and a neutrino). The indirect decay of a slepton into a chargino and its isospin partner was not considered, and the direct decay of charged slepton is not studied here.
When the charged leptons are τ , additional neutrinos are generated in the τ decay, producing more missing energy in the decay and leading to a smaller number of charged leptons in the final state.

Data samples
The total integrated luminosity collected by the DELPHI detector [11] during 1998 at centre-of-mass energies around 189 GeV was 158 pb −1 . An integrated luminosity of 153 pb −1 has been analysed, corresponding to high quality data, with the tracking detectors and the electromagnetic calorimeters in good working condition.
Signal events were generated with the SUSYGEN 2.20 program [19] followed by the full DELPHI simulation and reconstruction program (DELSIM). A faster simulation (SGV 1 ) was used to check that the efficiencies were stable at points without full simulation compared to their values at the nearest points determined with the full simulation. The R-parity violating couplings were set close to their experimental upper limit derived from the indirect R p searches (see section 1.2).
The χ 0 1 and χ ± 1 pair production was considered at several points in the MSSM parameter space, in order to scan neutralino masses from 15 to 80 GeV/c 2 and chargino masses from 45 to 95 GeV/c 2 . Moreover, for a given mass, several samples with different components and production processes were simulated. The pair production of heavier neutralinos and charginos has been taken into account since one can profit from the threefold increase in luminosity compared to the 1997 data.
For the study of slepton pair production, samples with sneutrino direct decay and samples with sneutrino or charged slepton indirect decay were generated for tanβ fixed at 1.5. Aν (l) mass range from 50 to 90 GeV/c 2 was covered; in the case of indirect decay, several ranges of mass difference between sleptons and neutralinos were considered.

Analysis descriptions 3.1 Analysis strategy and validity
Any of the possible R p signals produced via one of the λ ijk L i L jĒk couplings can be explored by the analyses described in this paper. In the analyses performed considering a dominant λ 133 coupling, the efficiencies and the rejection power are low, due to the presence of several taus in the final state. The highest efficiencies and background reduction are obtained if λ 122 is the dominant coupling. For final states produced by other λ ijk , the detection efficiencies lay between these two limiting cases. Analyses are then performed considering both the λ 122 and the λ 133 couplings. The weakest limits were derived considering the results of the analyses performed assuming a dominant λ 133 coupling. The studied final states are summarized in Table 1.
It was supposed that the Lightest Supersymmetric Particle (LSP) decays within a few centimeters of the production vertex. Since the mean LSP decay length depends on m −5 χ (if the LSP is a gaugino), and on λ −2 ijk , this assumption has two consequences on the analyses described here. First, they were not sensitive to light χ (Mχ LSP ≤ 10 GeV/c 2 ). processes final states with λ 122 final states with λ 133 Table 1: Pair production final states with λ 122 or λ 133 .
Second, analyses looking for neutralino decay products had a lower limit in the sensitivity of the λ coupling of the order of 10 −4 ; below this value, in some area in the MSSM space, the lightest neutralino has a non-negligible lifetime, and the corresponding event topology was not selected by the analyses. Inside the validity domain defined by the upper bound from indirect searches of R p effects and the lower bound due to the LSP flight, the coupling value has no influence on the efficiency of the analyses.

General analysis description
The applied selections were based on the criteria presented in [3], using mainly topological criteria, missing quantities, lepton identification and kinematic properties, and jet characteristics. Compared to the previous analyses, the electron identification has been improved at high energies and in the forward regions of the detector.
As already mentioned, indirect decays of gaugino pairs can add two or more jets to the leptons and missing energy final state, from the hadronic decay of W * and Z * . Moreover, in the case of the λ 133 coupling, thin jets are produced in τ decay. The jets were reconstructed with the DURHAM [20] algorithm. In order to cover the different topologies, the jet number was not fixed, and the jet charged multiplicity could be low (for instance thin jets with one charged particle), or could be zero in case of neutral jets. In the following, the transition value of the y cut in the DURHAM algorithm at which the event changes from a n-jet to a (n − 1)-jet configuration is noted y (n−1)n .
After a brief description of the λ 122 analyses, the selection procedures when λ 133 is the dominant coupling constant are detailed in the following sections.

Analyses applied in case of λ 122 coupling
As already mentioned, these analyses were based on the selection procedure described in [3]; they are not deeply detailed here.

Gaugino and slepton indirect decay searches
One analysis was designed to select leptonic channels with missing energy, with or without jets, in order to study gaugino decays and slepton indirect decays. Events with charged multiplicity greater than three and at least two charged particles with a polar angle between 40 • and 140 • were selected. The missing transverse momentum, p t , had to be greater than 5 GeV/c and the polar angle of the missing momentum to be between 20 • and 160 • . The missing energy had to be at least 0.2 √ s. This set of criteria reduced mainly the background coming from Bhabha scattering and two-photon processes. Then, requirements based on the lepton characteristics were applied: • at least two identified muons were required; • the energy of the most energetic identified lepton had to be greater than 0.1 √ s; • an isolation criterion was imposed for the identified leptons (no other charged particle in a half cone of seven degrees around the lepton); • at least two of the identified leptons had to be leading particles in the jets. One event remained in the data, compared to 1.1±0.3 expected from Standard Model processes contributing to the background (0.6, 0.3 and 0.2 from four-fermion, µ + µ − (γ) and γγ → µ + µ − processes respectively). The selection efficiencies were in the range 35-60% for gaugino decays, in the range 50-60% for sneutrino indirect decays, and in the range 20-50% for charged slepton indirect decays.

Sneutrino direct decay search
A different selection was used to search for final states with no missing energy and at least two muons, resulting from sneutrino direct decays via λ 122 . The thrust value had to be less than 0.95 and the polar angle of the thrust axis had to be between 25 • and 155 • . The total energy from charged particles had to be greater than 0.33 √ s, the missing transverse momentum had to be greater than 2 GeV/c and the missing energy to be less than 0.55 √ s. The charged multiplicity had to be four or six with the total event charge equal to 0. At least two muons were required and no other charged particle in a half cone of 20 • around each lepton was demanded. One event remained in the data after these criteria with 2.7±0.4 expected from standard background processes, mainly from the ℓ + ℓ − ℓ ′ + ℓ ′ − final states (1.4±0.2), and from the γγ → µ + µ − process (1.2±0.4). The efficiencies were from 62% to 51% in the explored sneutrino mass range of 60-90 GeV/c 2 . In the search for pair production of gauginos and sleptons in case of a dominant λ 133 coupling, the following criteria were required:

Analyses applied in case of
• at least one identified lepton; • more than three charged particles and at least two of them with a polar angle between 40 • and 140 • ; • the total energy and the energy from charged particles greater than 0.18 √ s and 0.16 √ s respectively; • the missing p t greater than 5 GeV/c; • the polar angle of the missing momentum between 27 • and 153 • . This was efficient in suppressing the background coming from Bhabha scattering and two-photon processes and in removing a large part of the ffγ contribution. After this preselection stage, 2114 events were selected compared with 1984±11 expected from the background sources (see Figure 1). There was an excess of data mostly concentrated in the low charged multiplicity events where γγ events contributed to the Standard Model background. A good agreement between data and the expected background was obtained when the contribution of γγ events was further reduced (see below).  Table 2: Selection criteria used in the search for neutralino and chargino decay via λ 133 . nj means n-jet topology, and a charged jet means a jet with at least one charged particle. The number of remaining data and Standard Model background events are reported; the quoted errors are statistical.

Neutralino and chargino search
Compared to the selection applied to 1997 data [3], it has been necessary to modify some criteria and to distinguish between low and high multiplicity cases in order to reach a higher purity. For events with a charged particle multiplicity from four to six (which corresponds to neutralino or chargino direct decay), the following criteria were applied: • the energy in a cone of 30 • around the beam axis was restricted to be less than 50% of the total visible energy; • the energy of the most energetic lepton had to be between 2 and 70 GeV; • there should be no other charged particle in a 10 • (6 • ) half cone around any identified lepton for a charged particle multiplicity equal to four (five or six); • the number of neutral particles had to be less than or equal to 10. For events with a charged particle multiplicity greater than six (corresponding to neutralino and chargino indirect decays), the criteria were: • the acollinearity 2 had to be greater than 7 • ; • the energy in a cone of 30 • around the beam axis was restricted to be less than 40% of the total visible energy; 2 the acollinearity is computed between the two vectors corresponding to the sum of the particle momenta in each event hemisphere.
• the energy of the most energetic lepton had to be between 5 and 60 GeV; • if there was only one identified lepton, no other charged particle in a 6 • half cone around it was allowed; and if there were more, there should not be any other charged particle in a 10 • half cone around at least two of them; • at least one well identified electron; • the number of neutral particles had to be less than or equal to 15. In both cases the missing energy had to be at least 30% of the available energy, and the polar angle of at least one lepton had to be between 40 • and 140 • . These criteria removed ffγ and hadronic ZZ and W + W − events.
The selection based on the jet characteristics and topologies was then applied. First, constraints have been imposed to y (n−1)n values to reduce, in particular the ffγ contribution. In events with more than six charged particles, at least one jet with low charged particle multiplicity was required. In four-or five-jet configurations, a minimum of four charged jets was required. In case of a four-jet topology, a cut was applied on the value of E j min × θ jaj b min where E j min is the energy of the least energetic jet, and θ jaj b min is the minimum angle between any pair of jets. These requirements significantly reduced the background from ffγ and W + W − production. The number of remaining real data and background events during the selection are reported in Table 2, and the contributions of the relevant Standard Model processes are detailed in Table 3. The main contribution comes from the W + W − production, with a semi-leptonic decay of the W pair, due to the specific design of the analysis to be efficient for channels with leptons (mainly taus) and jets in final states.  Table 3: Standard Model background contributions to the neutralino and chargino pair production analysis (λ 133 ). The results in the row labelled "Low" ("High") are obtained with the selection applied to the low (high) multiplicity events. The quoted errors are statistical.
Using the events produced with DELSIM, selection efficiencies have been studied on χ 0 1 χ 0 1 and χ + 1 χ − 1 signals. In order to benefit from the high luminosity, all e + e − → χ 0 i χ 0 j , e + e − → χ + k χ − l processes which contribute significantly have been simulated, at each MSSM point of this study. SUSYGEN followed by SGV was used for the scan. Then a global event selection efficiency was determined for each point, since the performed analyses were sensitive to many different topologies. The global selection efficiencies obtained with SGV simulated events have been cross-checked at several points with DELSIM simulated events. The efficiencies laid between 18% and 40%.

Sneutrino and charged slepton searches
Considering the λ 133 coupling, searches for sneutrino pair production and subsequent direct ( ν→ ℓ + ℓ − ) or indirect ( ν→ χ 0 1 ν) decay and searches for charged slepton pair production decaying indirectly (l→ χ 0 1 ℓ) have been performed. In these different searches, a large amount of energy is missing in the final states, due to neutrinos (from τ and/or χ 0 1 decays), except in the case of ν τ ν τ direct decay search (eeτ τ final state). Two different analyses were then performed, one applied to the channels with a large amount of missing energy, and the other one dedicated to the eeτ τ channel, with less missing energy.
• Analysis for channels with high value of missing energy The selection procedure was close to the one applied to low charged particle multiplicity events in the search for neutralino and chargino pair production. The same event characteristics were used. A large amount of missing energy was required, but only events with four to eight charged particles were selected. The criteria are listed in Table 4; the number of observed events and expected ones from the Standard Model background during the selection procedure is also given. At the end, one event remains in the data compared to 2.1±0.3 from the SM processes. The relevant contributions are listed in Table 5. The four-fermion contributions have been checked also with EXCALIBUR, and apart from the WW-like processes, the two other important background sources were the eeτ τ and eeµµ final states.

Selection criteria
Data  Table 4: Selection criteria used in the search for slepton pair production with R p decay via λ 133 . The number of remaining data and Standard Model background events are reported; the quoted errors are statistical.
• Analysis for channels with low value of missing energy In order to obtain higher efficiencies for the eeτ τ channel, the selection criteria were modified. In particular, the missing energy cut was reduced to 8% of the available energy.
The number of charged particles was restricted to be between four and six. Moreover, the lower limit on the energy of the most energetic lepton was increased to 20 GeV, and the isolation angle had to be greater than 10 • . Additional criteria were used: the acollinearity had to be greater than 2 • , and the presence of at least one identified electron was required. After the event selection 3 events remained, while 2.3±0.3 events were expected from SM processes. The main sources of background were the eeτ τ (57%) and the eeµµ (29%) four-fermion processes. The efficiencies were between 38% and 46%.

Interpretation of the results
The results of the searches presented in this paper, summarised in Table 6, were in agreement with the Standard Model expectation. They were used to extend the previously excluded part of the MSSM parameter space and to update limits obtained with the analysis of the 1997 data collected in DELPHI. In all the pair production processes

Coupling Process
Efficiency Selected events range in % Data MC λ 122  Table 6: LLĒ analyses: efficiency ranges in the different cases studied, and data and Monte Carlo events remaining after the applied selection.
studied, the weakest limits were derived from the results of the λ 133 analyses, and are hence valid for any choice of dominant λ ijk coupling, provided that the coupling is strong enough for the LSP to decay within a few centimetres.
In the searches for neutralino and chargino pair production, the number of expected events at each point of the explored MSSM parameter space was obtained by: } where L is the integrated luminosity, and ǫ g is the global efficiency determined as explained in section 3.4.2. This number has been compared to the number of signal events, N 95 , expected at a confidence level of 95% in presence of background [21]. All points which satisfied N exp >N 95 were excluded at 95% C.L. The excluded area in µ, M 2 planes obtained with the present searches are shown in Fig. 2, for m 0 = 90 GeV/c 2 (the t-channel contribution to the gaugino cross-sections has an important effect), m 0 = 300 GeV/c 2 (the t-channel contribution vanishes) and tanβ =1.5, 30. The smaller excluded area in the µ, M 2 planes for a given tanβ is obtained for high m 0 values.
For each tanβ, the highest value of neutralino mass which can be excluded has been determined in the µ, M 2 plane (-200 GeV/c 2 ≤ µ ≤ 200 GeV/c 2 , 5 < M 2 ≤ 400 GeV/c 2 ) for several m 0 values varying up to 500 GeV/c 2 . The smaller excluded area in the µ, M 2 plane is obtained for m 0 = 500 GeV/c 2 . The most conservative mass limit was obtained for high m 0 values, for which it reaches a plateau. The corresponding limit on neutralino mass as a function of tanβ is plotted in Fig. 3. From these studies, a neutralino lighter than 30 GeV/c 2 was excluded at 95% C.L. for 1 ≤tanβ ≤ 30. The same procedure was applied to determine the most conservative lower limit on the chargino masses. The result is less dependent on tanβ, almost reaching the kinematic limit for any value of tanβ: a chargino lighter than 94 GeV/c 2 was excluded at 95% C.L. Finally, using the same method, a lower limit of 50 GeV/c 2 for the χ 0 2 mass has been derived at 95% C.L.
The results obtained from the sneutrino pair production studies were used to derive limit on the sneutrino mass. In the case of the sneutrino direct decay, the results improved the upper limit on the sneutrino pair production cross-section. Taking into account the results of the two analyses and the efficiencies obtained when varying the sneutrino mass, the cross-section limits for 2e2τ and 4τ channels were derived and are reported in Fig. 4. The ν e ν e cross-section depends not only on the ν e mass but also on other MSSM parameters (due to the possible t−channel χ + 1 exchange contribution) and it is plotted for a specific MSSM point: M 2 = 100 GeV/c 2 and µ = -200 GeV/c 2 . The upper limit on the cross-section leads to a lower limit on the sneutrino mass of 78 GeV/c 2 .
In the case of theν indirect decay into ν χ 0 1 with the R p decay of the neutralino via λ 133 , the efficiencies depend on the sneutrino and neutralino masses. The search results allowed an area in the mχ0 versus mν plane to be excluded, as shown on Fig. 5. The same procedure has been followed for the charged slepton indirect decays. The indirect decay of aτ pair gives two taus and two neutralinos, and the final state selection was less efficient than for theẽ orμ pair; the results obtained for the τ R pair production gave the most conservative limits on the slepton mass for any flavour, assuming thatl R decays exclusively to ℓ χ 0 1 . The area excluded in the mχ0 versus ml R plane is plotted in Fig. 6. The region where ml R -mχ0 is less than 2-3 GeV/c 2 was not covered by the present analysis, since then the direct decay becomes the dominant mode, leading to two leptons and missing energy. Taking into account the limit on the neutralino mass at 30 GeV/c 2 , sneutrinos with mass lower than 76.5 GeV/c 2 and supersymmetric partners of the right-handed lepton, decaying indirectly, with mass lower than 83 GeV/c 2 were excluded at 95% C.L.

Summary
Searches for R-parity violating effects in e + e − collisions at √ s = 189 GeV have been performed with the DELPHI detector. The pair productions of neutralinos, charginos and sleptons have been studied under the assumption that the LLĒ term is responsible for the supersymmetric particle decays into standard particles. It was assumed that one λ ijk coupling is dominant at a time and that the λ ijk coupling is strong enough for the LSP to decay within a few centimetres. No evidence for R-parity violation has been observed, allowing to update the limits previously obtained at √ s =183 GeV. The present 95% C.L.
These limits are valid for tanβ ≥ 1 and m 0 < 500 GeV/c 2 and for all the generation indices i,j,k of the λ ijk coupling, and for any coupling value from 10 −4 up to the existing limits.