%================================================================
% 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}
\usepackage{color}
\usepackage{colordvi}
   \usepackage{epstopdf}
 \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}
\textheight22.8cm \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}
% 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{\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{\gevsq}{\ensuremath{\mathrm{GeV}^2} }
\newcommand{\et}{\ensuremath{E_t^*} }
\newcommand{\rap}{\ensuremath{\eta^*} }
\newcommand{\gp}{\ensuremath{\gamma^*}p }
\newcommand{\der}{{\mathrm d}}
\newcommand{\kd}{\textcolor{black}}
% Journal macro
\def\Journal#1#2#3#4{{#1} {\bf #2}, #4 (#3)}
\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\PR{{\em Phys. Rev.}    }
\def\ZPC{{\em Z. Phys.}      {\bf C}}
\def\ZP{{\em Z. Phys.}      }
\def\EJC{{\em Eur. Phys. J.} {\bf C}}
\def\EJA{{\em Eur. Phys. J.} {\bf A}}
\def\CPC{\em Comp. Phys. Commun.}

\begin{titlepage}

\begin{flushleft}
{\tt DESY 14-035    \hfill    ISSN 0418-9833} \\
{\tt March 2014}                              \\
\end{flushleft}
%\vspace{3mm}


\vspace*{10mm}

\rm

\begin{center}
\begin{Large}
{\boldmath \bf
   Measurement of Feynman-$x$ Spectra of Photons and Neutrons
   in the Very Forward Direction in  Deep-Inelastic Scattering at HERA
}

\vspace{2cm}
\rm
H1 Collaboration

\end{Large}
\end{center}

\vspace{2cm}

%
\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
in deep-inelastic $ep$ scattering at HERA are presented 
as a function of the Feynman variable $x_F$ and 
of the centre-of-mass energy of the
virtual photon-proton system $W$.
%
The data are taken with the H1 detector in the years 2006 and 2007 and
correspond to an integrated luminosity of $131~\mathrm{pb}^{-1}$.
The measurement is restricted to photons and neutrons in the
pseudorapidity range  $\eta>7.9$ and covers the range of
negative four momentum transfer squared at the positron vertex
$6<Q^2<100$~GeV$^2$, of inelasticity $0.05<y<0.6$ and of
$70<W<245~$GeV.
To test the Feynman scaling hypothesis the
$W$ dependence of the $x_F$ dependent cross sections is investigated.
Predictions of deep-inelastic scattering models and of models for hadronic 
interactions of high energy cosmic rays are compared to the measured cross
sections.

\end{abstract}

\vspace*{10mm}

\it \normalsize
\begin{center} submitted to {\it Eur.Phys.J.C}
\end{center}

\end{titlepage}
%----------------------------------------------------------------------------------
\begin{flushleft}
%-- H1AUTS Author list by names 
%-- Status: Wed Jan  8 13:54:33 CET 2014  Number of authors = 147 

\noindent
V.~Andreev$^{22}$,             %LPI -PD        8/88            Andreev             
A.~Baghdasaryan$^{34}$,        %YERE-PD        09/03           Baghdasaryana       
K.~Begzsuren$^{31}$,           %ULBA-PD        04/06           Begzsuren           
A.~Belousov$^{22}$,            %LPI -HON       12/12           Belousov            
P.~Belov$^{10}$,               %DESY-LEFT      08/13           Belov               
V.~Boudry$^{25}$,              %ECPL-PD        1/93            Boudry              
G.~Brandt$^{46}$,              %OXFU-PD        03/07           Brandt              
M.~Brinkmann$^{10}$,           %DESY-LEFT      07/13           Brinkmann           
V.~Brisson$^{24}$,             %ORSA-HON       01/12           Brisson             
D.~Britzger$^{10}$,            %DESY-PD        06/13           Britzger            
A.~Buniatyan$^{13}$,           %HDB1-PD        12/95           Buniatyan           
A.~Bylinkin$^{21,43}$,         %ITEP-ST        09/12           Bylinkin            
L.~Bystritskaya$^{21}$,        %ITEP-PD        05/99           Bystritskaya        
A.J.~Campbell$^{10}$,          %DESY-PD        10/84           Campbella           
K.B.~Cantun~Avila$^{20}$,      %MEX1-PD        05/13           Cantunavila         
F.~Ceccopieri$^{3}$,           %BRUX-LEFT      10/13           Ceccopieri          
K.~Cerny$^{28}$,               %PRG2-PD        11/08           Cernyk              
V.~Chekelian$^{23}$,           %MPIM-PD        01/90           Chekelian           
J.G.~Contreras$^{20}$,         %MEX1-PROF      03/98           Contreras           
J.B.~Dainton$^{17}$,           %LIVE-PROF      01/94           Dainton             
K.~Daum$^{33,38}$,             %WUPP-PROF      06/90           Daum                
C.~Diaconu$^{19}$,             %MARS-PD        09/96           Diaconu             
M.~Dobre$^{4}$,                %BUCH-PD        12/12           Dobre               
V.~Dodonov$^{10}$,             %DESY-PD        04/98           Dodonov             
A.~Dossanov$^{11,23}$,         %MPIM-LEFT      06/13           Dossanov            
G.~Eckerlin$^{10}$,            %DESY-PD        8/88            Eckerlin            
S.~Egli$^{32}$,                %PSI -PD        01/10           Egli                
E.~Elsen$^{10}$,               %DESY-PROF      01/06           Elsen               
L.~Favart$^{3}$,               %BRUX-PROF      10/99           Favart              
A.~Fedotov$^{21}$,             %ITEP-PD        8/88            Fedotov             
J.~Feltesse$^{9}$,             %SACL-HON       11/06           Feltesse            
J.~Ferencei$^{15}$,            %KOSI-PD        08/88           Ferencei            
M.~Fleischer$^{10}$,           %DESY-PD        07/95           Fleischer           
A.~Fomenko$^{22}$,             %LPI -PD        02/89           Fomenko             
E.~Gabathuler$^{17}$,          %LIVE-HON       01/12           Gabathulere         
J.~Gayler$^{10}$,              %DESY-HON       02/05           Gayler              
S.~Ghazaryan$^{10}$,           %DFLC-PD        09/09           Ghazaryan           
A.~Glazov$^{10}$,              %DESY-PD        01/04           Glazov              
L.~Goerlich$^{6}$,             %CRAC-PD        8/88            Goerlich            
N.~Gogitidze$^{22}$,           %LPI -PD        05/91           Gogitidze           
M.~Gouzevitch$^{10,39}$,       %DESY-PD        10/11           Gouzevitch          
C.~Grab$^{36}$,                %ZUTH-PROF      08/07           Grab                
A.~Grebenyuk$^{3}$,            %BRUX-PD        04/12           Grebenyuk           
T.~Greenshaw$^{17}$,           %LIVE-PD        8/88            Greenshaw           
G.~Grindhammer$^{23}$,         %MPIM-HON       09/11           Grindhammer         
D.~Haidt$^{10}$,               %DESY-HON       03/04           Haidt               
R.C.W.~Henderson$^{16}$,       %LANC-PROF      10/04           Henderson           
M.~Herbst$^{14}$,              %HDB2-PD        11/11           Herbst              
J.~Hladk\`y$^{27}$,            %PRAG-HON       12/02           Hladky              
D.~Hoffmann$^{19}$,            %MARS-PD        07/00           Hoffmann            
R.~Horisberger$^{32}$,         %PSI -PD        01/10           Horisberger         
T.~Hreus$^{3}$,                %BRUX-PD        10/08           Hreus               
F.~Huber$^{13}$,               %HDB1-PD        01/13           Huberf              
M.~Jacquet$^{24}$,             %ORSA-PD        09/96           Jacquet             
X.~Janssen$^{3}$,              %ANTW-PD        02/03           Janssenx            
H.~Jung$^{10,3}$,              %DESY-PROF      12/99           Jungh               
M.~Kapichine$^{8}$,            %JINR-PD        3/97            Kapichine           
C.~Kiesling$^{23}$,            %MPIM-HON       09/11           Kiesling            
M.~Klein$^{17}$,               %LIVE-PROF      12/06           Klein               
C.~Kleinwort$^{10}$,           %DESY-PD        8/88            Kleinwort           
R.~Kogler$^{11}$,              %HAM2-PD        12/10           Kogler              
P.~Kostka$^{35}$,              %ZEUT-PD        10/87           Kostka              
J.~Kretzschmar$^{17}$,         %LIVE-PD        03/08           Kretzschmar         
K.~Kr\"uger$^{10}$,            %DESY-PD        01/04           Kruegerk            
M.P.J.~Landon$^{18}$,          %QMWC-PD        10/83           Landon              
W.~Lange$^{35}$,               %ZEUT-PD        8/88            Lange               
P.~Laycock$^{17}$,             %LIVE-PD        11/03           Laycock             
A.~Lebedev$^{22}$,             %LPI -HON       12/12           Lebedev             
S.~Levonian$^{10}$,            %DESY-PD        08/88           Levonian            
K.~Lipka$^{10,42}$,            %DESY-PD        06/01           Lipka               
B.~List$^{10}$,                %DESY-PD        11/99           Listb               
J.~List$^{10}$,                %DFLC-PD        10/00           Listj               
B.~Lobodzinski$^{10}$,         %DESY-PD        01/06           Lobodzinski         
L.~Lytkin$^{8}$,               %JINR - ADD
E.~Malinovski$^{22}$,          %LPI -PD        01/89           Malinovskie         
H.-U.~Martyn$^{1}$,            %AAC1-HON       01/12           Martyn              
S.J.~Maxfield$^{17}$,          %LIVE-PD        8/88            Maxfield            
A.~Mehta$^{17}$,               %LIVE-PROF      06/00           Mehta               
A.B.~Meyer$^{10}$,             %DESY-PD        10/97           Meyeran             
H.~Meyer$^{33}$,               %WUPP-HON       12/00           Meyerhi             
J.~Meyer$^{10}$,               %DESY-PROF      10/84           Meyerj              
S.~Mikocki$^{6}$,              %CRAC-PROF      10/10           Mikocki             
A.~Morozov$^{8}$,              %JINR-PD        06/99           Morozova            
K.~M\"uller$^{37}$,            %ZUER-PD        11/94           Muellerk            
Th.~Naumann$^{35}$,            %ZEUT-PROF      01/05           Naumannt            
P.R.~Newman$^{2}$,             %BIRM-PROF      10/09           Newman              
C.~Niebuhr$^{10}$,             %DESY-PD        3/93            Niebuhr             
G.~Nowak$^{6}$,                %CRAC-PROF      02/12           Nowakg              
K.~Nowak$^{11}$,               %HAM2-LEFT      02/13           Nowakk              
J.E.~Olsson$^{10}$,            %DESY-HON       09/10           Olsson              
D.~Ozerov$^{10}$,              %DESY-PD        09/08           Ozerov              
P.~Pahl$^{10}$,                %DESY-ST        10/08           Pahl                
C.~Pascaud$^{24}$,             %ORSA-HON       09/05           Pascaud             
G.D.~Patel$^{17}$,             %LIVE-PD        8/88            Patel               
E.~Perez$^{9,40}$,             %SACL-PD        10/07           Perez               
A.~Petrukhin$^{10}$,           %DESY-PD        09/09           Petrukhin           
I.~Picuric$^{26}$,             %PODG-PROF      12/07           Picuric             
H.~Pirumov$^{10}$,             %DESY-PD        04/13           Pirumov             
D.~Pitzl$^{10}$,               %DESY-PD        8/88            Pitzl               
R.~Pla\v{c}akyt\.{e}$^{10,42}$, %DESY-PD        10/06           Placakyte           
B.~Pokorny$^{28}$,             %PRG2-ST        10/09           Pokorny             
R.~Polifka$^{28,44}$,          %PRG2-PD        05/11           Polifka             
B.~Povh$^{12}$,                %MPIH-PD - ADDED
V.~Radescu$^{10,42}$,          %DESY-PD        10/06           Radescu             
N.~Raicevic$^{26}$,            %PODG-PD        06/02           Raicevic            
T.~Ravdandorj$^{31}$,          %ULBA-PROF      09/00           Ravdandorj          
P.~Reimer$^{27}$,              %PRAG-PD        07/90           Reimer              
E.~Rizvi$^{18}$,               %QMWC-PROF      03/04           Rizvi               
P.~Robmann$^{37}$,             %ZUER-PD        8/88            Robmann             
R.~Roosen$^{3}$,               %BRUX-HON       03/12           Roosen              
A.~Rostovtsev$^{21}$,          %ITEP-PROF      02/11           Rostovtsev          
M.~Rotaru$^{4}$,               %BUCH-PD        06/11           Rotaru              
S.~Rusakov$^{22}$,             %LPI -HON       01/12           Rusakov             
D.~\v S\'alek$^{28}$,          %PRG2-PD        10/10           Salek               
D.P.C.~Sankey$^{5}$,           %RAL -PD        01/90           Sankey              
M.~Sauter$^{13}$,              %HDB1-PD        10/09           Sauter              
E.~Sauvan$^{19,45}$,           %MARS-PD        11/1            Sauvan              
S.~Schmitt$^{10}$,             %DESY-PD        09/99           Schmittst           
L.~Schoeffel$^{9}$,            %SACL-PROF      10/10           Schoeffel           
A.~Sch\"oning$^{13}$,          %HDB1-PROF      01/09           Schoening           
H.-C.~Schultz-Coulon$^{14}$,   %HDB2-PROF      01/04           Schultzcoulon       
F.~Sefkow$^{10}$,              %DFLC-PD        09/99           Sefkow              
S.~Shushkevich$^{10}$,         %DESY-PD        08/11           Shushkevich         
Y.~Soloviev$^{10,22}$,         %LPI -PD        08/89           Soloviev            
P.~Sopicki$^{6}$,              %CRAC-ST        09/07           Sopicki             
D.~South$^{10}$,               %DESY-PD        05/03           South               
V.~Spaskov$^{8}$,              %JINR-PD        12/97           Spaskov             
A.~Specka$^{25}$,              %ECPL-PROF      09/05           Specka              
M.~Steder$^{10}$,              %DESY-PD        09/08           Steder              
B.~Stella$^{29}$,              %ROME-HON       11/10           Stella              
U.~Straumann$^{37}$,           %ZUER-PROF      09/95           Straumann           
T.~Sykora$^{3,28}$,            %ANTW-PD        01/06           Sykora              
P.D.~Thompson$^{2}$,           %BIRM-PD        08/99           Thompsonp           
D.~Traynor$^{18}$,             %QMWC-PD        10/01           Traynor             
P.~Tru\"ol$^{37}$,             %ZUER-HON       10/06           Truoel              
I.~Tsakov$^{30}$,              %SOFI-PROF      03/12           Tsakov              
B.~Tseepeldorj$^{31,41}$,      %ULBA-PROF      06/06           Tseepeldorj         
J.~Turnau$^{6}$,               %CRAC-PROF      09/95           Turnau              
A.~Valk\'arov\'a$^{28}$,       %PRG2-PD        8/88            Valkarova           
C.~Vall\'ee$^{19}$,            %MARS-PD        06/86           Vallee              
P.~Van~Mechelen$^{3}$,         %ANTW-PROF      01/06           Vanmechelen         
Y.~Vazdik$^{22}$,              %LPI -PD        8/88            Vazdik              
D.~Wegener$^{7}$,              %DORT-PROF      01/06           Wegener             
E.~W\"unsch$^{10}$,            %DESY-PD        8/88            Wuensch             
J.~\v{Z}\'a\v{c}ek$^{28}$,     %PRG2-PROF      10/05           Zacek               
Z.~Zhang$^{24}$,               %ORSA-PD        10/92           Zhang               
R.~\v{Z}leb\v{c}\'{i}k$^{28}$, %PRG2-ST        03/12           Zlebcik             
H.~Zohrabyan$^{34}$,           %YERE-PD        11/02           Zohrabyan           
and
F.~Zomer$^{24}$                %ORSA-PROF      09/06           Zomer          

%\newpage

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

\bigskip\noindent
 $ ^{38}$ Also at Rechenzentrum, Universit\"at Wuppertal,
          Wuppertal, Germany \\
 $ ^{39}$ Also at IPNL, Universit\'e Claude Bernard Lyon 1, CNRS/IN2P3,
          Villeurbanne, France \\
 $ ^{40}$ Also at CERN, Geneva, Switzerland \\
 $ ^{41}$ Also at Ulaanbaatar University, Ulaanbaatar, Mongolia \\
 $ ^{42}$ Supported by the Initiative and Networking Fund of the
          Helmholtz Association (HGF) under the contract VH-NG-401 and S0-072 \\
 $ ^{43}$ Also at Moscow Institute of Physics and Technology, Moscow, Russia \\
 $ ^{44}$ Also at  Department of Physics, University of Toronto,
          Toronto, Ontario, Canada M5S 1A7 \\
 $ ^{45}$ Also at LAPP, Universit\'e de Savoie, CNRS/IN2P3,
          Annecy-le-Vieux, France \\
 $ ^{46}$ Department of Physics, Oxford University,
          Oxford, UK$^{ b}$ \\


\bigskip \noindent
 $ ^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 VEGA SR grant no. 2/7062/ 27 \\
 $ ^f$ Supported by the Ministry of Education of the Czech Republic
      under the projects  LC527, INGO-LA09042 and
      MSM0021620859 \\
 $ ^g$ Supported by the Swiss National Science Foundation \\
 $ ^h$ Supported by  CONACYT,
      M\'exico, grant 48778-F \\
 $ ^i$ Russian Foundation for Basic Research (RFBR), grant no 1329.2008.2
      and Rosatom \\
 $ ^j$ Supported by the Romanian National Authority for Scientific Research
      under the contract PN 09370101 \\
 $ ^k$ Partially Supported by Ministry of Science of Montenegro,
      no. 05-1/3-3352
}


\end{flushleft}

\newpage

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


\section{Introduction}

\noindent
%--------------------------------------------------------------------------------------
Measurements of particle production at very small polar
angles with respect to the proton beam direction ({forward direction})
in positron-proton collisions are
important inputs for the theoretical
understanding of proton fragmentation mechanisms.
%
Forward particle measurements are also valuable
for high energy cosmic ray experiments, as they 
provide important new constraints for high energy air shower models
\cite{Engel:1998hf,Bunyatyan:2009zz}.
%

The H1 and ZEUS experiments at the $ep$ collider HERA have
studied the production of forward baryons (protons and neutrons) and photons, 
which carry  a large fraction of the longitudinal
momentum of the incoming proton
\cite{Adloff:1998yg,Chekanov:2002pf,Chekanov:2007tv,Chekanov:2008tn,Aaron:2010ab,Aaron:2011pe}.
%
These analyses have demonstrated that models of deep-inelastic scattering (DIS) 
are able to reproduce the forward baryon measurements if 
contributions from different production mechanisms are considered,
such as string fragmentation, pion exchange, diffractive dissociation and elastic
scattering of the proton~\cite{Chekanov:2008tn,Aaron:2010ab}.
The forward photon production rate, however, is overestimated by the models
by $50$ to $70$\%~\cite{Aaron:2011pe}.
%
The measurements  % of forward particles 
also confirm the 
hypothesis of limiting fragmentation~\cite{Benecke:1969sh,Chou:1994dh},
according to which, in the high-energy limit, the cross
section for the inclusive production of particles in the target
fragmentation region is independent of the incident projectile energy.  

Measurements in the DIS regime provide a possibility to investigate
the process at different centre-of-mass (CM) energies
of the  virtual photon-proton system, $W$, within the same experiment.
The studies of the energy dependence of particle production 
allow a test of the Feynman scaling \cite{Feynman:1969ej} hypothesis,
according to which particle production 
is expected to show a scaling behaviour, 
i.e. independence of the CM energy in terms of the Feynman-$x$ variable,
$x_F=2p_{||}^*/W$.
Here  $p_{||}^*$ is the longitudinal momentum of the particle
in the virtual photon-proton CM frame with respect to the direction 
of the beam proton.
%
In several previous measurements Feynman scaling was found to be 
violated in the fragmentation process 
in the central rapidity region 
\cite{Arnison:1982ed,Arnison:1982rm,Banner:1984wh,Bernard:1985kh,Alner:1986xu,Alner:1987wb,Alexopoulos:1988na,Abe:1989td,Aaron:2007ds,Abramowicz:2010rz}.
%
On the contrary,  no  sizable violation of
Feynman scaling has been observed in the  target fragmentation region
in  $pp$ and $p\bar{p}$ collisions
by comparing the  $\pi^0$ production cross sections
at the SPS collider \cite{Pare:1989mr} with
$\pi^{\pm}$ 
measurements at the
ISR \cite{Albrow:1973kj,Albrow:1973ee,Capiluppi:1973fz,Capiluppi:1974rt}.
However, these conclusions are debated \cite{Ohsawa:1992dj}
and the scarcity of other experimental forward particle
production data motivates
further studies of forward particle production.

In this paper the production of forward neutrons and photons in
DIS is studied 
as a function of $x_F$\footnote{In the kinematic 
range of this measurement the variable $x_F$ is numerically almost 
equal to the longitudinal momentum fraction $x_L$ used in the 
previous publications \cite{Adloff:1998yg,Chekanov:2002pf,Chekanov:2007tv,Chekanov:2008tn,Aaron:2010ab,Aaron:2011pe}.
There, $x_L$ was defined as $x_L=E_{n,\gamma}/E_p$,
where $E_p$, $E_{n}$ and $E_{\gamma}$ are the energies of
the proton beam, the forward neutron and the forward photon in the
laboratory frame, respectively.} and $W$.
This is the first direct experimental test of Feynman scaling
for  photons and neutrons produced in the very forward direction.
%
Predictions from different DIS and different cosmic ray (CR) 
hadronic interaction 
Monte Carlo (MC) models are compared to the results.
The simultaneous measurement of  forward neutrons and photons 
provides a useful input for
MC model development also because of  the respective different 
production mechanisms:
forward photons almost exclusively originate from
decays of neutral mesons produced in the fragmentation of the proton remnant 
(Figure~\ref{diagram1}a),
while forward neutrons  are produced also
via a colour singlet exchange process
(Figure~\ref{diagram1}b).



\begin{figure}[h]
\epsfig{file=desy-14-035.fig1a.eps,height=62mm}
\hspace*{15mm}
\epsfig{file=desy-14-035.fig1b.eps,height=62mm}
\caption{ (a) Generic diagram for forward photon or neutron production 
$ep\rightarrow e'\gamma X$, $ep\rightarrow e'nX$ in deep-inelastic scattering.
(b) Diagram of forward neutron production via pion exchange.}

\vspace*{-79mm}
\hspace*{25mm}{\large \bf (a)} \hspace*{78mm} {\large \bf (b)}

\vspace*{77mm}  
\label{diagram1}
\end{figure}


The neutrons and photons studied here are produced at polar angles
below $0.75$~mrad and are measured in the Forward Neutron 
Calorimeter (FNC) of the H1 detector.
The data used in this analysis were collected with the H1 detector at
HERA in the years 2006 and 2007 and correspond to an integrated
luminosity of $131~ \rm pb^{-1}$. During this period 
HERA collided positrons and protons with energies of 
$E_e=27.6~\GeV$ and $E_p=920~\GeV$, respectively, corresponding to a
centre-of-mass energy of $\sqrt{s}=319~\GeV$.



\section{Experimental Procedure and Data Analysis}
%--------------------------------------------------------------------------------------
\subsection{H1 main detector}

A detailed description of the H1 detector can be found elsewhere
\cite{Abt:1996hi,Abt:1996xv,Appuhn:1996na,Pitzl:2000wz,Andrieu:1993kh,Laycock:2012xg}.
Only the detector components relevant to this analysis are briefly
described here.
 The origin of the right-handed
 H1 coordinate system is the nominal $e^+p$ interaction point. 
 The direction of the proton beam defines  the positive $z$--axis; 
 the polar angle $\theta$ is measured with respect to  this axis. 
 Transverse momenta are measured in the $x$--$y$ plane.
 The pseudorapidity is defined by  $\eta =  -\ln{[\tan(\theta/2)]}$ 
 and is measured in the laboratory frame.

The 
%$ep$ 
interaction region is surrounded by a two-layer silicon
strip detector and large concentric drift chambers,
operated inside a $1.16$~T solenoidal magnetic field.
Charged particle momenta are measured  in the
angular range $15\deg<\theta<165\deg$.
%
The forward tracking detector is used to supplement track reconstruction 
in the region $7\deg<\theta<30\deg$
and improves the hadronic final state 
reconstruction of forward going low momentum particles.
%
The tracking system is surrounded by a finely segmented
liquid argon (LAr) calorimeter, which covers the polar angle range
\mbox{$4\deg<\theta<154\deg$} with full azimuthal acceptance.  The LAr
calorimeter consists of an electromagnetic section with lead absorber
and a hadronic section with steel absorber.  The total depth of the
LAr calorimeter ranges from $4.5$ to $8$ hadronic interaction lengths $\lambda$.
 The absolute electromagnetic energy scale is
 known with a precision of $1\%$, while the absolute hadronic energy
 scale is known for the present data with a precision of $4\%$.

The backward region ($153\deg<\theta<177.8\deg$) is covered by a
lead/scintillating-fibre calorimeter called the SpaCal;
its main purpose is the detection of scattered positrons.  The polar angle of
the positron is measured with a precision of $1$~mrad.  The energy
resolution for
positrons is $\sigma(E)/E\approx
7.1\%/\sqrt{E[\GeV]}\oplus 1\%$ \cite{Nicholls:1995di} and the energy 
scale uncertainty is
less than $1\%$.  The hadronic energy scale in the SpaCal is known with
a precision of $7\%$.

 The luminosity is determined from the rate of the elastic QED Compton process 
 with the electron and the photon detected in the SpaCal calorimeter, 
 and the rate of DIS events measured in the SpaCal calorimeter 
 \cite{Aaron:2012kn}.


%--------------------------------------------------------------------------------------
\subsection{Forward detector for neutral particles}
\label{sec:fncdet}

Neutral particles produced at very small polar angles with respect to the proton
beam direction can be detected in the FNC, 
which is situated 
at a polar angle of $0\deg$ at $106$~m from the interaction point.
A detailed description
of the detector is given in \cite{Aaron:2010ab,Aaron:2011pe}.  
The FNC is a lead--scintillator sandwich calorimeter. It consists of two 
longitudinal sections: the Preshower Calorimeter 
with a  length corresponding to about $60$ radiation lengths
or $1.6 \lambda$
 and the
Main Calorimeter with a total length of $8.9 \lambda$.
%
The acceptance of the FNC is defined by the aperture of the HERA
beam-line magnets and is limited to scattering angles of
$\theta\lsim 0.8$~mrad with approximately $30\%$ azimuthal coverage. 

The longitudinal segmentation of the FNC allows an efficient
discrimination of photons from neutrons.  
Photons are absorbed completely in the Preshower Calorimeter, while
neutrons  have a significant fraction of their energy 
deposited in the Main Calorimeter.
Therefore, energy deposits in the FNC, which are contained in the
Preshower Calorimeter
with no energy deposits above the noise level in the Main Calorimeter,
are classified as electromagnetic clusters.
%
According to the Monte Carlo simulation about $98\%$ of all reconstructed 
photon and neutron candidates originate from  generated photons and neutrons,
respectively. 
%
Due to the relatively large size of the FNC readout modules in
combination with the small  geometrical acceptance window, 
two or more particles entering the FNC are reconstructed as 
a single cluster.
In the MC simulation about $1\%$ of all
hadronic clusters in the FNC associated with neutrons
are overlapped with  a photon, which was scattered
within the FNC acceptance together with the neutron.
%
At lower energies the electromagnetic clusters reconstructed in the FNC
mainly originate from single photons. At higher measured energies there is 
a significant contribution from two photons, with the 
fraction of two-photon events increasing from $0.5\%$ at $100$~GeV to 
$10\%$ at about $450$~GeV and to $80\%$ at $900$~GeV.
%
The two photons typically originate from the decay of the same meson.


The absolute electromagnetic and hadronic energy scales of the FNC are known to 
$5\%$~\cite{Aaron:2011pe} and $2\%$~\cite{Aaron:2010ab} precision, respectively.
%
The energy resolution of the FNC calorimeter for electromagnetic
showers is $\sigma(E)/E \approx 20\%/\sqrt{E~[\rm GeV]} \oplus 2\%$
and for hadronic showers $\sigma(E)/E \approx 63\%/\sqrt{E~[\rm GeV]} \oplus 3\%$,
as determined in test beam measurements~\cite{Uraev}.
%
The spatial resolution is
$\sigma(x,y)\approx 10\rm cm/\sqrt{E~[\rm GeV]} \oplus 0.6~\rm cm$ for 
hadronic showers starting in the Main Calorimeter.
A better spatial resolution of about $2\rm~ mm$ is achieved for the
electromagnetic showers and for those hadronic showers 
which start in the Preshower Calorimeter.


%------------------------------------------------------------------------
\subsection{Cross Section Definition}

The kinematics of semi-inclusive forward photon and neutron production are shown
in Figure~\ref{diagram1}a, where the four-vectors of the incoming
and outgoing particles and of the exchanged virtual photon $\gamma^*$
are indicated.
%
This measurement is restricted to the DIS kinematic range, determined by the 
photon virtuality $6<Q^2<100~\GeV^2$ and inelasticity $0.05<y<0.6$.
They are defined as
\begin{equation}
Q^2=-q^2\; , \hspace{3em}  
y=\frac{p \cdot q}{p \cdot k}\;,
\end{equation}
where $p$, $k$ and $q$ are the four-momenta of the incident proton,
the incident positron and the virtual photon, respectively.
The  CM  energy of the virtual photon-proton system, $W$, is
related to $Q^2$ and $y$ as $W \approx \sqrt{ys-Q^2}$, 
where $s$ is the squared
total CM energy of the positron-proton system.
The present analysis is restricted to the range $70<W<245~\rm GeV$.
 

 The analysis is performed in the pseudorapidity range $\eta>7.9$
 for forward neutrons and  photons.
 The pseudorapidity range $\eta>7.9$ corresponds to polar
 angles $\theta<0.75$~mrad\kd{, the geometrical acceptance limit of the FNC.}
%
 In the virtual photon-proton CM frame \kd{the $x_F$ of the neutron and photon are restricted to
  $0.1<x_F<0.94$ and   $0.1<x_F<0.7$, respectively.}
  \kd{These requirements limit the transverse momentum of neutrons to $0<p_T^*<0.6~\rm GeV$  and of photons to $0<p_T^*<0.4~\rm GeV$.}
  %the neutron transverse momentum $p_T^*$ and the neutron $x_F$ are restricted to the range $0<p_T^*<0.6~\rm GeV$ and $0.1<x_F<0.94$, respectively.
% For the forward photons measurement $p_T^*$ and $x_F$  are defined for the most energetic photon in the pseudorapidity range  $\eta>7.9$ and are  restricted to the range $0<p_T^*<0.4~\rm GeV$ and $0.1<x_F<0.7$.  
 The requirement that $x_F$ is below $0.7$ \kd{for photons} ensures that the 
 electromagnetic clusters reconstructed in the FNC mainly  
 originate from single photons, according to MC predictions.

 The kinematic phase space of the measurements is
 summarised in Table~\ref{tab:selection}.
%
  Cross sections of neutrons and photons produced in the forward
  direction, normalised to the inclusive DIS cross section, 
%  $\rm \frac{1}{\sigma_{DIS}}\frac{d\sigma}{dx_F}$, 
  $ 1/\sigma_{DIS}~{\rm d}\sigma/{\rm d}x_F$, 
  are
  determined differentially in $x_F$ in three ranges of $W$.
  In addition, the cross section ratios integrated over $x_F$,
%  $\rm \frac{\sigma^{\gamma,n}_{DIS}}{\sigma_{DIS}}$, 
  $ \sigma^{\gamma,n}_{DIS}/\sigma_{ DIS}$, 
  are measured  as a function of $W$.


\begin{table}[h]
        \centering
        \renewcommand{\arraystretch}{1.50}
  \begin{tabular}{| l | l |}
                        \hline 
%
        \multicolumn{2}{|c|}{\bf NC DIS Selection}  \\ \hline
        \multicolumn{2}{|c|}{$6 < Q^2 < 100$~GeV$^2$}  \\
        \multicolumn{2}{|c|}{$0.05 < y < 0.6$} \\ 
        \multicolumn{2}{|c|}{$70 < W < 245$~GeV} \\ \hline
        {\bf Forward photons} & {\bf Forward neutrons} \\ \hline
%       \multicolumn{2}{|c|}{$\eta_{lab}>7.9$}  \\ 
        $\eta>7.9$ & $\eta>7.9$ \\
        $0.1<x_F<0.7$ & $0.1<x_F<0.94$ \\
        $0<p_T^*<0.4$~GeV & $0<p_T^*<0.6$~GeV \\ \hline
        \multicolumn{2}{|c|}{\bf \boldmath $W$ ranges for cross sections 
        {\large $ \frac{1}{\sigma_{DIS}}\frac{{\rm d}\sigma}{{\rm d}x_F}$}}  \\
        \multicolumn{2}{|c|}{$~70 < W < 130$~GeV} \\ 
        \multicolumn{2}{|c|}{$130 < W < 190$~GeV} \\ 
        \multicolumn{2}{|c|}{$190 < W < 245$~GeV} \\ \hline
%
  \end{tabular}
   \vspace{\baselineskip} 
        \caption{Definition of the kinematic phase space of the measurements.}
  \label{tab:selection}
\end{table}



%--------------------------------------------------------------------------------------
\subsection{Event selection }

The  data selection and analysis procedures are similar to those 
described in previous publications using the FNC \cite{Aaron:2010ab,Aaron:2011pe}.
%
The data sample of this analysis was collected using
 triggers which require the scattered positron to be measured in the SpaCal.
 The trigger efficiency is about 96\% for the analysis phase space
 as determined from data using an independently triggered data sample.
The selection of DIS events is based on the identification
of the scattered positron as the most energetic, isolated compact calorimetric 
deposit
in the SpaCal  with an
energy $E_e'>11~\GeV$ and a polar angle $156\deg<\theta_e'<175\deg$.
The \mbox{$z$-coordinate} of the primary event vertex is required to be
within $\pm 35$~cm of the nominal interaction point.
%
The hadronic final state is reconstructed using an energy flow algorithm 
which combines charged particles measured in the trackers with information 
from the SpaCal and LAr calorimeters~\cite{Peez:2003zd,Hellwig:2004yp}.
%
To suppress events with hard initial state  QED
radiation, as well as
events originating  from  non-$ep$ interactions,
the quantity $\sum{E - p_z}$, summed over
all reconstructed particles including the positron, is required to lie
between  $35~\GeV$ and $70~\GeV$.
%
This cut also efficiently removes events from photoproduction processes, 
where the
positron is scattered into the backward beam-pipe and a particle from
the hadronic final state fakes the positron signature in the SpaCal.
%
The kinematic variables $Q^2$ and $y$ are reconstructed using a technique
which optimises the resolution throughout the
measured $y$ range, exploiting information from both the scattered
positron and the hadronic final state \cite{Adloff:1997sc}.
Events are restricted to the range
$6<Q^2<100~\GeV^2$ and $0.05<y<0.6$.
The DIS data sample contains about 9.3 million events.

A subsample of events containing forward photons or neutrons is 
selected by requiring 
either an electromagnetic or a hadronic cluster in
the FNC with a pseudorapidity above $7.9$ and an energy above $92~\rm GeV$. 
The data sample, called `FNC sample' in the following,
contains about $83 ,000$ events 
with photons and $230,000$ events with neutrons.


%---------------------------------------------------------------------------
\subsection{Monte Carlo simulations and corrections to the data}
\label{sec:mc}

Monte Carlo simulations are used to correct the data for the
effects of detector acceptance, inefficiencies, QED
radiation from the positron and migrations between
measurement bins
%measurement intervals
due to the finite detector resolution.
All generated events are passed through
a GEANT3~\cite{Brun:1978fy} based simulation of the H1 apparatus and
are subject to the same reconstruction and analysis chain as the data.

The DJANGOH~\cite{Charchula:1994kf}
%~\cite{Schuler:1991yg}
program is used to generate 
inclusive DIS events.  It is based on leading order
electroweak cross sections and takes into account QCD effects up to
order $\alpha_s$.
%
Higher order QCD effects are simulated using
leading log parton showers as implemented in LEPTO ~\cite{Ingelman:1996mq},
or using the Colour Dipole Model (CDM) as implemented in ARIADNE
\cite{Lonnblad:1992tz}.  Subsequent hadronisation effects are modelled
using the Lund string fragmentation model as implemented in JETSET
%%\cite{Sjostrand:1986hx,Bengtsson:1987kr,Sjostrand:1992mh}.
\cite{Andersson:1983ia,Sjostrand:1995iq}.  
Higher order electroweak processes are
simulated using an interface to HERACLES~\cite{Kwiatkowski:1990es}.
%
%
The LEPTO program optionally includes 
the simulation of soft colour interactions
(SCI)~\cite{Edin:1995gi}, in which the production of
diffraction-like configurations is enhanced via non-perturbative
colour rearrangements between the outgoing partons.
The SCI option in LEPTO is used for the simulation of forward photons
\kd{but not for neutrons because it has been shown~\cite{Aaron:2010ab} that the description of the forward neutron data is poor in this case.}
%
%
For the DJANGOH MC simulations the H1PDF~2009
parameterisation~\cite{Aaron:2009kv} of the parton distributions in the
proton is used.  In the following, the DJANGOH predictions based on LEPTO
and ARIADNE are denoted as LEPTO and CDM, respectively.
%
 In all DJANGOH simulations  forward particles originate exclusively
 from the
 hadronisation of the proton remnant and forward photons 
 are therefore mainly produced from the decay of $\pi^0$ mesons.

RAPGAP~\cite{Jung:1993gf} is a general purpose event generator for 
inclusive and diffractive $ep$ interactions.  
Higher order QCD effects are simulated
using parton showers and the final state hadrons are
obtained via the Lund string model.  As in DJANGOH higher order electroweak
processes are simulated using an
interface to HERACLES~\cite{Kwiatkowski:1990es}.
%
In the version denoted below as RAPGAP-$\pi$, the program
simulates exclusively the scattering of virtual or real photons
off an exchanged pion (Figure~\ref{diagram1}b). 
In this model the cross section for $ep$ scattering to the final state 
$nX$ takes the form
\begin{equation}
  \der\sigma (e p\rightarrow e'nX) =f_{\pi^+/p}(x_L,t)\cdot
  \der\sigma(e\pi^+\rightarrow e'X).
\label{crsec}
\end{equation}
Here $x_L$ is the longitudinal 
momentum fraction and $t$ is the squared four-momentum transfer
between the incident proton and the final state neutron; 
$f_{\pi^+/p}(x_L,t)$ represents the pion flux associated with the 
splitting of a proton into a $\pi^+ n$ system 
and $\der\sigma(e\pi^+\rightarrow e'X)$ is the cross section of
the positron-pion  interaction.
%
%
 There are several parameterisations of the pion flux
 \cite{Bishari:1972tx,Holtmann:1994rs,Kopeliovich:1996iw,Przybycien:1996zb,Szczurek:1997cw} and the one used here is taken from 
\cite{Holtmann:1994rs}. 
The details of the pion flux parameterisation
are described in \cite{Aaron:2010ab}. 
Using other parameterisations of the pion flux
affects mainly the absolute normalisation by up to $30~\%$.
%

As was shown in \cite{Aaron:2010ab},  the best description of the 
forward neutron data is  achieved 
by a combination of events with neutrons originating from pion exchange, 
as simulated by RAPGAP-$\pi$, and events with neutrons from 
proton remnant fragmentation, simulated by DJANGOH \kd{without SCI}. 
RAPGAP-$\pi$ mainly contributes at high neutron 
energies, while DJANGOH is significant at low energies.
%
In \cite{Aaron:2010ab} the contributions of 
 RAPGAP-$\pi$ and CDM were added using weighting factors
$0.65$ and $1.2$ for the respective predictions.
In the present analysis the best description of the neutron energy
distribution is obtained by the combination of RAPGAP-$\pi$ and CDM
using the respective weights $0.6$ and $1.4$,
or by the combination of RAPGAP-$\pi$ and LEPTO
using the respective weights $0.6$ and $0.7$.
The difference between the weighting factors for the 
combination of RAPGAP-$\pi$ and CDM in this analysis and 
in~\cite{Aaron:2010ab} is 
due to the different neutron energy range 
and the resulting different neutron energy dependence
in the two analyses.
Compared to  \cite{Aaron:2010ab} the current analysis is extended 
to much lower neutron energies, at which the contribution from the
fragmentation model dominates.

The measurements are also compared to predictions of several hadronic
interaction models which are commonly used for the simulation of cosmic ray
air shower
cascades: EPOS~LHC~\cite{Werner:2005jf,Pierog:2013ria}, 
QGSJET~01~\cite{Kalmykov:1993qe,Kalmykov:1997te},
QGSJET~II-04~
% \cite{Ostapchenko:2005nj,Ostapchenko:2007qb}  
\cite{Ostapchenko:2005nj,Ostapchenko:2010vb}  
and SIBYLL~2.1\cite{Engel:1992vf,Ahn:2009wx}.
% \cite{Engel:1992vf,Fletcher:1994bd,Ahn:2009wx}.
These models are based on Regge theory \cite{Collins:1977jy}, 
on the Reggeon calculus of Gribov \cite{Gribov:1968fc}
and on perturbative QCD. 
They use an unitarisation procedure to reconstruct amplitudes for
exclusive processes and to determine the total and elastic cross
sections.
Central elements of these models are the production of mini-jets and 
the formation of colour strings that fragment into hadrons.
Whereas the Regge-Gribov approximation is applied to hadrons as 
interacting objects in the case of QGSJET and SIBYLL, it is extended 
to include partonic constituents in EPOS~LHC. 
Compared to the earlier EPOS simulation~\cite{Werner:2005jf}, 
which was used for comparison with the 
previous H1 forward photon analysis \cite{Aaron:2011pe}, the new 
EPOS~LHC model~\cite{Pierog:2013ria} 
includes a  modified treatment of central diffraction and 
the diffractive remnant in order to reproduce rapidity gap 
measurements at the LHC.
%
%
The CR models also differ
in the treatment of saturation effects at high parton densities at small 
Bjorken-$x$
and in the treatment of the hadronic remnants in collisions.
%
 The programs are interfaced with the PHOJET program~\cite{Engel:1995yda}
 for the generation of the $ep$ scattering kinematics.
%
It was pointed out \cite{ostapchenko1} that the hadronic interaction models
have been developed for hadron-hadron interactions and therefore
the simulation of DIS events might
be affected by the superfluous contribution of 
multi-parton interactions. In order to investigate this assumption,
the QGSJET~01 model has been modified \cite{ostapchenko2} to exclude 
the multi-parton interactions. In the comparison with the  measurements
this modified model is denoted as `QGSJET~01~(no~mi)'.

The measured distributions may contain background arising from several
sources.  The background from photoproduction processes
is estimated using the PHOJET MC generator. It is found to
be about $1\%$ on average and is subtracted from the data distributions
bin-by-bin.
%
Background from misidentification of photons or
neutrons in the FNC is estimated
 from the DJANGOH ~MC simulation 
 to be $2\%$ on average
 and is subtracted from the data distributions bin-by-bin.
%
Background also arises from a random
coincidence of DIS events, causing activity in the central detector, 
with a beam-related background signal in the FNC, produced
from the interaction of another beam proton with a positron or
with residual gas in the beampipe. 
This contribution is estimated by combining DIS events with FNC clusters
originating from interactions in the bunch-crossings adjacent
to the bunch-crossings of the DIS events. 
It is found to be smaller than $1\%$ and is neglected.
%
%-------------------------------------------------------------------------------

The MC simulations are 
used to correct the distributions 
at the level of reconstructed particles
back to the hadron level on a bin-by-bin basis.
%
The correction factors are determined from  the 
MC simulations as the ratios of the normalised cross sections 
obtained from particles at hadron level without QED radiation
to the normalised cross sections calculated using reconstructed particles
and including QED radiation effects.
%
For the forward photon analysis 
the average of the correction factors determined from LEPTO and CDM 
is used. 
For the forward neutron analysis the correction factors are calculated
using the  combination of RAPGAP-$\pi$ and CDM  simulations,
with the weighting factors $0.6$ and $1.4$, as described above. 
%
The size of the 
correction  factors varies between $2$ and $3.5$ for the forward photon
and between $2$ and $6$ for the forward neutron $x_F$ distributions, 
and is about $3.2$ for the $W$  distribution in both cases.
The correction factors are large mainly  due to
 the non-uniform azimuthal acceptance of the FNC, 
which is about $30\%$ on average.
%
The bin purity, defined as the fraction
of events reconstructed in a particular bin that originate from
the same bin on hadron level,  varies between $75\%$ and $95\%$.

%--------------------------------------------------------------------------------------
\subsection{Systematic uncertainties}


Several sources of experimental uncertainties are considered and
their effect on the measured cross section is quantified. The systematic
uncertainties on the cross section measurements are determined using
MC simulations, by propagating the corresponding
uncertainty through the full analysis chain.

As the cross sections are normalised to the inclusive DIS cross
section measured in this analysis, some important systematic
uncertainties,
 such as those involving the trigger efficiency and the integrated luminosity
 and those related to the reconstruction of the scattered positron
 and the hadronic final state 
are largely reduced or cancel.
The following sources are considered for both the DIS sample and the 
FNC samples:
\begin{itemize}
\item
The uncertainty on the measurements of the scattered positron energy ($1\%$).
\item
The uncertainty on the measurements of the scattered positron angle ($1$~mrad).
\item
The uncertainty on the measurements of the energy of the hadronic final state ($4\%$). 
\item 
The uncertainty on the trigger efficiency ($1\%$).
\end{itemize}
%
These uncertainties are strongly correlated between the DIS and the 
forward photon and neutron samples. 
The resulting combined uncertainty of the cross section is 
about $2\%$ on average and
is considered as uncorrelated between the measurement points.


Several sources of uncertainties related to the reconstruction of the 
forward photons and neutrons in the FNC are considered:
\begin{itemize}
\item
The acceptance of the FNC calorimeter is
defined by the interaction point and the geometry of the HERA magnets
and is determined using MC simulations.
The uncertainty of the impact position of the particle on the FNC,
due to beam inclination and the uncertainty on the FNC position,
is estimated to be $5$~mm.  This results in uncertainties on the 
FNC acceptance determination of up to $15\%$ for the $x_F$ distributions. 
% which results up to $6.9\%$ uncertainty in the FNC acceptance.  
\item
The absolute electromagnetic energy scale of the FNC is known to a
precision of $5\%$~\cite{Aaron:2011pe}.
 This leads to an uncertainty of $1\%$ on the forward photon
cross section measurement at low energies, increasing to $11\%$ for the
largest $x_F$ values. 
\item
The uncertainty in the neutron detection efficiency and the $2\%$ uncertainty 
on the absolute hadronic energy scale of the FNC~\cite{Aaron:2010ab} 
lead to 
systematic errors 
on the cross section of  $2\%$ and up to $10\%$, respectively.
% These effects
\end{itemize}
These systematic uncertainties related to the FNC
are strongly correlated between measurement bins and
mainly contribute to the overall normalisation uncertainty.
For the normalised cross sections studied as a function of $W$
all above-mentioned FNC related systematic uncertainties 
contribute to a normalisation uncertainty of approximately $7\%$.

In the procedure of correcting the measurements to the hadron level, using
MC simulations, the following sources of systematic uncertainties are 
considered:
%
\begin{itemize}
\item
The systematic uncertainty arising from the radiative corrections 
and the model dependence of the data correction for the forward neutron cross 
section is estimated by varying
the DJANGOH and RAPGAP-$\pi$ scaling factors 
within values permitted by the data and by switching between
the CDM and LEPTO  models within DJANGOH.
 The resulting uncertainty on the cross section is typically $4\%$,
increasing to $5\%$ at lowest and highest $x_F$ values.
\item
For the forward photon cross section the systematic uncertainty,
taken as the difference of the acceptance corrections calculated
using the LEPTO and CDM models, increases from $1\%$ at low $x_F$
to $7\%$ at higher $x_F$.
\item
The use of different parton
distribution functions in the MC simulation results in a negligible
change of the correction factors.
\item
A $2\%$ uncertainty is attributed to the bin-by-bin subtraction
of background 
%In addition, a $2\%$ uncertainty is attributed to the contribution
%from background 
arising from the wrong identification of photon
and neutron candidates in the FNC and from photoproduction
processes.
\end{itemize}
These uncertainties
are assumed to be equally
shared between correlated and uncorrelated parts. 

The systematic uncertainties shown in the figures and tables 
are calculated using the quadratic sum of all contributions, 
which may vary from point to point. 
They are significantly larger than the statistical uncertainties 
in all measurement bins.
The total systematic uncertainty for the normalised cross section
measurements ranges between $8\%$ and $22\%$ for the $x_F$ dependent cross
sections and is about $8\%$ for the $W$ dependent cross sections.


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

\subsection{Normalised cross sections as a function of $W$}
\label{sec:resultW}

The ratios of the forward photon and forward neutron production
cross sections to the inclusive DIS cross section, 
%$\rm \frac{\sigma^{\gamma,n}_{DIS}}{\sigma_{DIS}}$, 
 $ \sigma^{\gamma,n}_{DIS}/\sigma_{DIS}$, 
in the kinematic range
$6<Q^2<100~\GeVS$ and $0.05<y<0.6$ and as a function of $W$
are presented in Tables~\ref{tab:table1} and \ref{tab:table2}
and are shown in  Figure~\ref{figW}.
Within uncertainties the measured ratios
 are consistent  with constant values of  about
$0.027$ (forwards photons) and $0.083$ (forward neutrons).
In other words, within uncertainties the $W$ dependence
   of the cross section is independent of the presence of a forward
   neutron or a forward photon, as predicted by the limiting
   fragmentation hypothesis ~\cite{Benecke:1969sh,Chou:1994dh}.

In Figure~\ref{figW} the MC model calculations are compared to the measurements.  
Both CDM and LEPTO predict 
a forward photon rate of about $70\%$ higher than observed.
A similar excess was observed earlier \cite{Aaron:2011pe}.
The photon production rate as a function of $W$ is rather flat in CDM
and shows a slight increase with $W$ in LEPTO. 
The shape of the $W$ distribution is in both models consistent with the data,
within errors.
%The shapes of distributions
%for both models are consistent with the data within errors.
%

The rate of forward neutron production predicted by LEPTO is consistent
with the data, while  CDM predicts a much lower rate. However, as was 
shown in the previous measurement~\cite{Aaron:2010ab}, 
the energy distribution
of forward neutrons can be described by MC simulation only if this
includes contributions both from standard fragmentation as simulated
in DJANGOH, and from a pion exchange mechanism as explicitly simulated in
RAPGAP-$\pi$ but not included in DJANGOH. 
In  Figure~\ref{figW}b the combinations of the RAPGAP-$\pi$ and DJANGOH
simulations, as described in section \ref{sec:mc},  
are compared to the measurement.
The weighting factors $1.4$ for the CDM, $0.7$ for the LEPTO and $0.6$ for the
RAPGAP-$\pi$
predictions are determined by fitting the observed neutron
energy distributions integrated over the full $W$ range.
The cross sections for inclusive DIS events, used for the normalisation 
of the forward neutron cross sections, $\sigma_{DIS}$, are taken from
the CDM and LEPTO simulations without additional weights. %ing.
The model combination describes the observed $W$ dependence well.
It is remarkable that 
the factors for
the CDM and LEPTO contributions differ by a factor two 
($1.4$ and $0.7$, respectively).
It is also notable that the CDM model, which overestimates
the rate of forward photons by about $70\%$, has to be scaled up 
in the combination to describe the forward neutron data.

In  Figure~\ref{figWcos} 
predictions of various cosmic ray hadronic interaction models (EPOS~LHC, 
SIBYLL~2.1 and the two versions of QGSJET) are compared to the measured
normalised cross sections as a function of $W$.
The CR model predictions show significant differences in
absolute values,
for both forward photons and forward neutrons. 
For photons all models predict too high rates by
$30$ to $40\%$, and these rates, with the exception of EPOS~LHC, show a slight
decrease with increasing $W$, not confirmed by  data.
For forward neutrons all CR predictions show a $W$ independent behaviour,
in accordance with the measured $W$ dependence.
The QGSJET~01 model predicts a much too high and SIBYLL~2.1 a much
too low neutron rate,
while the EPOS~LHC and QGSJET~II-04 models are closer to the measurement.


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

\subsection{Normalised cross sections as a function of \boldmath $x_F$ and test of Feynman scaling}
%
The measured normalised differential cross sections,
%  $\rm \frac{1}{\sigma_{DIS}}\frac{d\sigma}{dx_F}$,
  $ {1}/\sigma_{DIS}~  {\rm d}\sigma/{\rm d}x_F$,
 of the most energetic photon  
are presented as a function of $x_F$ 
in  Table~\ref{tab:table3} and in Figures~\ref{fig:xfgmc} and \ref{fig:xfgcos}
for the kinematic region defined in Table~\ref{tab:selection}.
In order to study the energy dependence of the $x_F$ distributions,
these cross sections are measured in three $W$ intervals.

The normalised differential cross sections as a function of $x_F$ 
are similar for the three $W$ ranges.
As shown in Figure~\ref{fig:xfgmc}
and already seen in the comparison of the $W$ dependence,
the LEPTO and CDM models predict a rate of forward photons about $70\%$ 
higher than measured.
The shapes of the measured distributions are well described by
LEPTO, while the CDM description is very poor
by showing a significantly harder spectrum than observed in data.
%
In Figure~\ref{fig:xfgcos} the predictions of the CR hadronic 
interaction models are compared to the same measurements.
Large differences between the CR models are observed, both in shape and
in normalisation.
All models tested here overestimate the forward photon 
rate by $30\%$ to $40\%$ at low $x_F$.
% , as was already shown for the measured $W$ dependent
% distributions in Figure~\ref{figWcos}.
The EPOS~LHC model describes the shapes of the photon $x_F$
distributions well. 
The SIBYLL~2.1 model
% (black line) 
predicts a harder $x_F$ dependence,
while the spectra obtained from the different variants of QGSJET are
softer than observed in the data.
%
% {\bf LHC} \\
Forward photon and neutral pion measurements at the LHC also revealed 
differences with respect to a similar selection of CR 
models~\cite{Adriani:2011nf,Adriani:2012ap,Adriani:2012fz}.


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

The normalised differential cross sections for forward neutrons
are presented in 
Table~\ref{tab:table4} and  in
Figures~\ref{fig:xfn} and \ref{fig:xfncos}
for the kinematic region defined in Table~\ref{tab:selection}.
%
The $x_F$ distributions %of forward neutrons  
are well reproduced %described
by a combination of CDM and RAPGAP-$\pi$, using the weighting factors
and normalisation as described in section \ref{sec:resultW}.
The individual contributions of the two models are shown in 
Figure~\ref{fig:xfn} as well.
Fragmentation, as simulated by CDM,  dominates the neutron 
production at lower $x_F$, while the contribution from 
pion exchange becomes significant at $x_F\gsim 0.7$.
The combination of LEPTO and RAPGAP-$\pi$ (not shown) also provides a 
good description of the measurements for the three $W$ ranges.
In Figure~\ref{fig:xfncos} the predictions of the CR hadronic 
interaction models are compared to the
forward neutron production cross sections.
The EPOS~LHC model provides a reasonable description
of the neutron $x_F$ distributions, except at the highest $x_F$ values.
The SIBYLL~2.1 model describes the shape of the $x_F$ spectra but
fails in the absolute rate.
The QGSJET~II-04 model 
shows a harder
$x_F$ dependence, 
and QGSJET~01 predicts a much too high neutron rate.

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

A modified version of the 
QGSJET~01 model, denoted 'QGSJET 01~(no~mi)' \cite{ostapchenko2}, 
in which the contribution of multi-parton interactions is excluded
(see section \ref{sec:mc}), is also compared to the measurements.
%As seen in Figures \ref{fig:xfgcos} and \ref{fig:xfncos} the effect of multi-parton interactions on the forward photon production  is  small, while 
%the $x_F$ spectrum of forward neutrons becomes significantly harder,
When  multi-parton interactions are switched off, the predicted $x_F$ spectra 
become harder without improving the data description.


The $W$ dependence of the $x_F$ distributions allows a test of the 
Feynman scaling hypothesis for particle production. For this test, the 
ratios of the normalised cross sections for different CM energy 
intervals are studied as a function of $x_F$. 
Figures~\ref{fig:xfmcratio}~and~\ref{fig:xfgcosratio} show the ratios of the second 
to the first and the third to the first $W$ range for photons.
The predictions from CDM, LEPTO and the CR models are also shown.
In Figure~\ref{fig:xfnratio} the same ratios are shown for forward neutrons and
the CR models are compared to the data.
For all data distributions the values of these ratios 
are %roughly 
consistent with unity and with being constant
within uncertainties, suggesting that Feynman scaling
in the target fragmentation region holds for photons and neutrons.
The LEPTO and CDM MC models, used for the comparison with forward
photon data, show a similar behaviour. 
All CR models indicate deviations from  scaling
for the forward photons, such that the production
rate decreases with increasing $W$. In particular,
this effect is strong for the SIBYLL~2.1 and QGSJET~01 models.
For forward neutrons the CR models
are consistent with Feynman scaling, with exception of SIBYLL~2.1.



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

\section{Summary}
\noindent
The production of high energy forward neutrons and 
photons has been studied at HERA
in deep-inelastic $ep$ scattering in the kinematic region
$6< Q^2 < 100~\GeV^2$ and $0.05<y<0.6$.
%
The normalised DIS cross sections $1/\sigma_{DIS}~{\rm d}\sigma/{\rm d}x_F$
for the production of  photons and neutrons at pseudorapidities 
$\eta>7.9$ and in the range of Feynman-$x$ of $0.1<x_F<0.7$ for the photons 
and $0.1<x_F<0.94$ for neutrons are presented.
%
The measured cross sections as a function of $x_F$ at 
different centre-of-mass
energies of the virtual photon-proton system agree within uncertainties,
confirming the validity of Feynman scaling in the energy range
of the virtual photon-proton system $70<W<245~\GeV$.
%

Different Monte Carlo models are compared to the measurements.
All these models overestimate the rate of photons by $30-70\%$.
The shapes of the measured forward photon 
cross sections are well described by the LEPTO
MC simulation, while predictions based on the colour dipole model fail,
especially at high $x_F$. 
The cross sections for forward neutrons are well described by 
a linear combination of the standard fragmentation model, as implemented in
DJANGOH, and 
the one-pion-exchange model RAPGAP-$\pi$.
%
Predictions of models,  which are commonly used for the simulation of
cosmic ray cascades, are also compared to the 
forward photon and neutron measurements.
None of the models describes the photon and neutron data 
simultaneously well. The best description of the shapes
of the photon and the neutron $x_F$ distributions is
provided by the EPOS~LHC model.
%
Within the kinematic range of the measurements, 
the relative rate of forward photons
and neutrons in DIS events is observed to be independent of the energy
of the virtual photon-proton CM, 
and therefore also consistent
% in agreement 
with the hypothesis of limiting fragmentation.
The present measurement provides new information to further improve
the understanding of proton fragmentation in collider and 
cosmic ray experiments.

%--------------------------------------------------------------------------------------
\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.
We would like to give credit to all partners contributing to the EGI 
computing infrastructure for their support for the H1 Collaboration.
We also wish to thank Tanguy Pierog, Ralph Engel and Sergey Ostapchenko
for providing the predictions of the cosmic ray models and for 
fruitful discussions.


%--------------------------------------------------------------------------------------
\newpage
% \bibliography{biblio}
% \bibliographystyle{f2lnart}
% \bibliographystyle{utphys}

\begin{thebibliography}{10}

\bibitem{Engel:1998hf}
R.~Engel,
\newblock Nucl. Phys. Proc. Suppl. {\bf 75A} (1999) 62 [astro-ph/9811225].
%%CITATION = ASTRO-PH/9811225;%%

\bibitem{Bunyatyan:2009zz}
A.~Buniatyan {\em et~al.},
{\it ''Cosmic Rays, HERA and the LHC: Working group summary."}
Proceedings of the Workshop on the Implications of HERA for LHC Physics,
Geneva, Switzerland, 26-30 May 2008, p.566, DESY-PROC-2009-02 
[arXiv:0903.3861].

\bibitem{Adloff:1998yg}
C.~Adloff {\em et~al.} [H1 Collaboration],
\newblock Eur. Phys. J. {\bf C6} (1999) 587 [hep-ex/9811013].
%%CITATION = HEP-EX/9811013;%%

\bibitem{Chekanov:2002pf}
S.~Chekanov {\em et~al.} [ZEUS Collaboration],
\newblock Nucl. Phys. {\bf B637} 3 (2002) 3  [hep-ex/0205076].
%%CITATION = HEP-EX/0205076;%%

\bibitem{Chekanov:2007tv}
S.~Chekanov {\em et~al.} [ZEUS Collaboration],
\newblock Nucl. Phys. {\bf B776} (2007)  1 [hep-ex/0702028].
%%CITATION = HEP-EX/0702028;%%

\bibitem{Chekanov:2008tn}
S.~Chekanov {\em et~al.} [ZEUS Collaboration],
\newblock JHEP {\bf 06} (2009) 074 [arXiv:0812.2416].
%%CITATION = 0812.2416;%%

\bibitem{Aaron:2010ab}
F.~Aaron {\em et~al.} [H1 Collaboration],
\newblock Eur. Phys. J. {\bf C68} (2010) 381 [arXiv:1001.0532].
%%CITATION = ARXIV:1001.0532;%%

\bibitem{Aaron:2011pe}
F.~Aaron {\em et~al.} [H1 Collaboration],
\newblock Eur. Phys. J. {\bf C71} (2011) 1771 [arXiv:1106.5944].
%%CITATION = ARXIV:1106.5944;%%

\bibitem{Benecke:1969sh}
J.~Benecke  {\em et~al.}, 
% J.~Benecke, T.~T. Chou, C.-N. Yang and E.~Yen,
\newblock Phys. Rev. {\bf 188} (1969) 2159.
%%CITATION = PHRVA,188,2159;%%

\bibitem{Chou:1994dh}
T.~T. Chou and C.-N. Yang,
\newblock Phys. Rev. {\bf D50} (1994) 590.
%%CITATION = PHRVA,D50,590;%%

\bibitem{Feynman:1969ej}
R.~P. Feynman,
\newblock Phys. Rev. Lett. {\bf 23} (1969) 1415.
%%CITATION = PRLTA,23,1415;%%

\bibitem{Arnison:1982ed}
G.~Arnison {\em et~al.} [UA1 Collaboration],
\newblock Phys. Lett. {\bf B118} (1982) 167.
%%CITATION = PHLTA,B118,167;%%

\bibitem{Arnison:1982rm}
G.~Arnison {\em et~al.} [UA1 Collaboration],
\newblock Phys. Lett. {\bf B123} (1983) 108.
%%CITATION = PHLTA,B123,108;%%

\bibitem{Banner:1984wh}
M.~Banner {\em et~al.} [UA2 Collaboration],
%  Bern-CERN-Copenhagen-Orsay-Pavia-Saclay Collaboration],
\newblock Z. Phys. {\bf C27} (1985) 329.
%%CITATION = ZEPYA,C27,329;%%

\bibitem{Bernard:1985kh}
D.~Bernard {\em et~al.} [UA4 Collaboration],
\newblock Phys. Lett. {\bf B166} (1986) 459.
%%CITATION = PHLTA,B166,459;%%

\bibitem{Alner:1986xu}
G.~Alner {\em et~al.} [UA5 Collaboration],
\newblock Z. Phys. {\bf C33} (1986) 1.
%%CITATION = ZEPYA,C33,1;%%

\bibitem{Alner:1987wb}
G.~Alner {\em et~al.} [UA5 Collaboration],
\newblock Phys. Rep. {\bf 154} (1987) 247.
%%CITATION = PRPLC,154,247;%%

\bibitem{Alexopoulos:1988na}
T.~Alexopoulos {\em et~al.},
\newblock Phys. Rev. Lett. {\bf 60} (1988) 1622.
%%CITATION = PRLTA,60,1622;%%

\bibitem{Abe:1989td}
F.~Abe {\em et~al.} [CDF Collaboration],
\newblock Phys. Rev. {\bf D41} (1990) 2330.
%%CITATION = PHRVA,D41,2330;%%

\bibitem{Aaron:2007ds}
F.~Aaron {\em et~al.} [H1 Collaboration],
\newblock Phys. Lett. {\bf B654} (2007) 148 [arXiv:0706.2456].
%%CITATION = ARXIV:0706.2456;%%

\bibitem{Abramowicz:2010rz}
H.~Abramowicz {\em et~al.} [ZEUS Collaboration],
\newblock JHEP {\bf 1006} (2010) 009 [arXiv:1001.4026].
%%CITATION = ARXIV:1001.4026;%%

\bibitem{Pare:1989mr}
E.~Pare {\em et~al.} [UA7 Collaboration],
\newblock Phys. Lett. {\bf B242} (1990) 531.
%%CITATION = PHLTA,B242,531;%%

\bibitem{Albrow:1973kj}
M.~Albrow {\em et~al.} [CHLM Collaboration],
\newblock Nucl. Phys. {\bf B56} (1973) 333.
%%CITATION = NUPHA,B56,333;%%

\bibitem{Albrow:1973ee}
M.~Albrow {\em et~al.} [CHLM Collaboration],
\newblock Nucl. Phys. {\bf B73} (1974) 40.
%%CITATION = NUPHA,B73,40;%%

\bibitem{Capiluppi:1973fz}
P.~Capiluppi {\em et~al.},
% P.~Capiluppi, G.~Giacomelli, A.~Rossi, G.~Vannini and A.~Bussiere,
\newblock Nucl. Phys. {\bf B70} (1974) 1.
%%CITATION = NUPHA,B70,1;%%

\bibitem{Capiluppi:1974rt}
P.~Capiluppi {\em et~al.},
\newblock Nucl. Phys. {\bf B79} (1974) 189.
%%CITATION = NUPHA,B79,189;%%

\bibitem{Ohsawa:1992dj}
A.~Ohsawa, 
\newblock Prog. Theor. Phys. {\bf 92} (1994) 1005.

\bibitem{Abt:1996hi}
I.~Abt {\em et~al.} [H1 Collaboration],
\newblock Nucl. Instrum. Meth. {\bf A386} (1997) 310.
%%CITATION = NUIMA,A386,310;%%

\bibitem{Abt:1996xv}
I.~Abt {\em et~al.} [H1 Collaboration],
\newblock Nucl. Instrum. Meth. {\bf A386} (1997) 348.
%%CITATION = NUIMA,A386,348;%%

\bibitem{Appuhn:1996na}
R.~D. Appuhn {\em et~al.} [H1 SPACAL Group],
\newblock Nucl. Instrum. Meth. {\bf A386} (1997) 397.
%%CITATION = NUIMA,A386,397;%%

\bibitem{Pitzl:2000wz}
D.~Pitzl {\em et~al.},
\newblock Nucl. Instrum. Meth. {\bf A454} (2000) 334 [hep-ex/0002044].
%%CITATION = HEP-EX/0002044;%%

\bibitem{Andrieu:1993kh}
B.~Andrieu {\em et~al.} [H1 Calorimeter Group],
\newblock Nucl. Instrum. Meth. {\bf A336} (1993) 460.
%%CITATION = NUIMA,A336,460;%%

\bibitem{Laycock:2012xg}
P.J.~Laycock {\em et~al.},
\newblock JINST {\bf 7} (2012) T08003 [arXiv:1206:4068].


\bibitem{Nicholls:1995di}
T.~Nicholls {\em et al.} [H1 SPACAL Group],
\newblock Nucl. Instrum. Meth. {\bf A374} (1996)  149.

\bibitem{Aaron:2012kn}
F.~Aaron {\em et al.} [H1 Collaboration], 
Eur. Phys. J. {\bf C72} (2012) 2163, Erratum-ibid. {\bf C74} (2014) 2733
 [arXiv:1205.2448].

\bibitem{Uraev}
A.~Uraev, 
  ``{\it A new Forward Neutron Calorimeter for the H1 Experiment at the HERA collider.}" 
  (In Russian), Diploma thesis, Moscow Engineering Physics Institute (2001) 
  (available at  http://www-h1.desy.de/publications/theses\_list.html).

\bibitem{Peez:2003zd}
M.~Peez,
``{\it Search for deviations from the standard model in high transverse energy
  processes at the electron proton collider HERA.}" (In French), PhD thesis,
  Univ.~Lyon (2003), CPPM-T-2003-04 
  (available at http://www-h1.desy.de/publications/theses\_list.html).
%%CITATION = DESY-THESIS-2003-023 ETC.;%%.

\bibitem{Hellwig:2004yp}
 S.~Hellwig, ``{\it Investigation of the $D^* - \pi_{slow}$ double tagging method in
  charm analyses.}" (In German), Diploma thesis, Univ.~Hamburg (2004) 
 (available at http://www-h1.desy.de/publications/theses\_list.html).

\bibitem{Adloff:1997sc}
C.~Adloff {\em et~al.} [H1 Collaboration],
\newblock Z. Phys. {\bf C76} (1997) 613 [hep-ex/9708016].
%%CITATION = HEP-EX/9708016;%%

\bibitem{Brun:1978fy}
R.~Brun  {\em et~al.}, GEANT3, CERN-DD/EE/84-1.
% R.~Brun, R.~Hagelberg, M.~Hansroul and J.~C. Lassalle,
%\newblock CERN-DD-78-2-REV.

\bibitem{Charchula:1994kf}
K.~Charchula, G.~A. Schuler and H.~Spiesberger,
\newblock DJANGOH~1.4,
\newblock Comput. Phys. Commun. {\bf 81} (1994) 381.
%%CITATION = CPHCB,81,381;%%

\bibitem{Ingelman:1996mq}
G.~Ingelman, A.~Edin and J.~Rathsman,
\newblock LEPTO 6.5,
\newblock Comput. Phys. Commun. {\bf 101} (1997) 108 [hep-ph/9605286].
%%CITATION = HEP-PH/9605286;%%

\bibitem{Lonnblad:1992tz}
L.~L\"onnblad,
\newblock ARIADNE~4.10,
\newblock Comput. Phys. Commun. {\bf 71} (1992) 15.
%%CITATION = CPHCB,71,15;%%

\bibitem{Andersson:1983ia}
B.~Andersson  {\em et~al.},
% B.~Andersson, G.~Gustafson, G.~Ingelman and T.~Sj\"ostrand,
\newblock Phys. Rep. {\bf 97} (1983) 31.
%%CITATION = PRPLC,97,31;%%

\bibitem{Sjostrand:1995iq}
T.~Sj\"ostrand,
\newblock PYTHIA 5.7 and JETSET 7.4,
\newblock hep-ph/9508391.
%%CITATION = HEP-PH/9508391;%%

\bibitem{Kwiatkowski:1990es}
A.~Kwiatkowski, H.~Spiesberger and H.~J. M\"ohring,
\newblock Comp. Phys. Commun. {\bf 69} (1992) 155.
%%CITATION = CPHCB,69,155;%%

\bibitem{Edin:1995gi}
A.~Edin, G.~Ingelman and J.~Rathsman,
\newblock Phys. Lett. {\bf B366} (1996) 371 [hep-ph/9508386].
%%CITATION = HEP-PH/9508386;%%

\bibitem{Aaron:2009kv}
F.~D. Aaron {\em et~al.} [H1 Collaboration],
\newblock Eur. Phys. J. {\bf C64} (2009) 561 [arXiv:0904.3513].
%%CITATION = 0904.3513;%%

\bibitem{Jung:1993gf}
H.~Jung,
\newblock RAPGAP 3.1,
\newblock Comp. Phys. Commun. {\bf 86} (1995) 147.
%%CITATION = CPHCB,86,147;%%

\bibitem{Bishari:1972tx}
M.~Bishari,
\newblock Phys. Lett. {\bf B38} (1972) 510.
%%CITATION = PHLTA,B38,510;%%

\bibitem{Holtmann:1994rs}
H.~Holtmann {\em et~al.},
\newblock Phys. Lett. {\bf B338} (1994) 363.
%%CITATION = PHLTA,B338,363;%%

\bibitem{Kopeliovich:1996iw}
B.~Kopeliovich, B.~Povh and I.~Potashnikova,
\newblock Z. Phys. {\bf C73} (1996) 125 [hep-ph/9601291].
%%CITATION = HEP-PH/9601291;%%

\bibitem{Przybycien:1996zb}
M.~Przybycie\'n, A.~Szczurek and G.~Ingelman,
\newblock Z. Phys. {\bf C74} (1997) 509 [hep-ph/9606294].
%%CITATION = HEP-PH/9606294;%%

\bibitem{Szczurek:1997cw}
A.~Szczurek, N.~N. Nikolaev and J.~Speth,
\newblock Phys. Lett. {\bf B428} (1998) 383 [hep-ph/9712261].
%%CITATION = HEP-PH/9712261;%%

\bibitem{Werner:2005jf}
K.~Werner, F.-M. Liu and T.~Pierog,
\newblock Phys. Rev. {\bf C74} (2006) 044902 [hep-ph/0506232].
%%CITATION = HEP-PH/0506232;%%

\bibitem{Pierog:2013ria}
T.~Pierog  {\em et~al.},
% T.~Pierog , I.~Karpenko, J.~Katzy, E.~Yatsenko and K.~Werner,
\newblock DESY-13-125, arXiv:1306.0121.
%%CITATION = ARXIV:1306.0121;%%

\bibitem{Kalmykov:1993qe}
N.~N. Kalmykov and S.~S. Ostapchenko,
\newblock Phys. Atom. Nucl. {\bf 56} (1993) 346.
%%CITATION = PANUE,56,346;%%

\bibitem{Kalmykov:1997te}
N.~N. Kalmykov, S.~S. Ostapchenko and A.~I. Pavlov,
\newblock Nucl. Phys. Proc. Suppl. {\bf 52B} (1997) 17.
%%CITATION = NUPHZ,52B,17;%%

\bibitem{Ostapchenko:2005nj}
S.~S. Ostapchenko,
\newblock Phys. Rev. {\bf D74} (2006) 014026 [hep-ph/0505259].
%%CITATION = HEP-PH/0505259;%%

\bibitem{Ostapchenko:2010vb}
S.~S.~Ostapchenko,
\newblock Phys. Rev. {\bf D83} (2011) 014018 [arXiv:1010.1869].
%%CITATION = ARXIV:1010.1869;%%

\bibitem{Engel:1992vf}
J.~Engel  {\em et~al.},
% J.~Engel, T.~K. Gaisser, T.~Stanev and P.~Lipari,
\newblock Phys. Rev. {\bf D46} (1992) 5013.
%%CITATION = PHRVA,D46,5013;%%

\bibitem{Ahn:2009wx}
E.-J. Ahn  {\em et~al.},
% E.-J. Ahn, R.~Engel, T.~K. Gaisser, P.~Lipari and T.~Stanev,
\newblock Phys. Rev. {\bf D80} (2009) 094003 [arXiv:0906.4113].
%%CITATION = 0906.4113;%%


\bibitem{Collins:1977jy}
P.~D.~B. Collins, ``{\it An Introduction to Regge Theory and High-Energy  Physics}.''
 Cambridge 1977, 445p.

\bibitem{Gribov:1968fc}
V.~N. Gribov,
\newblock Sov. Phys. JETP {\bf 26} (1968) 414.
%%CITATION = SPHJA,26,414;%%

\bibitem{Engel:1995yda}
R.~Engel and J.~Ranft,
\newblock Phys. Rev. {\bf D54} (1996) 4244 [hep-ph/9509373].
%%CITATION = HEP-PH/9509373;%%

\bibitem{ostapchenko1}
S.~S. Ostapchenko,
\newblock {\em {Private communication}}.

\bibitem{ostapchenko2}
S.~S. Ostapchenko and T.~Pierog,
\newblock {\em {Private communication}}.

\bibitem{Adriani:2011nf}
O.~Adriani {\em et~al.} [LHCf Collaboration],
\newblock Phys. Lett. {\bf B703} (2011) 128  [arXiv:1104.5294].
%%CITATION = ARXIV:1104.5294;%%

\bibitem{Adriani:2012ap}
O.~Adriani {\em et~al.} [LHCf Collaboration],
\newblock Phys. Rev. {\bf D86} (2012) 092001  [arXiv:1205.4578].
%%CITATION = ARXIV:1205.4578;%%

\bibitem{Adriani:2012fz}
O.~Adriani {\em et~al.} [LHCf Collaboration],
\newblock Phys. Lett. {\bf B715} (2012) 298  [arXiv:1207.7183].
%%CITATION = ARXIV:1207.7183;%%



\end{thebibliography}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\newpage


\begin{table}[hp]
\centering
{
\input{desy-14-035.table2.tex}
}
\caption{
The fraction of DIS events with forward photons
in the kinematic region given in Table~\ref{tab:selection}.
For each measurement, the  statistical ($\delta_{stat.}$), 
the total systematic  ($\delta_{total~sys.}$), the
uncorrelated  ($\delta_{uncorrel.sys.}$)
systematic uncertainties and the bin-to-bin correlated
systematic uncertainties due to the FNC absolute 
energy scale  ($\delta_{E_{FNC}}$), 
the measurement of the particle impact position 
in the FNC  ($\delta_{XY_{FNC}}$)
and the model dependence of the data  
correction ($\delta_{model}$) are given.
}
\label{tab:table1}
\end{table}



\begin{table}[hp]
\centering
{
\input{desy-14-035.table3.tex}
}
\caption{
The fraction of DIS events with forward neutrons
in the kinematic region given in Table~\ref{tab:selection}.
For each measurement, the  statistical ($\delta_{stat.}$), 
the total systematic  ($\delta_{total~sys.}$), the
uncorrelated  ($\delta_{uncorrel.sys.}$)
systematic uncertainties and the bin-to-bin correlated
systematic uncertainties due to the FNC absolute 
energy scale  ($\delta_{E_{FNC}}$), 
the measurement of the particle impact position 
in the FNC  ($\delta_{XY_{FNC}}$)
and the model dependence of the data
correction ($\delta_{model}$) are given.
}
\label{tab:table2}
\end{table}



\begin{table}[hp]
\centering
{
\input{desy-14-035.table4.tex}
}
\caption{
Normalised cross sections of forward photon production in DIS
as a function of $x_F$.
The kinematic phase space of the measurements is
given in Table~~\ref{tab:selection}.
For each measurement, the  statistical ($\delta_{stat.}$), 
the total systematic  ($\delta_{total~sys.}$), the
uncorrelated  ($\delta_{uncorrel.sys.}$)
systematic uncertainties and the bin-to-bin correlated
systematic uncertainties due to the FNC absolute 
energy scale  ($\delta_{E_{FNC}}$), 
the measurement of the particle impact position 
in the FNC  ($\delta_{XY_{FNC}}$)
and the model dependence of the data  
correction ($\delta_{model}$) are given.
}
\label{tab:table3}
\end{table}

\begin{table}[hp]
\centering
{
\input{desy-14-035.table5.tex}
}
\caption{
Normalised cross sections of forward neutron production in DIS
as a function of $x_F$.
The kinematic phase space of the measurements is
given in Table~~\ref{tab:selection}.
For each measurement, the  statistical ($\delta_{stat.}$), 
the total systematic  ($\delta_{total~sys.}$), the
uncorrelated  ($\delta_{uncorrel.sys.}$)
systematic uncertainties and the bin-to-bin correlated
systematic uncertainties due to the FNC absolute 
energy scale  ($\delta_{E_{FNC}}$), 
the measurement of the particle impact position 
in the FNC  ($\delta_{XY_{FNC}}$)
and the model dependence of the data 
correction ($\delta_{model}$) are given.
}
\label{tab:table4}
\end{table}




\newpage

%\input figsW.tex
%\input figsxf.mc.tex
%\input figsxf.cos.tex
%\input figsxfn.mc.tex
%\input figsxfn.cos.tex
%\input figsxf.ratio.tex
%\input figsxfn.ratio.tex

\begin{figure}[h]
\epsfig{file=desy-14-035.fig2a.eps,width=79mm}
\epsfig{file=desy-14-035.fig2b.eps,width=79mm}
\caption{
The fraction of DIS events with forward photons {\bf (a)} and
forward neutrons {\bf (b)} as a function of $W$
in the kinematic region given in Table~\ref{tab:selection}.
Also shown are the predictions of the LEPTO (solid line) and CDM (dashed line)
MC models. In the case of forward neutron production, the predictions of
RAPGAP-$\pi$ and the linear combinations of LEPTO and RAPGAP-$\pi$, as 
 well as  CDM and RAPGAP-$\pi$ are also shown.
}

\vspace*{-92mm}
{
\bf \large
\hspace*{67mm}(a) 
\hspace*{74mm}(b) 
		
}
\label{figW}
\end{figure}

\vspace*{92mm}



\begin{figure}[h]
\epsfig{file=desy-14-035.fig3a.eps,width=79mm}
\epsfig{file=desy-14-035.fig3b.eps,width=79mm}
\caption{
The fraction of DIS events with forward photons {\bf (a)} and
forward neutrons {\bf (b)} as a function of $W$
in the kinematic region given in Table~\ref{tab:selection}.
Also shown are the predictions of the cosmic ray hadronic interaction
models SIBYLL~2.1 (solid line), QGSJET~01 (dashed line), QGSJET~II-04 (dotted line)
and EPOS~LHC (dash-dotted line).
}
\vspace*{-85mm}
{
\bf \large
\hspace*{67mm}(a) 
\hspace*{74mm}(b) 
}
\label{figWcos}
\end{figure}
\vspace*{80mm}

\newpage



\begin{figure}[h]
\epsfig{file=desy-14-035.fig4a.eps,width=67mm}
\epsfig{file=desy-14-035.fig4b.eps,width=67mm}
\epsfig{file=desy-14-035.fig4c.eps,width=67mm}
\caption{
Normalised cross sections of forward photon production in DIS
as a function of $x_F$ in three $W$ intervals
in the kinematic region given in Table~\ref{tab:selection}.
The error bars show the total experimental uncertainty, calculated using the 
quadratic sum of the statistical and systematic uncertainties.
Also shown are the predictions of the LEPTO (solid line) and CDM (dashed line)
MC models.
}
\label{fig:xfgmc}
\end{figure}

\begin{figure}[h]
\centering
\epsfig{file=desy-14-035.fig5a.eps,width=65mm}
\epsfig{file=desy-14-035.fig5d.eps,width=65mm}

\epsfig{file=desy-14-035.fig5b.eps,width=65mm}
\epsfig{file=desy-14-035.fig5e.eps,width=65mm}

\epsfig{file=desy-14-035.fig5c.eps,width=65mm}
\epsfig{file=desy-14-035.fig5f.eps,width=65mm}

\caption{
Normalised cross sections of forward photon production in DIS
as a function of $x_F$ in three $W$ intervals
in the kinematic region given in Table~\ref{tab:selection}.
The inner error bars show the statistical uncertainty, while
the outer error bars show the total experimental uncertainty, 
calculated using the
quadratic sum of the statistical and systematic uncertainties.
Also shown are the predictions of the cosmic ray hadronic interaction
models SIBYLL~2.1 (solid line), QGSJET~01 (dashed line), 
QGSJET~01~(no~mi) (dash-double dotted line),  
QGSJET~II-04 (dotted line) and EPOS~LHC (dash-dotted line).
In the right column the ratios of the CR model predictions to the data are shown.
}
\label{fig:xfgcos}
\end{figure}


\begin{figure}[h]
\epsfig{file=desy-14-035.fig6a.eps,width=67mm}
\epsfig{file=desy-14-035.fig6b.eps,width=67mm}
\epsfig{file=desy-14-035.fig6c.eps,width=67mm}
%
\caption{
Normalised cross sections of forward neutron production in DIS
as a function of $x_F$ in three $W$ intervals
in the kinematic region given in Table~\ref{tab:selection}.
The inner error bars show the statistical uncertainty, while
the outer error bars show the total experimental uncertainty, 
calculated using the
quadratic sum of the statistical and systematic uncertainties.  
Also shown are the predictions of CDM (dotted line), RAPGAP-$\pi$ (dashed line)
and a linear combination of CDM and  RAPGAP-$\pi$ predictions (solid line).
}
\label{fig:xfn}
\end{figure}




\begin{figure}[h]
\centering
\epsfig{file=desy-14-035.fig7a.eps,width=65mm}
\epsfig{file=desy-14-035.fig7d.eps,width=65mm}

\epsfig{file=desy-14-035.fig7b.eps,width=65mm}
\epsfig{file=desy-14-035.fig7e.eps,width=65mm}

\epsfig{file=desy-14-035.fig7c.eps,width=65mm}
\epsfig{file=desy-14-035.fig7f.eps,width=65mm}

\caption{
Normalised cross sections of forward neutron production in DIS
as a function of $x_F$ in three $W$ intervals
in the kinematic region given in Table~\ref{tab:selection}.
The inner error bars show the statistical uncertainty, while
the outer error bars show the total experimental uncertainty, 
calculated using the
quadratic sum of the statistical and systematic uncertainties.  
Also shown are the predictions of the cosmic ray hadronic interaction
models SIBYLL~2.1 (solid line), QGSJET~01 (dashed line), 
QGSJET~01~(no~mi) (dash-double dotted line),  
QGSJET~II-04 (dotted line) and EPOS~LHC (dash-dotted line).
In the right column the ratios of the CR model predictions to the data are shown.
}
\label{fig:xfncos}
\end{figure}



\newpage

\begin{figure}[ht]
\epsfig{file=desy-14-035.fig8a.eps,width=74mm}
\hspace*{5mm}
\epsfig{file=desy-14-035.fig8b.eps,width=74mm}
\caption{
Ratios of normalised cross sections of forward photon production in DIS
corresponding to two different $W$ intervals, shown in Figure~\ref{fig:xfgmc},
 as a function of $x_F$:
{\bf (a)} ratio of the cross section in the $130<W<190$~GeV
interval to the cross section in the $70<W<130$~GeV interval;
{\bf (b)}   ratio of the cross section in the $190<W<245$~GeV
interval to the cross section in the $70<W<130$~GeV interval.
The kinematic phase space is defined in Table~\ref{tab:selection}.
The error bars show the total experimental uncertainty, calculated using the
quadratic sum of the statistical and systematic uncertainties.
Also shown are the predictions of the LEPTO (solid line) and 
CDM (dashed line) MC models.
}

\vspace*{-107mm}
{
\bf \large
\hspace*{17mm}(a)
\hspace*{74mm}(b)

}
\vspace*{102mm}
\label{fig:xfmcratio}
\end{figure}




\begin{figure}[h]
\epsfig{file=desy-14-035.fig9a.eps,width=74mm}
\hspace*{5mm}
\epsfig{file=desy-14-035.fig9b.eps,width=74mm}
\caption{
Ratios of normalised cross sections of forward photon production in DIS
corresponding to two different $W$ intervals, shown in Figure~\ref{fig:xfgcos}, 
as a function of $x_F$:
{\bf (a)} 
ratio of the cross section in the $130<W<190$~GeV
interval to the cross section in the $70<W<130$~GeV interval;
{\bf (b)} 
 ratio of the cross section in the $190<W<245$~GeV
interval to the cross section in the $70<W<130$~GeV interval.
The kinematic phase space is defined in Table~\ref{tab:selection}.
The error bars show the total experimental uncertainty, calculated using the
quadratic sum of the statistical and systematic uncertainties.
Also shown are the predictions of the cosmic ray hadronic interaction
models SIBYLL~2.1 (solid line), QGSJET~01 (dashed line),  
QGSJET~II-04 (dotted line) and EPOS~LHC (dash-dotted line).
}

\vspace*{-112mm}
{
\bf \large
\hspace*{17mm}(a)
\hspace*{74mm}(b)

}
\vspace*{110mm}
\label{fig:xfgcosratio}
\end{figure}


\begin{figure}[h]
\epsfig{file=desy-14-035.fig10a.eps,width=74mm}
\hspace*{5mm}
\epsfig{file=desy-14-035.fig10b.eps,width=74mm}
\caption{
Ratios of normalised cross sections of forward neutron production in DIS
corresponding to two different $W$ intervals, shown in 
Figure~\ref{fig:xfncos},
as a function of $x_F$:
{\bf (a)}
ratio of the cross section in the $130<W<190$~GeV
interval to the cross section in the $70<W<130$~GeV interval;
{\bf (b)}
 ratio of the cross section in the $190<W<245$~GeV
interval to the cross section in the $70<W<130$~GeV interval.
The kinematic phase space is defined in Table~\ref{tab:selection}.
The error bars show the total experimental uncertainty, calculated using the
quadratic sum of the statistical and systematic uncertainties.
Also shown are the predictions of the cosmic ray hadronic interaction
models SIBYLL~2.1 (solid line), QGSJET~01 (dashed line),
QGSJET~II-04 (dotted line) and EPOS~LHC (dash-dotted line).
}
\vspace*{-117mm}
{
\bf \large
\hspace*{17mm}(a)
\hspace*{74mm}(b)

}
\vspace*{110mm}
\label{fig:xfnratio}
\end{figure}



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