\documentclass[12pt]{article} \usepackage{epsfig} \usepackage{amsmath} \usepackage{hhline} \usepackage{amssymb} \usepackage{times} \usepackage{cite} \usepackage{psfrag} \usepackage{color} \definecolor{vlMJ}{rgb}{0.2,0.6,0.8} \definecolor{elMJ}{rgb}{0.53,0.4,0.34} \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============================= \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{\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}} \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^*} } \newcommand{\Rp}{\mbox{$\not \hspace{-0.15cm} R_p$}} %\newcommand{\lsim}{\raisebox{-1.5mm}{$\:\stackrel{\textstyle{<}}{\textstyle{\sim}}\:$}} %\newcommand{\gsim}{\raisebox{-0.5mm}{$\stackrel{>}{\scriptstyle{\sim}}$}} \newcommand{\GeV}{\rm GeV} \newcommand{\TeV}{\rm TeV} \newcommand{\pb}{\rm pb} \newcommand{\cm}{\rm cm} \newcommand{\hdick}{\noalign{\hrule height1.4pt}} % 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.} \begin{titlepage} \noindent \begin{flushleft} {\tt DESY 10-181 \hfill ISSN 0418-9833} \\ {\tt October 2010} \\ \end{flushleft} \vfill \begin{center} \begin{Large} {\bf Search for Squarks in \boldmath{$R$}--parity Violating Supersymmetry in \boldmath{$ep$} Collisions at HERA } \vspace{2cm} H1 Collaboration \end{Large} \end{center} \vspace{2cm} \begin{abstract} A search for squarks in $R$--parity violating supersymmetry is performed in $e^{\pm} p$ collisions at HERA using the H1 detector. The full data sample taken at a centre--of--mass energy $\sqrt{s}=319$~GeV is used for the analysis, corresponding to an integrated luminosity of $255$~pb$^{-1}$ of $e^{+}p$ and $183$~pb$^{-1}$ of $e^{-}p$ collision data. The resonant production of squarks via a Yukawa coupling $\lambda'$ is considered, taking into account direct and indirect $R$--parity violating decay modes. Final states with jets and leptons are investigated. No evidence for squark production is found and mass dependent limits on $\lambda'$ are obtained in the framework of the Minimal Supersymmetric Standard Model and in the Minimal Supergravity Model. In the considered part of the parameter space, for a Yukawa coupling of electromagnetic strength $\lambda'= 0.3$, squarks of all flavours are excluded up to masses of $275$~GeV at $95\%$ confidence level, with down--type squarks further excluded up to masses of $290$~GeV. \end{abstract} \vspace{1.5cm} \begin{center} Submitted to \EJC \end{center} \end{titlepage} \begin{flushleft} % \input{h1auts} %-- H1AUTS Author list by names %-- Status: Wed Oct 27 14:04:49 CEST 2010 Number of authors = 214 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 S.~Backovic$^{30}$, %PODG-PD 03/2 Backovic A.~Baghdasaryan$^{38}$, %YERE-PD 09/03 Baghdasaryana S.~Baghdasaryan$^{38}$, %YERE-ST 02/10 Baghdasaryans E.~Barrelet$^{29}$, %PARI-PD 11/99 Barrelet W.~Bartel$^{11}$, %DESY-PD 8/88 Bartel K.~Begzsuren$^{35}$, %ULBA-PD 04/06 Begzsuren A.~Belousov$^{25}$, %LPI -PD 8/88 Belousov J.C.~Bizot$^{27}$, %ORSA-PD 8/88 Bizot 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-PD 03/10 Brinkmann V.~Brisson$^{27}$, %ORSA-PD 8/88 Brisson D.~Britzger$^{11}$, %DESY-ST 10/09 Britzger D.~Bruncko$^{16}$, %KOSI-PD 8/88 Bruncko A.~Bunyatyan$^{13,38}$, %MPIH-PD 12/95 Bunyatyan G.~Buschhorn$^{26, \dagger}$, %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.~Ceccopieri$^{4}$, %BRUX-PD 10/09 Ceccopieri K.~Cerny$^{32}$, %PRG2-PD 09/08 Cernyk V.~Cerny$^{16,47}$, %KOSI-PD 06/04 Cernyv V.~Chekelian$^{26}$, %MPIM-PD 01/90 Chekelian A.~Cholewa$^{11}$, %DESY-ST 01/10 Cholewa J.G.~Contreras$^{22}$, %MEX1-PD 04/97 Contreras J.A.~Coughlan$^{6}$, %RAL -PD 8/88 Coughlan 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 B.~Delcourt$^{27}$, %ORSA-PD 8/88 Delcourt J.~Delvax$^{4}$, %BRUX-ST 10/06 Delvax E.A.~De~Wolf$^{4}$, %ANTW-PD 3/93 Dewolf C.~Diaconu$^{21}$, %MARS-PD 01/05 Diaconu M.~Dobre$^{12,51,52}$, %HAM2-ST 07/09 Dobre 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 S.~Egli$^{36}$, %PSI -PD 01/10 Egli A.~Eliseev$^{25}$, %LPI -PD 01/06 Eliseev E.~Elsen$^{11}$, %DESY-PD 8/88 Elsen 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 D.-J.~Fischer$^{11}$, %DESY-ST 03/08 Fischer 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$^{11}$, %DFLC-PD 09/09 Ghazaryan A.~Glazov$^{11}$, %DESY-PD 01/04 Glazov L.~Goerlich$^{7}$, %CRAC-PD 8/88 Goerlich N.~Gogitidze$^{25}$, %LPI -PD 8/88 Gogitidze M.~Gouzevitch$^{11,45}$, %DESY-PD 09/10 Gouzevitch C.~Grab$^{40}$, %ZUTH-PD 8/88 Grab A.~Grebenyuk$^{11}$, %DESY-ST 03/09 Grebenyuk T.~Greenshaw$^{18}$, %LIVE-PD 8/88 Greenshaw B.R.~Grell$^{11}$, %DESY-LEFT 10/10 Grell G.~Grindhammer$^{26}$, %MPIM-PD 8/88 Grindhammer S.~Habib$^{11}$, %DESY-PD 09/09 Habib D.~Haidt$^{11}$, %DESY-PD 8/88 Haidt C.~Helebrant$^{11}$, %DFLC-PD 09/09 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 M.~Herbst$^{15}$, %HDB2-ST 06/08 Herbst G.~Herrera$^{23}$, %MEX2-PD 07/98 Herrera M.~Hildebrandt$^{36}$, %PSI -PD 01/10 Hildebrandtm K.H.~Hiller$^{39}$, %ZEUT-PD 8/88 Hiller D.~Hoffmann$^{21}$, %MARS-PD 10/0 Hoffmann R.~Horisberger$^{36}$, %PSI -PD 01/10 Horisberger T.~Hreus$^{4,44}$, %BRUX-PD 10/08 Hreus F.~Huber$^{14}$, %HDB1-ST 09/09 Huberf M.~Jacquet$^{27}$, %ORSA-PD 09/96 Jacquet X.~Janssen$^{4}$, %ANTW-PD 02/03 Janssenx L.~J\"onsson$^{20}$, %LUND-PD 8/88 Joensson A.W.~Jung$^{15}$, %HDB2-LEFT 02/10 Junga H.~Jung$^{11,4}$, %DESY-PD 07/00 Jungh M.~Kapichine$^{9}$, %JINR-PD 3/97 Kapichine J.~Katzy$^{11}$, %DESY-LEFT 02/10 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.~Kluge$^{18}$, %LIVE-PD 05/04 Kluge A.~Knutsson$^{11}$, %DESY-LEFT 12/09 Knutsson R.~Kogler$^{26}$, %MPIM-ST 01/07 Kogler P.~Kostka$^{39}$, %ZEUT-PD 8/88 Kostka M.~Kraemer$^{11}$, %DESY-PD 10/09 Kraemer J.~Kretzschmar$^{18}$, %LIVE-PD 01/08 Kretzschmar K.~Kr\"uger$^{15}$, %HDB2-PD 01/04 Kruegerk K.~Kutak$^{11}$, %DESY-LEFT 01/10 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 V.~Lendermann$^{15}$, %HDB2-PD 01/2 Lendermann S.~Levonian$^{11}$, %DESY-PD 8/88 Levonian K.~Lipka$^{11,51}$, %DESY-PD 01/03 Lipka B.~List$^{12}$, %HAM2-PD 11/99 Listb J.~List$^{11}$, %DFLC-PD 01/05 Listj N.~Loktionova$^{25}$, %LPI -LEFT 01/10 Loktionova R.~Lopez-Fernandez$^{23}$, %MEX2-PD 03/2 Lopezfernandez V.~Lubimov$^{24}$, %ITEP-PD 01/95 Lubimov A.~Makankine$^{9}$, %JINR-PD 11/02 Makankine E.~Malinovski$^{25}$, %LPI -PD 01/89 Malinovskie P.~Marage$^{4}$, %BRUX-PD 8/88 Marage 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 A.B.~Meyer$^{11}$, %DESY-PD 01/00 Meyeran H.~Meyer$^{37}$, %WUPP-PD 8/88 Meyerhi J.~Meyer$^{11}$, %DESY-PD 8/88 Meyerj 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-LEFT 02/10 Mozer M.~Mudrinic$^{2}$, %BEOG-PD 01/07 Mudrinic K.~M\"uller$^{41}$, %ZUER-PD 8/88 Muellerk 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 D.~Nikitin$^{9}$, %JINR-PD 06/08 Nikitin G.~Nowak$^{7}$, %CRAC-PD 8/88 Nowakg K.~Nowak$^{11}$, %DESY-PD 10/09 Nowakk J.E.~Olsson$^{11}$, %DESY-PD 8/88 Olsson S.~Osman$^{20}$, %LUND-PD 06/09 Osman D.~Ozerov$^{24}$, %ITEP-PD 08/08 Ozerov P.~Pahl$^{11}$, %DESY-ST 10/08 Pahl 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 E.~Perez$^{10,45}$, %SACL-PD 10/07 Perez A.~Petrukhin$^{11}$, %DESY-PD 10/09 Petrukhin I.~Picuric$^{30}$, %PODG-PD 01/06 Picuric S.~Piec$^{11}$, %DESY-PD 11/09 Piec H.~Pirumov$^{14}$, %HDB1-ST 09/09 Pirumov D.~Pitzl$^{11}$, %DESY-PD 8/88 Pitzl R.~Pla\v{c}akyt\.{e}$^{12}$, %HAM2-PD 07/10 Placakyte B.~Pokorny$^{32}$, %PRG2-ST 10/09 Pokorny R.~Polifka$^{32}$, %PRG2-ST 10/06 Polifka B.~Povh$^{13}$, %MPIH-PD 8/88 Povh V.~Radescu$^{14}$, %HDB1-PD 10/06 Radescu N.~Raicevic$^{30}$, %PODG-PD 03/2 Raicevic 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 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 S.~Rusakov$^{25}$, %LPI -PD 8/88 Rusakov D.~\v S\'alek$^{32}$, %PRG2-PD 10/10 Salek D.P.C.~Sankey$^{6}$, %RAL -PD 8/88 Sankey M.~Sauter$^{14}$, %HDB1-PD 10/09 Sauter E.~Sauvan$^{21}$, %MARS-PD 11/1 Sauvan S.~Schmitt$^{11}$, %DESY-PD 09/07 Schmittst L.~Schoeffel$^{10}$, %SACL-PD 12/98 Schoeffel A.~Sch\"oning$^{14}$, %HDB1-PD 04/09 Schoening H.-C.~Schultz-Coulon$^{15}$, %HDB2-PD 01/04 Schultzcoulon F.~Sefkow$^{11}$, %DFLC-PD 09/99 Sefkow 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 Y.~Soloviev$^{25}$, %LPI -PD 8/88 Soloviev P.~Sopicki$^{7}$, %CRAC-ST 09/07 Sopicki D.~South$^{11}$, %DESY-PD 07/10 South V.~Spaskov$^{9}$, %JINR-PD 12/97 Spaskov A.~Specka$^{28}$, %ECPL-PD 3/95 Specka Z.~Staykova$^{11}$, %DESY-PD 10/10 Staykova M.~Steder$^{11}$, %DESY-PD 09/08 Steder B.~Stella$^{33}$, %ROME-PD 8/88 Stella G.~Stoicea$^{5}$, %BUCH-PD 02/08 Stoicea U.~Straumann$^{41}$, %ZUER-PD 8/88 Straumann D.~Sunar$^{4}$, %BRUX-LEFT 11/09 Sunar T.~Sykora$^{4,32}$, %ANTW-PD 01/06 Sykora G.~Thompson$^{19}$, %QMWC-LEFT 11/09 Thompsong P.D.~Thompson$^{3}$, %BIRM-PD 08/99 Thompsonp T.~Toll$^{11}$, %DESY-LEFT 05/10 Toll T.H.~Tran$^{27}$, %ORSA-PD 03/10 Tran D.~Traynor$^{19}$, %QMWC-PD 12/01 Traynor P.~Tru\"ol$^{41}$, %ZUER-PD 8/88 Truoel I.~Tsakov$^{34}$, %SOFI-PD 04/03 Tsakov B.~Tseepeldorj$^{35,50}$, %ULBA-PD 06/06 Tseepeldorj J.~Turnau$^{7}$, %CRAC-PD 8/88 Turnau K.~Urban$^{15}$, %HDB2-PD 06/09 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 Y.~Vazdik$^{25}$, %LPI -PD 8/88 Vazdik M.~von~den~Driesch$^{11}$, %DESY-ST 06/08 Vondendriesch D.~Wegener$^{8}$, %DORT-PD 8/88 Wegener E.~W\"unsch$^{11}$, %DESY-PD 8/88 Wuensch 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.~Zhokin$^{24}$, %ITEP-PD 04/99 Zhokine H.~Zohrabyan$^{38}$, %YERE-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 \\ $ ^{2}$ Vinca Institute of Nuclear Sciences, University of Belgrade, 1100 Belgrade, Serbia \\ $ ^{3}$ School of Physics and Astronomy, University of Birmingham, Birmingham, UK$^{ b}$ \\ $ ^{4}$ Inter-University Institute for High Energies ULB-VUB, Brussels and Universiteit Antwerpen, Antwerpen, Belgium$^{ c}$ \\ $ ^{5}$ National Institute for Physics and Nuclear Engineering (NIPNE) , Bucharest, Romania$^{ m}$ \\ $ ^{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/Irfu, 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, Aix-Marseille Universit\'e, CNRS/IN2P3, Marseille, France \\ $ ^{22}$ Departamento de Fisica Aplicada, CINVESTAV, M\'erida, Yucat\'an, M\'exico$^{ j}$ \\ $ ^{23}$ Departamento de Fisica, CINVESTAV IPN, M\'exico City, M\'exico$^{ j}$ \\ $ ^{24}$ Institute for Theoretical and Experimental Physics, Moscow, Russia$^{ k}$ \\ $ ^{25}$ Lebedev Physical Institute, Moscow, Russia$^{ e}$ \\ $ ^{26}$ Max-Planck-Institut f\"ur Physik, M\"unchen, Germany \\ $ ^{27}$ LAL, Universit\'e Paris-Sud, CNRS/IN2P3, Orsay, France \\ $ ^{28}$ LLR, Ecole Polytechnique, CNRS/IN2P3, Palaiseau, France \\ $ ^{29}$ LPNHE, Universit\'e Pierre et Marie Curie Paris 6, Universit\'e Denis Diderot Paris 7, CNRS/IN2P3, 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}$ Also at Ulaanbaatar University, Ulaanbaatar, Mongolia \\ $ ^{51}$ Supported by the Initiative and Networking Fund of the Helmholtz Association (HGF) under the contract VH-NG-401. \\ $ ^{52}$ Absent on leave from NIPHE-HH, Bucharest, Romania \\ \smallskip $ ^{\dagger}$ Deceased \\ \bigskip $ ^a$ Supported by the Bundesministerium f\"ur Bildung und Forschung, FRG, under contract numbers 05H09GUF, 05H09VHC, 05H09VHF, 05H16PEA \\ $ ^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 DPN/N168/DESY/2009 \\ $ ^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-LA09042 and MSM0021620859 \\ $ ^i$ Supported by the Swiss National Science Foundation \\ $ ^j$ Supported by CONACYT, M\'exico, grant 48778-F \\ $ ^k$ Russian Foundation for Basic Research (RFBR), grant no 1329.2008.2 \\ $ ^l$ This project is co-funded by the European Social Fund (75\%) and National Resources (25\%) - (EPEAEK II) - PYTHAGORAS II \\ $ ^m$ Supported by the Romanian National Authority for Scientific Research under the contract PN 09370101 \\ } \end{flushleft} \newpage \section{Introduction} The $ep$ collider HERA is ideally suited to search for new particles coupling to electron\footnote{In the following the generic term {\it electron} refers to both electron and positron unless explicitly stated otherwise. }--quark pairs. In supersymmetric (SUSY) models with $R$--parity violation (\Rp ), squarks can couple to electrons and quarks via Yukawa couplings $\lambda'$. At HERA, squarks with masses up to the electron--proton centre--of--mass energy, \mbox{$\sqrt{s}=319$~GeV}, could be produced resonantly via the fusion of the incoming $27.6$~GeV electron and a quark from the incoming $920$~GeV proton. Squark decays typically result in a number of high energetic particles in the final state, thus several complementary \mbox{multi--lepton} and \mbox{multi--jet} topologies are investigated. The data used in this analysis correspond to an integrated luminosity of \mbox{$255\,\mathrm{pb^{-1}}$} for $e^+p$ collisions and \mbox{$183\,\mathrm{pb^{-1}}$} for $e^-p$ collisions which represents the full data sample collected at \mbox{$\sqrt{s}=319$~GeV}. For the $e^-p$ sample, this represents an increase of a factor of thirteen compared to the previous H1 analysis~\cite{Aktas:2004ij}, while for the $e^+p$ sample this corresponds to a factor of four. The search presented here supersedes the results previously obtained by H1~\cite{Aktas:2004ij,PEREZ}. Complementary direct searches for \Rp\ SUSY have been carried out at the LEP $e^+e^-$ collider~\cite{LEP,L3RPV} and at the Tevatron $p\bar{p}$ collider~\cite{TEVATRON,D0RES}. Indirect constraints from low energy precision observables are also available~\cite{BETA0NU,APV,CCU,barbier}. \section{Phenomenology and Monte Carlo Simulation} \subsection{Production of squarks in \boldmath{$R$}--parity violating supersymmetry} Supersymmetric extensions of the Standard Model (SM) introduce new elementary particles which are the superpartners (sparticles) of SM particles but differ in spin by half a unit. A new quantum number $R_p=(-1)^{3B+L+2S}$ is defined, denoted $R$--parity, where $B$ is the baryon number, $L$ the lepton number and $S$ the spin of a particle. For particles $R_p=1$ and for their supersymmetric partners $R_p=-1$. Most of the collider searches focus on SUSY models that conserve $R$--parity, allowing only pair--production of sparticles. However, the most general supersymmetric theory that is renormalisable and gauge invariant with respect to the Standard Model gauge group does not impose $R$--parity conservation. Couplings between two SM fermions and a squark ($\tilde{q}$) or a slepton ($\tilde{l}$) are then possible, allowing the single production of sparticles. The \Rp\ Yukawa couplings responsible for squark production at HERA originate from a lepton number violating term $\lambda'_{ijk} L_{i}Q_{j}\overline{D}_k$ in the superpotential, where $i,j$ and $k$ are family indices. $L_i$, $Q_j$ and $\overline{D}_k$ are superfields, which contain the left--handed leptons, the left--handed up--type quarks and the right--handed down--type quarks, respectively, together with their SUSY partners. Non--vanishing couplings $\lambda'_{1jk}$ allow the resonant production of squarks at HERA via $eq$ fusion~\cite{RPVIOLATION}. Feynman diagrams of these processes are shown in figure~\ref{fig:directdecays}. The values of the couplings are not fixed by the theory but are required to be small to conform with present observations. For simplicity, it is assumed here that one of the $\lambda'_{1jk}$ couplings dominates over all the other trilinear \Rp\ couplings. At high Bjorken--$x$ the density of antiquarks in the proton is significantly smaller than that of the valence quarks. Hence $e^-p$ scattering gives sensitivity to the couplings $\lambda'_{11k}$ $(k=1,2,3)$ which dominate the production of $\tilde{d}_R$--type squarks (i.e. the superpartners $\tilde{d}_R$, $\tilde{s}_R$ and $\tilde{b}_R$ of down--type quarks). The dominant contribution to the production cross section is thus approximately proportional to $\lambda'^2_{11k} \cdot u(x)$ where $u(x)$ gives the probability to find a $u$ quark in the proton carrying the momentum fraction $x=M^2_{\tilde{q}}/s$, where $M^{2}_{\tilde{q}}$ is the squared mass of the produced squark. By contrast, $e^+p$ scattering provides sensitivity to the couplings $\lambda'_{1j1}$ $(j=1,2,3)$ which dominate the production of $\tilde{u}_L$--type squarks (i.e. the superpartners $\tilde{u}_L$, $\tilde{c}_L$ and $\tilde{t}_L$ of up--type quarks). Here the dominant contribution to the production cross section is approximately proportional to $\lambda'^2_{1j1} \cdot d(x)$. Due to the larger $u$ quark density in the proton at large $x$ with respect to the $d$ quark density, larger production cross sections are expected in $e^-p$ interactions for identical couplings and squark masses. Signal cross sections are obtained in the narrow width approximation by using the leading order amplitudes given in~\cite{BUCHMULL}, corrected to account for next--to--leading order QCD effects using multiplicative correction factors~\cite{SPIRANLO}. The parton densities are evaluated at the hard scale $M_{\tilde{q}}^2$. \subsection{Final states from squark decays} In \Rp\ SUSY all sparticles are unstable. Squarks can decay directly via the Yukawa coupling $\lambda'$ into SM fermions. The $\tilde{d}^k_R$--type ($k=1,2,3$) squarks can decay via the coupling $\lambda'_{11k}$ either into $e^- + {u}$ or $\nu_e + {d}$, while the $\tilde{u}^j_L$--type ($j=1,2,3$) squarks decay via the coupling $\lambda'_{1j1}$ into $e^+ + d$ only, as illustrated in figure~\ref{fig:directdecays}. Squarks may also decay via $R_p$ conserving gauge couplings as illustrated in figure~\ref{fig:gaugedecays}. The $\tilde{u}_L$--type squarks can undergo a gauge decay into states involving a neutralino $\chi^0_i$ ($i=1,2,3,4$), a chargino $\chi^{+}_i$ ($i=1,2$) or a gluino $\tilde{g}$. In contrast, $\tilde{d}_R$--type squarks mainly decay to $\chi^0_i$ or $\tilde{g}$ and decays into charginos are suppressed~\cite{RPVIOLATION}. \par The squark decay chains analysed in this paper are classified by event topology\cite{PEREZ,Aktas:2004ij}. This classification relies on the number of isolated electrons, muons and hadronic jets in the final state, and on the presence of missing energy (indicating undetected neutrinos). The channels labelled $eq$ and $\nu q$ are the squark decay modes that proceed directly via \Rp\ couplings resulting in event topologies with an isolated electron or neutrino and a single jet. The remaining channels result from the gauge decays of the squark and are characterised by multijet ($M\!J$) final states with additional leptons. The channels labelled $eM\!J$ and $\nu M\!J$ involve one or two gauginos ($\chi$ or $\tilde{g}$) in the decay cascade. In the $eM\!J$ channel $e^+$ and $e^-$ are possible in the final state, such that with respect to the incident lepton charge, a ``right'' (same sign) charge $eM\!J$(RC) and a ``wrong'' (opposite sign) charge $eM\!J$(WC) channel are distinguished. Channels with an electron and an additional charged lepton $\ell$ (where $\ell=e,\mu$) denoted by $ee M\!J$, $e\mu M\!J$ or neutrinos and a charged lepton $e \nu M\!J$, $\nu \mu M\!J$ (generically written as $e \ell M\!J$ and $\nu \ell M\!J$) necessarily involve two gauginos. Decay patterns involving more than two gauginos are kinematically suppressed and are therefore not explicitly studied here. Processes leading to final states with tau leptons are not expected to increase the sensitivity of this analysis and are not explicitly investigated. A dedicated search for isolated tau leptons in H1 shows good agreement with the SM~\cite{Aaron:2009wp}. \subsection{Event simulation} \label{sec:massgeneration}\label{sec:MC} For each of the signal topologies described above a dedicated Monte Carlo (MC) simulation is done. For the direct lepton--quark decay channels $eq$ and $\nu q$, as shown in figure~\ref{fig:directdecays}, events are generated using LEGO~\cite{LEGO}. For the gauge decays of squarks (figure~\ref{fig:gaugedecays}) events are generated using SUSYGEN3~\cite{SUSYGEN}. \par To allow a model independent interpretation of the results, the squark decay processes are simulated for a wide range of masses of the sparticles involved. The final states contain only SM fermions ($f$) considered as massless. The squark mass is varied from $100$~GeV to $290$~GeV in steps of typically $25$~GeV. For gauge decays of squarks involving a gaugino, which decays directly via \Rp, the process $\tilde{q} \rightarrow q \chi^0_1$ is generated for $\chi^0_1$ masses ranging between $30$~GeV and $M_{\tilde{q}}$. In order to study cascade gauge decays which involve two gauginos, the processes $\tilde{q} \rightarrow q \chi^+_1 \rightarrow q \chi^0_1 f \bar{f'}$ and $\tilde{q} \rightarrow q \chi^0_2 \rightarrow q \chi^0_1 f \bar{f'}$ are generated for $\chi^+_1$ and $\chi^0_2$ masses ranging between $40$~GeV and $M_{\tilde{q}}$, and for $\chi^0_1$ masses between $30$~GeV and $M_{\chi^+_1}$ or $M_{\chi^0_2}$. The masses of the gauginos are varied in steps of approximately $10$~GeV. The lower mass values for squarks and gauginos are motivated by the exclusion domains resulting from \Rp\ SUSY searches at LEP~\cite{LEP,L3RPV}. The mass intervals are sufficiently small to allow linear interpolation of signal detection efficiencies as a function of the masses of the sparticles involved. \par SM processes may mimic the characteristics of the final states of squark decays. This SM background is dominated by neutral current (NC) and charged current (CC) deep inelastic scattering (DIS), with additional small contributions from photoproduction, single $W$ boson production and lepton pair production. The RAPGAP~\cite{Jung:1993gf} event generator, which implements the Born level, QCD Compton and boson--gluon fusion matrix elements, is used to model inclusive NC DIS events. The QED radiative effects arising from real photon emission from both the incoming and outgoing electrons are simulated using the HERACLES~\cite{Kwiatkowski:1990es} program. Direct and resolved photoproduction of jets and prompt photon production are simulated using the PYTHIA~\cite{Sjostrand:2000wi} event generator. The simulation is based on Born level scattering matrix elements with radiative QED corrections. In RAPGAP and PYTHIA, jet production from higher order QCD radiation is simulated using leading logarithmic parton showers and hadronisation is modelled with Lund string fragmentation~\cite{lund}. Inclusive CC DIS events are simulated using the DJANGO~\cite{Schuler:yg} program, which includes first order leptonic QED radiative corrections based on HERACLES. The production of two or more jets in DJANGO is accounted for using the colour dipole model~\cite{Lonnblad:1992tz}. The leading order MC prediction of processes with two or more high transverse momentum jets in NC DIS, CC DIS and photoproduction is scaled by a factor of $1.2$ to account for the incomplete description of higher orders in the MC generators~\cite{Adloff:2002au,Aktas:2004pz}. Contributions arising from the production of single W bosons and multi--lepton events are modelled using the EPVEC~\cite{Baur:1991pp} and GRAPE~\cite{Abe:2000cv} event generators, respectively. \par Generated events are passed through a GEANT~\cite{Brun:1987ma} based simulation of the H1 apparatus, which takes into account the actual running conditions of the data taking. Simulated events are reconstructed and analysed using the same program chain as is used for the data. \section{Experimental Method} \subsection{H1 detector} A detailed description of the H1 experiment can be found elsewhere~\cite{Abt:h1}. Only the detector components relevant to this analysis are briefly described here. A right--handed cartesian coordinate system is used with the origin at the nominal primary $ep$ interaction vertex. The proton beam direction defines the positive $z$ axis (forward direction). The polar angle $\theta$ and the transverse momenta $P_T$ of all particles are defined with respect to this axis. The azimuthal angle $\phi$ defines the particle direction in the transverse plane. The pseudorapidity is defined as $\eta=-\ln {\tan {\frac{\theta}{2}}}$. The Liquid Argon (LAr) calorimeter~\cite{Andrieu:1993kh} covers the polar angle range \mbox{$4^\circ < \theta < 154^\circ$} with full azimuthal acceptance. The energies of electromagnetic showers are measured in the LAr with a precision of \mbox{$\sigma (E)/E \simeq 11\%/ \sqrt{E/\mbox{GeV}} \oplus 1\%$} and hadronic energy depositions with \mbox{$\sigma (E)/E \simeq 50\%/\sqrt{E/\mbox{GeV}} \oplus 2\%$}, as determined in test beam measurements~\cite{Andrieu:1993tz,Andrieu:1994yn}. A lead--scintillating fibre calorimeter (SpaCal)~\cite{SpaCal} covering the backward region \mbox{$153^\circ < \theta < 178^\circ$} completes the measurement of charged and neutral particles. For electrons a relative energy resolution of \mbox{$\sigma (E)/E \simeq 7\%/\sqrt{E/\mbox{GeV}} \oplus 1\%$} is reached, as determined in test beam measurements~\cite{SpaTestBeam}. The central \mbox{($20^\circ < \theta < 160^\circ$)} and forward \mbox{($7^\circ < \theta < 25^\circ$)} inner tracking detectors are used to measure charged particle trajectories and to reconstruct the interaction vertex. 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 \mbox{$\sigma_{P_T}/P_T = 0.005 P_T / \mbox{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 \mbox{($4^\circ < \theta < 171^\circ$)}. In the very forward region of the detector \mbox{($3^\circ < \theta < 17^\circ$)} a set of drift chambers detects muons and measures their momenta using an iron toroidal magnet. 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. \subsection{Particle identification and event reconstruction} Electromagnetic particle (electron and photon) candidates are identified as compact and isolated clusters of energy in the electromagnetic part of the LAr calorimeter. Electron candidates are identified as electromagnetic particle candidates with an associated track. Identification of muon candidates is based on a track in the inner tracking detectors, associated to a signal in the muon system. Tracks and calorimeter deposits not identified as originating from isolated electromagnetic particles or muons are combined into cluster--track objects to reconstruct the hadronic final state~\cite{benji}. Jets are reconstructed from these objects using an inclusive $k_T$ algorithm~\cite{kT1,kT2} with a minimum $P_T$ of $4$~GeV. The missing transverse momentum \mbox{${P}_T^{\rm{miss}}$}, which may indicate the presence of neutrinos in the final state, is derived from all reconstructed particles in the event. A neutrino four--vector $P_{\nu}$ is reconstructed by exploiting momentum and energy conservation. The transverse momentum of the neutrino is reconstructed by assuming one neutrino with significant energy in the event $ \vec{P}_T^{\nu} \equiv \vec{P}_T^{\rm{miss}} $. The energy of the neutrino is then reconstructed exploiting the energy and longitudinal momentum balance: $ \sum_{i} \left( E^i-P^i_z \right) + \left( E^{\nu} - P^{\nu}_z \right) = 2E^{0}_e = 55.2\,\mbox{GeV}, $ where the sum runs over all detected particles, $P_z$ is the momentum along the proton beam axis and $E_e^0$ denotes the energy of the incident electron. For further selection the following observables are used which in SM DIS events correspond to the Lorentz--invariant quantities: inelasticity $y$, negative four--momentum transfer squared $Q^2$ and Bjorken's scaling variable $x$. They can be reconstructed as: $$y_e = 1 - \frac{E_e (1 - \cos \theta_e) }{2E_e^0}, \;\;\;Q^2_e = \frac{P^2_{T,e}}{1-y_e},\;\;\;x_e = \frac{Q^2_e} {y_e s},$$ where the polar angle $\theta_e$, energy $E_e$ and transverse momentum $P_{T}^{e}$ of the electron with the highest $P_T$ found in the event are used. If no electron is reconstructed in the event similar quantities can be calculated using the Jacquet--Blondel method~\cite{BLONDEL} from the hadronic final state variables: $$y_h=\frac{\sum \left(E-P_z\right)_h}{2E_e^0}, \;\;\;Q^2_h= \frac{P^2_{T,h}}{1-y_h},\;\;\;x_h = \frac{Q^2_h} {y_h s},$$ where $P_{T,h}$ is the transverse momentum of the hadronic final system calculated from the before mentioned cluster--track objects. The sum ${\sum \left(E-P_z\right)_h}$ runs over energies and momenta of jets only. \subsection{Trigger and data quality} All data events used for this search are triggered by the LAr calorimeter~\cite{Adloff:2003uh}. Events with an electromagnetic deposit in the LAr with an energy greater than $10$~GeV are detected with an efficiency close to $100\%$~\cite{nikiforov}. Events are also triggered by hadronic jets, with a trigger efficiency above $95\%$ for a jet transverse momentum $P_T^{\rm{jet}} > 20$~GeV and almost $100\%$ for $P_T^{\rm{jet}} > 25$~GeV~\cite{matti}. For events with missing transverse energy of $20$~GeV, the trigger efficiency is about $90\%$ and increases to above $95\%$ for missing transverse energy above $30$~GeV~\cite{trinh}. \par In order to remove background events induced by cosmic showers and other non--$ep$ sources, the event vertex is required to be within $35$~cm in $z$ of the mean position for $ep$ collisions. In addition, topological filters and timing vetoes are applied~\cite{negri} in order to reject beam related and cosmic background. \section{Data Analysis} The event selection is carried out in several exclusive analysis channels. The resulting event rates for each channel and the range of efficiencies for the selection of signal events are given in table~\ref{tab:totnum}. \subsection{Electron--jet final state \boldmath{$eq$}} The final state of a squark decaying into an electron and a high $P_T$ jet is identical to the NC DIS signature at high $x$ and $Q^2$. For the signal the reconstructed invariant mass distribution \mbox{$M_e= \sqrt{x_e s}$} shows a resonance peak at the nominal squark mass with a resolution \mbox{$\delta M_e = 4-10$~GeV} depending on the mass of the squark. Differences in the $M_e$ and $y_e$ distributions of the two processes allow to discriminate them statistically. Squarks produced in the $s$--channel decay isotropically leading to a flat $d\sigma/dy$ distribution, whereas for NC DIS events a distribution proportional to $1/y^{2}$ is expected. \par Events are selected by requiring \mbox{${P}_T^{\rm{miss}}<15$~GeV} and \mbox{$40$~GeV $< \sum (E-P_z) < 70$~GeV} where the sum runs over all particles in the final state. An isolated electron with \mbox{$P_{T}^{e}>16$~GeV} in the region \mbox{$5^{\circ}<\theta_e<145^{\circ}$}, \mbox{$Q^2_e>2500$~GeV$^2$} and \mbox{$y_e<0.9$} is required. The high $y$ region is excluded to reduce background contributions arising from photoproduction events. An $M_e$ dependent cut on $y_e$ is determined by minimising the expected limit from signal and SM MC. The cut ranges from \mbox{$y_e>0.5$} for masses around $100$~GeV to \mbox{$y_e>0.2$} for masses around $290$~GeV. In order to remain exclusive with respect to other channels, events with muon candidates with \mbox{$P_T^{\mu}>5$~GeV} or two jets with \mbox{$P_{T}^{\rm{jet}}>15$~GeV} are rejected. \par The $M_e$ spectra after this selection for data and SM background are shown in figures~\ref{fig:masselec}a and~\ref{fig:massposi}a for $e^-p$ data and $e^{+}p$ data, respectively. For both data samples no significant deviation from the SM expectation is observed. In the $e^{-}p$ sample, $3121$ events are observed while the SM expectation yields \mbox{$3215\pm 336$}. In the $e^+p$ data sample $2946$ candidate events are found compared to \mbox{$2899\pm 302$} expected from SM processes. \subsection{Neutrino--jet final state \boldmath{$\nu q$}} Squarks decaying into a neutrino and a high $P_T$ jet lead to the same signature as CC DIS events with high missing transverse momentum. Similarly to the $eq$ channel, the resonant $s$--channel production and isotropic decay allows a statistical separation of signal and background. The resolution $\delta M_h$ of the reconstructed invariant mass $M_h= \sqrt{x_h s}$ varies between \mbox{$12$ and $22$~GeV} depending on the mass of the squark considered. \par The presence of a neutrino in the event is required by imposing \mbox{${P}_T^{\rm{miss}}>30$~GeV} and\linebreak \mbox{$\sum (E-P_z) < 50$~GeV}. The phase space is restricted to $Q^2_h>2500$~GeV$^2$ and $y_h<0.9$. Similar to the $eq$ channel, a cut on $y_h$ dependent on the reconstructed mass $M_h$ is applied. The cut ranges from \mbox{$y_h>0.3$} for masses around $100$~GeV to \mbox{$y_h>0.1$} for masses around $290$~GeV. In order to remain exclusive with respect to other channels, events with any electron or muon candidate with \mbox{$P_T^{e, \mu}>5$~GeV} or events containing two jets with \mbox{$P_{T}^{\rm{jet}}>15$~GeV} are rejected. \par Only $\tilde{d}_R$--type squarks, which are produced dominantly in $e^{-}p$ collisions, can undergo direct decay leading to a $\nu q$ final state. For selected events the $M_h$ spectrum of this data set and of the simulation of SM background events is shown in figure~\ref{fig:masselec}b. No significant deviation from the SM expectation is found. In data $2858$ events are observed while \mbox{$2983\pm 358$} are expected from SM processes. \subsection{Electron--multijet and electron--lepton--multijet final states} \subsubsection*{Common preselection for \boldmath{$eM\!J(RC),eM\!J(WC),eeM\!J,e\mu M\!J$ and $e\nu M\!J$}} Squarks decaying via neutralinos or charginos are expected to have a higher multiplicity of jets and leptons in the final state. The signatures correspond to final states detectable in higher order NC DIS processes. However, as heavy particles are boosted forward, in events with squarks the decay products are mainly emitted into the forward part of the detector. This feature is used to distinguish between the signal and SM background. \par A common preselection is applied for the $eM\!J(RC),eM\!J(WC),eeM\!J,e\mu M\!J$ and $e \nu M\!J$ channels: At least one isolated electron with \mbox{$P_{T}^{e}>6$~GeV} and \mbox{$E_e>11$~GeV} in the region \mbox{$5^{\circ}<\theta_e<110^{\circ}$} is required. In addition, the condition $y_e>0.3$ is used to reduce the background from NC DIS. Central electrons with \mbox{$30^\circ<\theta_e<110^\circ$} are required to have a well measured track associated to the cluster and the distance of closest approach between the track impact point and the centre--of--gravity of the cluster should not exceed $12$~cm. Furthermore the energy of the track is required to match the energy of the associated cluster according to \mbox{$E_{{\rm{cluster}}}/E_{{\rm{track}}}>0.5$}. At least two jets with \mbox{$P_{T}^{\rm{jet}}>15$~GeV} in the jet polar angle range \mbox{$7^{\circ}<\theta_{\rm{jet}}<145^{\circ}$} are also required. \par By requiring \mbox{$Q^2_e> 1000$~GeV$^2$} the steep decrease of the NC DIS cross section with increasing $Q^2$ is exploited. This corresponds to an implicit upper cut on the electron polar angle as $Q^2$ is strongly correlated with the polar angle of the scattered electron. For signal events at least one high $P_T$ particle is expected to be emitted in the forward direction, therefore the highest $P_T$ electron or one of the two highest $P_T$ jets has to fulfill \mbox{$\theta_{e,{\rm{jet}}}<40^{\circ}$}. Moreover, of the two highest $P_T$ jets, the one with the largest polar angle, $\theta_{\rm{backw}}$, must satisfy the condition \mbox{$\theta_{{\rm{backw}}} < (y_e - 0.3) \cdot 180^{\circ} $}, separating efficiently signal events from NC DIS background~\cite{JOHTHESIS}. For selected events an invariant mass $M_{\rm{rec}}$ is calculated as \mbox{$M_{\rm{rec}}= \sqrt{4E^0_e (\sum{E_i}-E^0_e)}$}, where the energies $E_i$ of electrons, muons and jets found in the event with \mbox{$P_T^{\rm{jet}}>5$~GeV} are included in the sum. The resolution $\delta M_{\rm{rec}}$ of this method ranges between \mbox{$6$ and $10$~GeV} depending on the mass of the squark considered. \subsubsection*{Electron--multijet final state \boldmath{$eM\!J(RC),eM\!J(WC)$}} Decays of squarks via gauginos are likely to produce a single isolated electron and multiple jets in the final state. Events produced in $e^{\pm}p$ collisions with exactly one electron and multiple jets in the final state are also expected in the SM, where in general the measured charge of the electron corresponds to that of the incident electron. A selection channel labelled ``right'' (same sign) charge \mbox{$eM\!J $ (RC)} is used for events fulfilling this criterion. A selection channel where the electron charge is identified as opposite to the incident electron, denoted ``wrong'' (opposite sign) charge \mbox{$eM\!J $ (WC)}, represents therefore a powerful test of the SM and is expected to be essentially background free. The distinction between RC and WC $eM\!J $ events is based on the curvature of the electron track measured in the central tracking system. Events are allocated to the WC channel if the electron is found in the central region \mbox{$30^{\circ} <\theta_e<110^{\circ}$} and its charge is measured to be opposite to that of the incident electron, with a charge significance greater than two standard deviations~\cite{JOHTHESIS}. Otherwise the event is assigned to the RC channel. Events are selected from the common preselection described above by requiring in addition \mbox{$P_{T}^{\rm{miss}} < 15$~GeV} and \mbox{$40$~GeV $ < \sum (E - P_{z}) < 70$~GeV} since no neutrinos are expected in these channels. To ensure that the selection channels are exclusive, no additional electron or muon candidate with \mbox{$P_{T}^{e, \mu} > 5$~GeV} may be present in the event. \par In the \mbox{$eM\!J $ (WC)} channel no event is observed in the $e^-p$ data while \mbox{$1.3\pm 0.3$} are expected from SM processes and one candidate event is observed in the $e^+p$ data while \mbox{$0.6\pm 0.4$} are expected. For $e^+p$ collisions the mass spectrum is shown in figure~\ref{fig:massposi}b. \par For events in the \mbox{$eM\!J $ (RC)} channel an $M_{\rm{rec}}$ dependent cut on $y_e$ is applied to increase the sensitivity for signal events. The cut ranges from \mbox{$y_e>0.7$} for masses around $100$~GeV to \mbox{$y_e>0.5$} for masses around $290$~GeV. The $M_{\rm{rec}}$ distributions for data and simulation are shown in figures~\ref{fig:masselec}c and~\ref{fig:massposi}c for $e^-p$ and $e^{+}p$ collisions, respectively. No significant deviation from the SM expectation is observed. In total $147$ events are observed in the $e^-p$ data while the SM simulation yields \mbox{$158.3\pm 23.9$} and in the $e^+p$ data $140$ events are observed for an expectation of \mbox{$146.0\pm 21.4$} from SM processes. \subsubsection*{Electron--lepton--multijet final states, \boldmath{$eeM\!J,e\mu M\!J,e \nu M\!J$}} Final states from squark decays may contain more than one isolated lepton if the decays proceed via cascades of gauginos. In addition to the common preselection the $eeM\!J$ and $e\mu M\!J$ channels require either an additional electron with the same criteria as described in the common preselection or an isolated muon with \mbox{$P_{T}^{\mu}>5$~GeV} in the polar angle range \mbox{$10^{\circ}<\theta_{\mu}<110^{\circ}$}. After applying this selection the SM background expectation is very low in these channels. For the \mbox{$e\mu M\!J$} channel there are no candidate events observed in the $e^-p$ and $e^+p$ collision data for a SM expectation of \mbox{$0.03\pm 0.02$} and \mbox{$0.03\pm 0.03$}, respectively. In the $e^-p$ collision data no candidate event for a SM expectation of \mbox{$1.5\pm 0.5$} is observed in the $eeM\!J$ channel and two events in the $e^+p$ data are observed compared to a SM expectation of \mbox{$1.7\pm 0.5$}. The mass spectrum for the $eeM\!J$ channel in the $e^+p$ data is shown in figure~\ref{fig:massposi}d. \par In the $ e\nu M\!J$ channel a neutrino is expected in the final state, therefore in addition to the common preselection, large missing transverse momentum \mbox{${P}_T^{\rm{miss}}>15$~GeV} is required. Due to the presence of the neutrino \mbox{$\sum (E-P_z)$} is significantly reduced causing $y_h$ to be substantially smaller than $y_e$, while in NC DIS events $y_h\approx y_e$ is expected. Thus a cut \mbox{$y_e(y_e-y_h)>0.04$} is used to discriminate the SUSY signal from background events~\cite{JOHTHESIS}. Exclusivity with respect to the $eeM\!J$ and $e\mu M\!J$ channels is achieved by rejecting events containing an additional electron or muon with \mbox{$P_T^{e, \mu}>5$~GeV}. The method used to reconstruct a mass $M_{\rm{rec}, \nu}$ for selected events taking the energy of the neutrino into account is explained in section~\ref{sec:vmj}. In the $e^-p$ collision data three events are observed in the $ e\nu M\!J$ channel while \mbox{$5.6\pm 1.2$} are expected and in $e^+p$ collision data five events are observed while \mbox{$8.2\pm 2.0$} are expected from SM processes. The mass spectra are shown in figures~\ref{fig:masselec}d and~\ref{fig:massposi}e for the $e^-p$ data and $e^{+}p$ data, respectively. \par \subsection{Neutrino--multijet and neutrino--muon--multijet final states} \label{sec:vmj} \subsubsection*{Common preselection for \boldmath{$\nu M\!J$} and \boldmath{$\nu \mu M\!J$}} Squark decays with single or multiple neutrinos produced via neutralino or chargino decays can result in final states similar to that of higher order CC DIS processes. A substantial missing transverse momentum \mbox{${P}_T^{\rm{miss}}>26$~GeV} is required and at least two jets must be found with \mbox{$P_{T}^{\rm{jet}}>15$~GeV} in the range \mbox{$7^{\circ}<\theta_{\rm{jet}}<145^{\circ}$}. No electron candidate with \mbox{$P_T^{e}>5$~GeV} is allowed. A cut \mbox{$\sum (E - P_{z})<50$~GeV} is used to ensure the neutrino energy is positive. For each selected event a squark mass $M_{\rm{rec},\nu}$ is calculated as \mbox{$M_{\rm{rec}, \nu}= \sqrt{4E^0_e (\sum{E_i}-E^0_e)}$} where the sum includes the energies of the reconstructed neutrino, electrons, muons and jets with \mbox{$P_T^{\rm{jet}}>5$~GeV} in the event. This method assumes that all missing energy is carried by a single neutrino and yields a resolution $\delta M_{\rm{rec}, \nu}$ of about \mbox{$15$ to $20$~GeV} depending on the squark mass. \subsubsection*{Neutrino--multijet final state \boldmath{$\nu M\!J$}} Squark decays via gauginos are likely to produce final states with multiple jets and a single neutrino. Events are selected in the $\nu M\!J$ channel if no muon candidate is found. \par A cut on $y_h$ dependent on the reconstructed mass $M_{\rm{rec}, \nu}$ is applied to enhance the signal. The cut ranges from \mbox{$y_h>0.5$} for masses around $100$~GeV to \mbox{$y_h>0.4$} for masses around $290$~GeV. The $M_{\rm{rec},\nu}$ spectra are shown in figures~\ref{fig:masselec}e and~\ref{fig:massposi}f for $e^-p$ and $e^{+}p$ collision data and SM background simulation. In the $e^-p$ data $204$ candidate events are selected while \mbox{$235.5\pm 63.3$} are expected and in the $e^+p$ data $113$ candidate events are selected while \mbox{$134.0 \pm 33.8$} are expected from SM processes. \subsubsection*{Neutrino--muon--multijet final state \boldmath{$\nu \mu M\!J$}} If an isolated muon with \mbox{$P_{T}^{\mu}>5$~GeV} in the polar angle range \mbox{$10^{\circ}<\theta_{\mu}<110^{\circ}$} is found in an event in the common $\nu M\!J$ selection, the event is classified as $\nu \mu M\!J$ candidate. No candidate events are found, in the $e^-p$ or in the $e^+p$ collision data, in agreement with the SM expectations of \mbox{$0.04\pm 0.02$} and \mbox{$0.06\pm 0.03$}, respectively. \subsection{Systematic uncertainties} 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 central region to $2\%$ in the forward region~\cite{trinh}. The polar angle measurement uncertainty of electromagnetic clusters is $3$~mrad. \item The jet energy scale is known within $2\%$~\cite{trinh}. The uncertainty on the jet polar angle determination is $10$ mrad. \item The luminosity measurement has an uncertainty of $3\%$. \end{itemize} The effects of the experimental systematic uncertainties on the SM expectation and signal efficiencies are determined by varying the corresponding experimental quantities within one standard deviation in the MC samples and propagating the variations to the final distributions. The resulting experimental uncertainties are determined for each analysis channel individually and added in quadrature. In the $eq$ channel the uncertainty on the overall SM event yield was found to be $3\%$, while in the ${\nu}q$ channel an uncertainty of $7\%$ is determined. In the ${e}M\!J$ (RC and WC) channels the resulting systematic uncertainty amounts to $4\%$, in the ${e}{\mu}M\!J$ and ${e\nu}M\!J$ channels to $7\%$, in the ${ee}M\!J$ channel to $10\%$ and in the ${\nu}M\!J$ and ${\nu}{\mu}M\!J$ channels to $20\%$. Additional model uncertainties are attributed to the SM MC event generators described in section~$\ref{sec:MC}$. A conservative error of $10\%$ is attributed to NC (RAPGAP) and CC (DJANGO) DIS processes with only one high $P_T$ jet. To account for the uncertainty on higher order QCD corrections, an uncertainty of $15\%$ is attributed to NC DIS and photoproduction processes (PYTHIA) with at least two high $P_T$ jets. The normalisation uncertainty of CC DIS processes with at least two high $P_T$ jets is estimated to be $20\%$~\cite{trinh}. A $5\%$ uncertainty is attributed to the contribution from multi--lepton events (GRAPE) and a $15\%$ uncertainty on the production of single W bosons (EPVEC). These uncertainties include contributions from the proton parton distribution functions and from missing higher order QCD corrections. % The total error on the SM prediction is determined by adding the effects of all model and experimental systematic uncertainties in quadrature. For the signal cross section further uncertainties arise from the determination of signal efficiencies ($10\%$ due to available MC statistics), the theoretical uncertainty on the squark production cross section ($7\%$ for low squark masses, up to $50\%$ for the highest masses from the PDF uncertainty) and an uncertainty due to the scale at which the PDFs are evaluated ($7\%$)~\cite{JOHTHESIS}. \section{Exclusion Limits} No significant deviation from the SM expectation is observed in any channel. Consequently the observations in all analysis channels are combined to set constraints on various supersymmetric models. Exclusion limits are obtained on the production of squarks parameterised by the strength of the \Rp\ couplings $\lambda'_{1j1}$ and $\lambda'_{11k}$ and dependent on the mass of the squark. \subsection{Procedure} For the interpretation of the results a version of the Minimal Supersymmetric Standard Model (MSSM) is considered where the masses of the neutralinos, charginos and gluinos are determined via the usual parameters: the ``Higgs--mass'' term $\mu$, which mixes the Higgs superfields; the SUSY soft--breaking mass parameter $M_2$; and the ratio of the vacuum expectation values of the two neutral scalar Higgs fields $\tan \beta$~\cite{MSSM}. The parameters are defined at the electroweak scale. A set of parameters $(\tan \beta,\mu,M_2)$ together with the sfermion masses and a coupling $\lambda'_{1jk}$ define a supersymmetric scenario where the masses of the gauginos and the branching ratios for squark decays into the different final state topologies are fixed and can be obtained using the SUSYGEN3~\cite{SUSYGEN} package. The branching ratios for the specific parameters of the model are taken into account in the combination~\cite{Junk:1999kv}. A sliding mass window technique is used in channels with high contributions of irreducible SM background ($eq$, $\nu q$, $eM\!J$ and $\nu M\!J$) to improve the signal to background ratio for the squark mass examined. The width of the mass window is determined by minimising the expected limit, and increases towards high squark masses, reflecting the corresponding mass reconstruction resolution. The small efficiency losses due to the finite mass window width are taken into account. A $95\%$ confidence level (CL) upper limit $\sigma_{\rm lim}$ on the squark production cross section compatible with the simultaneous observation in all channels is derived using a modified frequentist approach based on Likelihood ratios\cite{Junk:1999kv}. Sets of model parameters leading to signal cross sections above $\sigma_{\rm lim}$ are excluded. If the squark width is non--negligible, in particular for squark masses approaching the kinematic limit, the production cross section decreases at the resonance peak and contributions from the lower tail of the squark mass distribution become important, enhanced by the rapid increase of proton parton distributions at low Bjorken--$x$~\cite{JOHTHESIS}. This is taken into account by generating events for negligible squark widths to determine signal efficiencies at all masses. The selection efficiencies are then corrected for the actual squark width by reducing the efficiency for the signal selection accordingly~\cite{JOHTHESIS}. In the special case of stop and sbottom squark production, namely via \Rp\ couplings $\lambda'_{131}$ and $\lambda'_{113}$, the mixing of the weak eigenstates \mbox{$\tilde{t}_L,\tilde{t}_R$} (and \mbox{$\tilde{b}_L,\tilde{b}_R$}) to the mass eigenstates \mbox{$\tilde{t}_1,\tilde{t}_2$ ($\tilde{b}_1,\tilde{b}_2$)} via an angle $\theta_{\tilde{t}}$ ($\theta_{\tilde{b}}$) becomes important for the calculation of branching ratios and production cross sections. The gauge decay via a top quark would lead to decay products different from the first two generation squarks, for which the efficiencies are determined. Final states with top quarks are not considered explicitly. Since top signals would in any case be present in one of the selection topologies, this approach is conservative. \subsection{Constraints on a phenomenological MSSM} Constraints are set in a scenario of a phenomenological MSSM~\cite{MSSM} where the lightest supersymmetric particle (LSP) is the neutralino $\chi^0_1$. Slepton masses $M_{\tilde{l}}$ are fixed at $90$~GeV, close to the lowest mass bound from \Rp\ sfermion searches at LEP~\cite{LEP} and squark masses are treated as free parameters. For higher slepton masses only very small degradations in the derived constraints are expected~\cite{Aktas:2004ij}. For a single point in the parameter space, characterised by $\mu=-200$~GeV, $M_2=80$~GeV and $\tan\beta=2$, constraints on the strength of the \Rp\ couplings depending on the mass of the squark are derived for $\tilde{d}^k_R\,(k=1,2)$ (figure~\ref{fig:photinolikeelec}a) and $\tilde{u}^j_L\,(j=1,2)$ (figure~\ref{fig:photinolikeelec}b) production. The HERA sensitivity allows tests of $\lambda'$ values as low as $10^{-2}$ for squark masses of $100$~GeV. For high squark masses the sensitivity degrades since the production cross section decreases strongly. The limits from the previous H1 analysis on a smaller data sample~\cite{Aktas:2004ij} are also indicated. This choice of parameters leads to a dominant photino~($\tilde{\gamma}$) component to the neutralino's composition. As a consequence, gauge decays are likely to result in charged leptons in the final state. The branching ratios into the decay topologies are shown at the observed limit for $\tilde{d}^k_R\,(k=1,2)$ (figure~\ref{fig:photinolikeelec}c) and $\tilde{u}^j_L\,(j=1,2)$ (figure~\ref{fig:photinolikeelec}d) production. For $\tilde{d}^k_R\,(k=1,2)$ production the channels $eM\!J$ (RC) and $eM\!J$ (WC) each contribute about $40\%$ over a wide range of squark masses and only $10\%$ of squark decays appear in the $\nu M\!J$ channel. For squark masses approaching the kinematic limit of the centre--of--mass energy the lepton--quark channels $eq$ and $\nu q$ begin to dominate the decays of the squarks, because gauge decay modes become negligible at high values of the \Rp\ couplings. Over the whole mass range the sum of analysed branching ratios is close to $100\%$. For $\tilde{u}^j_L\,(j=1,2)$ production the channels $e \ell M\!J$ have the highest branching ratio over a wide mass range at the observed limit. \par A different point in the parameter space yields a complementary scenario with the choice of $\mu=200$~GeV, $M_2=150$~GeV and $\tan\beta=2$. Again the neutralino $\chi^0_1$ is the LSP but its composition is now dominated by a zino ($\tilde{Z}$) component. Squark decays are now more likely to produce neutrinos in the final state. The constraints on the couplings depending on the mass of the squark are shown for $\tilde{d}^k_R\,(k=1,2)$ (figure~\ref{fig:zinolikeelec}a) and $\tilde{u}^j_L\,(j=1,2)$ (figure~\ref{fig:zinolikeelec}b) production and are of the same order of magnitude as in the photino scenario. Branching ratios at the observed limits show dominant contributions from the $\nu M\!J$ and $\nu \ell M\!J$ channels~(figures~\ref{fig:zinolikeelec}c and~\ref{fig:zinolikeelec}d). \par These two scenarios illustrate the sensitivity for various model configurations achieved by the combination of the complementary search topologies. The sensitivity of the analysis is explored in a scan of the MSSM parameters. The parameters $M_2$ and $\mu$ are varied in the range $70$~GeV~$< M_2 < 350$~GeV and $-300$~GeV~$ < \mu < 300$~GeV for $\tan \beta = 6$. Parameter sets leading to a scalar LSP or to LSP masses below $30$~GeV are not considered. The latter restriction, as well as the lower boundary of the $M_2$ range, are motivated by the exclusion domains resulting from gaugino searches in \Rp\ SUSY at LEP~\cite{L3RPV}. Figures~\ref{fig:scan11k} and \ref{fig:scan1j1} show the resulting constraints on the couplings as a function of the squark mass. The region of values excluded at $95\%$~CL for the couplings in all scenarios and the best exclusion limit achieved in all scenarios are indicated for first and second generation squarks $\tilde{d}_R$, $\tilde{s}_R$ (figure~\ref{fig:scan11k}a) and $\tilde{u}_L$, $\tilde{c}_L$ (figure~\ref{fig:scan1j1}a) as well as for third generation squarks $\tilde{b}_R$ (figure~\ref{fig:scan11k}b) and $\tilde{t}_L$ (figure~\ref{fig:scan1j1}b). The resulting exclusion domains are compared to the previous H1 results~\cite{Aktas:2004ij}. Constraints on the \Rp\ couplings are also available as indirect limits from low energy experiments probing virtual squark contributions~\cite{barbier}. The production of up--type and down--type squarks via the $\lambda'_{111}$ coupling is strongly constrained by the non--observation of neutrinoless double beta decay ($\beta \beta 0 \nu$)~\cite{BETA0NU,barbier}. The best indirect limit on the couplings $\lambda'_{112}$ and $\lambda'_{113}$ results from tests of charged current universality (CCU)~\cite{CCU,barbier} and can be compared to the direct limits obtained in this analysis for $\lambda'_{11k}$ in figure~\ref{fig:scan11k}. The best indirect limit on the couplings $\lambda'_{121}$ and $\lambda'_{131}$ comes from atomic parity violation (APV) measurements~\cite{APV,barbier} and can be compared to the direct limits obtained for $\lambda'_{1j1}$ in figure~\ref{fig:scan1j1}. In the part of the parameter space considered here, Yukawa couplings of electromagnetic strength $\lambda'_{1j1}$ or $\lambda'_{11k}=\sqrt{4\pi\alpha_{\rm em}}= 0.3$, are excluded up to masses of $275$~GeV at $95\%$~CL for up--type squarks and up to masses of $290$~GeV for down--type squarks. \subsection{Constraints on the Minimal Supergravity Model} Constraints are also obtained on the Minimal Supergravity Model (mSUGRA)~\cite{MSUGRA} which is a complete SUSY model using the assumption of gauge coupling unification and radiative electroweak symmetry breaking (REWSB) with the choice of 5 parameters: the common mass of scalar sparticles $m_0$; the common mass of fermionic sparticles $m_{1/2}$; the common trilinear coupling $A_0$; the ratio of Higgs vacuum expectation values $\tan \beta$; and the sign of the Higgs mixing parameter $\mu$. % The masses of squarks, sleptons and gauginos as well as the branching ratios in the analysis channels are determined by the set ($m_0,m_{1/2},\tan\beta,sign(\mu),A_0$) for given values of the couplings $\lambda'_{11k}$ and $\lambda'_{1j1}$. The program SUSPECT~2.1~\cite{SUSPECT} is used to obtain the REWSB solution for $|\mu|$ and to calculate the full supersymmetric mass spectrum. $A_0$ enters only marginally in the interpretation and is set to zero. The parameter $\mu$ is taken with negative sign. Figures~\ref{fig:msugra11k} and~\ref{fig:msugra1j1} show constraints in the \mbox{$m_0,m_{1/2}$} plane when values of the couplings are assumed to be of the electromagnetic coupling strength $\lambda'_{11k}=0.3$ or $\lambda'_{1j1}=0.3$ for different values of $\tan \beta$. The excluded region typically covers masses of $m(\tilde{u})=275$~GeV, $m(\tilde{t})=270$~GeV and $m(\tilde{d})=280$~GeV, as indicated in the figures. Complementary constraints are obtained by the L3 experiment~\cite{L3RPV} at LEP and the \mbox{D\O} experiment~\cite{D0RES} at the Tevatron which exploit di--electron events. The LEP and Tevatron limits are independent of the Yukawa coupling. For $\tan \beta =2$, the parameter space is more strongly constrained by the searches for gauginos and sleptons at the L3 experiment at LEP, as shown in figures~\ref{fig:msugra11k} and \ref{fig:msugra1j1}. This is the only $\tan \beta$ value considered in the L3 analysis, although results for higher values are expected to be similar~\cite{L3RPV}. Compared to the \mbox{D\O} experiment, the H1 limits are more stringent only for low values of $m_0$ for $\tan \beta =2$ , whereas for $\tan \beta = 6$ the domain excluded by H1 is considerably larger. The exclusion limits in figure~\ref{fig:msugra11k} are very similar for all three flavours of down--type squarks. Significant differences are observed between the first two and the third generation of up--type squarks. The stronger limit for stop squark production results from strong mixing effects that occur for third generation squarks with increasing $\tan \beta$, leading to masses for stop squarks lower than for first and second generation up--type squarks. The $\tan \beta$ dependence of the mSUGRA exclusion limits is studied assuming a unified common mass $M=m_0=m_{1/2}$ and \Rp\ couplings of electromagnetic coupling strength. This is illustrated in figure~\ref{fig:tanbeta} for the individual flavours. For the first two generations of up--type and down--type squarks no dependence on $\tan \beta$ is observed and values of $M<105$~GeV and $M<110$~GeV, respectively, are excluded over the whole range. In the case of stop squark production, significantly higher values (up to $M<148$~GeV) are excluded due to the presence of a light stop squark state. This effect is also observed for sbottom production, where increased values of $\tan \beta$ allow higher values of $M$ to be excluded. While there is a steep increase for the limit from stop squark production at small $\tan \beta$ and a flat plateau over the remaining $\tan \beta$ range, the limit from sbottom production increases steadily over the range of $\tan \beta$ values. The sharp edge in the stop exclusion curve at $\tan\beta\approx38$ for stop production follows from mixing effects in the $\tilde {\tau}$ sector at high $\tan\beta$ leading to scenarios with strong contributions from events with $\tau$ leptons in the final state, which are not explicitely considered in the analysis. \section{Summary} A search for $R$--parity violating production of squarks in $255\,\mathrm{pb^{-1}}$ of $e^+ p$ and $183\,\mathrm{pb^{-1}}$ of $e^- p$ collisions at HERA is presented. No significant deviation from the Standard Model is observed in the study of final state topologies which may result from direct or indirect \Rp\ squark decays. Mass dependent limits on the \Rp\ couplings $\lambda'_{1j1}$ and $\lambda'_{11k}$ ($j,k=1,2,3$) are derived within a phenomenological version of the MSSM. The existence of $\tilde{u}_L$--type and $\tilde{d}_R$--type squarks of all three generations with masses up to $275$~GeV and $290$~GeV, respectively, is excluded at the $95\%$ CL for Yukawa couplings of electromagnetic strength. These mass limits set the most stringent direct bounds on $\lambda'_{1j1}$ and $\lambda'_{11k}$. For lower squark masses, the results improve the indirect bounds set by low--energy experiments. Exclusion limits are also derived in the mSUGRA model, and are competitive with and complementary to those derived at the LEP and Tevatron colliders. \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 support and for the hospitality which they extend to the non--DESY members of the collaboration. \vfill % \input{bib.sorted} % This file is automatically generated with Ytsen's script "sort.cites.pl" % DO NOT EDIT THIS FILE! % New bib items go in the unsorted bibliography file (bib.tex for isolepdraft) \begin{thebibliography}{99} \bibitem{Aktas:2004ij} A.~Aktas {\it et al.} [H1 Collaboration], % ``Search for squark production in R parity violating supersymmetry at HERA,'' Eur.\ Phys.\ J.\ {\bf C36} (2004) 425 [hep-ex/0403027]. %%CITATION = EPHJA,C36,425;%% \bibitem{PEREZ} C.~Adloff {\it et al.} [H1 Collaboration], Eur.\ Phys.\ J.\ {\bf C20} (2001) 639 [hep-ex/0102050]. \bibitem{LEP} A.~Heister {\it et al.} [ALEPH Collaboration], % ``Search for Supersymmetric Particles with R-Parity Violating Decays in e+e- Collisions at sqrt{s} up to 209GeV,'' Eur.\ Phys.\ J.\ {\bf C31} (2003) 1 [hep-ex/0210014]; \\ A.~Heister {\it et al.} [ALEPH Collaboration], % ``Search for R-Parity Violating Production of Single Sneutrinos in e+e- Collisions at sqrt(s) = 189-209GeV,'' Eur.\ Phys.\ J.\ {\bf C25} (2002) 1 [hep-ex/0201013]; \\ P.~Abreu {\it et al.} [DELPHI Collaboration], % ``Search for Spontaneous R-parity violation at sqrt(s) = 183 GeV and 189 GeV,'' Phys.\ Lett.\ {\bf B502} (2001) 24 [hep-ex/0102045]; \\ P.~Abreu {\it et al.} [DELPHI Collaboration], % ``Search for R-parity violation with a Ubar Dbar Dbar coupling at sqrt(s)=189 GeV,'' Phys.\ Lett.\ {\bf B500} (2001) 22 [hep-ex/0103015]; \\ P.~Abreu {\it et al.} [DELPHI Collaboration], % ``Search for SUSY with R-parity violating $LL\bar{E}$ couplings at sqrt(s) = 189 GeV,'' Phys.\ Lett.\ {\bf B487} (2000) 36 [hep-ex/0103006]; \\ P.~Achard {\it et al.} [L3 Collaboration], % ``Search for R-parity Violating Decays of Supersymmetric Particles in e+e- Collisions at LEP,'' Phys.\ Lett.\ {\bf B524} (2002) 65 [hep-ex/0110057]; \\ G.~Abbiendi {\it et al.} [OPAL Collaboration], % ``Search for R-Parity Violating Decays of Scalar Fermions at LEP,'' Eur.\ Phys.\ J.\ {\bf C33} (2004) 194 [hep-ex/0310054]; \\ % CERN-EP-2003-036, accepted by Eur.\ Phys.\ J.\ {\bf C} [hep-ex/0310054];\\ G.~Abbiendi {\it et al.} [OPAL Collaboration], % ``Searches for R-Parity Violating Decays of Gauginos at 183 GeV at LEP,'' Eur.\ Phys.\ J.\ {\bf C11} (1999) 619 [hep-ex/9901037]. %\vspace{-6mm} \bibitem{L3RPV} M.~Acciarri {\it et al.} [L3 Collaboration], % ``Search for R-parity Violating Decays of Supersymmetric Particles in e+e- Collisions at root(s)=189GeV,'' Eur.\ Phys.\ J.\ {\bf C19} (2001) 397 [hep-ex/0011087]. %\vspace{-2mm} \bibitem{TEVATRON} D.~Acosta {\it et al.} [CDF Collaboration], % ``Search For Pair Production of Scalar Top Quarks in R-parity Violating Decay Modes in ppbar Collisions at sqrt{s} = 1.8 TeV,'' Phys.\ Rev.\ Lett.\ {\bf 92} (2004) 051803 [hep-ex/0305010];\\ F.~Abe {\it et al.} [CDF Collaboration], % ``Search for R-parity Violating Supersymmetry using Like-Sign Dielectrons in ppbar Collisions at sqrt{s} = 1.8 TeV,'' Phys.\ Rev.\ Lett.\ {\bf 83} (1999) 2133 [hep-ex/9908063];\\ V.~M.~Abazov {\it et al.} [D0 Collaboration], % ``Search for the Production of Single Sleptons through R-parity Violation in p bar{p} Collisions at sqrt{s}=1.8 TeV,'' Phys.\ Rev.\ Lett.\ {\bf 89} (2002) 261801 [hep-ex/0207100];\\ V.~M.~Abazov {\it et al.} [D0 Collaboration], % ``Search for R-parity Violating Supersymmetry in Dimuon and Four-Jets Channel,'' Phys.\ Rev.\ Lett.\ {\bf 89} (2002) 171801 [hep-ex/0111053];\\ B.~Abbott {\it et al.} [D0 Collaboration], % ``Search for R-parity Violation in Multilepton Final States in p-barp Collisions at sqrt{s}=1.8 TeV,'' Phys.\ Rev.\ {\bf D62} (2000) 071701 [hep-ex/0005034]. %\vspace{-2mm} \bibitem{D0RES} B.~Abbott {\it et al.} [D0 Collaboration], % ``Search for R-parity Violating Supersymmetry in the Dielectron Channel,'' Phys.\ Rev.\ Lett.\ {\bf 83} (1999) 4476 [hep-ex/9907019]. %\vspace{-2mm} \bibitem{BETA0NU} R.~N.~Mohapatra, Phys.\ Rev.\ {\bf D34} (1986) 3457;\\ J.~D.~Vergados, Phys.\ Lett.\ {\bf B184} (1987) 55;\\ M.~Hirsch, H.~V.~Klapdor-Kleingrothaus and S.~G.~Kovalenko, Phys.\ Lett.\ {\bf B352} (1995) 1 [hep-ph/9502315];\\ M.~Hirsch, H.~V.~Klapdor-Kleingrothaus and S.~G.~Kovalenko, Phys.\ Rev.\ Lett.\ {\bf 75} (1995) 17;\\ M.~Hirsch, H.~V.~Klapdor-Kleingrothaus and S.~G.~Kovalenko, Phys.\ Rev.\ {\bf D53} (1996) 1329 [hep-ph/9502385]. %\vspace{-2mm} \bibitem{APV} P.~Langacker, Phys.\ Lett.\ {\bf B256} (1991) 277. %\vspace{-2mm} \bibitem{CCU} V.~Barger, G.~F.~Giudice and T.~Han, Phys.\ Rev.\ {\bf D40} (1989) 2987. %\vspace{-2mm} \bibitem{barbier} R.~Barbier {\it et al.}, % ``R-Parity-violating supersymmetry'', Phys.\ Rept.\ {\bf 420} (2005) 1. \bibitem{RPVIOLATION} J.~Butterworth and H.~Dreiner, Nucl.\ Phys.\ {\bf B397} (1993) 3 and references therein. \bibitem{BUCHMULL} W.~Buchm\"uller, R.~R\"uckl and D.~Wyler, % ``Leptoquarks in lepton-quark collisions,'' Phys.\ Lett.\ {\bf B191} (1987) 442 and Phys.\ Lett.\ {\bf B448} (1999) 320. %\vspace{-2mm} \bibitem{SPIRANLO} T.~Plehn, H.~Spiesberger, M.~Spira and P.~M.~Zerwas, % ``Formation and Decay of Scalar Leptoquarks/Squarks in ep collisions,'' Z.\ Phys.\ {\bf C74} (1997) 611 [hep-ph/9703433]. %\vspace{-2mm} \bibitem{Aaron:2009wp} F.~D.~Aaron {\it et al.} [H1 Collaboration], % ``Events with Isolated Leptons and Missing Transverse Momentum and % Measurement of $W$ Production at HERA,'' Eur.\ Phys.\ J.\ {\bf{C64}} (2009) 251 [arXiv:0901.0488]. %%CITATION = 0901.0488;%% \bibitem{LEGO} K.~Rosenbauer, LEGO~0.02, ``Suche nach Leptoquarks und Leptogluonen im H1-Experiment bei HERA'', Ph.D. thesis (in German), RWTH Aachen, PITHA 95/16 (1995). %\vspace{-2mm} \bibitem{SUSYGEN} S.~Katsanevas and P.~Morawitz, Comput.\ Phys.\ Commun.\ {\bf 112} (1998) 227 [hep-ph/9711417];\\ N.~Ghodbane, S.~Katsanevas, P.~Morawitz and E.~Perez, SUSYGEN\,3, hep-ph/9909499. %\vspace{-2mm} \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, available at \verb$http://www.desy.de/~jung/rapgap/$. %%CITATION = CPHCB,86,147;%% \bibitem{Kwiatkowski:1990es} A.~Kwiatkowski, H.~Spiesberger and H.~J.~M\"ohring, %``Heracles: An Event Generator For E P Interactions At Hera Energies Including Radiative Processes: Version 1.0,'' Comput.\ Phys.\ Commun.\ {\bf 69} (1992) 155. %%CITATION = CPHCB,69,155;%% \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], available at\\ \verb$http://home.thep.lu.se/~torbjorn/Pythia.html$. %%CITATION = HEP-PH 0010017;%% \bibitem{lund} B.~Andersson, G.~Gustafson, G.~Ingelman and T.~Sj\"ostrand, % ``Parton Fragmentation and String Dynamics'', Phys.\ Rept.\ {\bf 97} (1983) 31. \bibitem{Schuler:yg} G.~A.~Schuler and H.~Spiesberger, DJANGOH version 1.4, %``Django: The Interface for The Event Generators Heracles and Lepto'', Proceedings of the Workshop ``Physics at HERA'', eds. W.~Buchm\"uller and G.~Ingelman, DESY (1991), vol.~3, p.~1419. %href{http://www.slac.stanford.edu/spires/find/hep/www?irn=2712555}{SPIRES entry} \bibitem{Lonnblad:1992tz} L.~L\"onnblad, %``ARIADNE version 4: A Program for simulation of QCD cascades implementing the color dipole model,'' Comput.\ Phys.\ Commun.\ {\bf 71} (1992) 15. %%CITATION = CPHCB,71,15;%% \bibitem{Adloff:2002au} C.~Adloff {\it et al.} [H1 Collaboration], %``Measurement of dijet cross sections in photoproduction at HERA,'' Eur.\ Phys.\ J.\ {\bf C25} (2002) 13 [hep-ex/0201006]. %%CITATION = HEP-EX 0201006;%% \bibitem{Aktas:2004pz} A.~Aktas {\it et al.} [H1 Collaboration], % ``A general search for new phenomena in e p scattering at HERA,'' Phys.\ Lett.\ {\bf B602} (2004) 14 [hep-ex/0408044]. %%CITATION = PHLTA,B602,14;%% %\bibitem{Aktas:2004pz} %A.~Aktas {\it et al.} [H1 Collaboration], %``A general search for new phenomena in $ep$ scattering at HERA,'' %Phys.\ Lett.\ B {\bf 602} (2004) 14 %[hep-ex/0408044]. % %%CITATION = PHLTA,B602,14;%% \bibitem{Baur:1991pp} U.~Baur, J.~A.~Vermaseren and D.~Zeppenfeld, %``Electroweak vector boson production in high-energy e p collisions,'' Nucl.\ Phys.\ {\bf B375} (1992) 3. %%CITATION = NUPHA,B375,3;%% \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;%% \bibitem{Brun:1987ma} R.~Brun {\it et al.}, ``GEANT3'', CERN-DD/EE/84-1. \bibitem{Abt:h1} I.~Abt {\it et al.} [H1 Collaboration], % ``The H1 detector at HERA,'' % ``The Tracking, calorimeter and muon detectors of the H1 experiment at HERA ,'' Nucl.\ Instrum.\ Meth.\ {\bf A386} (1997) 310;\\ %%CITATION = NUIMA,A386,310;%% I.~Abt {\it et al.} [H1 Collaboration], % ``The H1 detector at HERA,'' % ``The Tracking, calorimeter and muon detectors of the H1 experiment at HERA,'' Nucl.\ Instrum.\ Meth.\ {\bf A386} (1997) 348. %%CITATION = NUIMA,A386,348;%% \bibitem{Andrieu:1993kh} B.~Andrieu {\it et al.} [H1 Calorimeter Group], %``The H1 liquid argon calorimeter system,'' Nucl.\ Instrum.\ Meth.\ {\bf A336} (1993) 460. %%CITATION = NUIMA,A336,460;%% \bibitem{Andrieu:1993tz} B.~Andrieu {\it et al.} [H1 Calorimeter Group], % ``Results From Pion Calibration Runs For The H1 Liquid Argon Calorimeter And % Comparisons With Simulations,'' Nucl.\ Instrum.\ Meth.\ {\bf A336} (1993) 499. %%CITATION = NUIMA,A336,499;%% \bibitem{Andrieu:1994yn} B.~Andrieu {\it et al.} [H1 Calorimeter Group], % ``Beam tests and calibration of the H1 liquid argon calorimeter with % electrons,'' Nucl.\ Instrum.\ Meth.\ {\bf A350} (1994) 57. %%CITATION = NUIMA,A350,57;%% \bibitem{SpaCal} R.-D.~Appuhn {\it et al.} [H1 SpaCal Group], % ``The H1 Lead/Scintillating-Fibre Calorimeter'', Nucl.\ Instr.\ and\ Meth. {\bf A386} (1997) 397. \bibitem{SpaTestBeam} T.~Nicholls {\it et al.} [H1 SpaCal Group], % ``Performance of an electromagnetic lead/ scintillating fiber calorimeter for the H1 detector'', Nucl.\ Instr.\ and\ Meth. {\bf A374} (1996) 149. \bibitem{Kleinwort:2006zz} C.~Kleinwort [H1 Collaboration], %``H1 alignment experience,'' Proceedings of the Workshop ``First LHC Detector Alignment Workshop'', eds. S.~Blusk {\it et al.}, CERN (2006), CERN-2007-04. \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 \verb$http://www-h1.desy.de/publications/theses_list.html$. \bibitem{kT1} S.D.~Ellis and D.E.~Soper, % ``Successive Combination Jet Algorithm For Hadron Collisions'', Phys.\ Rev.\ {\bf D48} (1993) 3160 [hep-ph/9305266]. \bibitem{kT2} S.~Catani {\it et al.}, % ``Longitudinally-invariant $k_T$-clustering algorithms for hadron-hadron collisions'', Nucl.\ Phys.\ {\bf B406} (1993) 187. \bibitem{BLONDEL} A.~Blondel and F.~Jacquet, Proceedings of ``The Study of an $ep$ Facility for Europe'', ed. U.~Amaldi, DESY 79/48 (1979) 391. \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.\ {\bf C30} (2003) 1 [hep-ex/0304003]. %%CITATION = HEP-EX 0304003;%% \bibitem{nikiforov} A.~Nikiforov, ``Measurements of the neutral current $e^{\pm}p$ cross sections using longitudinally polarised lepton beams at HERA II'', Ph.D. thesis, Ludwig-Maximilians University (2007), available at\\ \verb$http://www-h1.desy.de/publications/theses_list.html$. \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\\ \verb$http://www-h1.desy.de/publications/theses_list.html$. \bibitem{trinh} T.~N.~Trinh, ``Recherche de leptons excit\'es sur le collisionneur HERA avec le d\'etecteur H1'', Ph.D. thesis, Universit\'e de la M\'editerran\'ee Aix-Marseille II, (2008), CPPM-T-2008-02, available at\\ \verb$http://www-h1.desy.de/publications/theses_list.html$. \bibitem{negri} I.~Negri, ``Recherche de fermions excit\'es dans l'exp\'erience H1 aupr\`es du collisionneur positron-proton HERA'', Ph.D. thesis, Universit\'e de la M\'editerran\'ee Aix-Marseille II, (1998), CPPM-T-1998-02, DESY-THESIS-1998-022, available at\\ \verb$http://www-h1.desy.de/publications/theses_list.html$. \bibitem{JOHTHESIS} J.~Haller, ``Search for Squark Production in $R$-Parity Violating Supersymmetry at HERA'', Ph.D. thesis, University of Heidelberg, (2003), DESY-THESIS-2003-035, available at\\ \verb$http://www-h1.desy.de/publications/theses_list.html$. % (available at http://www-h1.desy.de/psfiles/theses/h1th-318.ps). %\vspace{-2mm} \bibitem{MSSM} H.~P.~Nilles, Phys.\ Rept.\ {\bf 110} (1984) 1;\\ H.~E.~Haber and G.~L.~Kane, Phys.\ Rept.\ {\bf 117} (1985) 75. %\vspace{-2mm} \bibitem{Junk:1999kv} T.~Junk, % ``Confidence level computation for combining searches with small statistics,'' Nucl.\ Instrum.\ Meth.\ {\bf A434} (1999) 435 [hep-ex/9902006]. %%CITATION = NUIMA,A434,435;%% \bibitem{MSUGRA} M.~Drees and M.~M.~Nojiri, Nucl.\ Phys.\ {\bf B369} (1992) 54;\\ H.~Baer and X.~Tata, Phys.\ Rev.\ {\bf D47} (1993) 2739;\\ G.~L.~Kane, C.~Kolda, L.~Roszkowski and J.~D.~Wells, Phys.\ Rev.\ {\bf D49} (1994) 6173. %\vspace{-2mm} \bibitem{SUSPECT} A.~Djouadi, J.~L.~Kneur and G.~Moultaka, SUSPECT 2.1 Comput. Phys. Commun. {\bf 176} (2007) 426 [hep-ph/0211331]. % , code written by A.~Djouadi, J.~L.~Kneur % and G.~Moultaka. %\vspace{-2mm} \end{thebibliography} \clearpage % \input{tables} \label{sec:limideri} \begin{table}[t] \renewcommand{\arraystretch}{1.2} \begin{center} \begin{tabular}{|c|cccc|cccc|ccc|} \multicolumn{12}{c}{{H1 Search for Squarks in \Rp\ SUSY}} \\ \hline \rule[-2mm]{0mm}{7mm} Selection& \multicolumn{4}{c|}{ $e^-p$ $(183$~pb$^{-1})$ } & \multicolumn{4}{c|}{ $e^+p$ $(255$~pb$^{-1})$ } & \multicolumn{3}{c|}{Range of Signal} \\ \multicolumn{1}{|c|}{ Channel} & \multicolumn{1}{c}{ Data} & \multicolumn{3}{c|}{SM Expectation} & { Data} & \multicolumn{3}{c|}{ SM Expectation} & \multicolumn{3}{c|}{ Efficiencies}\\ \hline \hline $eq$ & $ 3121 $ & $ 3215 $ & $\pm $ & $ 336 $ & $ 2946 $ & $ 2899 $ & $\pm $ & $ 302 $ & $30\%$&$-$&$40\%$\\ $\nu q$ & $ 2858 $ & $ 2983 $ & $\pm $ & $ 358 $ & -- & & -- & & $50\%$&$-$&$60\%$\\ $eMJ$ (RC) & $ 147 $ & $ 158.3 $ & $\pm $ & $ 23.9 $ & $ 140 $ & $ 146.0 $ & $\pm $ & $ 21.4 $ & $10\%$&$-$&$40\%$\\ $eMJ$ (WC) & $ 0 $ & $ 1.3 $ & $\pm $ & $ 0.3 $ & $ 1 $ & $ 0.6 $ & $\pm $ & $ 0.4 $ & $5\%$&$-$&$20\%$\\ $ee MJ$ & $ 0 $ & $ 1.5 $ & $\pm $ & $ 0.5 $ & $ 2 $ & $ 1.7 $ & $\pm $ & $ 0.5 $ & $5\%$&$-$&$35\%$\\ $e\mu MJ$ & $ 0 $ & $ 0.03 $ & $\pm $ & $ 0.02 $ & $ 0 $ & $ 0.03 $ & $\pm $ & $ 0.03 $ & $5\%$&$-$&$15\%$\\ $e \nu MJ$ & $ 3 $ & $ 5.6 $ & $\pm $ & $ 1.2 $ & $ 5 $ & $ 8.2 $ & $\pm $ & $ 2.0 $ & $5\%$&$-$&$40\%$\\ $\nu MJ$ & $ 204 $ & $ 235.5 $ & $\pm $ & $ 63.3 $ & $ 113 $ & $ 134.0 $ & $\pm $ & $ 33.8 $ & $5\%$&$-$&$15\%$\\ $\nu\mu MJ$ & $ 0 $ & $ 0.04 $ & $\pm $ & $ 0.02 $ & $ 0 $ & $ 0.06 $ & $\pm $ & $ 0.03 $ & $5\%$&$-$&$20\%$\\ \hline \end{tabular} \end{center} \caption[Selection summary: total event numbers] {Total numbers of selected events, SM expectations and ranges of signal efficiencies for the squark decay channels considered in $e^-p$ and in $e^+p$ collisions. The range of signal efficiencies gives the extreme values for squark masses ranging from $100$~GeV to $290$~GeV and gaugino masses ranging from $30$~GeV up to the squark mass. The $\nu q$ channel is not relevant for $e^+p$ data since the $\tilde{u}_L$--type squarks produced in $e^+p$ do not undergo this decay. Only $\tilde{d}_R$--type squarks, which are produced dominantly in $e^-p$ collisions, can undergo direct decay leading to a $\nu q$ final state. The total error on the SM prediction is determined by adding the effects of all model and experimental systematic uncertainties in quadrature.\\ } \label{tab:totnum} \end{table} \clearpage % \input{figures} \begin{figure}[t] \begin{center} \epsfig{file=d10-181f1a.eps,width=0.4\textwidth} \hspace{2cm} \epsfig{file=d10-181f1b.eps,width=0.4\textwidth} \end{center} \begin{picture}(1,2)(1,2) \put(20,5){(a)} \put(107,5){(b)} \end{picture} \caption{ Feynman diagrams for the single resonant s-channel production of right--handed down--type squarks in $e^-p$ collisions (a) and left--handed up--type squarks in $e^+p$ collisions (b) with subsequent decays into SM particles via Yukawa couplings $\lambda'_{11k}$ or $\lambda'_{1j1}$, respectively. The right--handed down--type squarks can decay either into $e^- + {u}$ or $\nu_e + {d}$, while the left-handed up--type squarks decay into $e^+ + d$ only. } \label{fig:directdecays} \end{figure} \begin{figure}[t] \begin{center} \epsfig{file=d10-181f2a.eps,width=0.49\textwidth} \epsfig{file=d10-181f2b.eps,width=0.49\textwidth} \end{center} \begin{picture}(1,2)(1,2) \put(20,5){(a)} \put(107,5){(b)} \end{picture} \caption{Feynman diagrams for squark decays proceeding via gauginos in the case of right--handed down--type squarks (a) and left--handed up--type squarks (b) with subsequent \Rp\ decay into SM fermions via Yukawa couplings $\lambda'_{11k}$ or $\lambda'_{1j1}$, respectively. The resulting final states may contain multi--leptons and multi--jets. The $\tilde{d}_R$--type squarks decay to \mbox{$\chi^0_i$ ($i=1,2,3,4$)} or $\tilde{g}$, decays into charginos are suppressed while $\tilde{u}_L$--type squarks couple also to charginos $\chi^{+}_i$ ($i=1,2$). } \label{fig:gaugedecays} \end{figure} \clearpage \begin{figure}[t] \begin{center} \epsfig{file=d10-181f3a.eps,width=0.45\textwidth} \epsfig{file=d10-181f3b.eps,width=0.45\textwidth}\\ \epsfig{file=d10-181f3c.eps,width=0.45\textwidth} \epsfig{file=d10-181f3d.eps,width=0.45\textwidth} \epsfig{file=d10-181f3e.eps,width=0.45\textwidth} \end{center} \begin{picture}(1,2)(1,2) \put(67,170){(a)} \put(140,170){(b)} \put(67,110){(c)} \put(140,110){(d)} \put(103,47){(e)} \end{picture} \caption{Reconstructed invariant mass distributions in all selection channels with data (points) events from $183$~pb$^{-1}$ of $e^{-}p$ collisions compared to SM MC predictions. The method used for the reconstruction ($M_e,M_h,M_{\rm{rec}},M_{\rm{rec},\nu}$) depends on the analysis channel. The error band gives all model and experimental systematic uncertainties on the SM prediction (solid histogram) added in quadrature. Error bars of data points show statistical uncertainties. The dashed histogram indicates the signal from a squark with $M_{\tilde q} = 150\,\text{GeV}$ with arbitrary normalisation.} \label{fig:masselec} \end{figure} \clearpage \begin{figure}[t] \begin{center} \epsfig{file=d10-181f4a.eps,width=0.45\textwidth} \epsfig{file=d10-181f4b.eps,width=0.45\textwidth}\\ \epsfig{file=d10-181f4c.eps,width=0.45\textwidth} \epsfig{file=d10-181f4d.eps,width=0.45\textwidth} \epsfig{file=d10-181f4e.eps,width=0.45\textwidth} \epsfig{file=d10-181f4f.eps,width=0.45\textwidth} \end{center} \begin{picture}(1,2)(1,2) \put(67,170){(a)} \put(140,170){(b)} \put(67,110){(c)} \put(140,110){(d)} \put(67,47){(e)} \put(140,47){(f)} \end{picture} \caption{Reconstructed invariant mass distributions in all selection channels with data (points) events from $255$~pb$^{-1}$ of $e^{+}p$ collisions compared to SM MC predictions. The method used for the reconstruction ($M_e,M_{\rm{rec}},M_{\rm{rec},\nu}$) depends on the analysis channel. The error band gives all model and experimental systematic uncertainties on the SM prediction (solid histogram) added in quadrature. Error bars of data points show statistical uncertainties. The dashed histogram indicates the signal from a squark with $M_{\tilde q} = 150\,\text{GeV}$ with arbitrary normalisation.} \label{fig:massposi} \end{figure} \clearpage \begin{figure}[t] \begin{center} \epsfig{file=d10-181f5a.eps,width=0.45\textwidth} \epsfig{file=d10-181f5b.eps,width=0.45\textwidth}\\ \vspace{1cm} \epsfig{file=d10-181f5c.eps,width=0.45\textwidth} \epsfig{file=d10-181f5d.eps,width=0.45\textwidth} \end{center} \begin{picture}(1,2)(1,2) \put(21,94){(a)} \put(95,94){(b)} \put(21,23){(c)} \put(95,23){(d)} %\put(12,150){{\sf MSSM scenario with lightest neutralino dominated by its photino component}} \put(27,143){\textsf{\textbf{MSSM scenario with photino-dominated lightest neutralino}}} \end{picture} \caption{Exclusion limits at $95\%$~CL on (a) $\lambda'_{11k}\,(k=1,2)$ and on (b) $\lambda'_{1j1}\,(j=1,2)$ in a phenomenological MSSM with a photino ($\tilde{\gamma}$) like neutralino ($\chi^0_1$). For comparison, the corresponding limit from the previous H1 analysis~\cite{Aktas:2004ij} is also indicated. Also shown are branching ratios to the decay channels considered in this analysis for (c) ${\lambda}'_{11k}$ and (d) ${\lambda}'_{1j1}$ values at the observed limit. } \label{fig:photinolikeelec} \end{figure} \clearpage \begin{figure}[t] \begin{center} \epsfig{file=d10-181f6a.eps,width=0.45\textwidth} \epsfig{file=d10-181f6b.eps,width=0.45\textwidth}\\ \vspace{1cm} \epsfig{file=d10-181f6c.eps,width=0.45\textwidth} \epsfig{file=d10-181f6d.eps,width=0.45\textwidth} \end{center} \begin{picture}(1,2)(1,2) \put(21,94){(a)} \put(95,94){(b)} \put(21,23){(c)} \put(95,23){(d)} %\put(12,150){{\sf MSSM scenario with lightest neutralino dominated by its zino component}} \put(30,143){\textsf{\textbf{MSSM scenario with zino-dominated lightest neutralino}}} \end{picture} \caption{Exclusion limits at $95\%$~CL on (a) $\lambda'_{11k}\,(k=1,2)$ and on (b) $\lambda'_{1j1}\,(j=1,2)$ in a phenomenological MSSM with a zino ($\tilde{Z}$) like neutralino ($\chi^0_1$). Also shown are branching ratios to the decay channels considered in this analysis for (c) ${\lambda}'_{11k}$ and (d) ${\lambda}'_{1j1}$ values at the observed limit. } \label{fig:zinolikeelec} \end{figure} \clearpage \begin{figure}[t] \begin{center} \epsfig{file=d10-181f7a.eps,width=0.65\textwidth}\\ \epsfig{file=d10-181f7b.eps,width=0.65\textwidth} \end{center} \begin{picture}(1,2)(1,2) \put(47,116){(a)} \put(47,28){(b)} \end{picture} \caption{Exclusion limits ($95\,\%$ CL) on $\lambda'_{11k}$ for (a) $k=1,2$ and (b) $k=3$ as a function of the squark mass from a scan of the MSSM parameter space. The dark filled region indicates values of the coupling $\lambda'_{11k}$ excluded in all investigated scenarios whereas the light filled region is excluded only in part of the scenarios. Indirect limits from neutrinoless double beta decay experiments ($\beta\beta 0\nu$)~\cite{BETA0NU,barbier} and tests of charged current universality (CCU)~\cite{CCU,barbier} are also shown. For comparison, the corresponding limits from the previous H1 analysis~\cite{Aktas:2004ij} are also indicated.} \label{fig:scan11k} \end{figure} \clearpage \begin{figure}[t] \begin{center} \epsfig{file=d10-181f8a.eps,width=0.65\textwidth}\\ \epsfig{file=d10-181f8b.eps,width=0.65\textwidth} \end{center} \begin{picture}(1,2)(1,2) \put(47,116){(a)} \put(47,28){(b)} \end{picture} \caption{Exclusion limits ($95\,\%$ CL) on $\lambda'_{1j1}$ for (a) $j=1,2$ and (b) $j=3$ as a function of the squark mass from a scan of the MSSM parameter space as indicated in the figures. The dark filled region indicates values of the coupling $\lambda'_{1j1}$ excluded in all investigated scenarios whereas the light filled region is excluded only in part of the scenarios. Indirect limits from neutrinoless double beta decay experiments ($\beta\beta 0\nu$)~\cite{BETA0NU,barbier} and atomic parity violation (APV)~\cite{APV,barbier} are also shown. For comparison, the corresponding limits from the previous H1 analysis~\cite{Aktas:2004ij} are also indicated.} \label{fig:scan1j1} \end{figure} \clearpage \begin{figure}[t] \begin{center} \epsfig{file=d10-181f9a.eps,width=0.65\textwidth}\\ \epsfig{file=d10-181f9b.eps,width=0.65\textwidth} \end{center} \begin{picture}(1,2)(1,2) \put(60,177){(a)} \put(60,90){(b)} \end{picture} \caption{Exclusion limits ($95\,\%$ CL) in the $m_0,m_{1/2}$ plane assuming $\lambda'_{11k}=0.3$ for (a) $\tan \beta=2$ and (b) $\tan\beta=6$ for $k=1,2$ (hatched region) and $k=3$ (light filled region). A curve of constant squark mass is illustrated for $m(\tilde{d})=280$~GeV. Also indicated are constraints obtained by the L3 experiment at LEP~\cite{L3RPV} and the D\O\ experiment at the Tevatron~\cite{D0RES}. The dark filled region labelled as ``not allowed'' indicates where no REWSB solution is possible or where the LSP is a sfermion. } \label{fig:msugra11k} \end{figure} \clearpage \begin{figure}[t] \begin{center} \epsfig{file=d10-181f10a.eps,width=0.65\textwidth}\\ \epsfig{file=d10-181f10b.eps,width=0.65\textwidth} \end{center} \begin{picture}(1,2)(1,2) \put(60,177){(a)} \put(60,90){(b)} \end{picture} \caption{Exclusion limits ($95\,\%$ CL) in the $m_0,m_{1/2}$ plane assuming $\lambda'_{1j1}=0.3$ for (a) $\tan \beta=2$ and (b) $\tan\beta=6$ for $j=1,2$ (hatched region) and $j=3$ (light filled region). Curves of constant squark mass are illustrated for $m(\tilde{u})=275$~GeV and $m(\tilde{t})=270$~GeV. Also indicated are constraints obtained by the L3 experiment at LEP~\cite{L3RPV} and the D\O\ experiment at the Tevatron~\cite{D0RES}. The dark filled region labelled as ``not allowed'' indicates where no REWSB solution is possible or where the LSP is a sfermion. } \label{fig:msugra1j1} \end{figure} \clearpage \begin{figure}[t] \begin{center} \epsfig{file=d10-181f11.eps ,width=0.8\textwidth} \end{center} \caption{ Exclusion limits for $M=m_0=m_{1/2}$ in mSUGRA as function of $\tan \beta$. Shown are the $95\,\%$ CL exclusion domains for the model parameters from the production of squarks of first and second generation ($\tilde{u},\tilde{c}$ and $\tilde{d},\tilde{s}$) and of third generation ($\tilde{t},\tilde{b}$) assuming a value of $\lambda'=0.3$ for the respective coupling. The area below the curves is excluded. The dark filled region region labelled as ``not allowed'' indicates where no REWSB solution is possible or where the LSP is a sfermion. } \label{fig:tanbeta} \end{figure} \clearpage \clearpage \end{document}