%================================================================
% LaTeX file with preferred layout for the contributed papers to
% the ICHEP Conference 98 in Vancouver
% process with:  latex hep98.tex
%                dvips -D600 hep98
%================================================================
\documentclass[12pt]{article}
\usepackage{epsfig}
\usepackage{amsmath}
\usepackage{hhline}
\usepackage{amssymb}
\usepackage{times}

\renewcommand{\topfraction}{1.0}
\renewcommand{\bottomfraction}{1.0}
\renewcommand{\textfraction}{0.0}
\newlength{\dinwidth}
\newlength{\dinmargin}
\setlength{\dinwidth}{21.0cm}
\textheight23.5cm \textwidth16.0cm
\setlength{\dinmargin}{\dinwidth}
\setlength{\unitlength}{1mm}
\addtolength{\dinmargin}{-\textwidth}
\setlength{\dinmargin}{0.5\dinmargin}
\oddsidemargin -1.0in
\addtolength{\oddsidemargin}{\dinmargin}
\setlength{\evensidemargin}{\oddsidemargin}
\setlength{\marginparwidth}{0.9\dinmargin}
\marginparsep 8pt \marginparpush 5pt
\topmargin -42pt
\headheight 12pt
\headsep 30pt \footskip 24pt
\parskip 3mm plus 2mm minus 2mm

%===============================title page=============================

% Some useful tex commands
%
\newcommand{\GeV}{\rm GeV}
\newcommand{\TeV}{\rm TeV}
\newcommand{\pb}{\rm pb}
\newcommand{\cm}{\rm cm}
\newcommand{\hdick}{\noalign{\hrule height1.4pt}}

\begin{document}

\pagestyle{empty}
\begin{titlepage}

\noindent
\begin{center}
%{\it {\large version of \today}} \\[.3em] 
\begin{small}
\begin{tabular}{llrr}
%Submitted to & \multicolumn{3}{r}{\footnotesize Electronic Access: {\it http://www-h1.desy.de/h1/www/publications/conf/conf\_list.html}} \\[.2em] \hline 
%Submitted to & \multicolumn{3}{r}{\footnotesize {\it www-h1.desy.de/h1/www/publications/conf/conf\_list.html}} \\[.2em] \hline 
Submitted to & & &
\epsfig{file=H1logo_bw_small.epsi
,width=2.cm} \\[.2em] \hline
\multicolumn{4}{l}{{\bf
                31st International Conference 
                on High Energy Physics, ICHEP02},
                July~24,~2002,~Amsterdam} \\
                 & Abstract:        & {\bf 1019}    &\\
                 & Parallel Session & {\bf 4,10}   &\\ \hline
 & \multicolumn{3}{r}{\footnotesize {\it
    www-h1.desy.de/h1/www/publications/conf/conf\_list.html}} \\[.2em]
\end{tabular}
\end{small}
\end{center}
\vspace*{2cm}

\begin{center}
  \Large
{\bf Multi-electron Production at High Transverse Momentum \\
in $ep$ collisions at HERA }

  \vspace*{1cm}
    {\Large H1 Collaboration} 
\end{center}

\begin{abstract}

\noindent
Multi-electron production has been measured at high transverse momentum 
in 115 pb$^{-1}$ of positron- and electron-proton collisions 
collected by the H1 experiment at HERA.
Good overall agreement is found with the Standard Model predictions,
dominated by photon-photon interactions. Six events are observed with
 a di-electron mass above 100 GeV, a domain where the Standard Model 
 prediction is low.
\end{abstract}

\end{titlepage}

\pagestyle{plain}
\section{Introduction}
The measurement of rare processes in 
electron\footnote{In this paper the term
``electron'' is used to describe generically electrons or positrons.}-proton 
interactions at HERA provides a unique method to 
search for new physics. In this paper we report on the first measurement 
of multi-electron production performed at high transverse momentum ($P_T$) 
at HERA. Within the Standard Model (SM) 
the production of multilepton events in $ep$ collisions is possible 
mainly through photon-photon interactions, where quasi-real photons 
radiated from the incident electron and proton interact for producing 
a pair of leptons 
$\gamma\gamma\rightarrow \ell^+ \ell^-$~\cite{verm}. 
An excess of events in a given kinematic domain compared to the SM 
might be a signature for new phenomena.
% or the production of supersymmetric particles\cite{SUSY}. 

\par
The analysis is based on the $e^{\pm}p$ collisions recorded 
by the H1 experiment in the 1994-2000 period.
The total integrated luminosity of 115.2 $pb^{-1}$ is shared between
36.5 $pb^{-1}$ and 65.1 $pb^{-1}$ of $e^+p$ collisions recorded at 
center of mass energies $\sqrt{s}$ of 300 GeV and 318 GeV respectively, and 
13.6 $pb^{-1}$ of $e^-p$ collisions recorded at $\sqrt{s}$ = 318 GeV. 

\section {Experimental conditions} 
 
A detailed description of the H1 detector can be found in \cite{H1detector}.
Only components essential for this analysis are shortly described
here\footnote{The origin of the H1 coordinate system is the nominal
$ep$ interaction point. The direction of the incoming proton defines
the z-axis for the measurement of polar and azimuthal angles as well
as for transverse momenta. It also defines the forward region of the
detector. 
%The z-axis is completed by a vertical y-axis, oriented
%to the top, and a horizontal x-axis to form a direct reference system.
}. 

A tracking system of central and forward drift chambers is used to measure the
charged particle trajectories and to determine the interaction vertex. 
Particle transverse momenta are determined from the curvature of
the trajectories in a magnetic field of 1.15 Tesla.
%yields  a measurement of the particle momenta with a nominal
%resolution  $\sigma_{P}/P^2=3\times10^{-3}$ GeV$^{-1}$.
 
Hadronic and electromagnetic final state particles
are absorbed in a highly segmented liquid argon calorimeter
\cite{h1cal} covering the polar angular range 
$4^\circ < \theta < 153^\circ$. The calorimeter is 5 to 8 hadronic interaction
lengths deep depending on the polar angle of the particle.
It includes an electromagnetic section which is 20 to 30 
radiation lengths deep.
Electromagnetic shower energies are measured with a precision of
$\sigma (E)/E = 12 \% / \sqrt{E/\mathrm{GeV}} \oplus 1\%$ and
 hadronic shower energies with
$\sigma (E)/E = 50 \% / \sqrt{E/\mathrm{GeV}} \oplus 2 \%$ as measured in test beams~\cite{h1calotestbeams}.
The electromagnetic energy scale is known to the level of 0.7 to 3\%
depending on the polar angle. 
The hadronic energy scale is understood at the 2\%
level.
% when comparing to the Monte Carlo expectations.

In the backward region the liquid argon calorimeter is completed
by a lead/scintillating-fibre calorimeter
\footnote{This device was installed in 1995 replacing a lead-scintillator
``sandwitch'' calorimeter~\cite{H1detector}.}\cite{SPACAL} covering
the range $155^\circ < \theta < 178^\circ$, and by electron and photon
taggers located downstream of the interaction point in the electron 
beam direction.
The calorimeter is surrounded by a superconducting coil
and an iron yoke instrumented with streamer
tubes. Leakage of hadronic showers outside the  calorimeter is measured
 from an analogue  charge sampling of the streamer tubes (Tail Catcher).
% with a resolution of $\sigma (E)/E = 100 \% / \sqrt{E/\mathrm{GeV}}$ .
%Tracks of
%penetrating charged particles, such as muons, escaping the calorimeter
% are reconstructed from
%their hit pattern in the streamer tubes with an efficiency higher than 90\%.
% 
%The trigger condition is based on calorimeter signals for multi-electron analysis and on muon detector and central  tracker measurements for multi-muon analysis.
% The trigger efficiency is higher than 95\% for events with a electron with an energy above 10 GeV, and typically between 70 and 90\% depending on running conditions for muon selection.

The analysis is restricted to periods  where the  central drift chambers, 
the liquid argon calorimeter and the luminosity system are  active 
and fully operational. For the analysed topologies trigger conditions are
derived from the liquid argon calorimeter signals. The trigger efficiency 
is higher than 95\% for events with an electron energy above 10 GeV.
%and is controlled with a precision of 3 \%.

\section{Standard Model Processes and their Simulation}

The main SM processes involved in multi-electron production at HERA
are summarized in figure~\ref{diagrams}. The dominant contribution is
the interaction of two photons radiated from the incident electron and proton,
either in the elastic or quasi-elastic mode (figure ~\ref{diagrams}a), 
or in the inelastic mode (figure ~\ref{diagrams}b).   
The Cabibbo-Parisi process (figure ~\ref{diagrams}c), which can also be elastic 
or inelastic, involves an $e^+e^-$
interaction where one of the electrons is issued from a photon radiated 
from the proton. Its contribution is one order of magnitude lower, 
except at high transverse
momentum where it contributes more significantly. The Drell-Yan process
(figure~\ref{diagrams}d) involves a quark-antiquark interaction where 
one of the quarks is issued from a photon radiated from the incident electron.
Its contribution is negligible in the entire measured phase 
space~\cite{romero}.

Multi-lepton event production is simulated with the GRAPE Monte Carlo 
generator~\cite{grape}. GRAPE is based on the symbolic graphs calculation 
system GRACE~\cite{GRACE}. It simulates the full $2\rightarrow 4$ electroweak 
matrix elements at tree level, except for the negligible Drell-Yan 
contribution, for interactions resulting in three leptons plus hadrons 
in the final state. 
GRAPE contains $\gamma\gamma$ graphs with interference effects between 
the scattered and produced electrons. It takes into account 
the Cabibbo-Parisi process, photon internal conversions (diagrams containing a radiated photon with subsequent $\gamma\rightarrow e^+e^-$ vertex) and $Z^0$ production or virtual 
exchange. Initial and final state radiation processes 
(QED and QCD - parton showers) 
are simulated in the leading log approximation. The proton dissociation is 
treated in three phase space regions: elastic, quasielastic and inelastic.
In order to simulate experimental distributions the GRAPE generator was 
interfaced to the H1 detector simulation. 
%Tau decay contributions to  
%multi-electron final states are included in the simulated samples.

GRAPE was cross-checked with the generator LPAIR~\cite{LPAIR}, 
which contains only the photon-photon process generated in the 
Weizs\"acker-Williams approximation without interference terms
or weak contributions. When restricted to the same processes,
GRAPE agrees with LPAIR at the percent level for both total 
and differential cross-sections. The additional diagrams in GRAPE
increase the predicted cross-section by 20 \% in average, and up to
40 \% for di-lepton masses either very low (photon internal conversions) 
or around 90 GeV ($Z^0$ production). At masses higher than the $Z^0$ mass,
the difference between GRAPE and LPAIR is dominated by the weak interference,
which has a 10\% negative effect for the dominant $e^+p$ sample.
 
The main experimental backgrounds to multi-electron production are 
processes where, in addition to a true electron, one or more fake 
electrons may be reconstructed from the final state particles. 
The dominant contribution is expected from
Neutral Current Deep Inelastic Scattering (NC-DIS) where a fake 
electron from the hadrons or some radiated photon is selected together 
with the scattered electron. Elastic Compton scattering can also 
contribute if the high $P_T$ photon is misidentified as an electron. 
The NC-DIS and elastic Compton processes are simulated with the 
DJANGO~\cite{DJANGO} and WABGEN~\cite{WABGEN} generators respectively.

The Monte-Carlo predictions are attributed theoretical
and experimental systematic uncertainties added in quadrature. 
The theoretical errors reflect the available freedom in
the proton parton distributions and in the choice of 
internal cuts between different phase 
space regions used in the generators. The main experimental uncertainty 
is due to the tracker efficiency in the electron identification procedure
(see below). Uncertainties on the energy scales of the calorimeters, the trigger 
efficiency and the luminosity measurement are also taken into account. 
The prediction of fake electron background contributions has been studied in control analyses described in the next section and the corresponding systematical error is taken into account.
%attributed 
%an error of 20 \% (see below).  

\section{Event Selection}

The multi-electron event selection is based on an electron identification
procedure designed to minimize the contribution of fake electrons,
while keeping a high efficiency for true electrons and a reliable
control of the overall selection performance.
 
Electrons with energies above 5 GeV are first preselected in the fiducial
volume of the liquid argon and SPACAL calorimeters, covering the polar
angular range $5^\circ < \theta < 175^\circ$. Isolated electromagnetic showers 
are sorted out using pattern recognition algorithms based on the 
geometric profiles expected from electrons. 
The remaining clusters~-~related to hadronic activity~-~are combined into hadronic jets using a
$k_T$ algorithm~\cite{ktalgo}. 
In the region of overlap between the liquid argon calorimeter and the
central drift chamber ($20^\circ < \theta < 150^\circ$), the calorimetric
identification is complemented by a tight tracking condition, requiring 
that an isolated high quality track matches the electromagnetic cluster 
both geometrically and in momentum. This additional constraint strongly
reduces the contribution of fake electrons issued from photons and hadrons.
Electrons selected in this polar angular 
range are called ``central electrons'' hereafter. 
In the forward ($5^\circ < \theta < 20^\circ$)
and backward ($150^\circ < \theta < 175^\circ$) regions no explicit track requirement 
is made. The forward electron energy threshold is raised to 10 GeV in 
order to minimize the contribution of NC-DIS forward jets to fake electrons. 
Once identified, only electron candidates isolated from each other
and from hadronic jets by more than 0.5 radians in the pseudorapidity-azimuth
plane are considered.
% ($D^{\mathrm e-e',e-jets}_{\eta\phi}>0.5$).
 
The performance and reliability of this harsh electron identification
procedure was controlled in several ways.~The tracking efficiency 
 was measured with a NC-DIS sample selected with no tracking conditions 
for the electron candidate, requiring a single electromagnetic cluster 
in the calorimeter with a transverse momentum above 10 GeV in the polar 
angle interval $20^\circ < \theta < 150^\circ$. 
Figure ~\ref{track_efficiency} shows
that the average efficiency is 90 \% and varies only slightly with
the track momentum and polar angle. The decrease observed in the 
forward region is due to the electron showering. The measured
efficiency is described by the simulation with an accuracy
ranging from 3 \% in the central region to 15 \% at the edge of the
angular acceptance of the central tracker.
Electron misidentification was studied by measuring the probability
to select a second electron in NC-DIS events when relaxing the 
track quality criteria. In case no tracking conditions are 
applied to the second electromagnetic cluster, fake electrons are dominated by 
photons of elastic Compton events (figure~\ref{fake_e2} left). 
This contribution is drastically reduced by a loose tracking requirement 
(figure~\ref{fake_e2} right). In both cases fake electron background
is described by the simulation to within 20 \%.
Photon conversions were investigated using a sample enriched with elastic 
Compton events, selected by requiring one central electron plus a second electromagnetic cluster (photon candidate) and no significant extra energy in the calorimeters.
The charged tracks associated to the photon candidate are studied 
in figure~\ref{photon_conversion}. 
Their number (figure ~\ref{photon_conversion} left) and 
starting radius in the transverse plane 
(figure~\ref{photon_conversion} right)
are well reproduced by the simulation.  
On figure~\ref{photon_conversion} right, the central tracker structure is 
visible as peaks in the distribution, corresponding to the photon conversions 
in the tracker walls. These conversions are well described.  

The selection of the multi-electron events requires
two central electron candidates,
out of which one must have $P_T > 10$ GeV and the second one
$P_T > 5$ GeV. Additional electron candidates are identified in
the forward, central or backward regions with no explicit $P_T$ cut.  
The electron candidates are indexed 
in decreasing $P_T$: $P_T^{e_i} > P_T^{e_{i+1}}$.

The selected events are classified as di-electrons (``2e'') in the case where
only the two central electron candidates are visible, and tri-electrons
(``3e'') in the case where exactly one additional electron candidate 
is identified. No event is observed with more that one additional electron.
In addition, a subsample of the ``2e'' sample, labelled ``$\gamma\gamma$'',
is selected in order to measure the photon-photon 
cross section in a well defined phase space with low background.
In this subsample, the two electrons must be of opposite charge,
and a significant deficit compared to the initial state 
must be observed in the difference $E - P_z$ of the energy and longitudinal momentum 
of all visible particles ($E-P_z<45$~GeV) 
\footnote{For events where only longitudinal momentum
along the proton direction is undetected, one expects \\
$E - P_z = 2E_e = 55$~GeV, where $E_e$ is the energy of the incident
electron. The threshold $E-P_z<45$~GeV corresponds to a cut on the fractional energy loss $y=\frac{E-P_z}{2E_e^0}<0.82$.}.
These two conditions ensure that the incident electron
is lost in the beam pipe after radiating a quasi-real photon of
squared momentum $Q^2$ lower than 1 GeV$^2$.

% ${ \gamma\gamma\rightarrow e^+e^-}$ 

\section{Results}

The event yields for each sample are summarized in table~\ref{metab1}.
The observed numbers of di-electron and tri-electron events are in agreement 
with the expectations. No event is seen with four or more identified electrons.

The distributions of the global event variables 
are presented in figure~\ref{global_variables}. 
The 3e events accumulate at $E - P_z$ values around 55 GeV, 
as expected since the scattered electron is visible in the detector.
The 2e events show a tail at lower $E - P_z$, 
due to the incident electron being lost in the beam pipe,
corresponding to the $\gamma\gamma$ topology described above. 
The transverse momentum is balanced within resolution as quantified by the 
vectorial sum $P_T^{miss}$ of the transverse momenta of all visible particles.
%No event shows evidence for undetected particles emitted with 
%a significant $P_T$.
The  $P_T^{miss}$ distribution is consistent with the expectation for no
emission of undetected particles with substantial transverse momentum.
The inelasticity spectrum, as measured from the transverse momentum $P_T^{hadrons}$ 
of all visible particles except identified electrons, is also well 
described by the SM.

The distributions of the individual electron transverse momenta $P_T^{e_i}$ 
are steeply falling as shown in figure~\ref{transverse_momenta}. 
The 2e and 3e samples are in good overall agreement with the SM.
Three 2e events 
are seen with $P_T^{e_1}$ above 50 GeV, a domain where the SM expectation 
is low. The electron polar angle $\theta^{e_i}$ distributions are rather 
uniform (figure ~\ref{polar_angles}) and in agreement with the SM. 

The distribution of the invariant mass $M_{12}$ of the two highest 
$P_T$ electrons in the event, and its correlation with respect to 
$P_T^{e_1} + P_T^{e_2}$, are shown in figure~\ref{m12}. 
Three 2e and three 3e events are seen with
invariant masses $M_{12}$ above 100 GeV, where the SM expectation is low.
The three 2e events are the same as those observed at high $P_T^{e1}$.
The comparison to the SM expectations of the observed event rates 
for masses above 100 GeV is presented in table~\ref{metab2}. 
The invariant masses of the other possible electron combinations 
in the 3e sample are shown in figure~\ref{m123}. 
No event is seen with an unexpected high $M_{13}$ or $M_{23}$.
The three high $M_{12}$ events have also a high 
tri-electron mass $M_{123}$.

%The six high-mass events observed in the data are discussed 
%in more detail in the following section.

The six events with $M_{12}>100$~GeV were recorded in the
positron-proton collisions. They are shown 
%together in figure~\ref{meevdall} and
 individually in figures~\ref{meevd1} to~\ref{meevd6} . 
For these events, all available detector signals support
the interpretation of the electron candidates as being electrons.
The electromagnetic shower shapes were checked individually 
to be as expected from the calorimeter response to electrons.  
All central tracks yield a specific ionization
in the central drift chamber as expected from  single electrons.
The forward electron candidates in events 4, 5 and 6 all have at
least one track pointing to the calorimetric energy cluster, though no such
requirement was set in their identification.  
The measurements of the central electron momenta by the 
tracker and the calorimeter are compatible with each other 
within errors.

Though classified as 2e, event 1 also contains a third
electron candidate with energy below the identification threshold.
Similarly event 3 has a compact electromagnetic energy deposit 
located at the forward boundary of the liquid argon calorimeter, 
outside the electron identification fiducial volume. This event also
shows a low energy converted photon radiated close to electron $e_2$.
%All events have a configuration of activity in the detector
%characteristic of an elastically scattered proton, except event 6 which
%has a significant forward hadronic jet. 
With the exception of event 6, which has a significant forward hadronic jet, the events show no hadronic activity in the detector.
%The electron and photon 
%taggers are empty for all six events. Within resolution, the events 
%are balanced in transverse momentum $P_T^{miss}$ and, apart for event 2, 
%show no significant deficit in $E - P_z$. 
It should be noted that the 
high mass electron pair topology differs for events classified 
as ``2e'' and ``3e'': in the former case, the high-mass pair 
is formed of two central high-$P_T$ electrons,
whereas in the latter case, it contains one forward 
and one central electron of intermediate $P_T$'s (figure~\ref{m12}). 

The photon-photon differential cross-sections measured with the 
$\gamma\gamma$ sample as function of $P_T^{hadrons}$, $P_T^{e_1}$
and $M_{12}$ are shown 
in figures~\ref{sigma_ptx},~\ref{sigma_pte}~and~\ref{sigma_m12} respectively.
These cross sections are extracted in the restricted kinematical domain
defined by $20^\circ < \theta^{e_{1,2}} < 150^\circ$,
$P_T^{e_1}>10$~GeV, $P_T^{e_2}>5$~GeV, incident electron 
fractional energy loss $y<0.82$, and $Q^2<1$~GeV$^{2}$.
A good agreement with the SM is found for all three variables
\footnote{It is worth noting that the agreement in $P_T^{hadrons}$ 
confirms our understanding of the SM expectation for the excess of
events with a high $P_T$ lepton
and missing transverse momentum $P_T^{miss}$ observed at 
high $P_T^{hadrons}$ by H1 in a separate analysis 
\cite{isolated_leptons}.}.


\section{Summary}

High-$P_T$ multi-electron production has been measured for
the first time in $ep$ scattering at HERA. 
No event was found with four or more identified electrons.
The di-electron and tri-electron event yields are in good overall agreement 
with the Standard Model predictions, dominated by photon-photon 
interactions. The differential cross-sections for the photon-photon
process have been extracted in a restricted phase space and
found to agree with the predictions.

For masses of the highest $P_T$ electron pair above 100 GeV, 
three events classified as di-electrons, and three events 
classified as tri-electrons are seen, 
compared to SM expectations of $0.25 \pm 0.05$ and
$0.23 \pm 0.04$ respectively. 
%The kinematic properties
%and electric charge structure of the six high-mass
%electron pairs discard their interpretation    
%in term of the decay of a single narrow resonance.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section*{Acknowledgements}

We are grateful to the HERA machine group whose outstanding
efforts have made and continue to make this experiment possible. 
We thank
the engineers and technicians for their work in constructing and now
maintaining the H1 detector, our funding agencies for 
financial support, the
DESY technical staff for continual assistance, 
and the DESY directorate for the
hospitality which they extend to the non DESY 
members of the collaboration.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\begin{thebibliography}{99}



%Vermaseren paper
\bibitem{verm} J.A.M.~Vermaseren, Nucl.\ Phys.\ {\bf B229} (1983) 347.

%Double charged Higgs paper
%\bibitem{Higgs} 
%E. Accomando and S. Petrarca, Phys. Lett. B323 (1994) 212.

%SUSY paper
%\bibitem{SUSY} 
%SUSY lepton pairs in ep collisions.

% H1 detector

\bibitem{H1detector} 
H1 Collaboration, I. Abt et al., Nucl. Instr. and Meth.
A386 (1997) 310 and 348. 

\bibitem{h1cal} 
H1 Calorimeter Group, B. Andrieu et al.,
Nucl. Instr. and Meth. A336 (1993) 460.

%
\bibitem{h1calotestbeams}
 H1 Calorimeter Group, B. Andrieu et al.,
 Nucl. Inst. and Meth. A344 (1994) 492;
 {\it idem}, Nucl. Instr. and Meth. A350 (1994) 57;
 {\it idem}, Nucl. Instr. and Meth. A336 (1993) 499.
%
%\bibitem{H1CALRES}
% H1 Calorimeter Group, B.~Andrieu et al.,
% Nucl. Instr. and Meth. A350 (1994) 57;
% {\it idem}, Nucl. Instr. and Meth. A336 (1993) 499.


\bibitem{SPACAL} 
H1 Spacal Group, R.D. Appuhn et al.,
Nucl. Instr. and Meth. A386 (1997) 397.

% Romero-Kessler
\bibitem{romero} 
N.~Artega-Romero, C.~Carimalo, P.~Kessler, Zeit.f.Phys {\bf C52} (1991) 289.

% MC generators

\bibitem{grape} 
T.~Abe, Comp.\ Phys.\ Comm. {\bf 136} (2001) 126,   
http://www.awa.tohoku.ac.jp/$\sim$tabe/grape/.

\bibitem{GRACE} 
T. Ishikawa et al., KEK preprint 92-19 (1993).

\bibitem{LPAIR} 
S. Baranov, O. D\"unger, H. Shooshtari, J.A.M. Vermaseren,
Proceedings of the workshop Physics at HERA, Vol. 3 Monte-Carlo Generators,
Eds. W. Buchm\"uller, G. Ingelmann (1991). 

\bibitem{DJANGO}
 DJANGO~2.1;
 G.A.~Schuler and H.~Spiesberger,
 Proceedings of the Workshop Physics at HERA,
 W.~Buchm\"uller and G.~Ingelman (Editors),
 (October 1991, DESY-Hamburg), vol. 3 p. 1419.

\bibitem{WABGEN} 
A new Generator for Wide Angle Bremsstrahlung: Ch. Berger and P. Kandel,
Proc. of the Monte Carlo Generators for HERA Physics Workshop, 
DESY-PROC-1999-02, p 596.

\bibitem{ktalgo}
%\bibitem{Ellis:tq}
S.~D.~Ellis and D.~E.~Soper,
%``Successive Combination Jet Algorithm For Hadron Collisions,''
Phys.\ Rev.\ D {\bf 48} (1993) 3160
[arXiv:hep-ph/9305266].\\
%%CITATION = HEP-PH 9305266;%%
%\cite{Adloff:1998ni}
%\bibitem{Adloff:1998ni}
C.~Adloff {\it et al.}  [H1 Collaboration],
%``Measurement of internal jet structure in dijet production in deep-inelastic scattering at HERA,''
Nucl.\ Phys.\ B {\bf 545} (1999) 3
[arXiv:hep-ex/9901010].
%%CITATION = HEP-EX 9901010;%%



\bibitem{isolated_leptons}
H1 Collaboration, C. Adloff et al., Eur. Phys. J. C5 (1998) 575.

\end{thebibliography}
\clearpage

\begin{table}[htb]
\begin{center}

\begin{tabular}{|l|c|c||c|c|}
\multicolumn{2}{l}{\Large \bf  H1 Preliminary $115$~pb$^{-1}$} &
\multicolumn{3}{r}{\Large \bf  Multi-electron Analysis} \\ [0.2cm] \hline
\hline
   Selection   & DATA  &  SM &  GRAPE &  NC-DIS + Compton \\ \hline \hline
  Visible 2e           &   105 &  $118.2\pm  12.8$ &  $ 93.3\pm  11.5$ &  $ 25.0\pm   5.5$ \\
  Visible 3e           &    16 &  $ 21.6\pm   3.0$ &  $ 21.5\pm   3.0$ &  $  0.1\pm   0.1$ \\
  Visible 4e or more   &     0 &  $  0.1\pm   0.0$ &  $  0.1\pm   0.0$ &  $  0.0\pm   0.0$ \\
   \hline
 $\gamma\gamma \rightarrow e^+e^-$ subsample  &    41 &  $ 48.3\pm   6.1$ &  $ 46.4\pm   6.1$ &  $  1.9\pm   0.9$ \\
   \hline
  
 \end{tabular}

\caption{
Observed and predicted multi-electron event rates for the samples described 
in the text. The prediction errors include model 
uncertainties and experimental systematical errors added in quadrature.
}
\label{metab1}
\end{center}
\end{table}


\begin{table}[htb]
\begin{center}
 
\begin{tabular}{|c|c|c||c|c|}
\multicolumn{2}{l}{\Large \bf  H1 Preliminary $115$~pb$^{-1}$} &
\multicolumn{3}{r}{\Large \bf  Multi-electron Analysis} \\ [0.2cm] \hline
   Selection   & DATA  &  SM &  GRAPE &  NC-DIS + Compton \\ \hline \hline
 Visible 2e $M_{12}>100$~GeV &   3 &  $ 0.25\pm  0.05$ &  $ 0.21\pm  0.04$ &  $ 0.04\pm  0.03$ \\
 Visible 3e $M_{12}>100$~GeV &   3 &  $ 0.23\pm  0.04$ &  $ 0.23\pm  0.04$ &  $ 0.00\pm  0.00$ \\
   \hline
 \end{tabular}
 
\caption{ Observed and predicted multi-electron event rates for 
masses $M_{12}>100$~GeV, for the samples described in the text.
The prediction errors include model uncertainties and experimental 
systematical errors added in quadrature.}
\label{metab2}
\end{center}
\end{table}


\newpage
% SM diagrams for lepton pairs
%
\begin{figure}[p] 
  \begin{center} \vspace*{-.5cm}
      \epsfig{file=H1prelim-02-052.diagrams.eps,width=14cm}
  \end{center} \vspace*{-.3cm}
  \caption{ Main processes involved in lepton pair production
in $ep$ collisions: a) elastic photon-photon; b) inelastic photon-photon;
c) Cabibbo-Parisi; d) Drell-Yan.}
  \label{diagrams}
\end{figure} 


%%%%%%%%%%%%%%%%%%%%%%%% CONTROL FIGURES


\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.control1.eps,width=14cm}%
  \end{center}
  \caption{Electron track association efficiency as function of the electron 
transverse momentum (left) and electron polar angle (right).}
  \label{track_efficiency}
\end{figure} 

\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.control2.eps,width=14cm}%
  \end{center}
  \caption{ Momentum spectrum of second electrons identified in NC-DIS events
with either no track requirement (left) or a loose track requirement (right),
compared to expectations.}
  \label{fake_e2}
\end{figure} 

\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.control3.eps,width=14cm}%
  \end{center}
  \caption{ Number (left) and starting radius (right)
of the tracks associated to photon candidates
in an elastic Compton sample.}
  \label{photon_conversion}
\end{figure} 


%%%%%%%%%%%%%%%%%%%%%%% DATA ANALYSIS CROSS SECTIONS
%
\begin{figure}[p] 
  \begin{center} \vspace*{-.5cm}
      \epsfig{file=H1prelim-02-052.fig1.eps,width=18cm}
  \end{center} \vspace*{-.3cm}
  \caption{ Global event variable distributions (see text) for 
events classified as di-electrons (top) and tri-electrons (bottom).}
  \label{global_variables}
\end{figure} 


\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.fig2.eps,width=18cm}%
  \end{center}
  \caption{ Electron transverse momentum distributions 
for events classified as di-electrons (top) and tri-electrons (bottom).}
  \label{transverse_momenta}
\end{figure} 

\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.fig3.eps,width=18cm}%
  \end{center}
  \caption{ Electron polar angle distributions for events classified 
as di-electrons (top) and tri-electrons (bottom).}
  \label{polar_angles}
\end{figure} 

\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.fig4.eps,width=16cm}%
  \end{center}
  \caption{ Distribution of the invariant mass $M_{12}$ of the two 
highest $P_T$ electrons (left) and its correlation with
the sum of the electron $P_T$'s (right), for events 
classified as di-electrons (top) and tri-electrons (bottom). 
The dots on the right plots represent the Standard Model prediction 
for a luminosity 500 times higher than in the data.}
  \label{m12}
\end{figure} 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.fig5.eps,width=16cm}%
  \end{center}
  \caption{ Invariant mass distributions of electron pairs 1-3
and 2-3 (top left and right) and of the tri-electron system (bottom right)
for events classified as tri-electrons.}
  \label{m123}
\end{figure} 




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
\clearpage
%\begin{figure}[htb]
%  \begin{center}
%     \epsfig{file=me100z.eps,width=17.5cm}%
%  \end{center}
%  \caption{Displays of the six multi-electron events with $M_{12}>100$~GeV
%in the $R - z$ view. Indicated are the reconstructed tracks and
%the energy depositions in the calorimeters. The HERA electrons and protons
%enter the detector from the left and right, respectively.}
%  \label{meevdall}
%\end{figure} 



\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.evd1.eps,width=11cm}%
  \end{center}
  \caption{Event Display  
%Run 83507   Event   16817 
of {\bf Event 1}, classified as 2e. Indicated are the reconstructed
tracks and the energy depositions in the calorimeters. The HERA electrons
and protons enter the detector from the left and right, respectively.
The third electron candidate visible in this event has an energy lower 
than the threshold of 5~GeV applied in the electron identification.
}
  \label{meevd1}
\end{figure} 

\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.evd2.eps,width=11cm}%
  \end{center}
  \caption{Event Display
%Run 89256 Event 224212 
of {\bf Event 2}, classified as 2e. 
}
  \label{meevd2}
\end{figure} 

\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.evd5.eps,width=11.5cm}%
  \end{center}
  \caption{Event Display
%Run 254959 Event 17892 
of {\bf Event 3}, classified as 2e. A third compact electromagnetic 
energy deposit 
is located at the forward boundary of the liquid argon calorimeter, 
outside the electron identification fiducial volume. 
This event also shows a low energy converted photon in the backward 
direction close to electron e2.}
  \label{meevd3}
\end{figure} 


\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.evd3.eps,width=11.5cm}%
  \end{center}
  \caption{Event Display
%Run 168058 Event 42123 
of {\bf Event 4}, classified as 3e.
}


  \label{meevd4}
\end{figure} 

\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.evd4.eps,width=11.5cm}%
  \end{center}
  \caption{Event Display 
%Run 192864 Event 123614 
of {\bf Event 5}, classified as 3e.}
  \label{meevd5}
\end{figure} 


\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.evd6.eps,width=11.5cm}%
  \end{center}
  \caption{Event Display 
%Run 267312 Event 203075 
of {\bf Event 6}, classified as 3e.    
}
  \label{meevd6}
\end{figure}


%%%%%%%%%%%%%%%%%%%%%%%%%%%% Cross sections

\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.fig6.eps,width=12cm}%
  \end{center}
  \caption{ Cross section measurement of the photon-photon process
as a function of the hadronic transverse momentum, compared to the
Standard Model prediction. }
  \label{sigma_ptx}
\end{figure} 


\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.fig7.eps,width=12cm}%
  \end{center}
  \caption{ Cross section measurement of the photon-photon process
as a function of the highest electron transverse momentum, compared to the
Standard Model prediction.}
  \label{sigma_pte}
\end{figure} 


\begin{figure}[htb]
  \begin{center}
     \epsfig{file=H1prelim-02-052.fig8.eps,width=12cm}%
  \end{center}
  \caption{ Cross section measurement of the photon-photon process
as a function of the invariant mass of the highest-$P_T$ electron pair, 
compared to the Standard Model prediction.}
  \label{sigma_m12}
\end{figure} 
\clearpage


\end{document}

