%================================================================
% LaTeX file with preferred layout for the contributed papers to
% the ICHEP Conference 98 in Vancouver
% process with:  latex hep98.tex
%                dvips -D600 hep98
%================================================================
\documentclass[12pt]{article}
\usepackage{epsfig}
\usepackage{amsmath}
\usepackage{hhline}
\usepackage{amssymb}
\usepackage{times}
\renewcommand{\topfraction}{1.0}
\renewcommand{\bottomfraction}{1.0}
\renewcommand{\textfraction}{0.0}
\newlength{\dinwidth}
\newlength{\dinmargin}
\setlength{\dinwidth}{21.0cm}
\textheight23.5cm \textwidth16.0cm
\setlength{\dinmargin}{\dinwidth}
\setlength{\unitlength}{1mm}
\addtolength{\dinmargin}{-\textwidth}
\setlength{\dinmargin}{0.5\dinmargin}
\oddsidemargin -1.0in
\addtolength{\oddsidemargin}{\dinmargin}
\setlength{\evensidemargin}{\oddsidemargin}
\setlength{\marginparwidth}{0.9\dinmargin}
\marginparsep 8pt \marginparpush 5pt
\topmargin -42pt
\headheight 12pt
\headsep 30pt \footskip 24pt
\parskip 3mm plus 2mm minus 2mm

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

% Some useful tex commands
%
\newcommand{\GeV}{\rm\,GeV}
\newcommand{\gev}{\rm\,GeV}
\newcommand{\mev}{\rm\,MeV}
\newcommand{\TeV}{\rm\,TeV}
\newcommand{\pb}{\rm\,pb}
\newcommand{\pbinv}{\rm\,pb^{-1}}
\newcommand{\cm}{\rm cm}
\newcommand{\hdick}{\noalign{\hrule height1.4pt}}
\newcommand{\PT}{p_{\perp}}
\newcommand{\pt}{p_{_{\rm T}}}
\newcommand{\dedxns}{${\rm d}E/{\rm d}x$}
\newcommand{\dedxf}{{\rm d}E/{\rm d}x}
\newcommand{\dedx}{${\rm d}E/{\rm d}x$~}
\newcommand{\thpl}{$\Theta^+$~}
\newcommand{\thplns}{$\Theta^+$}
\newcommand{\thplf}{\Theta^+}
\newcommand{\ksf}{K^0_s}
\newcommand{\knullf}{K^0}
\newcommand{\knull}{$K^0$~}
\newcommand{\GeVSq}{\rm\,GeV^2}
\newcommand{\ximm}{$\Xi^{--}_{5q}\,$}
\newcommand{\xipm}{$\Xi^{0}_{5q}\,$}
\newcommand{\xizero}{$\Xi(1530)^{0}\,$}


\newcommand{\coll}{Collaboration}
\newcommand{\etal}{{\it et al}}

\begin{document}

\pagestyle{empty}
\begin{titlepage}

\noindent
\begin{center}
%{\it {\large version of \today}} \\[.3em] 
\begin{small}
\begin{tabular}{llrr}
%Submitted to & \multicolumn{3}{r}{\footnotesize Electronic Access: {\it http://www-h1.desy.de/h1/www/publications/conf/conf\_list.html}} \\[.2em] \hline 
%Submitted to & \multicolumn{3}{r}{\footnotesize {\it www-h1.desy.de/h1/www/publications/conf/conf\_list.html}} \\[.2em] \hline 
Submitted to & & &
\epsfig{file=/h1/www/images/H1logo_bw_small.epsi
,width=2.cm} \\[.2em] \hline
\multicolumn{4}{l}{{\bf
                33rd International Conference 
                on High Energy Physics, ICHEP06},
                July-August,~2006,~Moscow} \\
                 & Abstract:        &     &\\
                 & Parallel Session & {\bf Hadron Spectroscopy and Exotics }   &\\ \hline
 & \multicolumn{3}{r}{\footnotesize {\it
    www-h1.desy.de/h1/www/publications/conf/conf\_list.html}} \\[.2em]
\end{tabular}
\end{small}
\end{center}
\vspace*{2cm}

\begin{center}
  \Large
  {\bf Search for Baryonic States Decaying to ${\boldmath \Xi\pi}$ \\
   in Deep Inelastic Scattering  at HERA}

  \vspace*{1cm}
    {\Large H1 Collaboration}
%  \footnote{For H1-preliminary-06-131; 
%    \quad {\it Contacts}: M.~Del Degan, C.~Grab} 
\end{center}

\begin{abstract}

\noindent
A search for narrow baryonic states decaying into $\Xi^{-}\pi^{-}$ 
or\ $\Xi^{-}\pi^{+}$ and their antiparticles is carried out
with the H1 detector using deep inelastic scattering events at HERA,
at centre-of-mass energies of 300 and 318$\gev$.
No signal is observed for a new baryonic state in 
the mass range 1600 - 2100$\mev$
 in neither  the $\Xi^{-}\pi^{-}$ nor the $\Xi^{-}\pi^{+}$ 
decay channels.
% that could be interpreted as a possible pentaquark state.
The standard baryon  \xizero is observed through the
decay  mode $\Xi(1530)^{0}\rightarrow\Xi^{-}\pi^{+}$.
Upper limits on the ratio of the  hypothetical pentaquark states 
\ximm and \xipm  to the \xizero signal are given.


\end{abstract}


\end{titlepage}

\pagestyle{plain}

\section{Introduction}


Various experiments have reported evidence for narrow signals with masses around 
1520$\mev$, decaying into  $nK^{+}$ and $pK^{0}_{S}$ final states \cite{kenhicks}.
Such states could be interpreted as an exotic strange pentaquark 
baryon state
with a minimal quark content of $uudd\bar{s}$. On the other hand, a number of other
experiments have reported negative evidence for the same state.

Such states are expected to occur as a flavour
antidecuplet within various theoretical approaches using QCD \cite{theory}.
The state reported above is referred to as the $\Theta^{+}$ \footnote{Unless explicitely
mentioned, the charge 
conjugate states are always implicitly included.}.
Searches for the other members of this decuplet are of interest, in particular
the doubly strange  \ximm $(ddss\bar{u})$  
and $\Xi^+_{5q}$ $(uuss\bar{d})$ states, which are manifestly exotic.
However, states from the standard isospin $3/2$ baryon multiplet 
such as the $\Xi^{0}$ and the $\Xi^{-}$ can also populate
states similar to these pentaquark states
%, but with normal flavour and 
%charge assignments,

The NA49 \cite{na49} collaboration 
reported the observation of the \ximm and \xipm  
members of the $\Xi_{3/2}$ multiplet, with a mass of $1862\pm 2\mev$ and 
a width of less than $18\mev$. 
%
These findings have not been confirmed by several other experiments  
\cite{ximinus-neg},
including the ZEUS collaboration \cite{Chekanov:2005at}.

This analysis describes results of a search for such new baryonic
states decaying into the $\Xi^{-}\pi^{\pm}$ and $\Xi^{+}\pi^{\pm}$
final states. 
The $\Xi^{-}$($\Xi^{+}$) particles were identified through
their decay channels into $\Lambda\pi^{-}$($\bar{\Lambda}\pi^{+}$). 
This search is performed
using 102.5$\pbinv$ of deep inelastic $ep$ scattering data taken with the H1 detector
in the years 1996 - 2000.

%--------------------------------------------------------------------
%
%%
\section{Experimental Procedure}
\label{method}
\subsection{H1 Apparatus}
\label{detector}
In the following only those detector components important 
for the present analysis are described.
A detailed description of the H1 detector can be found in \cite{h1det}. 

The tracks from charged particles used in this analysis are
reconstructed in the 
central tracker, whose main components are two cylindrical drift chambers,
the inner and outer central jet chamber (CJCs) \cite{Burger:eb}.
The inner and outer CJC are mounted concentrically around the beam-line,
covering the range of pseudorapidities\footnote{
The pseudorapidity is given by $\eta = - \ln \tan \theta / 2$,
where the polar angle $\theta$ is measured with respect to 
the \mbox{$z$ axis} given  by the proton beam direction.}
\mbox{$-1.75 < \eta<1.75$} for tracks coming from the 
nominal event vertex.
The CJCs lie within a homogeneous magnetic field of 
$1.15 \, {\rm T}$ which allows the transverse momentum, $\pt$, 
of charged particles to be measured.
Two additional drift chambers complement the CJCs by precisely
measuring the $z$ coordinates of track segments and hence assist in the
determination of polar angles. 
Two cylindrical multi-wire proportional chambers facilitate 
triggering on tracks.
The transverse momentum resolution of the central tracker is
$\sigma(\pt) / \pt \simeq 0.005 \, \pt \, /\,\GeV \, \oplus 0.015$.
Charge misidentification is negligible
for particles originating from the primary vertex and having 
transverse momenta in the range relevant to this analysis. 

The backward region
of H1 is equipped with a lead/scintillating fibre
``Spaghetti'' Calorimeter (SpaCal)~\cite{Appuhn:1996na},
which is optimized for the detection of the scattered positron in the
DIS kinematic range considered here.
It consists of an electromagnetic and a more coarsely
segmented hadronic section.
A planar drift chamber, positioned in front of the SpaCal, 
improves the measurement of the polar angle, and
is used to reject neutral particle background.
%%
The global properties of the hadronic final state are reconstructed 
using an algorithm which combines information from the
central tracker, the SpaCal and the Liquid Argon  
calorimeter, which surrounds the central tracker. 
The DIS events studied in this paper are triggered 
on the basis of an energy deposition detected in the SpaCal,
complemented by signals in the CJCs and the multi-wire proportional 
chambers in the central tracker. 

The luminosity measurement is based on the Bethe-Heitler process 
$ep \rightarrow ep\gamma$, where the photon is detected in a 
calorimeter located downstream of the interaction point.
%%

%--------------------------------------------------------------------
%
\subsection{Simulation of Baryonic states}
\label{models}
%%
To estimate the acceptance, the efficiency and the resolution for the detection of a 
hypothetical baryon state,
a Monte Carlo simulation based on the PYTHIA 6.2 \cite{pythia}
event generator is used, incorporating 
the Lund string model of fragmentation \cite{lund}.
%% 
The kinematic distributions of strange baryons in DIS data have been 
found to be reasonably well described \cite{strangenessdis} by the PYTHIA
simulation.
The \ximm state is introduced by changing the mass of the standard
$\Delta^{--}$ antibaryon to values in the required ranges from 1.6 
to 2.1$\gev$\ and
forcing it to decay into $\Xi \pi$ final states. 
%
Due to lack of knowledge of the pentaquark production mechanism, it 
is assumed that the production kinematics of the \ximm 
is similar to those of other baryons and that it decays isotropically.
In the simulation, the \ximm particle is produced on mass shell.
%%
The generated events are passed through the standard H1 detector 
simulation based on 
GEANT \cite{geant} and then through the identical
reconstruction and analysis chain as the data. 

%--------------------------------------------------------------------
%
\subsection{Selection of DIS Events}
%
The analysis is carried out using data corresponding to an integrated 
luminosity of ${\cal L}=102.5\pbinv$, taken in the years
1996 - 2000. During this time HERA collided 
electrons\footnote{The analysis uses data from periods when
the beam lepton was either a positron or an electron.}  
at an energy of $27.6\gev$ with protons at 
$820\gev$ (1996 - 1997) and $920\gev$ (1999 - 2000)\footnote{The
sample with a proton energy of $820 \, (920)\gev$\ corresponds to 
a luminosity of ${\cal L}= 26.8 \, (75.7) \, {\rm pb^{-1}}$.}.
%

Events are selected if the
$z$ coordinate of the event vertex, reconstructed using the 
central tracker, lies within $35{\rm\,cm}$ of
the mean position for $ep$ interactions. 
%
The scattered electron is required to be reconstructed in the SpaCal
with an energy, $E_e^{\prime}$, above 8$\gev$.
The negative four momentum transfer squared of the exchanged virtual 
photon, $Q^2$, is required to lie in the range 
$2\GeVSq < Q^2 < 120\GeVSq$, as
reconstructed from the energy and polar angle of the scattered electron.  
The \mbox{inelasticity $y$} of the event is reconstructed using
the scattered electron kinematics and is required to be in the range
$0.05 < y < 0.7$.
%The lower cut on $y$ ensures that the hadronic final state lies 
%in the central region of the detector, whilst the upper cut
%corresponds approximately to the cut on $E_e$.
The difference between the total energy $E$ and the longitudinal 
component of the total momentum $p_z$, calculated 
from the electron and the hadronic final state, is restricted
to $35\gev < E - p_z < 70\gev$. This requirement suppresses  
photoproduction background, in which the electron escapes detection and 
a hadron fakes the electron signature.
%


\subsection{Selection of the Baryon Candidates}
\label{lambda}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
The identification of the baryon states is done by the complete
reconstruction of their decay chain through $\Xi$ and
$\Lambda$ baryons into pions and protons:
\begin{equation}
\label{eq:decay}
\Xi^{--} \to \Xi^- \pi^- \to (\Lambda \pi^-) \pi^-  \to [(p \pi^-)\pi^- ] \pi^- 
\end{equation}
The decay daughter particles are fitted
in 3 dimensions to their respective decay vertices \cite{lg-fit}, 
which are referred to 
as the tertiary ($\Lambda$-decay) and the secondary ($\Xi$-decay) vertex.
The analysis is based on charged particles reconstructed in the 
central tracker. If not mentioned otherwise, tracks are accepted if they 
have transverse momenta $\pt > 0.12\gev$ and pseudorapidity 
$|\eta|<1.75$.

In the first step, the $\Lambda$ baryons  are identified by their
charged decay mode,  $\Lambda\to p\pi^{-}$, using pairs of oppositely charged
tracks, assuming a proton and pion hypothesis, respectively.
The track with the higher momentum is assigned the proton mass.
The particles are fitted to a tertiary vertex and required to have 
a vertex fit probability above 1\% and a 2-dim decay length (in the $r\,-\,\phi$ plane)
larger $7.5$\,mm.
The $\Lambda$ transverse momentum is required to be $p_T > 0.3\gev$.
The potential contamination from $K^0_S \to \pi^{+}\pi^{-}$  
is rejected by a cut on 
the $\pi \pi$ invariant mass of $\pm 10\mev$ around the nominal $K^0_S$ mass.
The resulting invariant-mass spectra
for $p\pi^{-}$ and $\bar{p}\pi^{+}$  are  shown in 
Fig.~\ref{lam_mass_spec}. 
%The total measured number of $\Lambda \to p\pi^{-}$ baryons
%is larger than the $\bar{\Lambda} \to \bar{p}\pi^{+}$ antibaryons.
%due to secondary interactions producing $\Lambda$ in the beampipe.    


In the second step, $\Xi^{-}$ candidates are formed by
fitting the $\Lambda$ candidates, taken within $\pm$ 8$\mev$ of the nominal
$\Lambda$ mass, with negatively charged tracks to a secondary vertex.
$\Xi^{-}$ candidates are retained if the vertex fit probability 
is above 1\% and the 2-dim distance of closest approach (DCA) 
of the $\Xi$ track 
was smaller than \mbox{2.5\,mm} with respect to the primary vertex.
%
To increase the significance of the $\Xi^{-}$ signal 
it is required that  
the angle between the secondary and the tertiary vertex vectors
be less than 0.6\,rad.
%
For the $\Xi$ candidate combinations, passing these criteria, the
invariant-mass spectra $\Lambda \pi^{-}$ and $\bar{\Lambda}\pi^{+}$  are 
shown in Fig.~\ref{xi_mass_spec}. 
%Because of the short decay length of the $\Xi$,
%the cuts eliminate the contributions from secondary interactions in the 
%beampipe; the number of $\Xi^{-}$ and $\bar{\Xi}^{+}$ are the same within
%their uncertainties.
%The measured width of the $\bar{\Xi}^{+}$ is somewhat larger than that of the
%$\Xi^{-}$ as expected from the different momentum resolution for 
%positive and negative particles.


In the third step, the $\Xi^{-}$ candidates within $\pm 15\mev$ 
of the nominal 
$\Xi$ mass are selected and combined with candidate  $\pi^{\pm}$ tracks.
The $\pi$ tracks are required to originate from the primary vertex and to have
a significance(DCA/$\sigma_{DCA}$) to the primary 
vertex of less than 4.
To reduce the background further, the transverse momentum of the $\Xi \pi$ combination is required to be $ > 1.0\gev$.
These $\Xi^{\pm} \pi^{\pm}$ and  $\Xi^{\pm} \pi^{\mp}$
 combinations are investigated further in section 3.
%with respect to containing the exotic states or standard baryon states:
%namely the \ximm state in the case of equal charges,
%and the \xipm or $\bar{\Xi}^{0}$ states in the case of opposite charges.


%--------------------------------------------------------------------

%%%%%%%%%%%%%%%%%%%%
\section{Resulting Spectra and Limit Calculation}
%%%%%%%%%%%%%%%%%%%%

The resulting invariant-mass spectra for each of the four charge combinations 
$\Xi^{-}\pi^{-}$, $\Xi^{-}\pi^{+}$, ${\Xi}^{+}\pi^{-}$, 
${\Xi}^{+}\pi^{+}$
are shown in Fig.~\ref{penta_mass_spec}a,b,d,e separately.
%
The sum of the equal charge ($++$ and $--$) 
and opposite charge ($+-$ and $-+$) combinations are also
shown in Fig.~\ref{penta_mass_spec}c and f.
In both the $\Xi^{-}\pi^{+}$ and ${\Xi}^{+}\pi^{-}$ 
(Fig.~\ref{penta_mass_spec}a,b spectra, and correspondingly in the
Fig.~\ref{penta_mass_spec}c),
the signal of the well known \xizero state 
\cite{pdg04} is clearly observed.
%
A fit of the data (Fig.~\ref{penta_mass_spec}c)
with a Gaussian for the signal and background is made using,
%
\begin{equation}
\label{bgeqn}
B(M)=P_{1}(M-m_{\Xi}-m_{\pi})^{P_{2}} \times (P_{3} +  P_{4}M + P_{5} M^2),
\end{equation}
%
where $M$ is the $\Xi\pi$ invariant mass, $m_{\Xi}$ and 
$m_{\pi}$ are the masses
of the $\Xi$ and the $\pi$, respectively, and $P_{i}$ are free parameters. 

The fit yields a total 
of $170\,\pm\,28\,(stat.)$ \xizero baryons. The reconstructed mass of $1530\,\pm\,2\,(stat.)\mev$
is consistent with the PDG value \cite{pdg04}.
The measured width of $8.9\,\pm\,1.3(stat.)\mev$ is 
consistent with the detector resolution, and is also correctly 
described by the simulations.

Figure \ref{penta_mass_spec_Q20} shows the sum of the invariant mass spectra for equal and opposite
charge ($+-$ and $-+$ and $--$ and $++$) combinations for events having $Q^2$ larger than 20$\GeVSq$.
The ZEUS collaboration reported to have seen a hint for a structure around 1690$\mev$  representing the $\Xi(1690)^0$ 
baryonic state \cite{Chekanov:2005at}. In this analysis there is no hint for such a structure observed.

The ratio of the equal to opposite charge mass spectra
(Fig.~\ref{penta_mass_spec}c divided by Fig.~\ref{penta_mass_spec}f),
is shown in Fig.\ref{penta_ratio}. It appears to be reasonably 
flat in the mass region above the 
\xizero. This supports
the use of the same background shape function for both charge combinations
and indicates that the reconstruction efficiencies
are very similar. 
%The differences are included in the systematic
%errors.

No signal is observed near $1860\mev$ in any of the equal charge 
mass spectra (Fig.~\ref{penta_mass_spec}a,b) nor 
in the opposite charge ones (Fig.~\ref{penta_mass_spec}d,e).
Therefore, upper limits are set on the production of new states
decaying to $\Xi^-\pi^-$ in the mass range 1.6 - 2.1$\gev$
in the kinematic region $0.05 < y < 0.7$ for $2\gev^2 < Q^2 < 120\gev^2$.
%
To quantify such limits
the ratio R(M) is determined using
\begin{equation}
\label{eq:ratio}
R_{u.l.}(M) = \frac{\frac{dN_{u.l.}}{dM}(\Xi^{\pm}\,\pi^{\pm})}{N(\Xi(1530)^{0})},
\end{equation}
where M is the $\Xi\pi$ invariant mass and $N(\Xi(1530)^{0})$ represents
the number of \xizero states, measured in the data with the
same kinematic conditions.
This ratio is presumed to be sensitive to the existence of new states,
because of the clean signal observed for the established \xizero.
The ratio also has the advantage that
the systematic effects of the acceptances and the reconstruction efficiencies
mostly cancel, making it insensitive to detector effects and thus provides
a robust method for the limit determination.

%
An upper limit at the 95\% confidence level (C.L.) on $R_{u.l.}(M)$ 
is obtained from the observed spectra
% from the background as modelled according to Eq.~(\ref{bgeqn}),
%and from the signal $M$ distributions (assumed Gaussian)  
using a modified 
frequentist approach based on likelihood ratios \cite{tjunk}. 
This method takes the 
statistical and systematic uncertainties in the 
number of signal and background combinations into account.
The M(\ximm) distribution for signal combinations is assumed to be  
a Gaussian with a mean $M$ and a width corresponding
to the experimental mass resolution as obtained from the
Monte Carlo simulations. 
This width $\sigma(M)$ 
varies from 7.5 to 18.4$\mev$ in the mass ranges
from 1.6 to 2.1$\gev$ considered here.
The background $M$ distribution is modelled according to 
the data given by equation (\ref{bgeqn}).
The limit is then calculated in a sliding mass window 
of mass-dependent width $2\sigma(M)$ of the 
hypothetical particle in the mass range.




%% %%%%%%%%%%%%%%%%%%%%%%%
%% systematic uncertainty 
%% %%%%%%%%%%%%%%%%%%%%%%%

The systematic uncertainties in the number $N(M)$ of observed
$\Xi^- \pi^-$ combinations due to \ximm production comprises the following
 contributions:
%%
\begin{itemize}
\item The uncertainty on the overall number of \xizero candidates.
%
\item The difference in reconstruction efficiency of the 
charged ($\Xi^-\pi^-$) and the
neutral ($\Xi^-\pi^+$) combinations.
%
\item The uncertainty on the width of the potential pentaquark signal.
%
\item 
The uncertainty on the background distribution.
%%
\item The uncertainties due to the inclusive DIS event selection
and the trigger efficiency are assumed to cancel in the ratio.
\end{itemize}
%

% the actual resulting limit plots:
The final result of the limit calculation is shown in 
Figure \ref{pentamm_limits},
which shows in the upper part again the 
$\Xi \pi$ invariant-mass spectrum summed for 
the two equal charge combinations (upper fig., solid dots)
as a function of the central value of the mass window $M$.
In the lower part of Figure \ref{pentamm_limits},
the $95\%$ C.L. upper limit on the ratio  $R_{u.l.}(M)$ is given; 
it is found to vary from 0.15 to 0.45 in the mass range 
1.6 - 2.1$\gev$.
In the NA49 signal region for the \ximm candidate particle at 1.86\gev, 
the ratio $R_{u.l.}(M)$ is found to be less than 0.2 at the $95\%$ C.L.

%
Furthermore, no structures are observed  in the opposite
charge mass spectra Fig.~\ref{penta_mass_spec}c, in the mass region
above the \xizero state. Therefore, separate 
$95\%$ C.L. upper limits are set
on the production of new states 
decaying to $\Xi^-\pi^+$ or $\Xi^+\pi^-$ in the mass range 1.6 - 2.1\gev.
Figure \ref{pentamp_limits} shows the resulting upper limit ratio $R_{u.l.}(M)$, which
is obtained in a procedure analogous to the $\Xi^-\pi^-$ case, excluding
the existence of an  $\Xi^{0}_{5q}$ pentaquark state at a level ranging
from 0.15 to 0.65 with respect to the \xizero normalisation.

Finally figure \ref{pentaall_limits} shows the mass spectra and $95\%$ C.L. upper limit for all 
 charge combinations (mass spectra in Fig. \ref{penta_mass_spec}c and \ref{penta_mass_spec}f 
summed up). In this case the upper limit $R_{u.l.}(M)$ ranges from 0.2 to 0.6.

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

A search for new baryonic states that decay to $\Xi^{-}\pi^{-}$ 
and $\Xi^{-}\pi^{+}$ and their conjugate states 
was performed with the H1 detector using a DIS 
data sample with $2 < Q^{2} < 120\gev^{2}$ and $0.05 < y < 0.7$,
corresponding to a total  integrated luminosity of $102.5\pbinv$. 
A clear signal from the established \xizero baryon state was seen.
No pentaquark candidate signal was found and 
upper limits at $95\%$ C.L. are set on the ratio of  the 
hypothetical \ximm and \xipm  signals to the  total number of
observed \xizero baryons in the mass range 1.6 to 2.1$\gev$.

The results reported here are fully compatible with limits measured
by the ZEUS collaboration.
The overall H1 statistics in the $\Xi$-samples is comparable 
with the NA49 data.
On the other hand, NA49 has a different production environment and kinematics 
and the limits obtained at HERA may therefore not contradict the 
NA49 findings.



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\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.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{thebibliography}{99}
%%
\bibitem{kenhicks} K.~Hicks, Prog Part. Nucl. Phys. {\bf 55} (2005) 647
[hep-ex/0504027].\\
%\bibitem{zeusspq} 
S.~Chekanov {\it et al.} [ZEUS Collaboration],
Phys. Lett. {\bf B591} (2004) 7 [hep-ex/0403051].
%
\bibitem{theory}
D.~Diakonov, V.~Petrov and M.~Polyakov, Z. Phys.  A {\bf 359}, 305 (1997)
[hep-ph/9703373];\\
R.~L.~Jaffe and F.~Wilczek, Phys. Rev. Lett.  {\bf 91} (2003) 23
[hep-ph/0307341];\\
M.~Karliner and H.~Lipkin, Phys. Lett. {\bf B575} (2003) 249 [hep-ph/0307243];\\
for a review on pentaquark phenomenology see
R.~L.~Jaffe, Phys.~Rept. {\bf 409} (2005) 1 [hep-ph/0409065].
%
\bibitem{na49}
NA49 \coll, C.~Alt {\it et al.},
Phys.\ Rev.\ Lett.{} 92~(2004)~042003 [hep-ex/0310014].
%
\bibitem{ximinus-neg}
WA89 \coll, M.I.~Adamovich {\it et al.},
Preprint \mbox{hep-ex/0405042}, 2004; \\
ALEPH \coll, S.~Schael {\it et al.},
Phys.\ Lett.{} B~599~(2004)~1; \\
BES \coll, J.Z.~Bai {\it et al.},
Phys.~Rev.{} D~70~(2004)~012004; \\
HERA-B \coll, I.~Abt {\it et al.},
Phys.\ Rev.\ Lett.{} 93~(2004)~212003; 
%ZEUS \coll. S.Chekanov \etal, DESY 05-18 (2005); \\
%
%\cite{Chekanov:2005at}
\bibitem{Chekanov:2005at}
S.~Chekanov {\it et al.}  [ZEUS Collaboration],
%``Search for pentaquarks decaying to Xi pi in deep inelastic scattering  at
%HERA,''
Phys.\ Lett.\ B {\bf 610} (2005) 212
[arXiv:hep-ex/0501069].
%
\bibitem{h1det} 
I. Abt {\it et al.} [H1 Collaboration], 
Nucl. Inst. Meth. A {\bf 386} (1997) 310; \\
I. Abt {\it et al.} [H1 Collaboration], Nucl. Inst. Meth. A {\bf 386} 
(1997) 348;
%
% CJC
%\cite{Burger:eb}
\bibitem{Burger:eb}
J.~B\"urger {\it et al.},
%``The Central Jet Chamber Of The H1 Experiment,''
Nucl.\ Instrum.\ Meth.\ A {\bf 279} (1989) 217.
%%CITATION = NUIMA,A279,217;%%

%%Spacal
\bibitem{Appuhn:1996na}
R.~D.~Appuhn {\it et al.} [H1 SpaCal Group],
Nucl. Inst. Meth. A {\bf 386} (1997) 397;

%Pythia
\bibitem{pythia} T. Sj\"{o}strand {\it et al.}, 
Comput. Phys. Commun. {\bf 135} (2001) 238 [hep-ph/0010017].
%%

\bibitem{lund} B. Andersson, G. Gustafson, G. Ingelman and T. Sj\"ostrand, 
Phys. Rept. {\bf 97} (1983) 31.
%%no hep-ph

\bibitem{strangenessdis}
M.~Derrick {\it et al.} [ZEUS Collaboration],
%  {\it Neutral strange particle production in deep inelastic
%    scattering at HERA},
Z. Phys. C {\bf 68} (1995) 29 [hep-ex/9505011];\\
S.~Aid {\it et al.} [H1 Collaboration],
%  {\it Strangeness Production in Deep Inelastic Positron-Proton
%    Scattering at HERA},
Nucl. Phys. B {\bf 480} (1996) 3 [hep-ex/9607010];\\
C. Risler, 
%`Die Produktion seltsamer neutraler Teilchen 
%in tiefinelastischer Streuung bei HERA`,
Ph.D. thesis, 2004, Universit\"{a}t
Hamburg (in German),
available from \\
http://www-h1.desy.de/publications/theses\_list.html.
%%
\bibitem{geant}  R. Brun {\it et al.}, GEANT3, Technical Report CERN-DD/EE/84-1,
CERN, 1987.
\bibitem{lg-fit} 
R.~Luchsinger and C.~Grab,
%``Vertex reconstruction by means of the method of Kalman filter,''
Comput.\ Phys.\ Commun.\  {\bf 76} (1993) 263.

%% 
\bibitem{pdg04} Particle Data Group, S. Eidelman {\it et al.}, 
Phys. Lett. {\bf B592} (2004) 1.
%%
\bibitem{tjunk} T.~Junk, Nucl. Inst. Meth. A {\bf 434} (1999) 435 [hep-ex/9902006].
%%
\end{thebibliography}
%\vspace{3cm}

\clearpage


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%% Figures
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Lambda signals: L + AL separately
\begin{figure}
\begin{center}
\includegraphics[width=77mm]{figs/LambdaMass_pm.eps} 
\includegraphics[width=77mm]{figs/LambdaMass_mp.eps} 
\caption{The invariant-mass spectra for (a) $p\pi^{-}$ and (b) $\bar{p}\pi^{+}$  
combinations in the DIS kinematic region defined in the text.
The solid line shows the result of a fit 
using a double Gaussian plus a second-order polynomial  
background function.} 
\label{lam_mass_spec}
\end{center}
\end{figure}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Xi- signal all ST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\newpage
\begin{figure}
\begin{center}
\includegraphics[width=77mm]{figs/XiMinusMass_0m.eps}
\includegraphics[width=77mm]{figs/XiMinusMass_0p.eps}
\caption{The invariant-mass spectra for 
(a) $\Lambda\pi^{-}$ and (b) $\bar{\Lambda}\pi^{+}$ 
particle combinations.
The solid line shows the result of a fit 
using a Gaussian plus a background function similar to that one defined by eq.~(\ref{bgeqn})
, while the dashed line shows the background only.
} 
\label{xi_mass_spec}
\end{center}
\end{figure}

%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Xi(1530)* signal   and Xi-- side-by-side
%%
\newpage
\begin{figure}
\begin{center}
\includegraphics[width=77mm]{figs/XiMass_mp.eps}
\includegraphics[width=77mm]{figs/XiMass_nn.eps} \\
\includegraphics[width=77mm]{figs/XiMass_pm.eps} 
\includegraphics[width=77mm]{figs/XiMass_pp.eps} \\
\includegraphics[width=77mm]{figs/XiMass_neutral.eps} 
\includegraphics[width=77mm]{figs/XiMass_charged.eps} \\
\caption{The  invariant mass spectra of
the $\Xi \pi$ particle combinations for:
(left column) the opposite charge combinations 
$-+, +-$ and their sum,
and (right column) equal charge combinations $--$,$++$ and their sum.
The solid line in (a) is the result of a fit to the data using  
a Gaussian plus the background function 
defined by eq.~(\ref{bgeqn}).  
The dashed line shows the background according to this fit. 
The same background function is also shown in figures in 
the right column.}
\label{penta_mass_spec}
\end{center}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\newpage
\begin{figure}
\begin{center}
\includegraphics[width=90mm]{figs/XiMass.AllComb.Q2gt20.eps}
\caption{The invariant mass spectra of the $\Xi \pi$ particle combinations for equal and opposite 
charges for $Q^2$ larger than 20$\gev^2$. The solid line shows the result of a fit using the background function 
defined by eq.~(\ref{bgeqn}).}
\label{penta_mass_spec_Q20}
\end{center}
\end{figure}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Ratio plot : Xi-- / Xi(1530)* signal 
%%
\begin{figure}
\begin{center}
\includegraphics[width=90mm]{figs/Ratio_neutral_charged.eps}
\caption{The ratio of the two invariant mass spectra for the 
opposite charge to equal charge $\Xi\pi$ combinations.
%i.e. the ratio of the $\Xi(1530)^0$ and the \ximm.
The solid line is the result of a fit to the ratio, 
yielding a constant value of 1.13.}
\label{penta_ratio}
\end{center}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%  Xi-- and limits overlaid:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\newpage
\begin{figure}
\begin{center}
\includegraphics[width=100mm]{figs/LimitPlots_charged.eps}%
\caption{The $\Xi \pi$ 
invariant-mass spectrum summed for
the two equal charge combinations (upper fig., solid dots).
The solid line is the result of a fit to the data using  
the background function defined in eq.~(\ref{bgeqn}).  
%
The lower fig. shows as a function of the invariant mass
the $95\%$ C.L. upper limit on the ratio $R$, as 
defined by eq.~(\ref{eq:ratio}).
}
\label{pentamm_limits}
\end{center}
\end{figure}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%  Xi^0_5q limits overlaid:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\newpage
\begin{figure}
\begin{center}
\includegraphics[width=100mm]{figs/LimitPlots_neutral.eps}%
\caption{The 
$\Xi \pi$ invariant-mass spectrum summed for
the two opposite charge combinations.
The solid line is the result of a fit to the data using  
a Gaussian plus the background function defined in eq.~(\ref{bgeqn}).  
Clearly visible at a mass of 1530\mev is the \xizero baryon. 
%
The lower fig. shows as a function of the invariant mass
the $95\%$ C.L. upper limit on the ratio $R$, as 
defined by eq.~(\ref{eq:ratio}).
}
\label{pentamp_limits}
\end{center}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%  Xi^0_5q limits overlaid:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\newpage
\begin{figure}
\begin{center}
\includegraphics[width=100mm]{figs/LimitPlots_allCombs.eps}%
\caption{The 
$\Xi \pi$ invariant-mass spectrum summed for
all charge combinations.
The solid line is the result of a fit to the data using  
a Gaussian plus the background function defined by eq.~(\ref{bgeqn}).  
Clearly visible at a mass of 1530$\mev$ is the \xizero baryon. 
%
The lower fig. shows as a function of the invariant mass
the $95\%$ C.L. upper limit on the ratio $R$, as 
defined by eq.~(\ref{eq:ratio}).
}
\label{pentaall_limits}
\end{center}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


\end{document}
%--------------------------------------------------------------------

