%================================================================
% LaTeX file with prefered layout for H1 paper drafts
% use: dvips -D600 file-name
%================================================================
%%%%%%%%%%%%% LATEX HEADER
%%%%%%%%%%%%% DO NOT DELET NOT CHANGE ANYTHING IN THE HEADER
%%%%%%%%%%%%% UNLESS CLEARLY RECOMMENDED
\documentclass[12pt]{article}
\usepackage{epsfig}
\usepackage{amsmath}
\usepackage{hhline}
\usepackage{amssymb}
\usepackage{times}
\usepackage{cite}
%%%%%%%%%%%% Comment the next two lines to remove the line numbering
%\usepackage[]{lineno}
%\linenumbers
%%%%%%%%%%%%
\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
%%%%%%%%%%%%%%%% END OF LATEX HEADER
%===============================title page=============================
\begin{document}  
%%%%%%%%%%%%%%%% Pre-defined commands, you can use for the most obvious notations
\newcommand{\pom}{{I\!\!P}}
\newcommand{\reg}{{I\!\!R}}
\newcommand{\slowpi}{\pi_{\mathit{slow}}}
%\newcommand{\gevsq}{\mathrm{GeV}^2}
\newcommand{\fiidiii}{F_2^{D(3)}}
\newcommand{\fiidiiiarg}{\fiidiii\,(\beta,\,Q^2,\,x)}
\newcommand{\n}{1.19\pm 0.06 (stat.) \pm0.07 (syst.)}
\newcommand{\nz}{1.30\pm 0.08 (stat.)^{+0.08}_{-0.14} (syst.)}
\newcommand{\fiidiiiful}{F_2^{D(4)}\,(\beta,\,Q^2,\,x,\,t)}
\newcommand{\fiipom}{\tilde F_2^D}
\newcommand{\ALPHA}{1.10\pm0.03 (stat.) \pm0.04 (syst.)}
\newcommand{\ALPHAZ}{1.15\pm0.04 (stat.)^{+0.04}_{-0.07} (syst.)}
\newcommand{\fiipomarg}{\fiipom\,(\beta,\,Q^2)}
\newcommand{\pomflux}{f_{\pom / p}}
\newcommand{\nxpom}{1.19\pm 0.06 (stat.) \pm0.07 (syst.)}
\newcommand {\gapprox}
   {\raisebox{-0.7ex}{$\stackrel {\textstyle>}{\sim}$}}
\newcommand {\lapprox}
   {\raisebox{-0.7ex}{$\stackrel {\textstyle<}{\sim}$}}
\def\gsim{\,\lower.25ex\hbox{$\scriptstyle\sim$}\kern-1.30ex%
\raise 0.55ex\hbox{$\scriptstyle >$}\,}
\def\lsim{\,\lower.25ex\hbox{$\scriptstyle\sim$}\kern-1.30ex%
\raise 0.55ex\hbox{$\scriptstyle <$}\,}
\newcommand{\pomfluxarg}{f_{\pom / p}\,(x_\pom)}
\newcommand{\dsf}{\mbox{$F_2^{D(3)}$}}
\newcommand{\dsfva}{\mbox{$F_2^{D(3)}(\beta,Q^2,x_{I\!\!P})$}}
\newcommand{\dsfvb}{\mbox{$F_2^{D(3)}(\beta,Q^2,x)$}}
\newcommand{\dsfpom}{$F_2^{I\!\!P}$}
\newcommand{\gap}{\stackrel{>}{\sim}}
\newcommand{\lap}{\stackrel{<}{\sim}}
\newcommand{\fem}{$F_2^{em}$}
\newcommand{\tsnmp}{$\tilde{\sigma}_{NC}(e^{\mp})$}
\newcommand{\tsnm}{$\tilde{\sigma}_{NC}(e^-)$}
\newcommand{\tsnp}{$\tilde{\sigma}_{NC}(e^+)$}
\newcommand{\st}{$\star$}
\newcommand{\sst}{$\star \star$}
\newcommand{\ssst}{$\star \star \star$}
\newcommand{\sssst}{$\star \star \star \star$}
\newcommand{\tw}{\theta_W}
\newcommand{\sw}{\sin{\theta_W}}
\newcommand{\cw}{\cos{\theta_W}}
\newcommand{\sww}{\sin^2{\theta_W}}
\newcommand{\cww}{\cos^2{\theta_W}}
\newcommand{\trm}{m_{\perp}}
\newcommand{\trp}{p_{\perp}}
\newcommand{\trmm}{m_{\perp}^2}
\newcommand{\trpp}{p_{\perp}^2}
\newcommand{\alp}{\alpha_s}

\newcommand{\alps}{\alpha_s}
\newcommand{\sqrts}{$\sqrt{s}$}
\newcommand{\LO}{$O(\alpha_s^0)$}
\newcommand{\Oa}{$O(\alpha_s)$}
\newcommand{\Oaa}{$O(\alpha_s^2)$}
\newcommand{\PT}{p_{\perp}}
\newcommand{\JPSI}{J/\psi}
\newcommand{\sh}{\hat{s}}
%\newcommand{\th}{\hat{t}}
\newcommand{\uh}{\hat{u}}
\newcommand{\MP}{m_{J/\psi}}
%\newcommand{\PO}{\mbox{l}\!\mbox{P}}
\newcommand{\PO}{I\!\!P}
\newcommand{\xbj}{x}
\newcommand{\xpom}{x_{\PO}}
\newcommand{\ttbs}{\char'134}
\newcommand{\xpomlo}{3\times10^{-4}}  
\newcommand{\xpomup}{0.05}  
\newcommand{\dgr}{^\circ}
\newcommand{\pbarnt}{\,\mbox{{\rm pb$^{-1}$}}}
\newcommand{\gev}{\,\mbox{GeV}}
\newcommand{\WBoson}{\mbox{$W$}}
\newcommand{\fbarn}{\,\mbox{{\rm fb}}}
\newcommand{\fbarnt}{\,\mbox{{\rm fb$^{-1}$}}}
\newcommand{\dsdx}[1]{$d\sigma\!/\!d #1\,$}
\newcommand{\eV}{\mbox{e\hspace{-0.08em}V}}
%
% Some useful tex commands
%
\newcommand{\qsq}{\ensuremath{Q^2} }
\newcommand{\gevsq}{\ensuremath{\mathrm{GeV}^2} }
\newcommand{\et}{\ensuremath{E_t^*} }
\newcommand{\rap}{\ensuremath{\eta^*} }
\newcommand{\gp}{\ensuremath{\gamma^*}p }
\newcommand{\dsiget}{\ensuremath{{\rm d}\sigma_{ep}/{\rm d}E_t^*} }
\newcommand{\dsigrap}{\ensuremath{{\rm d}\sigma_{ep}/{\rm d}\eta^*} }

%%% Dstar stuff
\newcommand{\dstar}{\ensuremath{D^*}}
\newcommand{\dstarp}{\ensuremath{D^{*+}}}
\newcommand{\dstarm}{\ensuremath{D^{*-}}}
\newcommand{\dstarpm}{\ensuremath{D^{*\pm}}}
\newcommand{\zDs}{\ensuremath{z(\dstar )}}
\newcommand{\Wgp}{\ensuremath{W_{\gamma p}}}
\newcommand{\ptds}{\ensuremath{p_t(\dstar )}}
\newcommand{\etads}{\ensuremath{\eta(\dstar )}}
\newcommand{\ptj}{\ensuremath{p_t(\mbox{jet})}}
\newcommand{\ptjn}[1]{\ensuremath{p_t(\mbox{jet$_{#1}$})}}
\newcommand{\etaj}{\ensuremath{\eta(\mbox{jet})}}
\newcommand{\detadsj}{\ensuremath{\eta(\dstar )\, \mbox{-}\, \etaj}}

% Journal macro
\def\Journal#1#2#3#4{{#1} {\bf #2} (#3) #4}
\def\NCA{\em Nuovo Cimento}
\def\NIM{\em Nucl. Instrum. Methods}
\def\NIMA{{\em Nucl. Instrum. Methods} {\bf A}}
\def\NPB{{\em Nucl. Phys.}   {\bf B}}
\def\PLB{{\em Phys. Lett.}   {\bf B}}
\def\PRL{\em Phys. Rev. Lett.}
\def\PRD{{\em Phys. Rev.}    {\bf D}}
\def\ZPC{{\em Z. Phys.}      {\bf C}}
\def\EJC{{\em Eur. Phys. J.} {\bf C}}
\def\CPC{\em Comp. Phys. Commun.}


%%%%%%%%%%%%% my new commands %%%%%%%%%%%%%%%%
\newcommand{\be}{\begin{equation}} 
\newcommand{\ee}{\end{equation}} 
\newcommand{\ba}{\begin{eqnarray}} 
\newcommand{\ea}{\end{eqnarray}} 


\begin{titlepage}

\noindent
\begin{flushleft}
{\tt DESY 08-065    \hfill    ISSN 0418-9833} \\
{\tt June 2008}                  \\
\end{flushleft}


\vspace{2cm}
\begin{center}
\begin{Large}

{\bf Multi--Lepton Production at High Transverse Momenta in \begin{boldmath}$ep$\end{boldmath} Collisions at HERA \\}

\vspace{2cm}

H1 Collaboration

\end{Large}
\end{center}

\vspace{2cm}

\begin{abstract}

Processes leading to a final state with at least two high transverse momentum leptons (electrons or muons) are studied using the full $e^\pm p$ data sample collected by the H1 experiment at HERA.
The data correspond to an integrated luminosity of $463$~pb$^{-1}$.
%
Di-lepton and tri-lepton event classes are investigated.
%
Cross sections of the production of $e^+e^-$ and $\mu^+\mu^-$ pairs are derived in a restricted phase space dominated by photon-photon collisions.
%
In general, good agreement is found with Standard Model predictions.
%
Events are observed with a total scalar sum of lepton transverse momenta above $100$~GeV where the Standard Model expectation is low.
In this region, combining di-lepton and tri-lepton classes, five events are observed in $e^+p$ collisions, compared to a Standard Model expectation of \mbox{$0.96 \pm 0.12$}, while no such event is observed in $e^-p$ data for $0.64 \pm 0.09$ expected.



\end{abstract}


\vspace{1cm}

\begin{center}
Dedicated to the memory of our dear friend and colleague, Beate Naroska.
\end{center}

\vspace{0.5cm}

\begin{center}
Submitted to \PLB
\end{center}

\end{titlepage}

%          THE PAPER DRAFTS HAVE NO AUTHORLIST
%
%          FOR PAPER ISSUED FOR THE FINAL READING 
%          COPY THE AUTHOR AND INSTITUTE LISTS 
%          INTO YOUR AREA
%
% from /h1/iww/ipublications/h1auts.tex 
%          AND UNCOMMENT THE NEXT THREE LINES 
%
\begin{flushleft}
%  \input{h1auts}
%-- H1AUTS Author list by names 
%-- Status: Thu May  8 12:48:44 CEST 2008  Number of authors = 273 

F.D.~Aaron$^{5,49}$,           %BUCH-PD        11/06           Aaron               
C.~Alexa$^{5}$,                %BUCH-PD        06/06           Alexa               
V.~Andreev$^{25}$,             %LPI -PD        8/88            Andreev             
B.~Antunovic$^{11}$,           %DESY-PD        05/07           Antunovic           
S.~Aplin$^{11}$,               %DESY-LEFT      01/08           Aplin               
A.~Asmone$^{33}$,              %ROME-ST        07/2            Asmone              
A.~Astvatsatourov$^{4}$,       %BRUX-LEFT      01/08           Astvatsatourov      
A.~Bacchetta$^{11}$,           %DESY-LEFT      02/08           Bacchetta           
S.~Backovic$^{30}$,            %PODG-PD        03/2            Backovic            
A.~Baghdasaryan$^{38}$,        %YERE-PD        09/03           Baghdasaryana       
P.~Baranov$^{25, \dagger}$,    %LPI -LEFT      05/07           Baranovp            
E.~Barrelet$^{29}$,            %PARI-PD        11/99           Barrelet            
W.~Bartel$^{11}$,              %DESY-PD        8/88            Bartel              
M.~Beckingham$^{11}$,          %DESY-LEFT      08/07           Beckingham          
K.~Begzsuren$^{35}$,           %ULBA-PD        04/06           Begzsuren           
O.~Behnke$^{14}$,              %HDB1-LEFT      12/07           Behnke              
A.~Belousov$^{25}$,            %LPI -PD        8/88            Belousov            
N.~Berger$^{40}$,              %ZUTH-LEFT      10/07           Bergern             
J.C.~Bizot$^{27}$,             %ORSA-PD        8/88            Bizot               
M.-O.~Boenig$^{8}$,            %DORT-LEFT      10/07           Boenig              
V.~Boudry$^{28}$,              %ECPL-PD        1/93            Boudry              
I.~Bozovic-Jelisavcic$^{2}$,   %BEOG-PD        03/06           Bozovicjelisavcic   
J.~Bracinik$^{3}$,             %BIRM-PD        01/2            Bracinik            
G.~Brandt$^{11}$,              %DESY-PD        01/20           Brandt              
M.~Brinkmann$^{11}$,           %DESY-ST        02/06           Brinkmann           
V.~Brisson$^{27}$,             %ORSA-PD        8/88            Brisson             
D.~Bruncko$^{16}$,             %KOSI-PD        8/88            Bruncko             
A.~Bunyatyan$^{13,38}$,        %MPIH-PD        12/95           Bunyatyan           
G.~Buschhorn$^{26}$,           %MPIM-PD        8/88            Buschhorn           
L.~Bystritskaya$^{24}$,        %ITEP-PD        05/99           Bystritskaya        
A.J.~Campbell$^{11}$,          %DESY-PD        8/88            Campbella           
K.B. ~Cantun~Avila$^{22}$,     %MEX1-ST        04/06           Cantunavila         
F.~Cassol-Brunner$^{21}$,      %MARS-PD        12/0            Cassolbrunner       
K.~Cerny$^{32}$,               %PRG2-ST        09/02           Cernyk              
V.~Cerny$^{16,47}$,            %KOSI-PD        06/04           Cernyv              
V.~Chekelian$^{26}$,           %MPIM-PD        01/90           Chekelian           
A.~Cholewa$^{11}$,             %DESY-ST        11/05           Cholewa             
J.G.~Contreras$^{22}$,         %MEX1-PD        04/97           Contreras           
J.A.~Coughlan$^{6}$,           %RAL -PD        8/88            Coughlan            
G.~Cozzika$^{10}$,             %SACL-PD        10/07           Cozzika             
J.~Cvach$^{31}$,               %PRAG-PD        8/88            Cvach               
J.B.~Dainton$^{18}$,           %LIVE-PD        8/88            Dainton             
K.~Daum$^{37,43}$,             %WUPP-PD        06/96           Daum                
M.~De\'{a}k$^{11}$,            %DESY-ST        08/06           Deak                
Y.~de~Boer$^{11}$,             %DESY-PD        01/08           Deboer              
B.~Delcourt$^{27}$,            %ORSA-PD        8/88            Delcourt            
M.~Del~Degan$^{40}$,           %ZUTH-ST        02/05           Deldegan            
J.~Delvax$^{4}$,               %BRUX-ST        10/06           Delvax              
A.~De~Roeck$^{11,45}$,         %DESY-PD        08/88           Deroeck             
E.A.~De~Wolf$^{4}$,            %ANTW-PD        3/93            Dewolf              
C.~Diaconu$^{21}$,             %MARS-PD        01/05           Diaconu             
V.~Dodonov$^{13}$,             %MPIH-PD        04/98           Dodonov             
A.~Dossanov$^{26}$,            %MPIM-ST        01/07           Dossanov            
A.~Dubak$^{30,46}$,            %PODG-PD        10/03           Dubak               
G.~Eckerlin$^{11}$,            %DESY-PD        8/88            Eckerlin            
V.~Efremenko$^{24}$,           %ITEP-PD        8/88            Efremenko           
S.~Egli$^{36}$,                %PSI -PD        01/01           Egli                
A.~Eliseev$^{25}$,             %LPI -PD        01/06           Eliseev             
E.~Elsen$^{11}$,               %DESY-PD        8/88            Elsen               
S.~Essenov$^{24}$,             %ITEP-PD        09/03           Essenov             
A.~Falkiewicz$^{7}$,           %CRAC-ST        07/04           Falkiewicz          
P.J.W.~Faulkner$^{3}$,         %BIRM-LEFT      03/08           Faulkner            
L.~Favart$^{4}$,               %BRUX-PD        8/88            Favart              
A.~Fedotov$^{24}$,             %ITEP-PD        8/88            Fedotov             
R.~Felst$^{11}$,               %DESY-PD        11/0            Felst               
J.~Feltesse$^{10,48}$,         %SACL-PD        03/05           Feltesse            
J.~Ferencei$^{16}$,            %KOSI-PD        01/05           Ferencei            
M.~Fleischer$^{11}$,           %DESY-PD        07/0            Fleischer           
A.~Fomenko$^{25}$,             %LPI -PD        8/88            Fomenko             
E.~Gabathuler$^{18}$,          %LIVE-PD        10/89           Gabathulere         
J.~Gayler$^{11}$,              %DESY-PD        8/88            Gayler              
S.~Ghazaryan$^{38}$,           %YERE-PD        8/88            Ghazaryan           
A.~Glazov$^{11}$,              %DESY-PD        01/04           Glazov              
I.~Glushkov$^{39}$,            %ZEUT-ST        11/03           Glushkov            
L.~Goerlich$^{7}$,             %CRAC-PD        8/88            Goerlich            
M.~Goettlich$^{12}$,           %HAM2-LEFT      11/07           Goettlich           
N.~Gogitidze$^{25}$,           %LPI -PD        8/88            Gogitidze           
M.~Gouzevitch$^{28}$,          %ECPL-ST        10/05           Gouzevitch          
C.~Grab$^{40}$,                %ZUTH-PD        8/88            Grab                
T.~Greenshaw$^{18}$,           %LIVE-PD        8/88            Greenshaw           
B.R.~Grell$^{11}$,             %DESY-ST        09/04           Grell               
G.~Grindhammer$^{26}$,         %MPIM-PD        8/88            Grindhammer         
S.~Habib$^{12,50}$,            %HAM2-ST        12/05           Habib               
D.~Haidt$^{11}$,               %DESY-PD        8/88            Haidt               
M.~Hansson$^{20}$,             %LUND-ST        04/04           Hansson             
C.~Helebrant$^{11}$,           %DFLC-ST        03/06           Helebrant           
R.C.W.~Henderson$^{17}$,       %LANC-PD        8/88            Henderson           
E.~Hennekemper$^{15}$,         %HDB2-ST        11/07           Hennekemper         
H.~Henschel$^{39}$,            %ZEUT-PD        06/99           Henschel            
G.~Herrera$^{23}$,             %MEX2-PD        07/98           Herrera             
M.~Hildebrandt$^{36}$,         %PSI -PD        10/99           Hildebrandtm        
K.H.~Hiller$^{39}$,            %ZEUT-PD        8/88            Hiller              
D.~Hoffmann$^{21}$,            %MARS-PD        10/0            Hoffmann            
R.~Horisberger$^{36}$,         %PSI -PD        8/88            Horisberger         
A.~Hovhannisyan$^{38}$,        %YERE-LEFT      09/07           Hovhannisyan        
T.~Hreus$^{4,44}$,             %BRUX-ST        10/04           Hreus               
M.~Jacquet$^{27}$,             %ORSA-PD        09/96           Jacquet             
M.E.~Janssen$^{11}$,           %DFLC-ST        06/06           Janssenm            
X.~Janssen$^{4}$,              %BRUX-PD        02/03           Janssenx            
V.~Jemanov$^{12}$,             %HAM2-LEFT      03/08           Jemanov             
L.~J\"onsson$^{20}$,           %LUND-PD        8/88            Joensson            
D.P.~Johnson$^{4, \dagger}$,   %BRUX-LEFT      05/07           Johnsond            
A.W.~Jung$^{15}$,              %HDB2-ST        11/04           Junga               
H.~Jung$^{11}$,                %DESY-PD        07/00           Jungh               
M.~Kapichine$^{9}$,            %JINR-PD        3/97            Kapichine           
J.~Katzy$^{11}$,               %DESY-PD        09/1            Katzy               
I.R.~Kenyon$^{3}$,             %BIRM-PD        8/88            Kenyon              
C.~Kiesling$^{26}$,            %MPIM-PD        8/88            Kiesling            
M.~Klein$^{18}$,               %LIVE-PD        8/88            Klein               
C.~Kleinwort$^{11}$,           %DESY-PD        8/88            Kleinwort           
T.~Klimkovich$^{}$,            %DFLC-PD        06/06           Klimkovich          
T.~Kluge$^{18}$,               %LIVE-PD        05/04           Kluge               
A.~Knutsson$^{11}$,            %DESY-PD        04/07           Knutsson            
R.~Kogler$^{26}$,              %MPIM-ST        01/07           Kogler              
V.~Korbel$^{11}$,              %DESY-LEFT      03/08           Korbel              
P.~Kostka$^{39}$,              %ZEUT-PD        8/88            Kostka              
M.~Kraemer$^{11}$,             %DESY-ST        02/06           Kraemer             
K.~Krastev$^{11}$,             %DESY-ST        02/05           Krastev             
J.~Kretzschmar$^{18}$,         %LIVE-PD        01/08           Kretzschmar         
A.~Kropivnitskaya$^{24}$,      %ITEP-ST        07/2            Kropivnitskaya      
K.~Kr\"uger$^{15}$,            %HDB2-PD        01/04           Kruegerk            
K.~Kutak$^{11}$,               %DESY-PD        01/07           Kutak               
M.P.J.~Landon$^{19}$,          %QMWC-PD        8/88            Landon              
W.~Lange$^{39}$,               %ZEUT-PD        8/88            Lange               
G.~La\v{s}tovi\v{c}ka-Medin$^{30}$, %PODG-PD        06/04           Lastovickamedin     
P.~Laycock$^{18}$,             %LIVE-PD        11/03           Laycock             
A.~Lebedev$^{25}$,             %LPI -PD        8/88            Lebedev             
G.~Leibenguth$^{40}$,          %ZUTH-PD        11/04           Leibenguth          
V.~Lendermann$^{15}$,          %HDB2-PD        01/2            Lendermann          
S.~Levonian$^{11}$,            %DESY-PD        8/88            Levonian            
G.~Li$^{27}$,                  %ORSA-PD        09/06           Li                  
K.~Lipka$^{12}$,               %HAM2-PD        01/03           Lipka               
A.~Liptaj$^{26}$,              %MPIM-ST        10/04           Liptaj              
B.~List$^{12}$,                %HAM2-PD        11/99           Listb               
J.~List$^{11}$,                %DFLC-PD        01/05           Listj               
N.~Loktionova$^{25}$,          %LPI -PD        03/99           Loktionova          
R.~Lopez-Fernandez$^{23}$,     %MEX2-PD        03/2            Lopezfernandez      
V.~Lubimov$^{24}$,             %ITEP-PD        01/95           Lubimov             
A.-I.~Lucaci-Timoce$^{11}$,    %DESY-LEFT      08/07           Lucacitimoce        
L.~Lytkin$^{13}$,              %MPIH-PD        8/88            Lytkine             
A.~Makankine$^{9}$,            %JINR-PD        11/02           Makankine           
E.~Malinovski$^{25}$,          %LPI -PD        01/89           Malinovskie         
P.~Marage$^{4}$,               %BRUX-PD        8/88            Marage              
Ll.~Marti$^{11}$,              %DESY-ST        09/05           Marti               
H.-U.~Martyn$^{1}$,            %AAC1-PD        8/88            Martyn              
S.J.~Maxfield$^{18}$,          %LIVE-PD        8/88            Maxfield            
A.~Mehta$^{18}$,               %LIVE-PD        8/88            Mehta               
K.~Meier$^{15}$,               %HDB2-PD        8/88            Meier               
A.B.~Meyer$^{11}$,             %DESY-PD        01/00           Meyeran             
H.~Meyer$^{11}$,               %DFLC-ST        06/06           Meyerhe             
H.~Meyer$^{37}$,               %WUPP-PD        8/88            Meyerhi             
J.~Meyer$^{11}$,               %DESY-PD        8/88            Meyerj              
V.~Michels$^{11}$,             %DESY-ST        03/05           Michels             
S.~Mikocki$^{7}$,              %CRAC-PD        8/88            Mikocki             
I.~Milcewicz-Mika$^{7}$,       %CRAC-ST        10/02           Milcewicz           
F.~Moreau$^{28}$,              %ECPL-PD        01/90           Moreau              
A.~Morozov$^{9}$,              %JINR-PD        06/99           Morozova            
J.V.~Morris$^{6}$,             %RAL -PD        8/88            Morris              
M.U.~Mozer$^{4}$,              %BRUX-PD        06/07           Mozer               
M.~Mudrinic$^{2}$,             %BEOG-PD        01/07           Mudrinic            
K.~M\"uller$^{41}$,            %ZUER-PD        8/88            Muellerk            
P.~Mur\'\i n$^{16,44}$,        %KOSI-PD        8/88            Murin               
K.~Nankov$^{34}$,              %SOFI-LEFT      10/07           Nankov              
B.~Naroska$^{12, \dagger}$,    %HAM2-PD        8/88            Naroska             
Th.~Naumann$^{39}$,            %ZEUT-PD        01/89           Naumannt            
P.R.~Newman$^{3}$,             %BIRM-PD        10/92           Newman              
C.~Niebuhr$^{11}$,             %DESY-PD        3/93            Niebuhr             
A.~Nikiforov$^{11}$,           %DESY-PD        05/07           Nikiforov           
G.~Nowak$^{7}$,                %CRAC-PD        8/88            Nowakg              
K.~Nowak$^{41}$,               %ZUER-ST        08/05           Nowakk              
M.~Nozicka$^{11}$,             %DESY-PD        11/06           Nozicka             
B.~Olivier$^{26}$,             %MPIM-PD        11/04           Olivier             
J.E.~Olsson$^{11}$,            %DESY-PD        8/88            Olsson              
S.~Osman$^{20}$,               %LUND-ST        02/04           Osman               
D.~Ozerov$^{24}$,              %ITEP-ST        08/98           Ozerov              
V.~Palichik$^{9}$,             %JINR-PD        01/04           Palichik            
I.~Panagoulias$^{l,}$$^{11,42}$, %DESY-ST        08/04           Panagoulias         
M.~Pandurovic$^{2}$,           %BEOG-ST        03/06           Pandurovic          
Th.~Papadopoulou$^{l,}$$^{11,42}$, %DESY-PD        06/04           Papadopoulou        
C.~Pascaud$^{27}$,             %ORSA-PD        8/88            Pascaud             
G.D.~Patel$^{18}$,             %LIVE-PD        8/88            Patel               
O.~Pejchal$^{32}$,             %PRG2-ST        12/06           Pejchal             
H.~Peng$^{11}$,                %DESY-LEFT      08/07           Peng                
E.~Perez$^{10,45}$,            %SACL-PD        10/07           Perez               
A.~Petrukhin$^{24}$,           %ITEP-ST        01/01           Petrukhin           
I.~Picuric$^{30}$,             %PODG-PD        01/06           Picuric             
S.~Piec$^{39}$,                %ZEUT-ST        01/06           Piec                
D.~Pitzl$^{11}$,               %DESY-PD        8/88            Pitzl               
R.~Pla\v{c}akyt\.{e}$^{11}$,   %DESY-PD        10/06           Placakyte           
R.~Polifka$^{32}$,             %PRG2-ST        10/06           Polifka             
B.~Povh$^{13}$,                %MPIH-PD        8/88            Povh                
T.~Preda$^{5}$,                %BUCH-PD        06/06           Preda               
V.~Radescu$^{11}$,             %DESY-PD        10/06           Radescu             
A.J.~Rahmat$^{18}$,            %LIVE-ST        01/05           Rahmat              
N.~Raicevic$^{30}$,            %PODG-PD        03/2            Raicevic            
A.~Raspiareza$^{26}$,          %MPIM-PD        12/06           Raspiareza          
T.~Ravdandorj$^{35}$,          %ULBA-PD        06/06           Ravdandorj          
P.~Reimer$^{31}$,              %PRAG-PD        8/88            Reimer              
E.~Rizvi$^{19}$,               %QMWC-PD        01/05           Rizvi               
P.~Robmann$^{41}$,             %ZUER-PD        8/88            Robmann             
B.~Roland$^{4}$,               %BRUX-ST        12/02           Roland              
R.~Roosen$^{4}$,               %BRUX-PD        8/88            Roosen              
A.~Rostovtsev$^{24}$,          %ITEP-PD        8/88            Rostovtsev          
M.~Rotaru$^{5}$,               %BUCH-ST        02/07           Rotaru              
J.E.~Ruiz~Tabasco$^{22}$,      %MEX1-ST        09/06           Ruiztabascojuliaelis
Z.~Rurikova$^{11}$,            %DESY-PD        05/06           Rurikova            
S.~Rusakov$^{25}$,             %LPI -PD        8/88            Rusakov             
D.~Salek$^{32}$,               %PRG2-ST        11/06           Salek               
F.~Salvaire$^{11}$,            %DESY-LEFT      09/07           Salvaire            
D.P.C.~Sankey$^{6}$,           %RAL -PD        8/88            Sankey              
M.~Sauter$^{40}$,              %ZUTH-ST        10/05           Sauter              
E.~Sauvan$^{21}$,              %MARS-PD        11/1            Sauvan              
S.~Schmidt$^{11}$,             %DFLC-LEFT      09/07           Schmidts            
S.~Schmitt$^{11}$,             %DESY-PD        09/07           Schmittst           
C.~Schmitz$^{41}$,             %ZUER-LEFT      04/08           Schmitz             
L.~Schoeffel$^{10}$,           %SACL-PD        12/98           Schoeffel           
A.~Sch\"oning$^{11,41}$,       %ZUER-PD        02/99           Schoening           
H.-C.~Schultz-Coulon$^{15}$,   %HDB2-PD        01/04           Schultzcoulon       
F.~Sefkow$^{11}$,              %DFLC-PD        09/99           Sefkow              
R.N.~Shaw-West$^{3}$,          %BIRM-ST        10/04           Shawwest            
I.~Sheviakov$^{25}$,           %LPI -LEFT      03/08           Sheviakov           
L.N.~Shtarkov$^{25}$,          %LPI -PD        8/88            Shtarkov            
S.~Shushkevich$^{26}$,         %MPIM-ST        08/07           Shushkevich         
T.~Sloan$^{17}$,               %LANC-PD        1/96            Sloan               
I.~Smiljanic$^{2}$,            %BEOG-PD        03/06           Smiljanic           
P.~Smirnov$^{25}$,             %LPI -LEFT      08/07           Smirnov             
Y.~Soloviev$^{25}$,            %LPI -PD        8/88            Soloviev            
P.~Sopicki$^{7}$,              %CRAC-ST        09/07           Sopicki             
D.~South$^{8}$,                %DORT-PD        06/03           South               
V.~Spaskov$^{9}$,              %JINR-PD        12/97           Spaskov             
A.~Specka$^{28}$,              %ECPL-PD        3/95            Specka              
Z.~Staykova$^{11}$,            %DESY-ST        08/06           Staykova            
M.~Steder$^{11}$,              %DESY-ST        05/05           Steder              
B.~Stella$^{33}$,              %ROME-PD        8/88            Stella              
U.~Straumann$^{41}$,           %ZUER-PD        8/88            Straumann           
D.~Sunar$^{4}$,                %ANTW-ST        03/05           Sunar               
T.~Sykora$^{4}$,               %ANTW-PD        01/06           Sykora              
V.~Tchoulakov$^{9}$,           %JINR-PD        05/03           Tchoulakov          
G.~Thompson$^{19}$,            %QMWC-PD        8/88            Thompsong           
P.D.~Thompson$^{3}$,           %BIRM-PD        08/99           Thompsonp           
T.~Toll$^{11}$,                %DESY-ST        07/05           Toll                
F.~Tomasz$^{16}$,              %KOSI-PD        07/05           Tomasz              
T.H.~Tran$^{27}$,              %ORSA-ST        10/06           Tran                
D.~Traynor$^{19}$,             %QMWC-PD        12/01           Traynor             
T.N.~Trinh$^{21}$,             %MARS-ST        11/05           Trinh               
P.~Tru\"ol$^{41}$,             %ZUER-PD        8/88            Truoel              
I.~Tsakov$^{34}$,              %SOFI-PD        04/03           Tsakov              
B.~Tseepeldorj$^{35,51}$,      %ULBA-PD        06/06           Tseepeldorj         
I.~Tsurin$^{39}$,              %ZEUT-LEFT      09/07           Tsurin              
J.~Turnau$^{7}$,               %CRAC-PD        8/88            Turnau              
E.~Tzamariudaki$^{26}$,        %MPIM-LEFT      08/07           Tzamariudaki        
K.~Urban$^{15}$,               %HDB2-ST        04/05           Urbank              
A.~Valk\'arov\'a$^{32}$,       %PRG2-PD        8/88            Valkarova           
C.~Vall\'ee$^{21}$,            %MARS-PD        8/88            Vallee              
P.~Van~Mechelen$^{4}$,         %ANTW-PD        12/98           Vanmechelen         
A.~Vargas Trevino$^{11}$,      %DFLC-PD        02/07           Vargastrevino       
Y.~Vazdik$^{25}$,              %LPI -PD        8/88            Vazdik              
S.~Vinokurova$^{11}$,          %DESY-ST        09/02           Vinokurova          
V.~Volchinski$^{38}$,          %YERE-PD        12/01           Volchinski          
D.~Wegener$^{8}$,              %DORT-PD        8/88            Wegener             
M.~Wessels$^{11}$,             %DESY-LEFT      10/07           Wessels             
Ch.~Wissing$^{11}$,            %DESY-PD        07/06           Wissing             
E.~W\"unsch$^{11}$,            %DESY-PD        8/88            Wuensch             
V.~Yeganov$^{38}$,             %YERE-LEFT      09/07           Yeganov             
J.~\v{Z}\'a\v{c}ek$^{32}$,     %PRG2-PD        8/88            Zacek               
J.~Z\'ale\v{s}\'ak$^{31}$,     %PRAG-PD        01/05           Zalesak             
Z.~Zhang$^{27}$,               %ORSA-PD        10/92           Zhang               
A.~Zhelezov$^{24}$,            %ITEP-LEFT      07/07           Zhelezov            
A.~Zhokin$^{24}$,              %ITEP-PD        04/99           Zhokine             
Y.C.~Zhu$^{11}$,               %DESY-LEFT      10/07           Zhu                 
T.~Zimmermann$^{40}$,          %ZUTH-ST        09/04           Zimmermannt         
H.~Zohrabyan$^{13}$,           %MPIH-PD        11/02           Zohrabyan           
and
F.~Zomer$^{27}$                %ORSA-PD        8/88            Zomer          


%-- H1 Institutes 
\bigskip{\it
 $ ^{1}$ I. Physikalisches Institut der RWTH, Aachen, Germany$^{ a}$ \\
 $ ^{2}$ Vinca  Institute of Nuclear Sciences, Belgrade, Serbia \\
 $ ^{3}$ School of Physics and Astronomy, University of Birmingham,
          Birmingham, UK$^{ b}$ \\
 $ ^{4}$ Inter-University Institute for High Energies ULB-VUB, Brussels;
          Universiteit Antwerpen, Antwerpen; Belgium$^{ c}$ \\
 $ ^{5}$ National Institute for Physics and Nuclear Engineering (NIPNE) ,
          Bucharest, Romania \\
 $ ^{6}$ Rutherford Appleton Laboratory, Chilton, Didcot, UK$^{ b}$ \\
 $ ^{7}$ Institute for Nuclear Physics, Cracow, Poland$^{ d}$ \\
 $ ^{8}$ Institut f\"ur Physik, TU Dortmund, Dortmund, Germany$^{ a}$ \\
 $ ^{9}$ Joint Institute for Nuclear Research, Dubna, Russia \\
 $ ^{10}$ CEA, DSM/DAPNIA, CE-Saclay, Gif-sur-Yvette, France \\
 $ ^{11}$ DESY, Hamburg, Germany \\
 $ ^{12}$ Institut f\"ur Experimentalphysik, Universit\"at Hamburg,
          Hamburg, Germany$^{ a}$ \\
 $ ^{13}$ Max-Planck-Institut f\"ur Kernphysik, Heidelberg, Germany \\
 $ ^{14}$ Physikalisches Institut, Universit\"at Heidelberg,
          Heidelberg, Germany$^{ a}$ \\
 $ ^{15}$ Kirchhoff-Institut f\"ur Physik, Universit\"at Heidelberg,
          Heidelberg, Germany$^{ a}$ \\
 $ ^{16}$ Institute of Experimental Physics, Slovak Academy of
          Sciences, Ko\v{s}ice, Slovak Republic$^{ f}$ \\
 $ ^{17}$ Department of Physics, University of Lancaster,
          Lancaster, UK$^{ b}$ \\
 $ ^{18}$ Department of Physics, University of Liverpool,
          Liverpool, UK$^{ b}$ \\
 $ ^{19}$ Queen Mary and Westfield College, London, UK$^{ b}$ \\
 $ ^{20}$ Physics Department, University of Lund,
          Lund, Sweden$^{ g}$ \\
 $ ^{21}$ CPPM, CNRS/IN2P3 - Univ. Mediterranee,
          Marseille - France \\
 $ ^{22}$ Departamento de Fisica Aplicada,
          CINVESTAV, M\'erida, Yucat\'an, M\'exico$^{ j}$ \\
 $ ^{23}$ Departamento de Fisica, CINVESTAV, M\'exico$^{ j}$ \\
 $ ^{24}$ Institute for Theoretical and Experimental Physics,
          Moscow, Russia \\
 $ ^{25}$ Lebedev Physical Institute, Moscow, Russia$^{ e}$ \\
 $ ^{26}$ Max-Planck-Institut f\"ur Physik, M\"unchen, Germany \\
 $ ^{27}$ LAL, Univ Paris-Sud, CNRS/IN2P3, Orsay, France \\
 $ ^{28}$ LLR, Ecole Polytechnique, IN2P3-CNRS, Palaiseau, France \\
 $ ^{29}$ LPNHE, Universit\'{e}s Paris VI and VII, IN2P3-CNRS,
          Paris, France \\
 $ ^{30}$ Faculty of Science, University of Montenegro,
          Podgorica, Montenegro$^{ e}$ \\
 $ ^{31}$ Institute of Physics, Academy of Sciences of the Czech Republic,
          Praha, Czech Republic$^{ h}$ \\
 $ ^{32}$ Faculty of Mathematics and Physics, Charles University,
          Praha, Czech Republic$^{ h}$ \\
 $ ^{33}$ Dipartimento di Fisica Universit\`a di Roma Tre
          and INFN Roma~3, Roma, Italy \\
 $ ^{34}$ Institute for Nuclear Research and Nuclear Energy,
          Sofia, Bulgaria$^{ e}$ \\
 $ ^{35}$ Institute of Physics and Technology of the Mongolian
          Academy of Sciences , Ulaanbaatar, Mongolia \\
 $ ^{36}$ Paul Scherrer Institut,
          Villigen, Switzerland \\
 $ ^{37}$ Fachbereich C, Universit\"at Wuppertal,
          Wuppertal, Germany \\
 $ ^{38}$ Yerevan Physics Institute, Yerevan, Armenia \\
 $ ^{39}$ DESY, Zeuthen, Germany \\
 $ ^{40}$ Institut f\"ur Teilchenphysik, ETH, Z\"urich, Switzerland$^{ i}$ \\
 $ ^{41}$ Physik-Institut der Universit\"at Z\"urich, Z\"urich, Switzerland$^{ i}$ \\

\bigskip
 $ ^{42}$ Also at Physics Department, National Technical University,
          Zografou Campus, GR-15773 Athens, Greece \\
 $ ^{43}$ Also at Rechenzentrum, Universit\"at Wuppertal,
          Wuppertal, Germany \\
 $ ^{44}$ Also at University of P.J. \v{S}af\'{a}rik,
          Ko\v{s}ice, Slovak Republic \\
 $ ^{45}$ Also at CERN, Geneva, Switzerland \\
 $ ^{46}$ Also at Max-Planck-Institut f\"ur Physik, M\"unchen, Germany \\
 $ ^{47}$ Also at Comenius University, Bratislava, Slovak Republic \\
 $ ^{48}$ Also at DESY and University Hamburg,
          Helmholtz Humboldt Research Award \\
 $ ^{49}$ Also at Faculty of Physics, University of Bucharest,
          Bucharest, Romania \\
 $ ^{50}$ Supported by a scholarship of the World
          Laboratory Bj\"orn Wiik Research
Project \\
 $ ^{51}$ Also at Ulaanbaatar University, Ulaanbaatar, Mongolia \\

\smallskip
 $ ^{\dagger}$ Deceased \\


\bigskip
 $ ^a$ Supported by the Bundesministerium f\"ur Bildung und Forschung, FRG,
      under contract numbers 05 H1 1GUA /1, 05 H1 1PAA /1, 05 H1 1PAB /9,
      05 H1 1PEA /6, 05 H1 1VHA /7 and 05 H1 1VHB /5 \\
 $ ^b$ Supported by the UK Science and Technology Facilities Council,
      and formerly by the UK Particle Physics and
      Astronomy Research Council \\
 $ ^c$ Supported by FNRS-FWO-Vlaanderen, IISN-IIKW and IWT
      and  by Interuniversity
Attraction Poles Programme,
      Belgian Science Policy \\
 $ ^d$ Partially Supported by Polish Ministry of Science and Higher
      Education, grant PBS/DESY/70/2006 \\
 $ ^e$ Supported by the Deutsche Forschungsgemeinschaft \\
 $ ^f$ Supported by VEGA SR grant no. 2/7062/ 27 \\
 $ ^g$ Supported by the Swedish Natural Science Research Council \\
 $ ^h$ Supported by the Ministry of Education of the Czech Republic
      under the projects  LC527, INGO-1P05LA259 and
      MSM0021620859 \\
 $ ^i$ Supported by the Swiss National Science Foundation \\
 $ ^j$ Supported by  CONACYT,
      M\'exico, grant 48778-F \\
 $ ^l$ This project is co-funded by the European Social Fund  (75\%) and
      National Resources (25\%) - (EPEAEK II) - PYTHAGORAS II \\
}\end{flushleft}
%
% Please not that the author list may need re-formatting.

\newpage

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Introduction}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


Within the Standard Model (SM) the production of multi-lepton events
in electron\footnote{In this paper the term ``electron'' is used generically to refer to both electrons and positrons, if not otherwise stated.}--proton collisions proceeds mainly via photon--photon interactions, $\gamma \gamma \rightarrow \ell^+ \ell^-$, where photons are radiated from the incident beam particles~\cite{Vermaseren:1982cz}. 
%
The clean experimental signature and the precise SM prediction of this process provide a high sensitivity to searches for new phenomena producing multi-lepton final states.
For example, the production of a doubly charged Higgs boson~\cite{Accomando:1993ar,Aktas:2006nu} or processes involving  generic bosons carrying two units of lepton number (bileptons)~\cite{Cuypers:1996ia}, could lead to multi-lepton events of large invariant mass.
%
Measurements of both multi-electron~\cite{Aktas:2003jg}  and muon pair~\cite{Aktas:2003sz} production at high transverse momentum ($P_T$) have already been performed by the H1 Collaboration using a data sample corresponding to an integrated luminosity of $115$~pb$^{-1}$.
%
Events with high invariant mass $M_{12}$ of the two highest $P_T$ electrons were measured. Three events with two electrons and three events with three electrons were observed in the region $M_{12}>100$~GeV, where the SM prediction is low~\cite{Aktas:2003jg}.

%
In the present paper a measurement of multi-lepton production at high $P_T$ using the complete $e^\pm p$ HERA collider data of the H1 experiment is presented.
%
The data are recorded at electron and proton beam energies of $27.6$~GeV and $820$~GeV or $920$~GeV, corresponding to centre-of-mass energies $\sqrt{s}$ of $301$~GeV or $319$~GeV, respectively.
The total integrated luminosity of the data is $463$~pb$^{-1}$, which represents a factor of four increase with respect to the previous published results.
%
The data comprise $178$~pb$^{-1}$ recorded in $e^-p$ collisions and $285$~pb$^{-1}$ in $e^+p$ collisions, of which $35$~pb$^{-1}$ were recorded at $\sqrt{s} = 301$~GeV.
%
While previous measurements were performed mainly using $e^+p$ data, an $e^-p$ data set with more than a ten-fold increase in integrated luminosity is now analysed.
%
The current analysis is extended with respect to those presented in~\cite{Aktas:2003jg,Aktas:2003sz}, such that all event final states with at least two high $P_T$ leptons, electrons ($e$) or muons  ($\mu$), are now investigated.
%
%
In addition, differential cross sections of the production of $e^+e^-$ and $\mu^+\mu^-$ pairs are measured in a restricted phase space dominated by photon-photon collisions.
%


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Standard Model Processes and their Simulation}\label{sec:MC}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Multi-lepton events are simulated with the GRAPE~\cite{Abe:2000cv}
generator, which includes all electroweak matrix elements at tree level.
The production mechanisms include $\gamma \gamma$, $\gamma Z$, $ZZ$
interactions, internal photon conversion and the decay of virtual or
real $Z$ bosons.
%
Different approaches for the cross section calculation are followed depending on the virtuality $Q^2$ of the photon coupled to the proton and on the mass $W$ of the hadronic final state.
%
In the elastic region, $ep \rightarrow e \: \ell^+ \ell^- p$, the proton vertex is described in terms of dipole form factors.
The quasi-elastic domain is defined by $W < 5$~GeV or $Q^2 < 1$~GeV$^2$.
In the region $W < 2$~GeV, a resonance parametrisation~\cite{Brasse:1976bf}
is used for the proton vertex.
In the remaining quasi-elastic phase space, a fit to photoproduction and deep inelastic scattering data is used~\cite{Abramowicz:1997ms}.
In the inelastic regime, corresponding to electron-quark interactions with $W > 5$~GeV and $Q^2 > 1$~GeV$^2$, proton parton densities taken from the CTEQ5L~\cite{Pumplin:2002vw} parametrisation are used.
%
%
Initial and final state QED radiation is included.  
%
The hadronic final state is simulated via
interfaces to PYTHIA~\cite{Sjostrand:2000wi} and SOPHIA~\cite{Mucke:1999yb} for the inelastic and quasi-elastic regimes, respectively.  
%
GRAPE predicts cross sections for
$ep \rightarrow e \: \mu^+ \mu^- X$ and $ep \rightarrow e \: e^+ e^- X$  processes, leading to $e\mu\mu$ and $eee$ final states.
 Events with only two leptons ($\mu\mu$, $e\mu$ or $ee$) are observed if the scattered electron or one lepton of the pair is not detected. 
%
The $ep \rightarrow e \: \tau^+ \tau^- X$ process with subsequent leptonic tau decays is also simulated with GRAPE and its 
contribution to the studied final states is found to be at most $4$\%.



Experimental backgrounds to $ee$ and $eee$ final states arise from events in which, in addition to a genuine electron, one or two fake electrons are reconstructed.
%
Neutral current (NC) deep-inelastic scattering (DIS) events ($ep \rightarrow e X$) in which hadrons or radiated photons are wrongly identified as electrons constitute the dominant background contribution.
%
QED Compton scattering $ep \rightarrow e \gamma X$ may also contribute if the photon is misidentified as an electron. 
%
Background to the $e\mu$ final state may arise from NC DIS events if hadrons are misidentified as muons.
%
%
The NC DIS and Compton processes are
simulated using the RAPGAP~\cite{Jung:1993gf} and WABGEN~\cite{Berger:kp}
Monte Carlo (MC) generators, respectively.


%
Generated events are passed through the full GEANT~\cite{Brun:1987ma} based simulation of the H1 apparatus, which takes into account the running conditions of the different data taking periods, and are reconstructed and analysed with the same program chain as is used for the data.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Experimental Conditions}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

A detailed description of the H1 experiment can be found in~\cite{Abt:h1}.
Only the detector components relevant to the
present analysis are briefly described here.  
%
The origin of the H1 coordinate system is the nominal $ep$ interaction point, with the direction of the proton beam defining the positive $z$-axis (forward region). Transverse momentum is measured in the $xy$ plane. The pseudorapidity $\eta$ is related to the polar angle $\theta$ by $\eta = -\ln \, \tan (\theta/2)$.
%
The Liquid Argon (LAr) calorimeter~\cite{Andrieu:1993kh} covers the polar angle range
$4^\circ < \theta < 154^\circ$ with full azimuthal acceptance.
Electromagnetic shower energies are measured with a precision of
$\sigma (E)/E = 12\%/ \sqrt{E/\mbox{GeV}} \oplus 1\%$ and hadronic energies
with $\sigma (E)/E = 50\%/\sqrt{E/\mbox{GeV}} \oplus 2\%$, as measured in test beams~\cite{Andrieu:1993tz,Andrieu:1994yn}.
%
In the backward region, energy measurements are provided by a lead/scintillating-fiber (SpaCal) calorimeter~\cite{Appuhn:1996na} covering the range $155^\circ < \theta < 178^\circ$.
%
The central ($20^\circ < \theta < 160^\circ$)  and forward ($7^\circ < \theta < 25^\circ$)  inner tracking detectors are used to
measure charged particle trajectories, to reconstruct the interaction
vertex and to complement the measurement of hadronic energies.
%
The LAr calorimeter and inner tracking detectors are enclosed in a super-conducting magnetic
coil with a field strength of $1.16$~T.
%
From the curvature of charged particle trajectories in the magnetic field, the central tracking system provides transverse momentum measurements with a resolution of $\sigma_{P_T}/P_T = 0.005 P_T / \rm{GeV} \oplus 0.015$~\cite{Kleinwort:2006zz}.
%
The return yoke of the magnetic coil is the outermost part of the detector and is
equipped with streamer tubes forming the central muon detector
($4^\circ < \theta < 171^\circ$).
%
%
%

The luminosity is determined from the rate of the Bethe-Heitler process $ep \rightarrow ep \gamma$,
measured using a photon detector located close to the beam pipe at $z=-103~{\rm m}$, in the backward direction.


Events having an electromagnetic deposit in the LAr with an energy greater than $10$~GeV are detected by the LAr trigger with an efficiency close to $100$\%~\cite{Adloff:2003uh}.
%
The muon trigger is based on single muon signatures from the central muon detector, combined with signals from the central tracking detector. The trigger efficiency is about $90$\% for di-muon events with muon transverse momenta larger than $10$ and $5$~GeV.
%
A combination of the LAr and muon triggers is used if an electron and a muon are both present in the final state, resulting in an efficiency close to $100$\% for electron and muon transverse momenta above $5$~GeV.




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Data Analysis}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


\subsection{Particle Identification}\label{sec:part_def}


Electron candidates are identified in the polar angle range $5^\circ < \theta < 175^\circ$ as a compact and isolated electromagnetic shower in either the LAr or SpaCal calorimeter.  
%
The electron energy and angular direction are measured by the calorimeters.
%
The electron energy threshold is $5$~GeV in the polar angle range $20^\circ < \theta < 175^\circ$ and is raised to $10$~GeV in the forward region ($5^\circ < \theta < 20^\circ$).
%
The calorimetric energy measured within a distance in the pseudorapidity-azimuth $(\eta - \phi)$ plane $R=\sqrt{\Delta \eta^2 + \Delta \phi^2} < 0.75$ around the electron is required to be below $2.5$\% of the electron energy.
%
%
In the region of angular overlap between the LAr and
the central tracking detectors ($20^\circ < \theta < 150^\circ$), hereafter referred to as the central region, the
calorimetric electron identification is complemented by tracking information.  
%
In this region it is required that a well measured track
geometrically matches the centre-of-gravity of the electromagnetic cluster within a distance of
closest approach of $12$~cm.  
%
The electron is required to be isolated from any other well measured track by a distance $R > 0.5$ to the electron direction.
%
%
Furthermore,
the distance from the first measured track point in the central drift
chambers to the beam axis is required to be below $30$~cm in order to
reject photons that convert late in the central tracker material. 
%
In addition, in the central region the transverse momentum of the associated electron track $P_T^{e_{tk}}$ is required to match the
calorimetric measurement $P_T^e$ such that $1/P_T^{e_{tk}} - 1/P_T^e <
0.02$~GeV$^{-1}$ in order to reject hadronic showers.  
%
Due to a lower track reconstruction efficiency and higher showering probability, no track conditions are required for electron candidates in the forward ($5^\circ < \theta < 20^\circ$) and backward ($150^\circ <\theta < 175^\circ$) regions. 
%
In these regions, no distinction between electrons and photons is made.
%
The resulting electron identification efficiency is $80$\% in the central region and larger than $95$\% in the forward and backward regions. 



Muon candidates are identified in the polar angle range $20^\circ <\theta < 160^\circ$ with a minimum transverse momentum of $2$~GeV.
%
Muon identification is based on the measurement of a track in the inner tracking systems associated with a track segment or an energy deposit in the central muon detector~\cite{Aktas:2003sz,Andreev:2003pm}.
%
The muon momentum is measured from the inner track curvature in the solenoidal magnetic field.
%
A muon candidate should have no more than $3.5$~GeV deposited in a
cylinder, centred on the muon track direction, of radius $25$~cm and $50$~cm in the electromagnetic and hadronic sections of the LAr calorimeter, respectively.
%
Misidentified hadrons are strongly suppressed by requiring that the muon candidate be separated from any jet and from any track by $R > 1$.
%
%
The efficiency to identify muons is $\sim 90$\%.


Calorimeter energy deposits and tracks not previously identified as electron or muon candidates are used to form combined cluster-track objects, from which the hadronic energy is reconstructed~\cite{matti,benji}.
%
Jets are reconstructed from these combined cluster-track objects using an inclusive $k_T$ algorithm~\cite{Ellis:1993tq,Catani:1993hr} with a minimum transverse momentum of $2.5$~GeV.



\subsection{Event Selection}


In order to remove background events induced by cosmic rays and other non-$ep$ sources, the event vertex is required to be reconstructed within $35$~cm in $z$ of the nominal interaction point. In addition, to remove non-$ep$ background topological filters and timing vetoes are applied.
%

Multi-lepton events are selected by requiring at least two
central ($20^\circ < \theta < 150^\circ$) electron or muon
candidates, of which one must have $P_T^\ell > 10$~GeV and the other
$P_T^{\ell} > 5$~GeV. 
%
Additional lepton candidates are identified in the
detector according to the criteria defined in section~\ref{sec:part_def}.
%
All lepton candidates are required to be isolated with respect to each other by a minimum distance in pseudorapidity-azimuth of $R > 0.5$.
%
Lepton candidates are
ordered according to decreasing transverse momentum, $P_T^{\ell_i} > P_T^{\ell_{i+1}}$.
%
Final states with all possible combinations of lepton candidates are investigated.
%
Selected events are classified into independent samples according to the flavour and the number of lepton candidates (e.g. $ee$, $e\mu$, $e\mu\mu$).
%
%

In order to measure the lepton pair production cross section in a well
defined region of phase space, sub-samples of $ee$ and $\mu\mu$ events dominated by
photon-photon collisions are selected, labelled $(\gamma\gamma)_{e}$ and $(\gamma\gamma)_{\mu}$, respectively.
%
%
In these subsamples the two leptons are required to be of opposite charge and a significant deficit compared to the initial state must be observed in the difference\footnote{For fully contained events or events where only longitudinal momentum along the proton direction ($+ z$) is undetected, one expects $E-P_z = 2E^0_e = 55.2$~GeV, where $E^0_e$ is the energy of the incident electron. If the scattered electron is undetected, the threshold $E-P_z < 45$~GeV corresponds to a cut on the fractional energy loss $y=(E-P_z)/2E^0_e < 0.82$.} $E -P_z$ of the energy and the longitudinal momentum of all visible particles, $E -P_z <45$~GeV.
%
These two conditions ensure that the incident electron is lost in the beam pipe after radiating a quasi-real photon of squared four-momentum $Q^2$ lower than $1$~GeV$^2$. 


\subsection{Systematic Uncertainties}
\label{sec:syst}

The following experimental systematic uncertainties are considered:

\begin{itemize}
\item The uncertainty on the electromagnetic energy scale varies depending on the polar angle from $0.7$\% in the backward and central region to $2$\% in the forward region. The polar angle measurement  uncertainty of electromagnetic clusters is $3$~mrad. The identification efficiency of electrons is known with an uncertainty of $3$ to $5$\%, depending on the polar angle. 
%
%
\item The scale uncertainty on the transverse momentum of high $P_T$ muons is $2.5$\%. The uncertainty on the reconstruction of the muon polar angle is $3$~mrad. The identification efficiency of muons is known with an uncertainty of $5$\%.
%
%
\item The hadronic energy scale is known within $2$\% at high transverse momentum and $5$\% for events with a total hadronic transverse momentum below $8$~GeV. 
%
%
\item The uncertainty on the trigger efficiency is estimated to be $3$\% if at least one electron candidate is detected, and $6$\% if only muons are present in the final state.
%
\item The luminosity measurement has an uncertainty of $3$\%.
\end{itemize}


The effect of the above systematic uncertainties on the SM expectation is determined by varying the experimental quantities by $\pm 1$ standard deviation in the MC samples and propagating these variations through the whole analysis chain.

Additional model systematic uncertainties are attributed to the SM Monte Carlo generators described in section~$\ref{sec:MC}$.
%
The theoretical uncertainty on the lepton pair production cross section calculated with GRAPE is $3$\%.
%
The uncertainty on the QED Compton and NC DIS background is $20$\%, as deduced from dedicated studies~\cite{Aktas:2003jg}.
%
The total error on the SM prediction is determined by adding the effects of all model and experimental systematic uncertainties in quadrature.




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Results}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\subsection{Multi--Lepton Event Samples}

The observed event yields are summarised in table~\ref{tab:mlepyields}.
Only classes for which at least one data event is selected are shown.
In all other classes, no event is observed and the SM prediction is also negligible.
% 
The observed numbers of events are in good agreement with the SM expectations.
% 
The $eee$ and $ee$ channels are dominated by electron pair production.
The $e\mu\mu$, $\mu\mu$ and $e\mu$ channels contain mainly events from muon pair production.
The $e\mu$ channel is populated if the scattered electron and only one muon of the pair is selected.
%
Four data events are classified as $ee\mu$ compared to a SM expectation of $1.43 \pm 0.26$, dominated by the production of muon pair events where one muon is lost and the second electron candidate is due to a radiated photon. 
%
One event with four electron candidates is observed compared to a prediction of $0.33 \pm 0.07$. %
According to MC simulations, this signature is due to 
tri-electron events with a radiated photon.
%

%
The distributions of the invariant mass $M_{12}$ of the two highest
$P_T$ leptons for the $eee$ and $e\mu\mu$ samples are shown in figures~\ref{fig:Masses}(a) and (b), respectively.
%
The distributions of the invariant mass $M_{12}$  of the two leptons in the
di-lepton event classes are presented in figures~\ref{fig:Masses}(c), (d) and (e). 
%
An overall agreement with the SM prediction is observed in all cases.
% 

High invariant mass events ($M_{12} > 100$~GeV) are observed in the data. 
The corresponding observed and predicted event yields are summarised for all channels in table~\ref{tab:mlepyieldsM100}.
%
%
The three $ee$, three $eee$ and one $\mu\mu$ high mass events have already been discussed extensively in previous H1 publications~\cite{Aktas:2003jg,Aktas:2003sz}. 
%
No additional events in these classes are observed in the new data.
%
One $e\mu$ and two $e\mu\mu$ high mass events are observed in the new data.
%
The $e\mu\mu$ event with the largest invariant mass $M_{12} = 127 \pm 10$~GeV is presented in figure~\ref{fig:display}.
In this event, $M_{12}$ is formed by the electron and the highest $P_T$ muon.
In the other $e\mu\mu$ event, $M_{12}$ is formed by the two muons. 
%
%
%
No $eeee$ or $ee\mu$ event is observed with a di-lepton invariant mass above $100$~GeV, in agreement with the corresponding SM expectations below $0.01$.


The results for $e^+p$ and $e^-p$ data are shown separately in table~\ref{tab:mlepyieldsM100}. 
The di-lepton ($ee$, $e\mu$ and $\mu\mu$) and tri-lepton ($eee$ and $e\mu\mu$) events  at high mass $M_{12}>100$~GeV are all observed in $e^+p$ collisions 
whereas no such event is observed in  the $e^-p$ data.
%
% 


The topology of the multi-lepton events can be further investigated using 
the scalar sum of the lepton transverse momenta $\sum P_T$. 
%
This variable indicates the ``hardness'' of the event and also offers a good sensitivity for new physics searches~\cite{Aktas:2004pz}. 
%
Figure~\ref{fig:2Dcorrelations} presents the correlation between $\sum P_T$ and the invariant mass $M_{12}$, separately for di-lepton and tri-lepton events.
%
It can be seen that tri-lepton events may have a large $M_{12}$ but only intermediate $\sum P_T$. 
In such cases, the high mass is formed by one forward and one backward lepton.  
%
However, for di-lepton event classes high $M_{12}$ values also imply  a large $\sum P_T$.

%
The quantity $\sum P_T$ can be used to select the most energetic events and allows the combination of different topologies  of di-lepton and tri-lepton events with electrons and muons.
%
Figure~\ref{fig:SumEt_All_lep} presents the distributions of $\sum P_T$ of the observed multi-lepton events compared to the SM expectation. 
%
A good overall agreement between the data and the SM prediction is observed.
%
For $\sum P_T >$~$100$~GeV, five events are
observed in total, compared to $1.60 \pm 0.20$ expected from the SM (see table~\ref{tab:mlepyieldsEt100}). 
%
These five events were all recorded in the $e^+p$ data, for which the SM expectation is $0.96$~$\pm$~$0.12$. 
%
Furthermore, the events correspond to the three $ee$ and the two $e\mu\mu$ events observed with $M_{12} > 100$~GeV.




\subsection{Cross Section Measurements}

Cross sections of the production of electron and muon pairs from photon-photon collisions are measured using the selected $(\gamma\gamma)_{e}$ and $(\gamma\gamma)_{\mu}$ samples.
The kinematic domain of the measurement is defined by $20^\circ < \theta^{\ell_{1,2}} < 150^\circ$, $P_T^{\ell_1} > 10$~GeV,  $P_T^{\ell_2} > 5$~GeV, $Q^2 < 1$~GeV$^2$ and $y < 0.82$.
%
%
The data samples collected at $\sqrt{s} = 301$~GeV and $319$~GeV are combined taking into account their respective luminosities.
Assuming a linear dependence of the cross section on the proton beam energy, as predicted by the SM, the resulting cross section corresponds to an effective $\sqrt{s} = 318$~GeV.
%
The total numbers of observed  $(\gamma\gamma)_{e}$ and $(\gamma\gamma)_{\mu}$ events are in agreement with the SM expectations, as summarised in table~\ref{tab:mlepyields}. 
%
In the $(\gamma\gamma)_{e}$ sample, the contamination from NC DIS and QED Compton background events is  $2$\%.
No significant background is present in the $(\gamma\gamma)_{\mu}$ sample.

The cross section is evaluated in each bin $i$ using the formula

\begin{equation}
\sigma_i = \frac{N_i^{\rm{data}}-N_i^{\rm{bgr}}}{ {\cal L} \cdot A_i},
\label{eq:xsection}
\end{equation}

\noindent where $N_i^{\rm{data}}$ is the number of observed events in bin $i$, $N_i^{\rm{bgr}}$ the expected contribution from background processes in bin $i$, ${\cal L}$ the integrated luminosity of the data and $A_i$ the signal acceptance in bin $i$. 
The signal acceptance is calculated using GRAPE MC events, as the ratio of the number of events reconstructed in bin $i$ divided by the number of events generated in the same bin. 
It accounts for detection efficiencies and migrations between bins.
%
%
The mean signal acceptance is $45$\% for $ep \rightarrow e \: e^+e^- X$ events and $60$\% for $ep \rightarrow e \: \mu^+\mu^- X$ events.
%
The systematic error of the measured cross section is determined by repeating the analysis after applying the appropriate variations to the MC for each systematic source, as described in section~\ref{sec:syst}.
%

The measured $ep \rightarrow e \: e^+e^- X$ cross section integrated over the phase space defined above is $\sigma = 0.67 \pm 0.06 \pm 0.05$~pb, where the first error is statistical and the second systematic. 
%
The measured cross section for muon pair production, $ep \rightarrow e \: \mu^+\mu^- X$, in the same phase space is $\sigma = 0.63 \pm 0.05 \pm 0.09$~pb.
%
The results are in agreement with the SM expectation of $0.63 \pm 0.02$~pb calculated using the GRAPE generator.
%
Combining the di-electron and di-muon samples, an average lepton pair production cross section of $\sigma = 0.65 \pm 0.04 \pm 0.06$~pb is measured.


The differential cross sections of lepton pair production as a function of the transverse momentum of the leading lepton $P_T^{\ell_1}$, the invariant mass of the lepton pair $M_{\ell\ell}$ and the hadronic transverse momentum $P_T^X$ are listed for each sample in table~\ref{tab:xsection}  and shown in figure~\ref{fig:XSec_comb} for the combined electron and muon samples.
%
The measurements are in good agreement with the SM cross sections.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Conclusion}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

The production of multi-lepton (electron or muon) events at high transverse momenta is studied in $e^+p$ and $e^-p$ scattering.
%
%
The full $e^\pm p$ data sample collected by the H1 experiment at HERA with an integrated luminosity of $463$~pb$^{-1}$ is analysed.
%
The yields of di-lepton and tri-lepton events are in good agreement with the SM predictions.
%
% 
%
In each sample distributions of the invariant mass $M_{12}$ of the two highest $P_T$ leptons and of the scalar sum of the lepton transverse momenta $\sum P_T$ are studied and found to be in good overall agreement with the SM expectation. 

Events are observed in the di-lepton and tri-lepton channels with high invariant masses $M_{12}$ above $100$~GeV. All such events are observed in $e^+p$ collisions.
%
%
Five of them have a $\sum P_T > 100$~GeV, whereas the corresponding SM expectation for $e^+p$ collisions is $0.96 \pm 0.12$. 
%


%
Differential cross sections for electron and muon pair production are measured in a restricted phase space dominated by photon-photon interactions.
%
The measured cross sections are in agreement with the SM expectations. 
%



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

We are grateful to the HERA machine group whose outstanding
efforts have made this experiment possible. 
We thank the engineers and technicians for their work in constructing 
and 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}



%%%%%%% gamma-gamma processes
%\cite{Vermaseren:1982cz}
\bibitem{Vermaseren:1982cz}
J.~A.~M.~Vermaseren,
%``Two Photon Processes At Very High-Energies,''
Nucl.\ Phys.\ B {\bf 229} (1983) 347.
%%CITATION = NUPHA,B229,347;%%

%%%%%% H++, bileptons
%\cite{Accomando:1993ar}
\bibitem{Accomando:1993ar}
  E.~Accomando and S.~Petrarca,
%  ``Searching A Doubly Charged Higgs Boson At Hera,''
  Phys.\ Lett.\  B {\bf 323} (1994) 212
  [hep-ph/9401242].
  %%CITATION = PHLTA,B323,212;%%

%\cite{Aktas:2006nu}
\bibitem{Aktas:2006nu}
  A.~Aktas {\it et al.}  [H1 Collaboration],
%  ``Search for doubly-charged Higgs boson production at HERA,''
  Phys.\ Lett.\  B {\bf 638} (2006) 432
  [hep-ex/0604027].
  %%CITATION = PHLTA,B638,432;%%
%\cite{Cuypers:1996ia}

\bibitem{Cuypers:1996ia}
  F.~Cuypers and S.~Davidson,
%  ``Bileptons: Present limits and future prospects,''
  Eur.\ Phys.\ J.\  C {\bf 2} (1998) 503
  [hep-ph/9609487].
  %%CITATION = EPHJA,C2,503;%%

%%%%%%% previous H1 analyses
%\cite{Aktas:2003jg}
\bibitem{Aktas:2003jg}
A.~Aktas {\it et al.}  [H1 Collaboration],
%``Multi-electron production at high transverse momenta in e p collisions at HERA,''
Eur.\ Phys.\ J.\ C {\bf 31} (2003) 17
[hep-ex/0307015].
%%CITATION = HEP-EX 0307015;%%


%\cite{Aktas:2003sz}
\bibitem{Aktas:2003sz}
  A.~Aktas {\it et al.}  [H1 Collaboration],
%``Muon pair production in e p collisions at HERA,''
Phys.\ Lett.\ B {\bf 583} (2004) 28
[hep-ex/0311015].
%%CITATION = HEP-EX 0311015;%%


%\cite{Abe:2000cv}
\bibitem{Abe:2000cv}
T.~Abe, GRAPE-Dilepton version 1.1,
%``GRAPE-Dilepton (Version 1.1): A generator for dilepton production in e  p collisions,''
Comput.\ Phys.\ Commun.\  {\bf 136} (2001) 126
[hep-ph/0012029].
%%CITATION = HEP-PH 0012029;%%

%\cite{Brasse:1976bf}
\bibitem{Brasse:1976bf}
  F.~W.~Brasse {\it et al.},
%  ``Parametrization Of The Q**2 Dependence Of Virtual Gamma P Total
%  Cross-Sections In The Resonance Region,''
  Nucl.\ Phys.\  B {\bf 110} (1976) 413.
  %%CITATION = NUPHA,B110,413;%%

%\cite{Abramowicz:1997ms}
\bibitem{Abramowicz:1997ms}
  H.~Abramowicz and A.~Levy,
  %``The ALLM parameterization of sigma(tot)(gamma* p): An update,''
  hep-ph/9712415.
  %%CITATION = HEP-PH/9712415;%%


%\cite{Pumplin:2002vw}
\bibitem{Pumplin:2002vw}
  J.~Pumplin {\it et al.},
%  ``New generation of parton distributions with uncertainties from global  QCD analysis,''
  JHEP {\bf 0207} (2002) 012
  [hep-ph/0201195].
  %%CITATION = JHEPA,0207,012;%%


%\cite{Sjostrand:2000wi}
\bibitem{Sjostrand:2000wi}
T.~Sj\"ostrand {\it et al.}, PYTHIA version 6.1,
%``High-energy-physics event generation with PYTHIA 6.1,''
Comput.\ Phys.\ Commun.\  {\bf 135} (2001) 238
[hep-ph/0010017].
%%CITATION = HEP-PH 0010017;%%

%\cite{Mucke:1999yb}
\bibitem{Mucke:1999yb}
A.~Mucke {\it et al.},
%``Monte Carlo simulations of photohadronic processes in astrophysics,''
Comput.\ Phys.\ Commun.\  {\bf 124} (2000) 290
[astro-ph/9903478].
%%CITATION = ASTRO-PH 9903478;%%

%\cite{Jung:1993gf}
\bibitem{Jung:1993gf}
H.~Jung, RAPGAP version 3.1,
%``Hard diffractive scattering in high-energy e p collisions and the Monte Carlo generation RAPGAP,''
Comput.\ Phys.\ Commun.\  {\bf 86} (1995) 147.
%%CITATION = CPHCB,86,147;%%


%\cite{Berger:kp}
\bibitem{Berger:kp}
C.~Berger and P.~Kandel,
``A New Generator For Wide Angle Bremsstrahlung,''
%\href{http://www.slac.stanford.edu/spires/find/hep/www?irn=4270703}{SPIRES entry}
Proceedings of the Workshop ``Monte Carlo Generators for HERA Physics'', eds. A.~T.~Doyle {\it et al.}, DESY (1998).

%\cite{Brun:1987ma}
\bibitem{Brun:1987ma}
R.~Brun {\it et al.},
``GEANT3,''
CERN-DD/EE/84-1.

  
%%%%%%%%%%%%%%%%%%%%%%% H1 detector %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
\bibitem{Abt:h1}
  I.~Abt {\it et al.}  [H1 Collaboration],
%   ``The H1 detector at HERA,''
  Nucl.\ Instrum.\ Meth.\ A {\bf 386} (1997) 310;\\
  %%CITATION = NUIMA,A386,310;%%
 I.~Abt {\it et al.}  [H1 Collaboration],
%   ``The Tracking, calorimeter and muon detectors of the H1 experiment at HERA,''
  Nucl.\ Instrum.\ Meth.\ A {\bf 386} (1997) 348.
  %%CITATION = NUIMA,A386,348;%%
 

%\cite{Andrieu:1993kh}
\bibitem{Andrieu:1993kh}
B.~Andrieu {\it et al.}  [H1 Calorimeter Group Collaboration],
%``The H1 liquid argon calorimeter system,''
Nucl.\ Instrum.\ Meth.\ A {\bf 336} (1993) 460.
%%CITATION = NUIMA,A336,460;%%



%%%%%%%%%% test beam LAr calibrations %%%%%%%%
%\cite{Andrieu:1993tz}
\bibitem{Andrieu:1993tz}
  B.~Andrieu {\it et al.}  [H1 Calorimeter Group Collaboration],
%  ``Results From Pion Calibration Runs For The H1 Liquid Argon Calorimeter And
%  Comparisons With Simulations,''
  Nucl.\ Instrum.\ Meth.\  A {\bf 336} (1993) 499.
  %%CITATION = NUIMA,A336,499;%%

 %\cite{Andrieu:1994yn}
\bibitem{Andrieu:1994yn}
  B.~Andrieu {\it et al.}  [H1 Calorimeter Group Collaboration],
%  ``Beam tests and calibration of the H1 liquid argon calorimeter with
%  electrons,''
  Nucl.\ Instrum.\ Meth.\  A {\bf 350} (1994) 57.
  %%CITATION = NUIMA,A350,57;%%
  
  
  
  
%\cite{Appuhn:1996na}
\bibitem{Appuhn:1996na}
R.~D.~Appuhn {\it et al.}  [H1 SPACAL Group Collaboration],
%``The H1 lead/scintillating-fibre calorimeter,''
Nucl.\ Instrum.\ Meth.\ A {\bf 386} (1997) 397.
%%CITATION = NUIMA,A386,397;%%


%%%%%%%%% CJC resolution
%\cite{Kleinwort:2006zz}
\bibitem{Kleinwort:2006zz}
  C.~Kleinwort  [H1 Collaboration],
  ``H1 Alignment Experience,''
%\href{http://www.slac.stanford.edu/spires/find/hep/www?irn=7543131}{SPIRES entry}
%{\it Prepared for 1^{st} LHC Detection Alignment Workshop, Geneva, Switzerland, 4-6 Sep 2006}
Proceedings of the Workshop ``First LHC Detector Alignment Workshop'', eds. 
S.~Blusk {\it et al.}, CERN (2006), CERN-2007-04.

%%%%%%%%%%%%%%%% trigger efficiencies %%%%%%%%%%%%%%%%%%%
%\cite{Adloff:2003uh}
\bibitem{Adloff:2003uh}
C.~Adloff {\it et al.}  [H1 Collaboration],
%``Measurement and QCD analysis of neutral and charged current cross  sections at HERA,''
Eur.\ Phys.\ J.\ C {\bf 30} (2003) 1 [hep-ex/0304003].
%%CITATION = HEP-EX 0304003;%%


%%%%%%%%%%%%%%%% mu ID : isolated leptons %%%%%%%%%%%%%%%%%%%
%\cite{Andreev:2003pm}
\bibitem{Andreev:2003pm}
V.~Andreev {\it et al.}  [H1 Collaboration],
%``Isolated electrons and muons in events with missing transverse momentum  at HERA,''
Phys.\ Lett.\ B {\bf 561} (2003) 241 [hep-ex/0301030].
%%CITATION = HEP-EX 0301030;%%


%%%%%%%%%%%%%%%% HFS finder %%%%%%%%%%%%%%%%%%%
\bibitem{matti}
M.~Peez, ``Recherche de d\'eviations au Mod\`ele Standard dans les processus de grande \'energie 
transverse sur le collisionneur \'electron - proton HERA'', Ph.D. thesis, Universit\'e de Lyon (2003), DESY-THESIS-2003-023.
(available at http://www-h1.desy.de/psfiles/theses/).

\bibitem{benji}
B.~Portheault, ``Premi\`ere mesure des sections efficaces de courant charg\'e et neutre avec le faisceau de positrons polaris\'e \`a HERA II et analyses QCD-\'electrofaibles'', Ph.D. thesis, Universit\'e Paris XI (2005), LAL-05-05.
(available at http://www-h1.desy.de/psfiles/theses/).


%%%%%%%%%%%%%%%% kT jets %%%%%%%%%%%%%%%%%%%
%\cite{Ellis:1993tq}
\bibitem{Ellis:1993tq}
S.~D.~Ellis and D.~E.~Soper,
%``Successive combination jet algorithm for hadron collisions,''
Phys.\ Rev.\ D {\bf 48} (1993) 3160
[hep-ph/9305266].
%%CITATION = HEP-PH 9305266;%%
 
%\cite{Catani:1993hr}
\bibitem{Catani:1993hr}
S.~Catani {\it et al.},
%, Y.~L.~Dokshitzer, M.~H.~Seymour and B.~R.~Webber,
%``Longitudinally invariant K(t) clustering algorithms for hadron-hadron collisions,''
Nucl.\ Phys.\ B {\bf 406} (1993) 187.
%%CITATION = NUPHA,B406,187;%%


%%%%%%%%%% general search HERA-1 %%%%%%%%%%%%
%\cite{Aktas:2004pz}
\bibitem{Aktas:2004pz}
  A.~Aktas {\it et al.}  [H1 Collaboration],
%  ``A general search for new phenomena in e p scattering at HERA,''
  Phys.\ Lett.\  B {\bf 602} (2004) 14
  [hep-ex/0408044].
  %%CITATION = PHLTA,B602,14;%%

 
  
\end{thebibliography}




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\clearpage

\begin{table}[]
\begin{center}
\begin{tabular}{ c c c c c }
\multicolumn{5}{c}{Multi-Leptons at HERA ($463$ pb$^{-1}$)}\\
\hline
Selection & Data & SM & Pair Production (GRAPE) & NC DIS + Compton \\
\hline                                        
$ee$ & $368$ & $390 \pm 46$ & $332 \pm 26$  & $58 \pm 30$ \\ 
%
$\mu\mu$ & $201$   & $211 \pm 32$ & $211 \pm 32$ &  $< 0.005$ \\
% 
$e\mu$ & $132$   & $128 \pm 9$~~ & $118 \pm 8$~~ & $\!\!\!10.0 \pm 2.5$ \\
% 
%\hline
$eee$ & $73$ & $70 \pm 7$ & $69.8 \pm 7.0$ & $0.2 \pm 0.1$ \\    
%
$e\mu\mu$ & $97$ & $102 \pm 14$  & $102 \pm 14$  &  $< 0.005$    \\   
%\hline
$ee\mu$   & $4$  & ~~$1.43 \pm 0.26$ & ~~$1.18 \pm 0.20$ & $0.25 \pm 0.14$ \\ 
%
$eeee$   &  $1$  & ~~$0.33 \pm 0.07$ & ~~$0.33 \pm 0.07$ & $< 0.005$  \\ 
 
\hline
$(\gamma\gamma)_{e}$ & $146$  & $138 \pm 12$ & $135 \pm 11$ & $3.0 \pm 1.0$\\ 
%   
$(\gamma\gamma)_{\mu}$ & $163$ & $162 \pm 24$ & $162 \pm 24$ & $< 0.005$  \\  
%  
\hline
\end{tabular}
\end{center}
\caption{Observed and predicted event yields for the different event classes.
  The errors on the predictions include model uncertainties and experimental systematic errors added in quadrature. The limits on the background estimations correspond to the selection of no event in the simulated sample and are quoted at $95$\% confidence level.}
\label{tab:mlepyields}
\end{table}



%%%%%%%%%%%%%%%%%%% V08C %%%%%%%%%%%%%%%%%%%%%%
\begin{table}[]
\begin{center}
\begin{tabular}{ c c c c c }
\multicolumn{5}{c}{Multi-Leptons at HERA ($463$ pb$^{-1}$)}\\
\hline
\multicolumn{5}{c}{$M_{12}>$$100$~GeV}\\
\hline
Selection & Data & SM & Pair Production (GRAPE) & NC DIS + Compton \\
\hline
\multicolumn{5}{c}{All data ($463$ pb$^{-1}$)}\\
\hline                                        
%
$ee$  & $3$ & $1.34 \pm 0.20$ & $0.83 \pm 0.11$  & $0.51 \pm 0.13$ \\ 
%
$\mu\mu$  & $1$  & $0.17 \pm 0.07$ & $0.17 \pm 0.07$ &  $< 0.005$  \\ 
%
$e\mu$  & $1$   & $0.59 \pm 0.06$ & $0.59 \pm 0.06$ & $< 0.005$  \\ 
%
$eee$   & $3$ & $0.66 \pm 0.09$ & $0.66 \pm 0.09$ &  $< 0.005$    \\    
%
$e\mu\mu$   & $2$ & $0.16 \pm 0.05$ & $0.16 \pm 0.05$ &  $< 0.005$    \\
\hline
\multicolumn{5}{c}{$e^+p$ collisions ($285$ pb$^{-1}$)}\\
\hline                                        
$ee$  & $3$ & $0.76 \pm 0.11$ & $0.49 \pm 0.07$  & $0.27 \pm 0.07$ \\ 
%
$\mu\mu$  & $1$  & $0.10 \pm 0.04$ & $0.10 \pm 0.04$ &  $< 0.005$  \\ 
%
$e\mu$  & $1$   & $0.35 \pm 0.04$ & $0.35 \pm 0.04$ & $< 0.005$  \\ 
%
%
$eee$   & $3$ & $0.39 \pm 0.05$ & $0.39 \pm 0.05$ &  $< 0.005$    \\    
%
$e\mu\mu$ & $2$ & $0.09 \pm 0.03$ & $0.09 \pm 0.03$ &  $< 0.005$    \\
%
\hline                                        
\multicolumn{5}{c}{$e^-p$ collisions ($178$ pb$^{-1}$)}\\
\hline                                        
$ee$  & $0$ & $0.58 \pm 0.09$ & $0.34 \pm 0.04$ & $0.24 \pm 0.07$ \\ 
%
$\mu\mu$  & $0$  & $0.07 \pm 0.03$ & $0.07 \pm 0.03$ &  $< 0.005$  \\ 
%
$e\mu$  & $0$   & $0.24 \pm 0.03$ & $0.24 \pm 0.03$  & $< 0.005$  \\ 
%
%
$eee$   & $0$ & $0.27 \pm 0.04$ & $0.27 \pm 0.04$ &  $< 0.005$    \\    
%
$e\mu\mu$  & $0$ & $0.07 \pm 0.03$ & $0.07 \pm 0.03$ &  $< 0.005$    \\
\hline                                        
%
% \hline
\end{tabular}
\end{center}
\caption{Observed and predicted multi-lepton event yields for masses $M_{12} > 100$~GeV for the different event classes in all analysed samples. 
  The errors on the predictions include model uncertainties and experimental systematic errors added in quadrature. The limits on the background estimations correspond to the selection of no event in the simulated sample and are quoted at $95$\% confidence level.}
\label{tab:mlepyieldsM100}
\end{table}



\begin{table}[]
\begin{center}
\begin{tabular}{ c c c c c }
\multicolumn{5}{c}{Multi-Leptons at HERA ($463$ pb$^{-1}$)}\\
\hline
\multicolumn{5}{c}{$\sum P_T>$$100$ GeV}\\
\hline
Data sample & Data & SM & Pair Production (GRAPE) & NC DIS + Compton \\
\hline                                        
e$^{+}$p ($285$ pb$^{-1}$)  & $5$ & $0.96 \pm 0.12$ & $0.78 \pm 0.09$  & $0.18 \pm 0.05$ \\ 
%
e$^{-}$p ($178$ pb$^{-1}$)  & $0$ & $0.64 \pm 0.09$ & $0.51 \pm 0.07$  & $0.13 \pm 0.04$ \\ 
%
All      ($463$ pb$^{-1}$)  & $5$ & $1.60 \pm 0.20$ & $1.29 \pm 0.15$  & $0.31 \pm 0.09$ \\ 
\hline                                        
\end{tabular}
\end{center}
\caption{Observed and predicted multi-lepton event yields for $\sum P_T >$ $100$~GeV. Di-lepton and tri-lepton events are combined.
  The errors on the predictions include model uncertainties and experimental systematic errors added in quadrature.}
\label{tab:mlepyieldsEt100}
\end{table}



% 
%%%%%%%%  Cross section table %%%%%%%%%%%%
\begin{center}
 \renewcommand{\arraystretch}{1.15} 
 \begin{table}[]
%\footnotesize{
\begin{center}
%\rotatebox{90}{
\begin{tabular}{ l l l l l l l l l l l l  }
\multicolumn{12}{c}{Multi-Leptons at HERA ($463$~pb$^{-1}$)}\\
\hline
\multicolumn{1}{c}{Variable}  &  \multicolumn{3}{c}{Measured} & \multicolumn{3}{c}{Measured}        &    \multicolumn{3}{c}{Measured}  & \multicolumn{2}{c}{Pair Production} \\
\multicolumn{1}{c}{ range}    &   \multicolumn{3}{c}{  ($e^+e^-$) }  &   \multicolumn{3}{c}{ ($\mu^+\mu^-$)}  &    \multicolumn{3}{c}{(average)}  &  \multicolumn{2}{c}{ (GRAPE)}     \\
\multicolumn{1}{c}{\mbox{[GeV]}} & \multicolumn{3}{c}{ [fb/GeV] }  &  \multicolumn{3}{c}{[fb/GeV] }  &   \multicolumn{3}{c}{[fb/GeV]}    &   \multicolumn{2}{c}{[fb/GeV]}   \\
\hline
\multicolumn{1}{c}{$P_T^{\ell_1}$} & \multicolumn{10}{c}{ $d\sigma/dP_T^{\ell_1}$}\\
\hline
$[10,15]$ & ~~~~$108$&$\!\!\!\!\!\pm \, 10$&$\!\!\!\!\!\pm \, 8$&~~~~$95$&$\!\!\!\!\!\pm \, 9$&$\!\!\!\!\!\pm \, 14$&~~~~$101$&$\!\!\!\!\!\pm \, 7$&$\!\!\!\!\!\pm \, 9$&~~~~$92$&$\!\!\!\!\!\pm \, 3$ \\
$[15,20]$ & ~~~~$16$&$\!\!\!\!\!\pm \, 4$&$\!\!\!\!\!\pm \, 2$&~~~~$23$&$\!\!\!\!\!\pm \, 4$&$\!\!\!\!\!\pm \, 4$&~~~~$20$&$\!\!\!\!\!\pm \, 3$&$\!\!\!\!\!\pm \, 2$&~~~~$22$&$\!\!\!\!\!\pm \, 1$ \\
$[20,25]$ & ~~~~$7.4$&$\!\!\!\!\!\pm \, 2.5$&$\!\!\!\!\!\pm \, 0.6$&~~~~$3.7$&$\!\!\!\!\!\pm \, 1.7$&$\!\!\!\!\!\pm \, 0.5$&~~~~$5.5$&$\!\!\!\!\!\pm \, 1.5$&$\!\!\!\!\!\pm \, 0.5$&~~~~$6.6$&$\!\!\!\!\!\pm \, 0.2$ \\ 
$[25,50]$ & ~~~~$0.60$&$\!\!\!\!\!\pm \, 0.30$&$\!\!\!\!\!\pm \, 0.06$&~~~~$0.70$&$\!\!\!\!\!\pm \, 0.30$&$\!\!\!\!\!\pm \, 0.10$&~~~~$0.70$&$\!\!\!\!\!\pm \, 0.20$&$\!\!\!\!\!\pm \, 0.07$&~~~~$0.86$&$\!\!\!\!\!\pm \, 0.03$ \\ %
\hline
\multicolumn{1}{c}{$M_{\ell\ell}$}  & \multicolumn{10}{c}{ $d\sigma/dM_{\ell\ell}$}\\
\hline
$[15,25]$   & ~~~~$34$&$\!\!\!\!\!\pm \, 4$&$\!\!\!\!\!\pm \, 3$&~~~~$29$&$\!\!\!\!\!\pm \, 3$&$\!\!\!\!\!\pm \, 4$&~~~~$32$&$\!\!\!\!\!\pm \, 3$&$\!\!\!\!\!\pm  \,3$&~~~~$28$&$\!\!\!\!\!\pm  \,1$ \\
$[25,40]$   & ~~~~$16$&$\!\!\!\!\!\pm  \,2$&$\!\!\!\!\!\pm \, 2$&~~~~$17$&$\!\!\!\!\!\pm \, 2$&$\!\!\!\!\!\pm  \,3$&~~~~$17$&$\!\!\!\!\!\pm \, 2$&$\!\!\!\!\!\pm  \,2$&~~~~$18.0$&$\!\!\!\!\!\pm \, 0.5$ \\
$[40,60]$   & ~~~~$2.7$&$\!\!\!\!\!\pm \, 0.8$&$\!\!\!\!\!\pm \, 0.2$&~~~~$1.7$&$\!\!\!\!\!\pm \, 0.6$&$\!\!\!\!\!\pm \, 0.3$&~~~~$2.2$&$\!\!\!\!\!\pm \, 0.5$&$\!\!\!\!\!\pm \, 0.2$&~~~~$2.9$&$\!\!\!\!\!\pm \, 0.1$ \\ 
$[60,100]$  & ~~~~$0.20$&$\!\!\!\!\!\pm \, 0.15$&$\!\!\!\!\!\pm \, 0.02$&~~~~$0.40$&$\!\!\!\!\!\pm \, 0.20$&$\!\!\!\!\!\pm \, 0.07$&~~~~$0.30$&$\!\!\!\!\!\pm \, 0.13$&$\!\!\!\!\!\pm \, 0.04$&~~~~$0.24$&$\!\!\!\!\!\pm  \,0.01$ \\ 
%
%
\hline
\multicolumn{1}{c}{$P_T^X$}  & \multicolumn{10}{c}{ $d\sigma/dP_T^X$}\\
\hline
$[0,5]$   & ~~~~$95$&$\!\!\!\!\!\pm \, 9$&$\!\!\!\!\!\pm \, 8$&~~~~$94$&$\!\!\!\!\!\pm \, 8$&$\!\!\!\!\!\pm \, 14$&~~~~$94$&$\!\!\!\!\!\pm \, 6$&$\!\!\!\!\!\pm \, 9$&~~~~$93$&$\!\!\!\!\!\pm \, 3$ \\
$[5,12]$  & ~~~~$25$&$\!\!\!\!\!\pm \, 6$&$\!\!\!\!\!\pm \, 2$&~~~~$13$&$\!\!\!\!\!\pm \, 4$&$\!\!\!\!\!\pm \, 2$&~~~~$18$&$\!\!\!\!\!\pm \, 3$&$\!\!\!\!\!\pm \, 2$&~~~~$15.4$&$\!\!\!\!\!\pm \, 0.5$ \\
$[12,25]$ & ~~~~$2.8$&$\!\!\!\!\!\pm \, 1.3$&$\!\!\!\!\!\pm \, 0.3$&~~~~$4.2$&$\!\!\!\!\!\pm \, 1.2$&$\!\!\!\!\!\pm \, 0.6$&~~~~$3.8$&$\!\!\!\!\!\pm \, 0.9$&$\!\!\!\!\!\pm \, 0.4$&~~~~$3.8$&$\!\!\!\!\!\pm \, 0.1$ \\ 
$[25,50]$ & ~~~~$0.20$&$\!\!\!\!\!\pm  \,0.30$&$\!\!\!\!\!\pm \, 0.03$&~~~~$0.50$&$\!\!\!\!\!\pm \, 0.30$&$\!\!\!\!\!\pm \, 0.09$&~~~~$0.40$&$\!\!\!\!\!\pm \, 0.20$&$\!\!\!\!\!\pm \, 0.06$&~~~~$0.20$&$\!\!\!\!\!\pm \, 0.01$ \\ 
%
\hline                                        
\end{tabular}
%}
\end{center}
\caption{Differential cross sections $d\sigma/dP_T^{\ell_1}$, $d\sigma/dM_{\ell\ell}$ and $d\sigma/dP_T^X$ averaged for each quoted interval for the process $ep \rightarrow e \ell^+ \ell^- X$ in a restricted phase space dominated by the photon-photon process (see text for details). 
Cross sections are measured for $e^+e^-$ or $\mu^+\mu^-$ pairs.
The cross section obtained from the combination of $e^+e^-$ and $\mu^+\mu^-$ channels is also presented.
The first error is statistical and the second is systematic.
Theoretical predictions, calculated with GRAPE, are shown in the last column.}
\label{tab:xsection}
%}
\end{table}
\end{center}





%=========================================================================
\vfill
\newpage


\begin{figure}[!htbp] 
  \begin{center}
  \begin{center}
 \includegraphics[width=.5\textwidth]{d08-065f1a.eps}\put(-11,33){{(a)}}
 \includegraphics[width=.5\textwidth]{d08-065f1b.eps}\put(-11,33){{(b)}}\\
 \includegraphics[width=.5\textwidth]{d08-065f1c.eps}\put(-11,33){{(c)}}
 \includegraphics[width=.5\textwidth]{d08-065f1d.eps}\put(-11,33){{(d)}}\\
\hspace{-.5\textwidth} \includegraphics[totalheight=6.4cm]{d08-065f1e.eps}\put(-11,33){{(e)}}
 \end{center}
  \end{center}
  \caption{The distribution of the invariant mass of the two highest $P_T$ leptons for events classified as $eee$ (a), $e\mu\mu$ (b) and $ee$ (c), $\mu\mu$ (d) and $e\mu$ (e). 
  The points correspond to the observed data events and the open histogram to the SM expectation. The total error on the SM expectation is given by the shaded band. The component of the SM expectation arising from lepton pair production is given by the hatched histogram.
   }
\label{fig:Masses}  
\end{figure} 




\begin{figure}[!htbp] 
  \begin{center}
 \includegraphics[width=.65\textwidth]{d08-065f2.eps}
  \end{center}
  \caption{Display (side view) of the $e\mu\mu$ event observed with the highest $M_{12}$ invariant mass. Indicated are the reconstructed tracks and the energy depositions in the calorimeters. The beam positrons enter the detector from the left and the protons from the right.
   }
\label{fig:display}  
\end{figure} 



%%%%%%%%%%% 2D plots SPt vs M12
\begin{figure}[!htbp] 
  \begin{center}
\includegraphics[width=.5\textwidth]{d08-065f3a.eps}\put(-11,15){{(a)}}
 \includegraphics[width=.5\textwidth]{d08-065f3b.eps}\put(-11,15){{(b)}}\\
  \end{center}
  \caption{Correlation of the invariant mass $M_{12}$ with the scalar sum of the transverse momenta $\sum P_T$ for di-lepton (a) and tri-lepton (b) events. The bold dots represent the data while the small points represent the pair production (GRAPE) prediction for a luminosity $\sim 1000$ times higher than that of data.
   }
\label{fig:2Dcorrelations}  
\end{figure} 



%=========================================================================
\vfill
\newpage



\begin{figure}[htbp] 
\begin{center}
 \includegraphics[width=.5\textwidth]{d08-065f4a.eps}\put(-11,33){{(a)}}\\
 \includegraphics[width=.5\textwidth]{d08-065f4b.eps}\put(-11,33){{(b)}}\\
 \includegraphics[width=.5\textwidth]{d08-065f4c.eps}\put(-11,33){{(c)}}
\end{center}
\vspace{-0.5cm}
  \caption{The distribution of the scalar sum of the transverse momenta $\sum P_T$ for combined di-lepton and tri-lepton event samples for all data (a) as well as for $e^+p$ (b) and $e^-p$ (c).
  The points correspond to the observed data events and the open histogram to the SM expectation. The total error on the SM expectation is given by the shaded band. The component of the SM expectation arising from lepton pair production is given by the hatched histogram.
}
\label{fig:SumEt_All_lep}
\end{figure}






\begin{figure}[htbp] 
\begin{center}
\includegraphics[totalheight=13cm]{d08-065f5.eps}\put(-87,47){{(c)}}\put(-87,112){{(a)}}\put(-12,112){{(b)}}
\end{center}
  \caption{The measured cross section for lepton pair production in a restricted phase space dominated by the photon-photon process as a function of the leading lepton transverse momentum $P_T^{\ell_1}$ (a), the invariant mass of the lepton pair $M_{\ell\ell}$ (b) and the hadronic transverse momentum $P_T^X$ (c).
The differential cross section is averaged over the intervals shown.
The inner error bars represent the statistical errors, the outer error bars the statistical and systematic errors added in quadrature.
The bands represent the SM prediction with its one standard deviation uncertainty.  }
\label{fig:XSec_comb}
\end{figure}




\end{document}




