%================================================================
% LaTeX file with prefered layout for H1 paper drafts
% use: dvips -D600 file-name
%================================================================
\documentclass[12pt]{article}
\usepackage{epsfig}
\usepackage{amsmath}
\usepackage{hhline}
\usepackage{amssymb}
\usepackage{times}
\usepackage{cite}
\usepackage{graphicx}
\usepackage{ifpdf}
\ifpdf
   \usepackage{epstopdf}
   \usepackage{pdfsync}
\fi
%%%%%%%%%%%% Comment the next two lines to remove the line numbering
%\usepackage[]{lineno}
%\linenumbers
%\usepackage{hyperref}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\renewcommand{\topfraction}{1.0}
\renewcommand{\bottomfraction}{1.0}
\renewcommand{\textfraction}{0.0}
\newlength{\dinwidth}
\newlength{\dinmargin}
\setlength{\dinwidth}{21.0cm}
\textheight21.0cm \textwidth19.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}
% The rest
\newcommand{\pom}{{I\!\!P}}
\newcommand{\reg}{{I\!\!R}}
\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{\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{\PT}{p_{\perp}}
\newcommand{\JPSI}{J/\psi}
\newcommand{\PO}{I\!\!P}
\newcommand{\xbj}{x}
\newcommand{\xpom}{x_{\PO}}
\newcommand{\dgr}{^\circ}
\newcommand{\gev}{\,\mbox{GeV}}
\newcommand{\GeV}{\mathrm{GeV}}
\newcommand{\mrad}{\mathrm{mrad}}
\newcommand{\GeVS}{\rm GeV^2}
\newcommand{\xp}{x_p}
\newcommand{\xpi}{x_\pi}
\newcommand{\xg}{x_\gamma}
\newcommand{\xgj}{x_\gamma^{jet}}
\newcommand{\xpj}{x_p^{jet}}
\newcommand{\xpij}{x_\pi^{jet}}
\renewcommand{\deg}{^\circ}
\newcommand{\qsq}{\ensuremath{Q^2} }
\newcommand{\et}{\ensuremath{E_t^*} }
\newcommand{\rap}{\ensuremath{\eta^*} }
\newcommand{\gp}{\ensuremath{\gamma^*}p }
\newcommand{\der}{{\mathrm d}}% Journal macro
 \normalsize



\begin{center}
\Large
{\boldmath \bf
Analysis of Feynman scaling of Photon and Neutron Production in the Very Forward Direction
in Deep-Inelastic Scattering at HERA
}
\end{center}

\begin{abstract}
\noindent
\rm
Measurements of normalised cross sections for the production of
photons and neutrons at very small angles with respect to the proton
beam direction as a function of Feynman-$x$ in deep-inelastic
positron-proton scattering at HERA are presented.
The data are taken with the H1 detector in the years 2006 and 2007 and
correspond to an integrated luminosity of $126~\mathrm{pb}^{-1}$.
The analysis covers the range of negative four momentum transfer squared
at the positron vertex $6<Q^2<100~ \GeV^2$, inelasticity $0.05<y<0.6$ and
the centre-of-mass energy of the virtual photon-proton system
$70<W<245~\GeV$.
The dependence of cross sections on $W$ is investigated.
The cross sections are compared to the predictions of DIS models and
models of the hadronic interactions of high energy cosmic rays.
\end{abstract}

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

\newpage
\begin{figure}[h]
\epsfig{file=H1prelim-13-012.fig1a.eps,width=75mm}
\epsfig{file=H1prelim-13-012.fig1b.eps,width=75mm}
\caption{
The fraction of DIS events with forward photons
in the kinematic region $6< Q^2 < 100~$GeV$^2$ and $0.05<y<0.6$
and the pseudorapidity of the photon $\eta>7.9$.
as a function of $W$.
The predictions of MC models  are compared to the measurements.
}
\label{fig:1}
\end{figure}

\begin{figure}[h]
\epsfig{file=H1prelim-13-012.fig2a.eps,width=75mm}
\epsfig{file=H1prelim-13-012.fig2b.eps,width=75mm}
\caption{
The fraction of DIS events with forward neutrons
in the kinematic region $6< Q^2 < 100~$GeV$^2$ and $0.05<y<0.6$
and the pseudorapidity of the photon $\eta>7.9$.
as a function of $W$.
The predictions of MC models  are compared to the measurements.
}
\label{fig:1n}
\end{figure}


\newpage

\begin{figure}[h]
\hspace*{-5mm}
\epsfig{file=H1prelim-13-012.fig3a.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig3b.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig3c.eps,width=62mm}
\end{figure}
\begin{figure}[h]
\hspace*{-5mm}
\epsfig{file=H1prelim-13-012.fig4a.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig4b.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig4c.eps,width=62mm}
\caption{
Normalised cross sections of forward photons production in DIS
as a function of $x_F$ in the region $\eta>7.9$, $6< Q^2 < 100~$GeV$^2$ 
and $0.05<y<0.6$, and in three $W$ intervals.
The data are compared to two predictions of the DJANGOH Monte Carlo simulation, 
using LEPTO and CDM to simulate higher orders.
The lower row shows the ratios of the MC predictions to the data.
%
% Also shown are models of hadronic interactions, QGSJET,
% EPOS and SIBYLL. 
The error bars show the total experimental uncertainty, defined as the 
quadratic sum of the statistical and systematic uncertainties.
}
\label{fig:2}
\end{figure}

\newpage

\begin{figure}[h]
\epsfig{file=H1prelim-13-012.fig5a.eps,width=80mm}
\epsfig{file=H1prelim-13-012.fig5b.eps,width=80mm}
\caption{
Ratios of normalised cross sections of forward photons production in DIS
corresponding to two different $W$ intervals, shown in Fig.2, as a function 
of $x_F$.
The left side figure shows the ratio of the cross section in $130<W<190$~GeV
interval to the cross section at $70<W<130$~GeV interval.
The right side figure shows the ratio of the cross section in $190<W<250$~GeV
interval to the cross section at $70<W<130$~GeV interval.
Predictions of LEPTO and CDM MC models are compared to the measurement.
}
\label{fig:3}
\end{figure}

\newpage

\begin{figure}[h]
\hspace*{-5mm}
\epsfig{file=H1prelim-13-012.fig6a.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig6b.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig6c.eps,width=62mm}
\end{figure}
\begin{figure}[h]
\hspace*{-5mm}
\epsfig{file=H1prelim-13-012.fig7a.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig7b.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig7c.eps,width=62mm}
\caption{
Normalised cross sections of forward photons production in DIS
as a function of $x_F$ in the region $\eta>7.9$, $6< Q^2 < 100~$GeV$^2$ 
and $0.05<y<0.6$, and in three $W$ intervals.
The data are compared to predictions os hadronic interactions, QGSJET,
EPOS and SIBYLL.
The lower row shows the ratios of the MC predictions to the data.
The error bars show the total experimental uncertainty, defined as the 
quadratic sum of the statistical and systematic uncertainties.
}
\label{fig:4}
\end{figure}

\newpage


\begin{figure}[h]
\epsfig{file=H1prelim-13-012.fig8a.eps,width=80mm}
\epsfig{file=H1prelim-13-012.fig8b.eps,width=80mm}
\caption{
Ratios of normalised cross sections of forward photons production in DIS
corresponding to two different $W$ intervals, shown in Fig.4, as a function 
of $x_F$.
The left side figure shows the ratio of the cross section in $130<W<190$~GeV
interval to the cross section at $70<W<130$~GeV interval.
The right side figure shows the ratio of the cross section in $190<W<250$~GeV
interval to the cross section at $70<W<130$~GeV interval.
Predictions of LEPTO and CDM MC models are compared to the measurement.
Predictions of CR models are compared to the measurement.
}
\label{fig:5}
\end{figure}

\newpage


\begin{figure}[h]
\hspace*{-5mm}
\epsfig{file=H1prelim-13-012.fig9a.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig9b.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig9c.eps,width=62mm}
\end{figure}
\begin{figure}[h]
\hspace*{-5mm}
\epsfig{file=H1prelim-13-012.fig10a.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig10b.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig10c.eps,width=62mm}
\caption{
Normalised cross sections of forward neutron production in DIS
as a function of $x_F$ in the region $\eta>7.9$, $6< Q^2 < 100~$GeV$^2$ 
and $0.05<y<0.6$, and in three $W$ intervals.
The data are compared to two predictions of the Monte Carlo simulations
CDM (Djangoh) and RAPGAP-$\pi$, 
and the combination of CDM and RAPGAP-$\pi$, multipled by the normalisation factors.
The lower row shows the ratios of the MC predictions to the data.
The error bars show the total experimental uncertainty, defined as the 
quadratic sum of the statistical and systematic uncertainties.
}
\label{fig:2n}
\end{figure}

\newpage

\begin{figure}[h]
\hspace*{-5mm}
\epsfig{file=H1prelim-13-012.fig11a.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig11b.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig11c.eps,width=62mm}
\end{figure}
\begin{figure}[h]
\hspace*{-5mm}
\epsfig{file=H1prelim-13-012.fig12a.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig12b.eps,width=62mm}
\epsfig{file=H1prelim-13-012.fig12c.eps,width=62mm}
\caption{
Normalised cross sections of forward neutron production in DIS
as a function of $x_F$ in the region $\eta>7.9$, $6< Q^2 < 100~$GeV$^2$ 
and $0.05<y<0.6$, and in three $W$ intervals.
The data are compared to predictions os hadronic interactions, QGSJET,
EPOS and SIBYLL.
The lower row shows the ratios of the MC predictions to the data.
The error bars show the total experimental uncertainty, defined as the 
quadratic sum of the statistical and systematic uncertainties.
}
\label{fig:4n}
\end{figure}


\newpage


\begin{figure}[h]
\epsfig{file=H1prelim-13-012.fig13a.eps,width=80mm}
\epsfig{file=H1prelim-13-012.fig13b.eps,width=80mm}
\caption{
Ratios of normalised cross sections of forward neutron production in DIS
corresponding to two different $W$ intervals, shown in Fig.8, as a function 
of $x_F$.
The left side figure shows the ratio of the cross section in $130<W<190$~GeV
interval to the cross section at $70<W<130$~GeV interval.
The right side figure shows the ratio of the cross section in $190<W<245$~GeV
interval to the cross section at $70<W<130$~GeV interval.
Predictions of CR models are compared to the measurement.
}
\label{fig:5n}
\end{figure}

\end{document}
