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


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

% Some useful tex commands
%


\newcommand{\CHECK}{{\bf CHECK}\xspace}
%\renewcommand{\mu}{\upmu}
\renewcommand{\Theta}{\theta}
%\renewcommand{\theta}{\uptheta}
%\renewcommand{\Phi}{\Upphi}
%\renewcommand{\phi}{\upphi}

\newcommand{\Cascade}{{\sc Cascade}\xspace}
\newcommand{\Pythia}{{\sc Pythia}\xspace}
\newcommand{\DiffVM}{{\sc DIFFVM}\xspace}

\newcommand{\be}{\begin{eqnarray}}
\newcommand{\ee}{\end{eqnarray}}


\newcommand{\mr}[1]{{{#1}}}

\newcommand{\kT}{\ensuremath{{k_{T}}}\xspace}

\newcommand{\GeV}{\rm GeV}
\newcommand{\TeV}{\rm TeV}
\newcommand{\pb}{\rm pb}
\newcommand{\invpb}{\rm pb^{-1}}
\newcommand{\cm}{\rm cm}
\newcommand{\hdick}{\noalign{\hrule height1.4pt}}

\newcommand{\xgobs}{$x_{\gamma}^{obs}$}
\newcommand{\xgobsm}{x_{\gamma}^{obs}}
\newcommand{\ptjet}{$p_{t}^{jet_{1}}$}
\newcommand{\ptjetm}{p_{t}^{jet_{1}}}
\newcommand{\ptjettwo}{$p_{t}^{jet_{2}}$}
\newcommand{\ptjets}{$p_{t}^{jet_{1(2)}}$}
\newcommand{\ptjetsm}{p_{t}^{jet_{1(2)}}}
\newcommand{\pt}{$p_{t}$}
\newcommand{\ptm}{p_{t}}
\newcommand{\fbbb}{$f^{b\bar{b}}$}
\newcommand{\fccb}{$f^{c\bar{c}}$}
\newcommand{\bbb}{\ensuremath{b\bar{b}}\xspace}
\newcommand{\ccb}{\ensuremath{c\bar{c}}\xspace}
\newcommand{\ra}{\rightarrow}

\newcommand{\mum}{\ensuremath{\mu m \xspace}}

\newcommand{\Mmumu}{\ensuremath{{m_{\mu\mu}}}\xspace}
\newcommand{\Mee}{\ensuremath{{m_{ee}}\xspace}}
\newcommand{\Mll}{\ensuremath{{m_{\ell\ell}}}\xspace}

\newcommand{\JPsi}{\ensuremath{{J/\psi}}\xspace}
\newcommand{\iPsi}{\ensuremath{{\psi}}\xspace}
\newcommand{\PsiPrime}{\ensuremath{{\psi(2S)}}\xspace}

\newcommand{\Pt}{\ensuremath{{P_{T}}}\xspace}
\newcommand{\PtMu}{\ensuremath{{P_{T,\mu}}}\xspace}
\newcommand{\PtLepton}{\ensuremath{{P_{T,\ell}}}\xspace}
\newcommand{\ThetaLepton}{\ensuremath{{\Theta_{\ell}}}\xspace}
\newcommand{\ThetaElectron}{\ensuremath{{\Theta_{\mr{e}}}}\xspace}
\newcommand{\ThetaMuon}{\ensuremath{{\Theta_{\mu}}}\xspace}
\newcommand{\PtTrks}{\ensuremath{{P_{T,Trks}}}\xspace}
\newcommand{\ThetaTrks}{\ensuremath{\Theta_{Trks}}\xspace}

\newcommand{\Qsquared}{\ensuremath{{Q^{2}}}\xspace}
\newcommand{\Wgp}{\ensuremath{{W_{\gamma p}}}\xspace}

\newcommand{\ThetaTrk}{\ensuremath{\Theta}\xspace}
\newcommand{\ThetaJPsi}{\ensuremath{\Theta_{\iPsi}}\xspace}
\newcommand{\ZJPsi}{\ensuremath{{z}}\xspace}
\newcommand{\PtJPsi}{\ensuremath{{P_{T,\iPsi}}}\xspace}
\newcommand{\PtJPsiSquare}{\ensuremath{{P^2_{T,\iPsi}}}\xspace}

\newcommand{\PtStarJPsi}{\ensuremath{{P^{*}_{T,\iPsi}}}\xspace}
\newcommand{\PtStarJPsiSquare}{\ensuremath{{P^{*2}_{T,\iPsi}}}\xspace}

\newcommand{\PtPtStarJPsi}{\ensuremath{{P_{T,\iPsi}(P^{*}	_{T,\iPsi})}}\xspace}

\newcommand{\CosThetaStar}{\ensuremath{\cos(\Theta^{*})}\xspace}
\newcommand{\SinThetaStar}{\ensuremath{\sin(\Theta^{*})}\xspace}

\newcommand{\CosThetaStarSquared}{\ensuremath{{cos^2(\Theta^{*})}}\xspace}

\newcommand{\PhiStar}{\ensuremath{{\phi^{*}}}\xspace}
\newcommand{\SinPhiStar}{\ensuremath{\sin(\phi^{*})}\xspace}
\newcommand{\CosPhiStar}{\ensuremath{\cos(\phi^{*})}\xspace}


\newcommand{\Yesigma}{\ensuremath{{Y_{e\Sigma}}}\xspace}
\newcommand{\Epz}{\ensuremath{{\Sigma E-p_z}}\xspace}
	

\DeclareMathAlphabet{\mathitbf}{OML}{cmm}{b}{it}
\DeclareMathSymbol{\psi}{\mathalpha}{letters}{"20}

%\xdefinecolor

\begin{document}

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



\pagestyle{empty}

\begin{titlepage}


\noindent
\begin{flushleft}
{\tt DESY 09-225    \hfill    ISSN 0418-9833} \\
{\tt December 2009}                  \\
\end{flushleft}


\noindent
%Version:     3.0 (after final reading, 17.12.2009) \\
%Editors:     M.~Steder (michael.steder@desy.de), A.~Meyer (andreas.meyer@desy.de) \\
%Referees:    H.~Jung (hannes.jung@desy.de), D.~Wegener (wegener@physik.uni-dortmund.de) \\

\vspace*{2cm}


\begin{center}
  \Large
  {\bf 
   Inelastic Production of $\mathitbf  J\mathitbf/\mathitbf\psi$ Mesons\\
in Photoproduction and Deep Inelastic Scattering\\
at HERA
  }	

  \vspace*{2cm}
    {\Large H1 Collaboration} 
\end{center}
\vspace*{2cm}

%\linenumbers

\begin{abstract}
\noindent
A measurement is presented of inelastic photo- and electroproduction of \JPsi 
mesons in
$ep$ scattering at HERA. The data were recorded with the H1 detector 
in the period from 2004 to 2007. 
Single and double differential cross sections are determined and the 
\mbox{helicity} distributions of the \JPsi mesons are analysed. 
The results are compared to theoretical predictions in the colour singlet model and in the framework of non-relativistic QCD.
Calculations in the colour singlet model using a \kT factorisation ansatz are able to give a good description of the data, 
while colour singlet model calculations to next-to-leading order in collinear factorisation underestimate the data.


\end{abstract}

\vspace{1cm}
\centerline{Submitted to \EJC}

\end{titlepage}


\renewcommand{\thepage}{\arabic{page}}
\setcounter{page}{2}
\noindent
%\include{h1auts}
%-- H1AUTS Author list by names 
%-- Status: Tue Dec  1 09:22:13 CET 2009  Number of authors = 228 

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


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

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


\bigskip
 $ ^a$ Supported by the Bundesministerium f\"ur Bildung und Forschung, FRG,
      under contract numbers 05H09GUF, 05H09VHC, 05H09VHF,  05H16PEA \\
 $ ^b$ Supported by the UK Science and Technology Facilities Council,
      and formerly by the UK Particle Physics and
      Astronomy Research Council \\
 $ ^c$ Supported by FNRS-FWO-Vlaanderen, IISN-IIKW and IWT
      and  by Interuniversity Attraction Poles Programme,
      Belgian Science Policy \\
 $ ^d$ Partially Supported by Polish Ministry of Science and Higher
      Education, grant PBS/DESY/70/2006 \\
 $ ^e$ Supported by the Deutsche Forschungsgemeinschaft \\
 $ ^f$ Supported by VEGA SR grant no. 2/7062/ 27 \\
 $ ^g$ Supported by the Swedish Natural Science Research Council \\
 $ ^h$ Supported by the Ministry of Education of the Czech Republic
      under the projects  LC527, INGO-1P05LA259 and
      MSM0021620859 \\
 $ ^i$ Supported by the Swiss National Science Foundation \\
 $ ^j$ Supported by  CONACYT,
      M\'exico, grant 48778-F \\
 $ ^k$ Russian Foundation for Basic Research (RFBR), grant no 1329.2008.2 \\
 $ ^l$ This project is co-funded by the European Social Fund  (75\%) and
      National Resources (25\%) - (EPEAEK II) - PYTHAGORAS II \\
}

\newpage

%\linenumbers

\pagestyle{plain}
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Introduction}
\label{intro}

The description of the process 
of charmonium production in interactions of photons and hadrons is a challenge to theory,
since it involves both the production of the heavy quark system and the formation of the bound state. 
Charmonium production in electron\footnote{In this paper "electron" is used to denote both electron and positron.}-proton collisions at HERA is dominated by photon-gluon fusion: 
a photon emitted from the incoming electron interacts with a gluon from 
the proton to produce a $c \bar c$ pair that evolves into a charmonium state.
In the colour singlet model, only those states with the same quantum numbers as the 
resulting charmonium contribute to the formation of a bound \ccb state.
This is achieved by radiating a hard gluon in a perturbative process.
In the factorisation ansatz of non-relativistic quantum chromodynamics, also colour octet \ccb states
contribute to the charmonium production cross section via soft gluon radiation.	


Previous measurements in electroproduction ($ep$) and photoproduction ($\gamma p$) at HERA~\cite{Aid:1996dn,Adloff:1999zs,Adloff:2002ex,Adloff:2002ey,Breitweg:1997we,Chekanov:2002at,ZEUS05} are not described by predictions in the colour singlet model to leading order.
In contrast, the calculation of photoproduction cross sections to next-to-leading order (NLO)~\cite{Kraemer} showed a reasonable description of the photoproduction cross sections.
The calculation proved that the corrections with respect to leading order results are very large, 
increasing towards large transverse momentum of the \JPsi meson. 
The same calculation, repeated recently with an up-to-date set of theoretical parameters~\cite{Maltoni}, results in a prediction which is about a factor of three below the measured cross sections, 
indicating that corrections beyond NLO are needed and/or that contributions 
from colour octet states may be sizable.

In this paper a measurement is presented of inelastic \JPsi meson production at HERA.
The measurement uses a larger data sample than previous results~\cite{Aid:1996dn,Adloff:1999zs,Adloff:2002ex,Adloff:2002ey} and benefits from improved systematics.
The data sets were collected in the years 2004 to 2007 with the H1 detector.
The \JPsi meson candidates are identified by the leptonic decay into two 
muons or electrons.
The cross sections are measured for both electroproduction and photoproduction. 
For the photoproduction sample \JPsi meson polarisation variables are determined. 
The data samples are restricted to the region of phase space where contributions from 
diffractive charmonium production are suppressed.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Theoretical Models}
\label{theory}

In order to describe inelastic charmonium production in the framework of 
perturbative QCD different models have been proposed, such as the colour-eva\-po\-ration model~\cite{Halzen:1977rs,Eboli:1998xx}, the colour-singlet model 
(CSM)~\cite{CSM_Chang,Berger:1980ni,Baier:1981uk,Baier:1981zz,Baier:1983va}, the factorisation ansatz in non-relativistic quantum chromodynamics (NRQCD)~\cite{Caswell:1985ui,Thacker:1990bm,Bodwin:1994jh} and soft 
colour interactions~\cite{Edin:1997zb}. 
In this paper the most recent calculations using the CSM or NRQCD are compared to the data.

In the CSM, only charm quark pairs in a colour singlet state with the same quantum numbers as the resulting charmonium contribute to the formation of a bound \ccb state.
This is achieved by radiating a hard gluon in the perturbative process.
The factorisation ansatz in NRQCD includes also colour octet $c\bar{c}$ states 
in the charmonium production cross section. 
The size of these colour octet contributions, described by long distance 
matrix elements (LDME), is defined by additional free parameters which
were determined in fits to the Tevatron data~\cite{Braaten:1999qk}.
The NRQCD factorisation approach contains also the colour singlet model
which is recovered in the limit in which the colour-octet LDME tend to zero.

The following calculations are compared to the measurements presented in this paper:
%
\begin{itemize}
\item A calculation of \JPsi meson photoproduction via a colour singlet mechanism \cite{Maltoni} 
provides predictions for both cross sections and helicity distributions to next-to-leading order. 
The uncertainty of this calculation is estimated by variations of the charm quark mass  
and the factorisation and renormalisation scales.
%
\item A calculation at NLO for photoproduction cross sections includes the full framework 
of NRQCD~\cite{Kniehl}. 
The uncertainty of this calculation is dominated by the limited knowledge of 
the LDMEs.
%
\item CSM predictions in the \kT factorisation approach are employed as implemented in the MC generator \Cascade \cite{Jung:2000hk}.
Higher order parton emissions based on the CCFM evolution equations~\cite{ccfm} are matched to ${\cal O}(\alpha_s)$ matrix elements in which the incoming parton can be off-shell.
The uncertainty on the calculation is estimated by varying the renormalisation scale by a factor of two.
In addition polarisation variables in the \kT factorisation approach are calculated analytically~\cite{Baranov}. 
%
\end{itemize}
Parameters and variations used in the theoretical calculations are given in table~\ref{tab:theoryparam}.

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

\begin{table}[tb]
\begin{center}
\begin{tabular}{ l p{8cm}}
%
%
%% \hline
%%  CSM (NLO), M.~Kr\"amer~\cite{Kraemer} &  \\
%% \hline
%% PDF                                                        & MRST 
%% renormalisation and factorisation scale                    & $\sqrt{2}m_c, \mr{max}\left[\sqrt{2}m_c,1/2\cdot\sqrt{m_c^2 + \PtJPsiSquare}~\right] (\mathrm{for}~\PtJPsiSquare) $\\
%% CS LDME  & $\left<\mathcal{O}\left[\underline{1},^3S_1\right]\right> = \unit[1.16]{GeV^3}$\\
%% $m_c$                                                      & $1.3 < m_c < \unit[1.5]{GeV}$\\
%% $\alpha_s(M_Z)$                                            & $0.1200 \pm 0.0025$\\

%
%
\toprule
 CSM (NLO), P.~Artoisenet et al.~\cite{Maltoni} &  \\
\midrule
PDF                                                        & CTEQ6M~\cite{CTEQ6M}\\
renormalisation and factorisation scale                    & $\mu_0 = 4 m_c$\\
scale variation                                            & $0.5 \mu_0 < \mu_\mr{f}$, $\mu_\mr{r} < 2 \mu_0$ and $0.5 < \mu_\mr{r}/\mu_\mr{f} < 2$\\
CS LDME  & $\left<\mathcal{O}\left[\underline{1},^3S_1\right]\right> = \unit[1.16]{GeV^3}$\\
$m_c$                                                      & $1.4 < m_c < \unit[1.6]{GeV}$\\
$\alpha_s(M_Z)$                                            & $0.118$ (+ running at 2 loops)\\
%
%
\midrule
 NRQCD (NLO), M.~Butensch\"on et al.~\cite{Kniehl} &  \\
\midrule
PDF                                                        & CTEQ6M~\cite{CTEQ6M}\\
renormalisation and factorisation scale                    & $\mu_0 = \sqrt{4m_c^2 + \PtJPsiSquare}$\\
NRQCD scale                                                & $\mu_\Lambda = \mr{m_c}$\\
%CS LDME  & $\left<\mathcal{O}\left[\underline{1},^3S_1\right]\right> = \unit[1.16]{GeV^3}$\\
$m_c$                                                      & $m_\JPsi / 2 \approx \unit[1.55]{GeV}$\\
$\alpha_s(M_Z)$                                            & $0.1176 \pm 0.002$\\
%
%
\midrule
 CSM (\kT factorisation), \Cascade~\cite{Jung:2000hk} &  \\
\midrule
PDF                                                        & CCFM set A0~\cite{SetA0}\\
                                                           &  (\lq set A0$^\pm$\rq\xspace for $\mu_r$ uncertainties)\\
renormalisation scale                                      & $\mu_0 = \sqrt{m_\psi^2 + \PtJPsiSquare}$\\
renormalisation scale variation                            & $0.5 \mu_0 < \mu_\mr{r} < 2 \mu_0$\\
factorisation scale                                        & $\sqrt{\hat{s} + Q^2_\perp}$\\
%CS LDME  & $$\left<\mathcal{O}\left[\underline{1},^3S_1\right]\right>$ = 1.16 \unit[]{GeV^3}$\\
$m_c$                                                      & $\unit[1.5]{GeV}$\\
$\Lambda_\mathrm{QCD}^{(3)}$                                  & $\unit[200]{MeV}$\\
%
%
\midrule
 CSM (\kT factorisation), S.~Baranov~\cite{Baranov} &  \\
\midrule
PDF                                                        & CCFM set A0~\cite{SetA0}\\
renormalisation and factorisation scale                    & $\mu_0 = \sqrt{m_\psi^2 + \PtJPsiSquare}$\\
%$\left<\mathcal{O}\left[\underline{1},^3S_1\right]\right>$ & $\unit[]{GeV^3}$\\
$m_c$                                                      & $\unit[1.5]{GeV}$\\
$\Lambda_\mathrm{QCD}^{(3)}$                                  & $\unit[200]{MeV}$\\
\bottomrule
\end{tabular}
\end{center}
\caption{Summary of the parameters employed in the CSM and NRQCD calculations 
used to compare to the measurements in this paper. In this table PDF means 
parton distribution function of the proton, $\hat{s}$ denotes the invariant mass square of the hard subprocess and $Q_\perp$ the initial transverse momentum of the partonic system ($\gamma$g).}
\label{tab:theoryparam}
\end{table}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{H1 Detector}
\label{sec:h1}

The H1 detector is described in detail elsewhere~\cite{Abt:1997xv}. Here only
the components essential to the present analysis are briefly described.
A right handed Cartesian  coordinate system is used with the origin at the nominal primary $ep$ interaction vertex. 
The proton beam direction defines the $z$ axis. The polar angles $\theta$ and transverse momenta $P_T$ of all particles 
are defined with respect to this axis.  The azimuthal angle $\phi$ defines the particle direction in the transverse plane. The pseudorapidity  is defined as $\eta=-\ln {\tan {\frac{\theta}{2}}}$.

Charged particles emerging from the $ep$ interaction region
are measured by the central tracking detector (CTD) in the pseudo-rapidity
range $|\eta| < 1.74$.
The CTD consists of two large cylindrical central jet drift chambers (CJC) which
are interleaved by a $z$-chamber and arranged concentrically 
around the beam-line in a magnetic field of $1.16~\mathrm{T}$. 
The CTD provides triggering information based on track segments from the CJC~\cite{FTT, FTT2},
and on the \ZJPsi-position of the vertex from the 5-layer multi-wire proportional chamber~\cite{CIP} which
is situated inside the inner CJC.
To provide the best possible spatial track reconstruction, CTD tracks are linked to 
hits in the vertex detector, the central silicon tracker CST~\cite{CST}.
The CST is installed close to the interaction point, surrounding the beam pipe in the 
pseudo-rapidity range $|\eta| < 1.3$ and consists of two layers of double sided silicon 
strip sensors. 

Charged and neutral particles are measured in the liquid argon calorimeter (LAr)~\cite{LAr} which surrounds the tracking chambers and covers the range 
$-1.5 < \eta < 3.4$ and a lead/scintillating-fibre calorimeter SpaCal~\cite{Nicholls:1996di}, 
covering the backward region $-4.0 < \eta < -1.4$.
The calorimeters are surrounded by the solenoidal 
magnet and the iron return yoke. The yoke is instrumented with
16 layers of limited streamer tubes, forming the central muon detector (CMD) 
in the range $-2.5 < \eta < 3.4$.

The luminosity determination is based on the measurement of the 
Bethe-Heitler process $ep \to ep\gamma$, where the photon is detected 
in a calorimeter located downstream of the interaction 
point in the electron beam direction at $z=-\unit[104]{m}$.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Data Analysis}
\label{sec:anal}

The kinematics of inelastic charmonium production at HERA are  
described using the following variables: the square of the $ep$ centre 
of mass energy $s = (p+k)^2$, where $p$ and $k$ denote the four vectors of electron and proton respectively; 
the negative squared four momentum transfer $\Qsquared = -q^2$, where $q$ is the
four vector of the virtual photon; and the mass of the hadronic final state $\Wgp = \sqrt{(p+q)^2}$. 
\Wgp is related to the scaled energy transfer $y = (p\cdot q) / (p \cdot k)$ via $\Wgp^2 = ys-\Qsquared$.  
In addition, the elasticity of the \JPsi meson production process is defined as
$\ZJPsi = (p_\psi\cdot p)/(q\cdot p)$, where $p_\psi$ is the four momentum of the \JPsi meson.
The elasticity denotes the fractional energy of the photon 
transferred to the \JPsi meson in the proton rest system.

Events are selected separately in the photoproduction and electroproduction regimes.
Photoproduction events are selected 
by requiring that no isolated high energy electromagnetic cluster, 
consistent with a signal from a scattered electron, is detected in 
the calorimeters. This limits the virtuality to values of 
$\Qsquared \lesssim \unit[2.5]{GeV^2}$, resulting in a mean value of $\langle \Qsquared\rangle \approx \unit[0.085]{GeV^2}$. 
Conversely, for the electroproduction sample,
a scattered electron with energy of more than $\unit[10]{GeV}$
 is required to be 
reconstructed in the backward calorimeter (SpaCal), corresponding to a range
in photon virtuality $3.6 < \Qsquared < \unit[100]{GeV^2}$.

In this analysis the photon virtuality \Qsquared is reconstructed from the
scattered electron energy $E_e'$ and polar angle $\Theta_e'$ as $\Qsquared=4E_eE_e'\cos^2(\Theta_e'/2)$, where $E_e$ denotes the energy of the beam electron.
The variable $y$ is reconstructed using 
the relation $y=\sum_h (E-p_z)/2E_e$ for photoproduction~\cite{jb} 
and $y=\sum_h (E-p_z)/\sum (E-p_z)$ for electroproduction~\cite{esigma}.
The sums in the numerator include all particles of the hadronic final state without the scattered electron, which is only included in the sum of the denominator for electroproduction.
The elasticity \ZJPsi is then obtained from 
$z = (E-p_z)_{\JPsi}/\sum_h (E-p_z)$, where $(E-p_z)_{\JPsi}$ is calculated from the decay particles of the \JPsi meson.
The kinematics of the final state particles are obtained from charged particle
tracks reconstructed in the CTD and energy depositions in the LAr and 
SpaCal calorimeters~\cite{hadrooII,peez}. 
\begin{figure}[t]
\unitlength1cm
%\frame{
\begin{picture}(16.,4.6)
\put( 0.,0.2){\includegraphics*[width=5.cm]{d09-225f1a.eps}}
\put( 5.5,0.2){\includegraphics*[width=5.cm]{d09-225f1b.eps}}
\put(11.,0.2){\includegraphics*[width=5.cm]{d09-225f1c.eps}}
\put(1.,4.2){a)}
\put(6.5,4.2){b)}
\put(12.,4.2){c)}
\end{picture}
%}
\caption{Invariant mass spectra of two oppositely charged leptons after all selection
cuts for a) the photoproduction ($\gamma p$) sample and b,c) the 
electroproduction ($ep$) samples as listed in table~\ref{tab:selection}. 
The lines show the results of fits to signals and backgrounds.}
\label{fig:ctrl:mass} 
\end{figure}

\begin{table}[t]
\begin{center}
\begin{tabular}{ rccc }
\toprule
 & Photoproduction & \multicolumn{2}{c}{Electroproduction} \\
 & $\JPsi \rightarrow \mu^+\mu^-$ & $\JPsi \rightarrow \mu^+\mu^-$ & $\JPsi \rightarrow e^+e^-$ \\
 \midrule
 \multicolumn{4}{c}{kinematic range} \\
\midrule
              & $ \Qsquared< \unit[2.5]{GeV^2}$      & \multicolumn{2}{c}{$3.6 < \Qsquared < \unit[100]{GeV^2}$} \\
                & \multicolumn{3}{c}{$ \unit[60]{GeV} < \Wgp < \unit[240]{GeV}$} \\ 
                & \multicolumn{1}{c}{$\PtJPsi > \unit[1]{GeV}$} &\multicolumn{2}{c}{$\PtStarJPsi > \unit[1]{GeV}$} \\ 
                & \multicolumn{3}{c}{$0.3 < \ZJPsi < 0.9$} \\ 
 \midrule
 \multicolumn{4}{c}{event selection} \\
 \midrule
           & \multicolumn{3}{c}{$\PtLepton >\unit[800]{MeV}$} \\ 
        & $20^{\circ} < \ThetaMuon < 160^{\circ}$ & $20^{\circ} < \ThetaMuon < 160^{\circ}$ & $20^{\circ} < \ThetaElectron < 150^{\circ}$ \\
      & \multicolumn{3}{c}{$\mathrm{N_{Trk}} \geq 5$ ~~(in the range $20^\circ < \ThetaTrk < 160^\circ$)} \\ 
 \midrule
 \multicolumn{4}{c}{event samples} \\
 \midrule
  $N_{\JPsi}$ & $2320 \pm 54$ & $501\pm34$  &  $290\pm24$ \\ 
%  $\mathcal{BR}$\cite{PDG} & $(5.93 \pm 0.06)\%$ & $(5.93 \pm 0.06)\%$ & $(5.94 \pm 0.06)\%$  \\
  $\mathcal{L}_\mr{int}$     & $\unit[165]{pb^{-1}}$ & $\unit[315]{pb^{-1}}$ & $\unit[315]{pb^{-1}}$ \\
\bottomrule
\end{tabular}
\end{center}
\caption{List of selection cuts and event yields for each of the three data samples.}
\label{tab:selection}
\end{table}

The \JPsi meson candidates are reconstructed through their decays
into two oppositely charged muons or electrons. These decay leptons
are reconstructed as charged particles in the CTD with a transverse momentum of at least $\unit[800]{MeV}$.
Muon candidates are identified as minimum ionising particles in the 
LAr calorimeter or through track segments in the CMD ($20^\circ < \ThetaMuon < 160^\circ$)~\cite{msteder}. 
Electron candidates are identified through their energy deposit 
in the central calorimeter ($20^\circ < \ThetaElectron < 150^\circ$)~\cite{msauter}.
For trigger reasons the photoproduction sample of \JPsi meson events is  
restricted to decays into $\mu^+\mu^-$, while the electroproduction sample includes
both leptonic decay channels. 
The photoproduction sample was recorded in the years 2006 and 2007 and
corresponds to an integrated luminosity of $\mathcal{L} = \unit[165]{pb^{-1}}$, while the electroproduction sample was recorded in the years 2004 to 2007 and corresponds to an integrated luminosity of $\mathcal{L} = \unit[315]{pb^{-1}}$.



The measurement is performed in the kinematic range 
$60 < \Wgp < \unit[240]{GeV}$, $0.3 < \ZJPsi < 0.9$ and $\PtPtStarJPsi > \unit[1]{GeV}$.
In photoproduction the transverse momentum \PtJPsi is measured in the lab frame, while in electroproduction the transverse momentum \PtStarJPsi is calculated in the $\gamma^
*p$ rest frame.
To suppress contributions from diffractive 
production of \JPsi and \PsiPrime mesons, selected events are required to contain at least 
five reconstructed tracks in the central region of the detector 
($20^\circ < \ThetaTrk < 160^\circ$). 
The reconstruction efficiency accounts for this experimental cut and the measured cross sections are corrected for this track multiplicity cut.


Figure~\ref{fig:ctrl:mass} shows the invariant mass spectra of the leptons in the selected event samples. 
The number of signal events, $N_\JPsi$, is obtained in all bins of the cross section measurements from a fit to the mass distributions in the interval $2 < \Mll < 6$ GeV.
For the decay into muons the signal peak is described using a modified Gaussian~\cite{ZEUS05}. 
In the case of a decay into two electrons an exponential is added to the lower mass region of the signal Gaussian in order to account for the radiative tail~\cite{msteder}.
For the differential cross section measurements, the width and asymmetry term of
the mass peak in each bin are fixed to the values obtained from the full samples. 
For both decay channels, the background is parametrised by a polynomial of third order. 
At $m_{\ell\ell} \approx \unit[3.7]{GeV}$, the nominal mass of \PsiPrime mesons, an additional Gaussian with fixed position and width is allowed in all analysis bins.

The selection criteria and the obtained event samples are summarised in table~\ref{tab:selection}.



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Monte Carlo Simulations}
\label{sec:mc}

Cross sections and polarisation parameters are derived by correcting the measured number of events
and angular distributions for detector effects, such as detector 
resolutions and inefficiencies.
Several Monte Carlo generator programs are used to determine
the corrections. 
All samples are passed through a detailed simulation 
of the H1 detector response based on the GEANT program~\cite{geant} and through the
same reconstruction and analysis algorithms as used for the data.

Signal events are generated using the Monte 
Carlo generator \Cascade \cite{Jung:2000hk}. 
Elastic and proton-dissociative production of \PsiPrime mesons is simulated using  
\DiffVM~\cite{diffvm} with parameters tuned to describe the results of previous H1 measurements~\cite{psi2s97,psi2s02}.
% Psi(2s) 2002: pd: 0.59+/-0.13+/-0.12 / elas: 4.31+-0.57+-0.46
% \DiffVM:    : pd: 1.07                       4.57
The Monte Carlo generator \Pythia~\cite{pythia} is used for the description of 
the contribution from $b$ hadron decays as described in section~\ref{sec:bg}. 
All generators use the JETSET part of the \Pythia program~\cite{pythia} 
to simulate the hadronisation and decay processes.


Signal events as simulated with the Monte Carlo generator \Cascade are compared with the 
data after final selection in figures~\ref{fig:ctrl:muon} and \ref{fig:ctrl:dis}.
All data distributions in these figures are corrected for contributions from non-resonant 
background events using a sideband method described in~\cite{msteder}.

Corrections as a function of \Wgp and \PtJPsi in bins of the elasticity \ZJPsi 
are applied to the \Cascade Monte Carlo simulation in order to describe the data. 
Details of the procedure are described in~\cite{msteder}. 
In figure~\ref{fig:ctrl:muon} distributions for the photoproduction sample are compared to \Cascade Monte Carlo predictions before and after correction for the observables \PtMu, \ThetaMuon, \PtJPsi, \ThetaJPsi, \Wgp and \ZJPsi.
Similarly, in figure~\ref{fig:ctrl:dis}, the summed distributions for the two 
electroproduction samples ($\JPsi \ra \mu^+\mu^-$ and $\JPsi \ra e^+e^-$) are shown for the observables \PtStarJPsi, \ThetaJPsi, 
\Qsquared, $\Sigma P_{T,\rm charged}$, \Wgp and \ZJPsi. Here, $\Sigma P_{T,\rm charged}$ 
is the scalar sum over the transverse momenta of all measured charged particles except for  
the scattered electron and the \JPsi meson decay leptons.
The corrected \Cascade simulation gives a good description of all aspects of the data and
is used to correct the data for losses due to limited acceptance and efficiency
of the detector.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Backgrounds}
\label{sec:bg}
Remaining backgrounds to prompt \JPsi meson production in the selected sample originate 
from feed-down processes, i.e.\,\JPsi mesons produced in decays of diffractively or 
inelastically produced \PsiPrime mesons and $\chi_c$ mesons or of $b$ hadrons. 

Inelastic production of \PsiPrime mesons with a subsequent decay  into \JPsi mesons is expected
to contribute about $15-\unit[20]{\%}$ to the selected \JPsi meson samples~\cite{ZEUSPsi2S, Kraemer}. Since the production processes are the same, the inelastic \PsiPrime mesons show similar dependences on the kinematic variables.

Diffractive production of \PsiPrime mesons contributes at large values of \ZJPsi by decays into a \JPsi meson and two charged pions. 
These events typically contain three or four reconstructed
charged tracks in the central detector ($20^\circ < \ThetaTrk < 160^\circ$). 
In figure~\ref{fig:ctrl:bgandcuts}a) the distribution of the charged track 
multiplicity measured in the central detector is shown for the photoproduction 
sample selected using all selection criteria given in table~\ref{tab:selection} except for the track multiplicity cut, which is relaxed to $\mathrm{N_{Trk}} \geq 3$.
The data are described by the sum 
of the \Cascade simulation and the prediction for diffractive \PsiPrime production, 
as simulated using the \DiffVM Monte Carlo generator. 
In the final selection remaining contributions from diffractive \PsiPrime meson production amount
to about 1.3\% in the total sample and to about 5\% in the highest elasticity bin, $0.75 < \ZJPsi < 0.9$.

The fraction of events arising from $b$ hadrons decaying into $\JPsi+X$ is 
estimated using the \Pythia simulation. The \Pythia prediction is scaled by
a factor of 2, based on results from previous measurements of beauty 
production at HERA~\cite{btomujj,btojj}. 
This scaled prediction by \Pythia amounts to 5\% in the total sample and 
about 20\% in the lowest \ZJPsi bin. It is confirmed within uncertainties by the following 
determination using data.
The fraction of events in the photoproduction sample containing $b$ hadrons
is estimated using the impact parameter of the decay muons to exploit the 
lifetime signature of $b$ hadrons. 
The impact parameter, $\delta$, of 
the decay muon tracks is defined as the distance of closest approach 
in the transverse plane to the reconstructed
primary vertex. The sign of the impact parameter is defined as positive if the 
angle between the decay muon and the \JPsi meson momentum direction is less than 
$\unit[90]{^\circ}$, and is defined as negative otherwise.
A signed significance $\mathcal{S} = \delta/\sigma(\delta)$ is reconstructed by weighting the reconstructed signed impact parameter with its uncertainty~\cite{CSTresol}.
Figure \ref{fig:ctrl:bgandcuts}b) shows the distribution of the signed 
significance for events in the interval $0.3 < \ZJPsi < 0.4$.
The histogram is filled with the signed significance of the decay muons 
for all events where both muon tracks have at least one hit in the CST.
The fraction of events coming from the decay of $b$ hadrons
is obtained from a fit of the significance distribution of \Cascade 
(simulating prompt \JPsi meson production) plus \Pythia 
(simulating \bbb events with subsequent decays into \JPsi + $X$) to that of the data.
The fit results are dominated by the region of small signed significances, $\mathcal{S} < 3$, due to large statistical uncertainties at larger values of $\mathcal{S}$.
The distribution of the data is corrected for non-resonant contributions using the side bands~\cite{msteder}.
The relative contribution from $b$ hadrons as resulting from the fits are 
shown in figure \ref{fig:ctrl:bgandcuts}c) for three bins of \ZJPsi. 
The scaled predictions from \Pythia are found to be in good agreement with the 
measured fractions, indicating that the background from $b$ hadrons is under control.

The contribution from $\chi_c$ production and decay was studied~\cite{Adloff:2002ex}
and found negligibly small in the present kinematic region, $0.3 < \ZJPsi < 0.9$.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Systematic Uncertainties}
\label{sec:syserr}

\begin{table}[tb] \centering
\begin{tabular}{lccc} 
\toprule
Source & \multicolumn{3}{c}{Uncertainty [$\%$]}\\
 \midrule  
 & Photoproduction & \multicolumn{2}{c}{Electroproduction} \\
 & $\JPsi \rightarrow \mu^+\mu^-$ & $\JPsi \rightarrow \mu^+\mu^-$ & $\JPsi \rightarrow e^+e^-$ \\
\midrule
Decay leptons reconstruction       &  $1$  & $1$ & $2$ \\
Decay leptons identification       &  $3$  & $3$ & $3$ \\
Number of signal events            &  $2$  & $2$ & $4$ \\
Trigger                            &  $3$  & $2$ & $2$ \\
Scattered electron energy scale    & ---   & $2$ & $2$ \\
Hadronic final state energy scale  &  $4$  & $3$ & $3$ \\
Integrated luminosity              &  $4$  & $3.2$ & $3.2$ \\
Model uncertainties                &  $5$  & $5$ & $5$ \\
Decay branching ratio              &  $1$  & $1$ & $1$ \\
\midrule
Sum                               &  $9.0$  & $8.2$ & $9.1$ \\
\bottomrule


\end{tabular}
\caption{Systematic uncertainties of the \JPsi meson production cross section. 
The total systematic uncertainty is the sum of the contributions 
added in quadrature.}
\label{tab:sys}
\end{table}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

The sources of systematic uncertainties of the cross section measurement are listed in table~\ref{tab:sys} and are detailed in the following:

\begin{itemize}

\item 
The uncertainty on the cross section due to the track and vertex reconstruction efficiency has been 
determined to be $\unit[1]{\%}$ for $\JPsi \ra \mu\mu$ and $\unit[2]{\%}$ for $\JPsi \ra ee$.

\item
The efficiency for the identification of the leptons is determined using a high statistics sample of events of elastically produced \JPsi mesons~\cite{msteder}.
The detector simulation is reweighted to match the efficiency measured in the data as necessary. 
Remaining differences are smaller than $\unit[3]{\%}$ everywhere and are taken as systematic uncertainty. 

\item 
The systematic uncertainty on the determination of the number of 
signal events, obtained by a fit to the mass distributions in every analysis bin, is determined by a variation of the extraction method.
Comparing the number of signal events for binned and unbinned 
log-likelihood fits yields a systematic uncertainty of $\unit[0.5]{\%}$.
In addition, the result from the fit to background and signal is compared 
to the number of signal events above the fitted background function 
in the mass window between $2.95$ and $\unit[3.2]{GeV}$.
An uncertainty of $\unit[2]{\%}$ for the decay into muons and $\unit[4]{\%}$ for the electrons is found. 
The uncertainty for the electron is larger due to an additional uncertainty originating from the description of the radiative tail.

\item 
The trigger efficiencies are determined using independent trigger channels. 
For the electroproduction sample the trigger efficiency is measured to be $\unit[(97 \pm 2)]{\%}$.
In the photoproduction sample the trigger efficiency depends mainly on the identification of the decay muons in the central muon system.
The efficiency amounts to about $\unit[70]{\%}$ with a systematic uncertainty of $\unit[3]{\%}$.
A detailed description of the determination of the trigger efficiencies can be found in~\cite{msteder}.


\item For the electroproduction sample the measurement of the scattered 
electron energy is known with a scale uncertainty of $\unit[1]{\%}$. The uncertainty of the scattering angle is 1 mrad.
Both uncertainties combined lead to an uncertainty of the cross section 
measurement of $\unit[2]{\%}$ on average.

\item 
The hadronic energy scale uncertainty is $\unit[4]{\%}$ in the LAr and  $\unit[7]{\%}$ in the SpaCal.
This leads to an uncertainty on the cross sections measurement of  $\unit[3]{\%}$ for the electroproduction sample and $\unit[4]{\%}$ for the photoproduction sample.

\item 
The integrated luminosity is known to a precision of $\unit[3.2]{\%}$ for the electroproduction sample and $\unit[4.0]{\%}$ for the photoproduction sample. 

\item 
The dependence of the result on model assumptions made in the \Cascade
Monte Carlo simulation were investigated and found to amount to $\unit[5]{\%}$ in total.
The model uncertainty arising from the knowledge of the decay angular 
distributions, explained in section~\ref{sec:angles}, is determined 
by variation of the parameter $\alpha$  in the simulation by $\pm 0.3$.
This variation results in a change of the cross section of up to $\unit[4]{\%}$.
The systematic uncertainty originating from the uncertainty of the slope of the \PtJPsi(\PtStarJPsi) distribution 
in the simulation is determined by a variation of this distribution as described in~\cite{msteder}.
This variation results in a change of the cross section of up to $\unit[4]{\%}$.


\item
The branching ratios of the leptonic decay channels of the \JPsi meson are known with an accuracy of $\unit[1]{\%}$ \cite{PDG}.


\end{itemize}

The total systematic uncertainty is obtained by adding all the above contributions 
in quadrature. 
A total systematic uncertainty of $\unit[9]{\%}$ is determined for the photoproduction sample.
For the combined electroproduction cross section the total systematic uncertainty is $\unit[8.5]{\%}$
The same uncertainties are attributed to all bins of the cross section measurement. 
For the measurement of the helicity distributions only the uncorrelated systematic uncertainties
are taken into account. 
They amount to about $\unit[3.5]{\%}$ and are negligible compared to the statistical uncertainties.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Cross Section Measurements}
\label{sec:xsec}

The cross section measurement is performed in the kinematic range $60 < \Wgp < \unit[240]{GeV}$, $0.3 < \ZJPsi < 0.9$ and $\PtPtStarJPsi > \unit[1]{GeV}$. The photon virtuality \Qsquared is limited in the electroproduction analysis to $3.6 < \Qsquared < \unit[100]{GeV^2}$ and for the photoproduction sample to $\Qsquared < \unit[2.5]{GeV^2}$.
 
For the measurement of differential cross sections the number of signal 
events in each bin is corrected for detector inefficiencies and acceptance
and normalised to integrated luminosity and branching ratio. 
They are not corrected for QED radiative effects.
The electroproduction cross sections, measured from $\JPsi \ra \mu\mu$ and $\JPsi \ra ee$, are combined~\cite{msteder}.
The differential cross sections are bin-centre corrected using MC simulations.
In order to avoid model dependencies, the measured cross sections are not corrected for contributions from backgrounds as described in section~\ref{sec:bg}.
All measured cross sections are listed in tables~\ref{tab:xsecs:gammap:pt2} -- \ref{tab:xsecs:dis:Z_PtStar} together with statistical and systematic uncertainties.

For the photoproduction sample the measured $ep$ cross sections are transformed to $\gamma p$ 
cross sections using the photon flux factors presented in
table~\ref{tab:fluxes}, calculated in the Weizs\"acker Williams approximation~\cite{WWA}.
The differential \JPsi meson photoproduction cross section is measured as function
of the elasticity \ZJPsi and the squared transverse momentum \PtJPsiSquare 
of the \JPsi meson. The total $\gamma p$ cross section is measured in bins of the 
photon proton centre of mass energy \Wgp.
The results are displayed in figure~\ref{fig:res:gp:xsecs:1d} and 
show a reasonable agreement with
the prediction from the \Cascade MC generator.
A variation of the renormalisation scale by a factor of two ($0.5 \mu_0 < \mu_\mr{r} < 2 \mu_0$) has little effect as shown by the band in the figures.
In addition to the \Cascade prediction, the remaining contributions from
diffractive \PsiPrime mesons and from $b$ hadrons are shown.
The distributions in \PtJPsiSquare and \ZJPsi are further investigated by 
dividing the sample into bins of \PtJPsiSquare and \ZJPsi, respectively as shown in figure~\ref{fig:res:gp:xsecs:2d}.
The \ZJPsi distribution tends to flatten off towards larger
values of \PtJPsi presented in figure~\ref{fig:res:gp:xsecs:2d}a).
It can be seen that differences between the data and the \Cascade prediction are localised at low elasticities and low transverse momenta of the \JPsi mesons, where \Cascade overshoots the data, and at large elasticities and large transverse momenta, where \Cascade is below the data.
Taking into account that the measured cross section in the lowest elasticity bin
includes a significant fraction of about $\unit[20]{\%}$ of events originating
from $b$ hadron decays, the difference to the \Cascade prediction is even more
significant.

Results for electroproduction are shown in figure~\ref{fig:dis:xsecs2} and figure~\ref{fig:res:xsecsptz}.
Differential $ep$ cross sections are measured as functions of the photon virtuality \Qsquared, the squared transverse momentum of the \JPsi meson in the photon proton rest frame \PtStarJPsiSquare, the energy \Wgp and the elasticity \ZJPsi. 
Figure~\ref{fig:res:xsecsptz} shows differential cross sections as a function of the elasticity \ZJPsi in bins of \PtStarJPsi and as a function of \PtStarJPsiSquare in bins of \ZJPsi.
A comparison of the electroproduction data with predictions from the Monte 
Carlo generator \Cascade reveals in general a reasonable agreement with the data. 
Differences in shape can be seen in the differential cross section as a function of \PtStarJPsiSquare.

For photoproduction, several theory calculations to next-to-leading order have 
been performed and are compared with the data in figure~\ref{fig:res:gp:xsecs:1d:CSMNLO}.
A calculation in the CSM at NLO~\cite{Kraemer}
	was repeated using up-to-date sets of scale parameters~\cite{Maltoni, Kniehl}, 
yielding predictions as shown in figure~\ref{fig:res:gp:xsecs:1d:CSMNLO}a)-b). 
The shapes of the data are reasonably described, whereas the normalisation of the prediction 
is about a factor three below the data, with large uncertainties, 
indicating that corrections beyond next-to-leading order are necessary in order to describe the data.
Estimates of the NNLO contribution for charmonium production at the Tevatron~\cite{maltoni1,maltoni2} indicate that these contributions can be large indeed. 

The calculation to next-to-leading order has been extended to include 
colour octet contributions resulting in a larger cross section~\cite{Kniehl}.
A comparison of this prediction with the data is shown in figure~\ref{fig:res:gp:xsecs:1d:CSMNLO}c)-d).
The dominant uncertainty arises from the difference in the predicted cross section when using LO colour octet LDMEs or higher order improved LDMEs~\cite{Kniehl}.
The NRQCD prediction fails however in describing the shape of the differential cross section as a function of the elasticity \ZJPsi,
even within the presently large uncertainties of the calculation.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Polarisation Measurement}
\label{sec:angles}

The measurement of the \JPsi meson helicity distributions provides an independent method to distinguish between different production mechanisms.
The measurement is performed for the photoproduction data sample.
The \JPsi meson polarisation is measured by analysing the decay angle
distributions of the \JPsi meson, and their dependence on \PtJPsi and \ZJPsi, in
two complementary 
frames~\cite{polarization}: the helicity frame and the Collins-Soper frame.
In the helicity frame the polarisation axis $z$ in the \JPsi meson rest frame 
is defined by the flight direction of the \JPsi meson in the $\gamma p$ 
rest frame, whereas the polarisation in the Collins-Soper frame is 
measured with respect to the bisector of proton ($-\vec{p_p}$) and 
photon ($\vec{p_\gamma}$) in the \JPsi meson rest frame \cite{beneke}.
Subsequently, the frame-dependent polarisation axis is taken as $z$
axis of a right handed coordinate system, where the $x$ and $z$ axis 
lie in a plane spanned by the photon and proton directions.
The $y$ axis is perpendicular to this plane and is the same in both reference frames.
The polar ($\Theta^*$) and azimuthal ($\phi^*$) angles of the positive decay muons are used.

The parametrisation of the measured decay angle
distributions as function of \CosThetaStar and \PhiStar is given by~\cite{beneke}:
\be
\frac{\mr{d}\sigma}{\mr{d}\cos{\Theta^*}}&\propto&1+\alpha\cos^2{\Theta^*}\, ;
\label{eq:theta}\\[.3em]
\frac{\mr{d}\sigma}{\mr{d}\phi^*}&\propto&1+\frac{\alpha}{3}+\frac{\nu}{3}\,\cos{2\phi^*}.
\label{eq:phi}
\ee
The polarisation variables $\alpha$ and $\nu$ can be related to elements of the spin density matrix for the \JPsi meson. Moreover, $\alpha = +1$ and $-1$  corresponds to fully transverse and longitudinal polarisation of the \JPsi meson, respectively. 

A $\chi^2$ fit is performed in each bin of the polarisation measurement, 
comparing data to Monte Carlo samples on reconstruction level probing values for $\alpha$ and $\nu$ between $-1$ and $+1$.
Systematic uncertainties on this measurement are negligible compared to rather large statistical uncertainties.
The results for $\alpha$ and $\nu$ as a function of \PtJPsi and \ZJPsi are presented for the helicity frame  in figure~\ref{fig:gp:Polarisation} and in figure~\ref{fig:gp:Polarisation:CS} for the Collins-Soper frame.
The values for the polarisation parameters in both frames are listed in
table~\ref{tab:pol:heli}.

Within uncertainties the \JPsi mesons produced inelastically at HERA are unpolarised.
The measurements are compared to predictions using 
a \kT factorisation ansatz~\cite{Baranov} and to calculations in the CSM in collinear 
factorisation at leading order~\cite{Baranov} and next-to-leading order~\cite{Maltoni}.
The predictions in the \kT factorisation ansatz describe the data. 
The NLO calculations show a similar
trend within large uncertainties.
In contrast, the leading order CSM calculation predicts larger values for the polarisation variables than the measured ones for many bins and is disfavoured by the measurement.
A similar measurement was published by the ZEUS collaboration in a different kinematic range~\cite{ZEUS_Pol}.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Conclusions}
\label{sec:conclusions}
A measurement of inelastic \JPsi meson production is performed. 
Differential cross sections with improved statistical and systematic uncertainties are presented for both electroproduction
and photoproduction. Polarisation parameters for the photoproduction of \JPsi mesons are measured in
two different reference frames, the helicity frame and the Collins-Soper frame.  

The data are compared to a number of recent theory predictions.
It is found that predictions based on \kT factorisation in the colour singlet model are able to describe the cross sections and the helicity distributions well. 
Calculations based on collinear factorisation in the colour singlet model at next-to-leading order
produce a reasonable description of the shape of the measured cross sections, but are lower in normalisation.
They give an acceptable description of the polarisation parameter measurements
within the large uncertainties.
The failure to describe the cross section measurements and the strong sensitivity to scale variations indicate that calculations beyond next-to-leading order are necessary.
Moreover contributions from colour octet states may be significant. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\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 thank 
Pierre Artoisenet, Sergey Baranov, Mathias Butensch\"on, Bernd Kniehl, Michael Kr\"amer and Fabio Maltoni 
for providing theory calculations for this paper as well as for helpful discussions.

\clearpage



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{thebibliography}{99}
% measurements of charm
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


\bibitem{Aid:1996dn} S.~Aid {\it et al.}  [H1 Collaboration],
%remove for draft
%``Elastic and Inelastic Photoproduction of \JPsi Mesons at HERA''
Nucl.\ Phys.\ B {\bf 472} (1996) 3 [hep-ex/9603005].

\bibitem{Adloff:1999zs} C.~Adloff {\it et al.}  [H1 Collaboration],
%remove for draft
%``Charmonium production in deep inelastic scattering at HERA''
Eur.\ Phys.\ J.\ C {\bf 10} (1999) 373 [hep-ex/9903008].

\bibitem{Adloff:2002ex} C.~Adloff {\it et al.}  [H1 Collaboration],
%remove for draft
%``Inelastic photoproduction of J/psi mesons at HERA''
Eur.\ Phys.\ J.\ C {\bf 25} (2002) 25 [hep-ex/0205064].

\bibitem{Adloff:2002ey} C.~Adloff {\it et al.}  [H1 Collaboration],
%remove for draft
%``Inelastic leptoproduction of J/psi mesons at HERA''
Eur.\ Phys.\ J.\ C {\bf 25} (2002) 41 [hep-ex/0205065].

\bibitem{Breitweg:1997we} J.~Breitweg {\it et al.}  [ZEUS Collaboration],
%remove for draft
%``Measurement of inelastic J/psi photoproduction at HERA''
Z.\ Phys.\ C {\bf 76} (1997) 599 [hep-ex/9708010].

\bibitem{Chekanov:2002at} S.~Chekanov {\it et al.}  [ZEUS Collaboration],
%remove for draft
%``Measurements of inelastic J/psi and psi' photoproduction at HERA''
Eur.\ Phys.\ J.\ C {\bf 27} (2003) 173 [hep-ex/0211011].



\bibitem{ZEUS05}
  S.~Chekanov {\it et al.}  [ZEUS Collaboration],
%  ``Measurement of inelastic $J/\psi$ production in deep inelastic scattering at HERA,''
  Eur.\ Phys.\ J.\  C {\bf 44} (2005) 13
  [hep-ex/0505008].
  %%CITATION = EPHJA,C44,13;%%


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

\bibitem{Kraemer} M.~Kr\"amer, 
%remove for draft
%``QCD Corrections to Inelastic \JPsi Photoproduction,''
Nucl.\ Phys.\  B {\bf 459} (1996) 3 [hep-ph/9508409].
  %%CITATION = NUPHA,B459,3;%%

\bibitem{Maltoni} 
  P.~Artoisenet, J.~M.~Campbell, F.~Maltoni and F.~Tramontano,
%  ``J/psi production at HERA,''
  Phys.\ Rev.\ Lett.\  {\bf 102} (2009) 142001
  [arXiv:0901.4352].


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

\bibitem{Halzen:1977rs} F.~Halzen,
%remove for draft
%``CVC For Gluons And Hadroproduction Of Quark Flavors''
Phys.\ Lett.\ B {\bf 69} (1977) 105.

\bibitem{Eboli:1998xx} O.~J.~P.~Eboli, E.~M.~Gregores and F.~Halzen,
%remove for draft
%``Inelastic photoproduction at HERA: A second charmonium crisis?''
Phys.\ Lett.\ B {\bf 451} (1999) 241 [hep-ph/9802421].

\bibitem{CSM_Chang} 
  C.~H.~Chang,
%  ``Hadronic Production Of J / Psi Associated With A Gluon,''
  Nucl.\ Phys.\  B {\bf 172} (1980) 425.
  %%CITATION = NUPHA,B172,425;%%

\bibitem{Berger:1980ni} E.~L.~Berger and D.~L.~Jones,
%remove for draft
%``Inelastic Photoproduction Of J / Psi And Upsilon By Gluons''
Phys.\ Rev.\ D {\bf 23} (1981) 1521.

\bibitem{Baier:1981zz} R.~Baier and R.~R\"uckl,
%remove for draft
%``On Inelastic Leptoproduction Of Heavy Quarkonium States''
Nucl.\ Phys.\ B {\bf 201} (1982) 1.

\bibitem{Baier:1981uk} R.~Baier and R.~R\"uckl,
%remove for draft
%``Hadronic Production Of J / Psi And Upsilon: Transverse Momentum Distributions''
Phys.\ Lett.\ B {\bf 102} (1981) 364.

\bibitem{Baier:1983va} R.~Baier and R.~R\"uckl,
%remove for draft
%``Hadronic Collisions: A Quarkonium Factory''
Z.\ Phys.\ C {\bf 19} (1983) 251.



\bibitem{Caswell:1985ui} W.~E.~Caswell and G.~P.~Lepage,
%remove for draft
%``Effective Lagrangians For Bound State Problems In QED, QCD, And Other Field Theories''
Phys.\ Lett.\ B {\bf 167} (1986) 437.

\bibitem{Thacker:1990bm} B.~A.~Thacker and G.~P.~Lepage,
%remove for draft
%``Heavy Quark Bound States In Lattice QCD''
Phys.\ Rev.\ D {\bf 43} (1991) 196.

\bibitem{Bodwin:1994jh} G.~T.~Bodwin, E.~Braaten and G.~P.~Lepage,
%remove for draft
%``Rigorous QCD analysis of inclusive annihilation and production of heavy quarkonium''
Phys.\ Rev.\ D {\bf 51} (1995) 1125 [Erratum-ibid.\ D {\bf 55} (1997) 5853] [hep-ph/9407339].

\bibitem{Edin:1997zb} A.~Edin, G.~Ingelman and J.~Rathsman,
%remove for draft
%``Quarkonium production at the Tevatron through soft colour interactions''
Phys.\ Rev.\ D {\bf 56} (1997) 7317 [hep-ph/9705311].

\bibitem{Braaten:1999qk} E.~Braaten, B.~A.~Kniehl and J.~Lee,
%remove for draft
%``Polarization of prompt J/psi at the Tevatron''
Phys.\ Rev.\ D {\bf 62} (2000) 094005 [hep-ph/9911436].


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

  
\bibitem{Kniehl} 
  M.~Butensch\"on and B.~A.~Kniehl,
%  ``Complete next-to-leading-order corrections to J/psi photoproduction in nonrelativistic quantum chromodynamics,''
  arXiv:0909.2798.
  %%CITATION = ARXIV:0909.2798;%%


% CASCADE
\bibitem{Jung:2000hk} H.~Jung and G.P.~Salam,
%remove for draft
%``Hadronic final state predictions from CCFM: The hadron-level Monte  Carlo generator CASCADE,''
Eur.\ Phys.\ J.\ C {\bf 19} (2001) 351 [hep-ph/0012143];\\
H.~Jung,
%remove for draft
%``The CCFM Monte Carlo generator CASCADE,''
Comput.\ Phys.\ Commun.\  {\bf 143} (2002) 100 [hep-ph/0109102].

\bibitem{ccfm} M.~Ciafaloni,
%remove for draft
%``Coherence Effects In Initial Jets At Small Q**2 / S,''
Nucl.\ Phys.\ B {\bf 296} (1988) 49;  \\
S.~Catani, F.~Fiorani and G.~Marchesini,
%remove for draft
%``QCD Coherence In Initial State Radiation,''
Phys.\ Lett.\ B {\bf 234} (1990) 339; 
%\cite{Catani:1989sg}
%S.~Catani, F.~Fiorani and G.~Marchesini,
%``Small X Behavior Of Initial State Radiation In Perturbative QCD,''
{\it idem}, Nucl.\ Phys.\ B {\bf 336} (1990) 18;\\
G.~Marchesini,
%remove for draft
%``QCD coherence in the structure function and associated distributions at small x,''
Nucl.\ Phys.\ B {\bf 445} (1995) 49.


%\bibitem{kTME}


\bibitem{Baranov}
  S.~P.~Baranov,
%  ``Polarization Variables In Inclusive J/Psi Production At Hera In The K(T)-Factorization Approach'', 
private communication, 2009, based on: JETP Lett.\  {\bf 88} (2008) 471.

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

\bibitem{SetA0}
  H.~Jung,
  ``Un-integrated PDFs in CCFM,''
  hep-ph/0411287.
  %%CITATION = HEP-PH/0411287;%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% detectors
\bibitem{Abt:1997xv} I.~Abt {\it et al.}  [H1 Collaboration],
%remove for draft
%``The Tracking, calorimeter and muon detectors of the H1 experiment at HERA ,''
Nucl.\ Instrum.\ Meth.\ A {\bf 386} (1997) 310 and 348.

\bibitem{FTT}
  A.~Baird {\it et al.},
%  ``A fast high resolution track trigger for the H1 experiment,''
  IEEE Trans.\ Nucl.\ Sci.\  {\bf 48} (2001) 1276
  [hep-ex/0104010].
  %%CITATION = IETNA,48,1276;%%

\bibitem{FTT2}
D.~Meer {\it et al.},
%``A Multifunctional Processing Board for the Fast Track Trigger of the H1 Experiment,''
IEEE Trans.\ Nucl.\ Sci.\  {\bf 49} (2002) 357
[hep-ex/0107010]. 

\bibitem{CIP}
  J.~Becker {\it et al.},
%  ``A vertex trigger based on cylindrical multiwire proportional chambers,''
  Nucl.\ Instrum.\ Meth.\  A {\bf 586} (2008) 190
  [physics/0701002].
  %%CITATION = NUIMA,A586,190;%%

\bibitem{CST} 
  D.~Pitzl {\it et al.},
%  ``The H1 silicon vertex detector,''
  Nucl.\ Instrum.\ Meth.\  A {\bf 454} (2000) 334
  [hep-ex/0002044];\\
  B.~List,
  %``The H1 central silicon tracker,''
  Nucl.\ Instrum.\ Meth.\  A {\bf 501} (2001) 49.
  %%CITATION = NUIMA,A501,49;%%


\bibitem{LAr}
  B.~Andrieu {\it et al.}  [H1 Calorimeter Group],
%  ``The H1 liquid argon calorimeter system,''
  Nucl.\ Instrum.\ Meth.\  A {\bf 336} (1993) 460.

\bibitem{Nicholls:1996di} T.~Nicholls {\it et al.}  [H1 SPACAL Group],
%remove for draft
%``Performance of an electromagnetic lead / scintillating fiber calorimeter for the H1 detector,''
Nucl.\ Instrum.\ Meth.\ A {\bf 374} (1996) 149.

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

\bibitem{jb} A. Blondel and F. Jacquet, in: Proc. Study of an ep Facility for Europe (Ed. U. Amaldi), 
DESY 79/48, (1979) 391. 

\bibitem{esigma} 
  U.~Bassler and G.~Bernardi,
%  ``On the kinematic reconstruction of deep inelastic scattering at HERA: The Sigma method,''
  Nucl.\ Instrum.\ Meth.\  A {\bf 361} (1995) 197
  [hep-ex/9412004].

\bibitem{hadrooII}
  S.~Hellwig,
  ``Investigation of the D* - pi(slow) double tagging method in charm analyses. (In German),''
  diploma thesis, University of Hamburg (2004), H1 thesis 341 
  (available at http://www-h1.desy.de/psfiles/theses/).	

\bibitem{peez}
  M.~Peez,
  ``Search for deviations from the standard model in high transverse energy processes at the electron proton collider HERA. (In French)'', 
  PhD thesis, University of Lyon (2003), H1 thesis 317
  (available at http://www-h1.desy.de/psfiles/theses/).	


\bibitem{msteder}
  M.~Steder,
  ``Measurement of inelastic charmonium production at HERA'', 
  PhD thesis, University of Hamburg (2008), H1 thesis 488
  (available at http://www-h1.desy.de/psfiles/theses/).	

\bibitem{msauter}
  M.~Sauter,
  ``Measurement of Beauty Photoproduction at Threshold using Di-Electron Events with the H1 Detector at HERA'', 
  PhD thesis, ETH Z\"urich (2009), H1 thesis 517
  (available at http://www-h1.desy.de/psfiles/theses/).	


\bibitem{DStar_paper}
  F.~D.~Aaron {\it et al.}  [H1 Collaboration],
%  ``Study of Charm Fragmentation into $\mr{D^{*\pm}}$ Mesons in Deep-Inelastic Scattering at HERA,''
  Eur.\ Phys.\ J.\  C {\bf 59} (2009) 589
  [arXiv:0808.1003].

\bibitem{PDG}
  C.~Amsler {\it et al.}  [Particle Data Group],
%  ``Review of particle physics,''
  Phys.\ Lett.\  B {\bf 667} (2008) 1 and 2009 partial update for the 2010 edition.

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

\bibitem{geant} R. Brun {\it et al.}, CERN-DD/EE-84-1 (1987).

\bibitem{diffvm}  
B. ~List and A. ~Mastroberardino,
%remove for draft
%``DIFFVM: A Monte Carlo generator for diffractive processes in e p
%scattering'' 
in: Proc.~of the Workshop on Monte Carlo generators for HERA physics,
ed.~A.~T.~Doyle {\it et al.,} 
% G.~Grindhammer, G.~Ingelman and H.~Jung, 
DESY-PROC-1999-02 (1999) 396.

%\cite{Adloff:1997yv}
\bibitem{psi2s97}
  C.~Adloff {\it et al.}  [H1 Collaboration],
  %remove for draft
%``Photo-production of psi(2S) mesons at HERA,''
  Phys.\ Lett.\  B {\bf 421} (1998) 385
  [hep-ex/9711012].
  %%CITATION = PHLTA,B421,385;%%

\bibitem{psi2s02}
  C.~Adloff {\it et al.}  [H1 Collaboration],
  Phys.\ Lett.\  B {\bf 541} (2002) 251 [hep-ex/0205107].


\bibitem{pythia}
  T.~Sj\"ostrand,
  %  ``The Lund Monte Carlo For Jet Fragmentation And E+ E- Physics: Jetset Version 6.2,''
  Comput.\ Phys.\ Commun.\  {\bf 39} (1986) 347;\\
  T.~Sj\"ostrand and M.~Bengtsson,
%  ``The Lund Monte Carlo For Jet Fragmentation And E+ E- Physics. Jetset Version 6.3: An Update,''
  Comput.\ Phys.\ Commun.\  {\bf 43} (1987) 367;\\
  T.~Sj\"ostrand,
%  ``High-energy physics event generation with PYTHIA 5.7 and JETSET 7.4,''
  Comput.\ Phys.\ Commun.\  {\bf 82} (1994) 74;\\
  T.~Sj\"ostrand {\it et al.},
%  ``High-energy physics event generation with PYTHIA 6.1,''
  Comput.\ Phys.\ Commun.\  {\bf 135} (2001) 238
  [hep-ph/0010017].

\bibitem{ZEUSPsi2S}
  S.~Chekanov {\it et al.}  [ZEUS Collaboration],
  %``Measurements of inelastic J/psi and psi' photoproduction at HERA,''
  Eur.\ Phys.\ J.\  C {\bf 27} (2003) 173
  [hep-ex/0211011].
  %%CITATION = EPHJA,C27,173;%%

\bibitem{btomujj}
  A.~Aktas {\it et al.}  [H1 Collaboration],
%  ``Measurement of beauty production at HERA using events with muons and jets,''
  Eur.\ Phys.\ J.\  C {\bf 41} (2005) 453
  [hep-ex/0502010].
  %%CITATION = EPHJA,C41,453;%%

\bibitem{btojj}
  A.~Aktas {\it et al.}  [H1 Collaboration],
%  ``Measurement of charm and beauty dijet cross sections in photoproduction  at HERA using the H1 vertex detector,''
  Eur.\ Phys.\ J.\  C {\bf 47} (2006) 597
  [hep-ex/0605016].
  %%CITATION = EPHJA,C47,597;%%

\bibitem{CSTresol}
  F.~D.~Aaron {\it et al.}  [H1 Collaboration],
  %``Measurement of the Charm and Beauty Structure Functions using the H1 Vertex
  %Detector at HERA,''
  Eur.\ Phys.\ J.\  C {\bf 65} (2010) 89
  [arXiv:0907.2643].
  %%CITATION = EPHJA,C65,89;%%


%\bibitem{H1:charmbeauty_beamspot}
%  F.D.~Aaron {\it et al.},[H1~Collaboration],
%  %remove for draft
%``Measurement of the Charm and Beauty Structure Functions using the H1 Vertex Detector at HERA,''
%  0907.2643 [hep-ex].
%  %%CITATION = ARXIV:0907.2643;%%


\bibitem{WWA}
  V.~M.~Budnev, I.~F.~Ginzburg, G.~V.~Meledin and V.~G.~Serbo,
%  ``The Two photon particle production mechanism. Physical problems.Applications. Equivalent photon approximation,''
  Phys.\ Rept.\  {\bf 15} (1974) 181.


\bibitem{maltoni1} 
  J.~M.~Campbell, F.~Maltoni and F.~Tramontano,
%  ``QCD corrections to J/psi and Upsilon production at hadron colliders,''
  Phys.\ Rev.\ Lett.\  {\bf 98} (2007) 252002
  [hep-ph/0703113].

\bibitem{maltoni2} 
  P.~Artoisenet, J.~P.~Lansberg and F.~Maltoni,
%  ``Hadroproduction of $J/\psi$ and $\upsilon$ in association with a heavy-quark pair,''
  Phys.\ Lett.\  B {\bf 653} (2007) 60
  [hep-ph/0703129].



\bibitem{polarization} 
  P.~Faccioli, C.~Lourenco, J.~Seixas and H.~K.~W\"ohri,
%  ``J/psi polarization from fixed-target to collider energies'',
  Phys.\ Rev.\ Lett.\  {\bf 102} (2009) 151802
  [arXiv:0902.4462].
  %%CITATION = PRLTA,102,151802;%%

\bibitem{beneke}
M.~Beneke, M.~Kr\"amer and M.~Vanttinen,
%remove for draft
%``Inelastic photoproduction of polarised J/psi,''
Phys.\ Rev.\ D {\bf 57} (1998) 4258 [hep-ph/9709376].

\bibitem{ZEUS_Pol}
  S.~Chekanov {\it et al.}  [ZEUS Collaboration],
  %``Measurement of $J/\psi$ helicity distributions in inelastic photoproduction at HERA,''
  JHEP {\bf 0912} (2009) 007
  [arXiv:0906.1424].
  %%CITATION = JHEPA,0912,007;%%


%
%\bibitem{Kniehl:2001tk} B.~A.~Kniehl and L.~Zwirner,
%%``J/psi inclusive production in e p deep-inelastic scattering at DESY HERA,''
%Nucl.\ Phys.\  B {\bf 621} (2002) 337 [hep-ph/0112199].
%%
%\bibitem{PhDLueders}
%  S.~L\"uders, PhD thesis (2001),
%%  ``A measurement of the beauty production cross-section via $B \to J/\psi X$ at HERA,''
%  ETHZ-IPP-INTERNAL-REPORT-2001-10;%%
%
%\bibitem{Martin:1998sq} A.D.~Martin, R.G.~Roberts, W.J.~Stirling and R.S.~Thorne,
%%``Parton distributions: A New global analysis,''
%Eur. Phys. J. {\bf C4}, 463 (1998) hep-ph/9803445.
%
%\bibitem{Lai:1999wy} H.L.~Lai {\it et al.}
%[CTEQ Collaboration], hep-ph/9903282.
%
%\bibitem{epjpsi} H.~Jung, in: Proc.~of the Workshop on Physics at HERA, 
%  ed.~W.~Buchm\"uller and G.~Ingelman, Hamburg (1992), Vol.~3, p.~1488; \\
%H.~Jung, D.~Kr\"ucker, C.~Greub and D.~Wyler,
%%``Relativistic corrections to photoproduction of J / psi,''
%Z.\ Phys.\ C {\bf 60} (1993) 721.
%
%\bibitem{kniehl}
%B.A.~Kniehl, Private Communication.
%
%\bibitem{Gribov:ri} V.N.~Gribov and L.N.~Lipatov,
%%``Deep Inelastic E P Scattering In Perturbation Theory,''
%Yad.\ Fiz.\  {\bf 15} (1972) 781 [Sov.\ J.\ Nucl.\ Phys.\  {\bf 15} (1972) 438]; 
%%\bibitem{Altarelli:1977zs}
%G.~Altarelli and G.~Parisi,
%%``Asymptotic Freedom In Parton Language,''
%Nucl.\ Phys.\ B {\bf 126} (1977) 298; 
%%\bibitem{Dokshitzer:1977sg}
%Y.L.~Dokshitzer,
%% ``Calculation Of The Structure Functions For Deep Inelastic Scattering And E+
%%E- Annihilation By Perturbation Theory In Quantum Chromodynamics. (In
%%Russian),''
%Sov.\ Phys.\ JETP {\bf 46} (1977) 641 [Zh.\ Eksp.\ Teor.\ Fiz.\  {\bf 73} (1977) 1216].
%
%\bibitem{kruecker}
%K. ~Kr\"uger and S. ~Mohrdieck, in: Proc.~of the Workshop on Monte Carlo
%generators for HERA physics, ed.~A.~T.~Doyle  {\it et al.,} 
%%G.~Grindhammer, G.~Ingelman and H.~Jung, Hamburg DESY - 
%DESY-PROC-1999-02 (1999) 497;
%D. ~Kr{\"u}cker, ``Models for Inelastic \JPsi Production
%at HERA'', PhD thesis (1995) RWTH Aachen (in German, 
%http://www-h1.desy.de/publications/theses\_list.html).

% theoretical calculations
% cross sections


\clearpage
\end{thebibliography}







%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[p]
\centering
\unitlength1cm
\begin{picture}(16,19.5)
\put( 0. ,13.){\includegraphics*[width=7.5cm]{d09-225f2a.eps}}%fin_control_PtTracks.eps}}
\put( 8.5,13.){\includegraphics*[width=7.5cm]{d09-225f2b.eps}}%fin_control_ThetaTracks.eps}} 
\put( 0. ,6.5){\includegraphics*[width=7.5cm]{d09-225f2c.eps}}%fin_control_PtJPsi.eps}} 
\put( 8.5,6.5){\includegraphics*[width=7.5cm]{d09-225f2d.eps}}%fin_control_ThetaJPsi.eps}}
\put( 0. , 0.){\includegraphics*[width=7.5cm]{d09-225f2e.eps}}%fin_control_ZJPsi.eps}}
\put( 8.5, 0.){\includegraphics*[width=7.5cm]{d09-225f2f.eps}}%fin_control_Wgp.eps}}
\put(4.7,19.5){\large \sf Inelastic \JPsi Photoproduction}

\put(1.1,18.9){a)}
\put(9.6,18.9){b)}
\put(1.1,12.4){c)}
\put(9.6,12.4){d)}
\put(1.1,5.9){e)}
\put(9.6,5.9){f)}
\end{picture}
\caption{Control distributions of the photoproduction sample: 
a) the transverse momentum \PtMu of the muon tracks, 
b) the polar angle \ThetaMuon of the muon tracks, 
c) the transverse momentum \PtJPsi of the \JPsi meson, 
d) the polar angle \ThetaJPsi of the \JPsi meson, 
e) the elasticity \ZJPsi and
f) the photon proton centre-of-mass energy \Wgp.
The data are compared with predictions
from the corrected \Cascade Monte Carlo simulation (solid lines), normalised to
the number of entries in the data. The uncorrected \Cascade Monte Carlo prediction 
is shown as dashed line.}
\label{fig:ctrl:muon} 
\end{figure}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[p]
\centering
\unitlength1cm
\begin{picture}(16,19.5)
\put( 0. ,13.){\includegraphics*[width=7.5cm]{d09-225f3a.eps}}%fin_control_PtStar2.eps}}
\put( 8.5,13.){\includegraphics*[width=7.5cm]{d09-225f3b.eps}}%fin_control_ThetaJPsi.eps}} 
\put( 0. ,6.5){\includegraphics*[width=7.5cm]{d09-225f3c.eps}}%fin_control_Q2e.eps}}
\put( 8.5,6.5){\includegraphics*[width=7.5cm]{d09-225f3d.eps}}%fin_control_PtSum_Other.eps}}
\put( 0. , 0.){\includegraphics*[width=7.5cm]{d09-225f3e.eps}}%fin_control_Z.eps}} 
\put( 8.5, 0.){\includegraphics*[width=7.5cm]{d09-225f3f.eps}}%fin_control_Wgp.eps}}
\put(4.55,19.5){\large \sf Inelastic \JPsi Electroproduction}

\put(1.1,18.9){a)}
\put(9.6,18.9){b)}
\put(1.1,12.4){c)}
\put(9.6,12.4){d)}
\put(1.1,5.9){e)}
\put(9.6,5.9){f)}
\end{picture}
\caption{Control distributions of the electroproduction sample:
a) The squared transverse momentum of the \JPsi meson in the photon proton rest 
frame \PtStarJPsiSquare, 
b) the polar angle of the \JPsi meson \ThetaJPsi,
c) the photon virtuality \Qsquared, 
d) the scalar transverse sum $\Sigma P_{T,\rm charged}$,
e) the elasticity \ZJPsi and
f) the photon proton centre of mass energy \Wgp.
The data are compared with predictions from the corrected \Cascade Monte Carlo simulation (solid lines), normalised to the number of entries in the data. The uncorrected \Cascade Monte Carlo prediction 
is shown as dashed line.}
\label{fig:ctrl:dis} 
\end{figure}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[p]
\centering
\unitlength1cm
\begin{picture}(16,13)
\put(3.8 ,6.5){\includegraphics*[width=7.5cm]{d09-225f4a.eps}}%fin_control_Inel_3Trk.eps}}
\put(0. ,0.){\includegraphics*[width=7.5cm]{d09-225f4b.eps}}%fin_lowZ_SIG_Z_3.eps}}
\put(8.5 ,0.){\includegraphics*[width=7.5cm]{d09-225f4c.eps}}%fin_lowZ_BFrac_Summary.eps}}

\put(4.7,12.8){\large \sf Inelastic \JPsi Photoproduction}

\put(4.9,12.35){a)}
\put(1.1,5.85){b)}
\put(9.6,5.85){c)}

\end{picture}

\caption{
a) Distribution of the multiplicity of tracks, $N_{Trk}$, in the central region of the detector 
($20^\circ < \ThetaTrk < 160^\circ$) for the photoproduction sample,
b) signed significance distribution $\mathcal{S}$ for the photoproduction sample at low elasticities 
($0.3 < \ZJPsi < 0.4$) and
c) measured contribution from $b$ hadron decays for three bins of the elasticity $z$
in comparison with the prediction based on \Pythia (scaled up by a factor of two) and \Cascade. }
\label{fig:ctrl:bgandcuts} 
\end{figure}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[p]
\centering
\unitlength1cm
\begin{picture}(16,13)
\put(0. , 6.5){\includegraphics*[width=7.5cm]{d09-225f5a.eps}}%fin_xs_Z_Cascade.eps}} 
\put(8.5, 6.5){\includegraphics*[width=7.5cm]{d09-225f5b.eps}}%fin_xs_Wgp_Cascade.eps}} 
\put(4.25 , 0.){\includegraphics*[width=7.5cm]{d09-225f5c.eps}}%fin_xs_Pt2JPsi_Cascade.eps}} 
		
\put(4.7,13){\large \sf Inelastic \JPsi Photoproduction}

\put(1.4,12.3){a)}
\put(9.9,12.3){b)}
\put(5.6,5.8){c)}
\end{picture}
\caption{Differential \JPsi meson photoproduction cross sections for the kinematic range
$60 < \Wgp < \unit[240]{GeV}$, $0.3 < \ZJPsi < 0.9$ and $\PtJPsi > \unit[1]{GeV}$,
as functions of 
a) the elasticity \ZJPsi, 
b) the photon proton centre of mass energy \Wgp  and 
c) the squared transverse momentum of the \JPsi meson \PtJPsiSquare.
The inner error bar represents the statistical uncertainty and the outer error bar indicates the statistical and systematic uncertainties added in quadrature.
The data are compared to the predictions from \Cascade (solid line). 
The uncertainty band of the \Cascade prediction arises from a scale variation by a factor of two.
The dashed and dotted lines indicate the remaining background from diffractive \PsiPrime or $b$ hadron decays respectively as estimated using MC simulations.}
\label{fig:res:gp:xsecs:1d}
\end{figure}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[p]
\centering
\unitlength1cm
\begin{picture}(16,13)
\put(3.8,6.5){\includegraphics*[width=7.5cm]{d09-225f6a.eps}}%fin_xs_Z_PtBINS_Cascade.eps}}
\put(3.8,0.){\includegraphics*[width=7.5cm]{d09-225f6b.eps}}%fin_xs_Pt2_ZBINS_Cascade.eps}}
\put(5.1,12.3){a)}
\put(5.1,5.8){b)}
\put(4.7,13.0){\large \sf Inelastic \JPsi Photoproduction}
\end{picture}
\caption{a) Differential \JPsi meson cross sections as a function 
of \ZJPsi in four bins of \PtJPsi and
b) differential \JPsi meson cross sections as a function 
of $\PtJPsiSquare$ in four bins of \ZJPsi. 
The inner error bar represents the statistical uncertainty and the outer error bar indicates the statistical and systematic uncertainties added in quadrature.
For visibility, the measured cross sections are scaled by the factors indicated in the figures. 
The data are compared to the predictions from \Cascade (lines).}
\label{fig:res:gp:xsecs:2d} 
\end{figure}




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[p]
\centering
\unitlength1cm
\begin{picture}(16,13)
\put(0. ,6.5){\includegraphics*[width=7.5cm]{d09-225f7a.eps}}%fin_xs_Q2e_Cascade.eps}}
\put(8.5,6.5){\includegraphics*[width=7.5cm]{d09-225f7b.eps}}%fin_xs_PtStar2_Cascade.eps}} 
\put(0., 0.){\includegraphics*[width=7.5cm]{d09-225f7c.eps}}%fin_xs_Wgp_Cascade.eps}} 
\put(8.5,0.){\includegraphics*[width=7.5cm]{d09-225f7d.eps}}%fin_xs_Z_Cascade.eps}} 
\put(4.55,13){\large \sf Inelastic \JPsi Electroproduction}
\put(1.3,12.3){a)}
\put(9.8,12.3){b)}
\put(1.3,5.8){c)}
\put(9.8,5.8){d)}
\end{picture}
\caption{Differential \JPsi meson cross sections for the kinematic range 
$3.6 < \Qsquared < \unit[100]{GeV^2}$, $60 < \Wgp < \unit[240]{GeV}$, $0.3 < \ZJPsi < 0.9$ and 
$\PtStarJPsiSquare > \unit[1]{GeV}$, as functions of 
a) the photon virtuality \Qsquared, 
b) the squared transverse momentum of the \JPsi meson in the photon proton rest 
frame \PtStarJPsiSquare, 
c) the energy in the photon proton rest frame \Wgp and
d) the elasticity \ZJPsi.
The inner error bar represents the statistical uncertainty and the outer error bar indicates the statistical and systematic uncertainties added in quadrature.
The data are compared to the predictions from \Cascade (solid line).
The uncertainty band of the \Cascade prediction arises from a scale variation by a factor of two.
}
\label{fig:dis:xsecs2}
\end{figure}



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[p]
\centering
\unitlength1cm
\begin{picture}(16,13.)
\put(3.8,6.5){\includegraphics*[width=7.5cm]{d09-225f8a.eps}}%fin_xs_Z_PtStarBINS_Cascade.eps}}
\put(3.8,0.){\includegraphics*[width=7.5cm]{d09-225f8b.eps}}%fin_xs_PtStar2_ZBINS_Cascade.eps}}
\put(5.15,12.3){a)}
\put(5.15,5.8){b)}
\put(4.55,13.0){\large \sf Inelastic \JPsi Electroproduction}
\end{picture}
\caption{a) Differential \JPsi meson cross sections as a function of \ZJPsi in three bins of \PtStarJPsi
and b) differential cross section as a function of \PtStarJPsiSquare in three bins of \ZJPsi.
For visibility, the measured cross sections are scaled by the factors indicated in the figures. 
Predictions from \Cascade are shown as solid line.}
\label{fig:res:xsecsptz} 
\end{figure}




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[p]
\centering
\unitlength1cm
\begin{picture}(16,13.)
\put(0.,6.5){\includegraphics*[width=7.5cm]{d09-225f9a.eps}}%fin_xs_Pt2JPsi_CSMNLO.eps}}
\put(8.5,6.5){\includegraphics*[width=7.5cm]{d09-225f9b.eps}}%fin_xs_Z_CSMNLO.eps}}
\put(0.,0.){\includegraphics*[width=7.5cm]{d09-225f9c.eps}}%fin_xs_Pt2JPsi_NRQCD.eps}}
\put(8.5,0.){\includegraphics*[width=7.5cm]{d09-225f9d.eps}}%fin_xs_Z_NRQCD.eps}}
\put(1.2,12.3){a)}
\put(9.8,12.3){b)}
\put(1.2,5.8){c)}
\put(9.8,5.8){d)}
\put(4.7,13.){\large \sf Inelastic \JPsi Photoproduction}
\end{picture}
\caption{Differential \JPsi meson photoproduction cross sections for the kinematic range
$60 < \Wgp < \unit[240]{GeV}$, $0.3 < \ZJPsi < 0.9$ and $\PtJPsi > \unit[1]{GeV}$
as functions of the squared transverse momentum of the \JPsi meson \PtJPsiSquare (a) and c)) and the elasticity \ZJPsi (b) and d)).
The inner error bar represents the statistical uncertainty and the outer error bar indicates the statistical and systematic uncertainties added in quadrature.
The data are compared with calculations to next-to-leading order:
a,b) a colour singlet model (CSM) calculation~\cite{Maltoni} and 
c,d) a NRQCD calculation including contributions from colour octet states (CS + CO)~\cite{Kniehl}. 
The colour singlet component (CS) of the latter calculation is shown separately in addition.}
\label{fig:res:gp:xsecs:1d:CSMNLO}
\end{figure}




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[p]
\centering
\unitlength1cm
\begin{picture}(16,13)
\put(3.8 , 6.5){\includegraphics*[width=7.5cm]{d09-225f10a.eps}}%fin_Polarisation_con_CosTheta_Z_3.eps}}
\put(3.8 , 0.){\includegraphics*[width=7.5cm]{d09-225f10b.eps}}%fin_Polarisation_con_Phi_Pt_1.eps}}
	
\put(5.1,12.35){a)}
\put(5.1,5.85){b)}
\put(5.2,13.5){\large \sf Inelastic \JPsi Polarisation }
\put(6.425,12.9){\large \sf Helicity Frame}

\end{picture}
\caption{Distributions for a) \CosThetaStar in the range $0.6 < \ZJPsi < 0.75$ and b) \PhiStar in the range $1 < \PtJPsi < \unit[2]{GeV}$. The data are compared to the corrected \Cascade simulation with three different assumptions for the polarisation variables $\alpha$ or $\nu$.}
\label{fig:gp:con:Polarisation}
\end{figure}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[p]
\centering
\unitlength1cm
\begin{picture}(16,13)
\put(0. , 6.5){\includegraphics*[width=7.5cm]{d09-225f11a.eps}}%fin_Polarisation_Alpha_Pt_kT_CSMNLO.eps}}
\put(8.5, 6.5){\includegraphics*[width=7.5cm]{d09-225f11b.eps}}%fin_Polarisation_Alpha_Z_kT_CSMNLO.eps}}
\put(0. , 0.){\includegraphics*[width=7.5cm]{d09-225f11c.eps}}%fin_Polarisation_Nu_Pt_kT_CSMNLO.eps}}
\put(8.5, 0.){\includegraphics*[width=7.5cm]{d09-225f11d.eps}}%fin_Polarisation_Nu_Z_kT_CSMNLO.eps}}

\put(5.2,14){\large \sf Inelastic \JPsi Polarisation }
\put(6.425,13.4){\large \sf Helicity Frame}

\put(1.1,12.8){a)}
\put(9.6,12.8){b)}
\put(1.1,6.3){c)}
\put(9.6,6.3){d)}
\end{picture}
\caption{Polarisation parameters $\alpha$ and $\nu$ measured in the helicity frame 
for the kinematic range $60 < \Wgp < \unit[240]{GeV}$, $0.3 < \ZJPsi < 0.9$ and $\PtJPsi > \unit[1]{GeV}$, 
as a function of \ZJPsi and \PtJPsi. The measurement is compared with predictions calculated 
in a \kT factorisation ansatz ~\cite{Baranov} 
and with calculations in CSM (collinear 
factorisation) at leading~\cite{Baranov} and next-to-leading order~\cite{Maltoni}.}
\label{fig:gp:Polarisation}
\end{figure}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[p]
\centering
\unitlength1cm
\begin{picture}(16,13)
\put(0. , 6.5){\includegraphics*[width=7.5cm]{d09-225f12a.eps}}%fin_Polarisation_Alpha_Pt_kT_CSMNLO_CS.eps}}
\put(8.5, 6.5){\includegraphics*[width=7.5cm]{d09-225f12b.eps}}%fin_Polarisation_Alpha_Z_kT_CSMNLO_CS.eps}}
\put(0. , 0.){\includegraphics*[width=7.5cm]{d09-225f12c.eps}}%fin_Polarisation_Nu_Pt_kT_CSMNLO_CS.eps}}
\put(8.5, 0.){\includegraphics*[width=7.5cm]{d09-225f12d.eps}}%fin_Polarisation_Nu_Z_kT_CSMNLO_CS.eps}}

\put(5.2,14){\large \sf Inelastic \JPsi Polarisation }
\put(5.7,13.4){\large \sf Collins-Soper Frame}

\put(1.1,12.8){a)}
\put(9.6,12.8){b)}
\put(1.1,6.3){c)}
\put(9.6,6.3){d)}
\end{picture}
\caption{Polarisation parameters $\alpha$ and $\nu$ in the Collins-Soper frame 
for the kinematic range $60 < \Wgp < \unit[240]{GeV}$, $0.3 < \ZJPsi < 0.9$ and $\PtJPsi > \unit[1]{GeV}$, 
as a function of 
\ZJPsi and \PtJPsi. The measurement is compared with predictions calculated in a \kT factorisation 
ansatz~\cite{Baranov} 
and with calculations in CSM (collinear factorisation) at leading~\cite{Baranov} and 
next-to-leading order~\cite{Maltoni}.}
\label{fig:gp:Polarisation:CS}
\end{figure}



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

\clearpage





%\renewcommand{\arraystretch}{1.1}
\begin{table}[ptb]
\begin{center}
%\begin{tabular}{D{.}{.}{4.0} l D{.}{.}{4.1} D{.}{.}{2.1} l c l c l}
\begin{tabular}{r c r c l p{0.1cm} l p{0.1cm} l}
%\hline 
\multicolumn{9}{ c }{Inelastic \JPsi Photoproduction} \\
%\multicolumn{5}{ c }{ $0.3 < \ZJPsi < 0.9, 60 < \Wgp < \unit[240]{GeV}$} \\
\toprule
\multicolumn{3}{ c }{$\PtJPsiSquare~[\mathrm{GeV}^2]$} & \multicolumn{1}{ c }{$\left< \PtJPsiSquare \right>~[\mathrm{GeV}^2]$} & \multicolumn{5}{ c }{$\mathrm{d}\sigma_\mathrm{\gamma p}/\mathrm{d}\PtJPsiSquare~\mathrm{[nb/GeV^2]}$}  \\
\midrule
 $1.0$ &$\div$& $  2.1$ & $~~1.5$  & $7.75$   & $\pm$ & $ 0.82    $ & $\pm$ & $ 0.70    $ \\
 $2.1$ &$\div$& $  3.5$ & $~~2.7$  & $4.43$   & $\pm$ & $ 0.48    $ & $\pm$ & $ 0.40    $ \\
 $3.5$ &$\div$& $  5.4$ & $~~4.3$  & $2.55$   & $\pm$ & $ 0.28    $ & $\pm$ & $ 0.23    $ \\
 $5.4$ &$\div$& $  7.6$ & $~~6.3$  & $1.06$   & $\pm$ & $ 0.13    $ & $\pm$ & $ 0.10   $ \\
 $7.6$ &$\div$& $ 10.0$ & $~~8.6$  & $0.677$  & $\pm$ & $ 0.084   $ & $\pm$ & $ 0.061   $ \\
$10.0$ &$\div$& $ 13.5$ & $11.4$  & $0.391$   & $\pm$ & $ 0.048   $ & $\pm$ & $ 0.035   $ \\
$13.5$ &$\div$& $ 20.0$ & $15.6$  & $0.156$   & $\pm$ & $ 0.020   $ & $\pm$ & $ 0.014   $ \\
$20.0$ &$\div$& $ 26.5$ & $22.1$  & $0.0509$  & $\pm$ & $ 0.0078  $ & $\pm$ & $ 0.0046  $ \\
$26.5$ &$\div$& $ 40.0$ & $30.0$  & $0.0175$  & $\pm$ & $ 0.0029  $ & $\pm$ & $ 0.0015  $ \\
$40.0$ &$\div$& $ 60.0$ & $46.0$  & $0.0049 $ & $\pm$ & $ 0.0012  $ & $\pm$ & $ 0.0004 $ \\
$60.0$ &$\div$& $100.0$ & $70.0$  & $0.00090$ & $\pm$ & $ 0.00035 $ & $\pm$ & $ 0.00008 $ \\
\specialrule{\heavyrulewidth}{\cmidrulesep}{\cmidrulesep}
\multicolumn{3}{ c }{\ZJPsi} & \multicolumn{1}{ c }{$\left< \ZJPsi \right>$} & \multicolumn{5}{ c }{$ \mathrm{d}\sigma_\mathrm{\gamma p}/\mathrm{d}\ZJPsi~\mathrm{[nb]}$}  \\
\midrule
$0.30$ &$\div$& $0.45$ & $0.375$  & $ 23.4 $ & $\pm$ & $ 2.6 $ & $\pm$ & $ 2.1 $ \\
$0.45$ &$\div$& $0.60$ & $0.525$  & $ 47.6 $ & $\pm$ & $ 4.7 $ & $\pm$ & $ 4.3 $ \\
$0.60$ &$\div$& $0.75$ & $0.675$  & $ 51.3 $ & $\pm$ & $ 5.0 $ & $\pm$ & $ 4.6 $ \\
$0.75$ &$\div$& $0.90$ & $0.825$  & $ 54.2 $ & $\pm$ & $ 5.6 $ & $\pm$ & $ 4.9 $ \\
\bottomrule
\end{tabular}
\end{center}
\caption{Measured differential photoproduction cross sections in the kinematic range $0.3 < \ZJPsi < 0.9, \PtJPsi > \unit[1]{GeV}$ and $60 < \Wgp < \unit[240]{GeV}$ as function of the squared transverse momentum \PtJPsiSquare and the elasticity \ZJPsi of the \JPsi meson. The bin centre values, $\left< \PtJPsiSquare \right>$ and $\left< \ZJPsi \right>$, are also given in the table.}
\label{tab:xsecs:gammap:z}
\label{tab:xsecs:gammap:pt2}
\end{table}

\begin{table}[ptb]
\begin{center}
\begin{tabular}{ r c r c l l p{0.1cm} l p{0.1cm} l }
\multicolumn{10}{ c }{Inelastic \JPsi Photoproduction} \\
\toprule
\multicolumn{3}{ c }{$\Wgp~\mathrm{[GeV]}$} & \multicolumn{1}{ c }{$\left< \Wgp \right>~\mathrm{[GeV]}$} & \multicolumn{1}{ c }{$\Phi_\gamma$} & \multicolumn{5}{ c }{$ \sigma_\mathrm{\gamma p}~[nb]$}  \\
\midrule
 $60$ &$\div$&  $80$ & $~~69$  & $0.0269 $  & $ 22.9 $ & $\pm $ & $ 4.1 $ & $\pm $ & $ 2.1 $ \\
 $80$ &$\div$& $100$ & $~~89$  & $0.0192 $  & $ 24.1 $ & $\pm $ & $ 3.3 $ & $\pm $ & $ 2.2 $ \\
$100$ &$\div$& $120$ & $110$   & $0.0145 $  & $ 24.0 $ & $\pm $ & $ 3.0 $ & $\pm $ & $ 2.2 $ \\
$120$ &$\div$& $140$ & $130$   & $0.0112 $  & $ 30.3 $ & $\pm $ & $ 3.6 $ & $\pm $ & $ 2.7 $ \\
$140$ &$\div$& $160$ & $150$   & $0.00891$  & $ 35.7 $ & $\pm $ & $ 4.3 $ & $\pm $ & $ 3.2 $ \\
$160$ &$\div$& $180$ & $170$   & $0.00716$  & $ 30.4 $ & $\pm $ & $ 3.9 $ & $\pm $ & $ 2.7 $ \\
$180$ &$\div$& $210$ & $194$   & $0.00832$  & $ 31.7 $ & $\pm $ & $ 4.2 $ & $\pm $ & $ 2.9 $ \\
$210$ &$\div$& $240$ & $224$   & $0.00621$  & $ 33.8 $ & $\pm $ & $ 5.6 $ & $\pm $ & $ 3.0 $ \\
\bottomrule
\end{tabular}
\end{center}
\caption{Measured photoproduction cross sections in the kinematic range $\PtJPsi > \unit[1]{GeV}$ and $0.3 < \ZJPsi < 0.9$ in bins of the photon proton centre-of-mass energy \Wgp. The bin centre values $\left< \Wgp \right>$ are also given in the table. $\Phi_\gamma$ denotes the photon flux factors~\cite{WWA} employed in the photoproduction analysis using an upper \Qsquared boundary of $\Qsquared = \unit[2.5]{GeV^2}$. For the range $60 < \Wgp < \unit[240]{GeV}$ a photon flux factor of $\Phi_\gamma = 0.1024$ is calculated.}
\label{tab:xsecs:gammap:wgp}
\label{tab:fluxes}
\end{table}



\begin{table}[ptb]
\begin{center}
\begin{tabular}{r c r c l p{0.1cm} l p{0.1cm} l}
\multicolumn{9}{ c }{Inelastic \JPsi Photoproduction} \\
\toprule
\multicolumn{3}{ c }{$\PtJPsiSquare~[\mathrm{GeV}^2]$} & \multicolumn{1}{ c }{$\left< \PtJPsiSquare \right>~[\mathrm{GeV}^2]$} & \multicolumn{5}{ c }{$ \mathrm{d}\sigma_\mathrm{\gamma p}/\mathrm{d}\PtJPsiSquare~\mathrm{[nb/GeV^2]}$ } \\
\specialrule{\lightrulewidth}{\cmidrulesep}{0pt}
\multicolumn{9}{ c }{ $0.30 < \ZJPsi < 0.45$} \\
\specialrule{\lightrulewidth}{0pt}{\cmidrulesep}
 $1.0$ &$\div$& $2.0 $  & $~~1.4$  & $ 1.02    $ & $\pm$ & $ 0.20    $ & $\pm$ & $ 0.09 $ \\
 $2.0$ &$\div$& $3.0 $  & $~~2.5$  & $ 0.64    $ & $\pm$ & $ 0.13    $ & $\pm$ & $ 0.06 $ \\
 $3.0$ &$\div$& $4.5 $  & $~~3.6$  & $ 0.402   $ & $\pm$ & $ 0.077   $ & $\pm$ & $ 0.036 $ \\
 $4.5$ &$\div$& $7.0 $  & $~~5.5$  & $ 0.180   $ & $\pm$ & $ 0.036   $ & $\pm$ & $ 0.016 $ \\
 $7.0$ &$\div$& $10.0$  & $~~8.2$  & $ 0.093   $ & $\pm$ & $ 0.021   $ & $\pm$ & $ 0.008 $ \\
$10.0$ &$\div$& $14.0$  & $11.6$   & $ 0.047   $ & $\pm$ & $ 0.011   $ & $\pm$ & $ 0.004 $ \\
$14.0$ &$\div$& $20.0$  & $16.2$   & $ 0.0210  $ & $\pm$ & $ 0.0052  $ & $\pm$ & $ 0.0019 $ \\
$20.0$ &$\div$& $40.0$  & $25.0$   & $ 0.0065  $ & $\pm$ & $ 0.0018  $ & $\pm$ & $ 0.0006 $ \\
$40.0$ &$\div$& $100.0$ & $49.0$   & $ 0.00065 $ & $\pm$ & $ 0.00032 $ & $\pm$ & $ 0.00006 $ \\
\specialrule{\lightrulewidth}{\cmidrulesep}{0pt}
\multicolumn{9}{ c }{ $0.45 < \ZJPsi < 0.60$} \\
\specialrule{\lightrulewidth}{0pt}{\cmidrulesep}
$1.0 $ &$\div$& $2.0  $ & $~~1.4$   & $ 2.17    $ & $\pm$ & $ 0.29    $ & $\pm$ & $ 0.19 $ \\
$2.0 $ &$\div$& $3.0  $ & $~~2.5$   & $ 1.21    $ & $\pm$ & $ 0.18    $ & $\pm$ & $ 0.11 $ \\
$3.0 $ &$\div$& $4.5  $ & $~~3.6$   & $ 0.74    $ & $\pm$ & $ 0.11    $ & $\pm$ & $ 0.07 $ \\
$4.5 $ &$\div$& $7.0  $ & $~~5.5$   & $ 0.392   $ & $\pm$ & $ 0.057   $ & $\pm$ & $ 0.035 $ \\
$7.0 $ &$\div$& $10.0 $ & $~~8.2$   & $ 0.219   $ & $\pm$ & $ 0.033   $ & $\pm$ & $ 0.020 $ \\
$10.0$ &$\div$& $14.0 $ & $11.6$    & $ 0.107   $ & $\pm$ & $ 0.014   $ & $\pm$ & $ 0.010 $ \\
$14.0$ &$\div$& $20.0 $ & $16.2$    & $ 0.0497  $ & $\pm$ & $ 0.0084  $ & $\pm$ & $ 0.0045 $ \\
$20.0$ &$\div$& $40.0 $ & $25.0$    & $ 0.0072  $ & $\pm$ & $ 0.0015  $ & $\pm$ & $ 0.0007 $ \\
$40.0$ &$\div$& $100.0$ & $49.0$    & $ 0.00072 $ & $\pm$ & $ 0.00030 $ & $\pm$ & $ 0.00007 $ \\
\specialrule{\lightrulewidth}{\cmidrulesep}{0pt}
\multicolumn{9}{ c }{ $0.60 < \ZJPsi < 0.75$} \\
\specialrule{\lightrulewidth}{0pt}{\cmidrulesep}
$1.0 $ &$\div$& $2.0 $ & $~~1.4$   & $ 2.40   $ & $\pm$ & $ 0.31   $ & $\pm$ & $ 0.22 $ \\
$2.0 $ &$\div$& $3.0 $ & $~~2.5$   & $ 1.79   $ & $\pm$ & $ 0.18   $ & $\pm$ & $ 0.11 $ \\
$3.0 $ &$\div$& $4.5 $ & $~~3.6$   & $ 1.01   $ & $\pm$ & $ 0.13   $ & $\pm$ & $ 0.09 $ \\
$4.5 $ &$\div$& $7.0 $ & $~~5.5$   & $ 0.506  $ & $\pm$ & $ 0.070  $ & $\pm$ & $ 0.046 $ \\
$7.0 $ &$\div$& $10.0$ & $~~8.2$   & $ 0.200  $ & $\pm$ & $ 0.032  $ & $\pm$ & $ 0.018 $ \\
$10.0$ &$\div$& $14.0$ & $11.6$    & $ 0.112  $ & $\pm$ & $ 0.018  $ & $\pm$ & $ 0.010 $ \\
$14.0$ &$\div$& $20.0$ & $16.2$    & $ 0.0413 $ & $\pm$ & $ 0.0076 $ & $\pm$ & $ 0.0037 $ \\
$20.0$ &$\div$& $40.0$ & $25.0$    & $ 0.0068 $ & $\pm$ & $ 0.0014 $ & $\pm$ & $ 0.0006 $ \\
\specialrule{\lightrulewidth}{\cmidrulesep}{0pt}
\multicolumn{9}{ c }{ $0.75 < \ZJPsi < 0.90$} \\
\specialrule{\lightrulewidth}{0pt}{\cmidrulesep}
$1.0 $ &$\div$& $2.0 $ & $~~1.4$  & $ 2.40   $ & $\pm$ & $ 0.36   $ & $\pm$ & $ 0.22 $ \\
$2.0 $ &$\div$& $3.0 $ & $~~2.5$  & $ 1.69   $ & $\pm$ & $ 0.27   $ & $\pm$ & $ 0.15 $ \\
$3.0 $ &$\div$& $4.5 $ & $~~3.6$  & $ 0.86   $ & $\pm$ & $ 0.15   $ & $\pm$ & $ 0.08 $ \\
$4.5 $ &$\div$& $7.0 $ & $~~5.5$  & $ 0.437  $ & $\pm$ & $ 0.076  $ & $\pm$ & $ 0.039 $ \\
$7.0 $ &$\div$& $10.0$ & $~~8.2$  & $ 0.226  $ & $\pm$ & $ 0.042  $ & $\pm$ & $ 0.020 $ \\
$10.0$ &$\div$&$ 14.0$ & $11.6$   & $ 0.099  $ & $\pm$ & $ 0.022  $ & $\pm$ & $ 0.009 $ \\
$14.0$ &$\div$&$ 20.0$ & $16.2$   & $ 0.0428 $ & $\pm$ & $ 0.0098 $ & $\pm$ & $ 0.0039 $ \\
$20.0$ &$\div$&$ 40.0$ & $25.0$   & $ 0.0076 $ & $\pm$ & $ 0.0021 $ & $\pm$ & $ 0.0007 $ \\
\bottomrule
\end{tabular}
\end{center}
\caption{Measured differential photoproduction cross sections in the kinematic
range $0.3 < \ZJPsi < 0.9$ and $60 < \Wgp < \unit[240]{GeV}$ as a function of the squared transverse momentum of the \JPsi meson in bins of the elasticity \ZJPsi. The bin centre values $\left< \PtJPsiSquare \right>$ are also given in the table.} 
\label{tab:xsecs:gammap:pt2_Z}
\end{table}




\begin{table}[ptb]
\begin{center}
\begin{tabular}{r c r p{0.5cm} c p{0.5cm} l p{0.1cm} l p{0.1cm} l}
\multicolumn{11}{ c }{Inelastic \JPsi Photoproduction} \\
\toprule
\multicolumn{3}{ c }{$\ZJPsi$} & & \multicolumn{1}{ c }{$\left< \ZJPsi \right>$} & & \multicolumn{5}{ c }{$ \mathrm{d}\sigma_\mathrm{\gamma p}/\mathrm{d}\ZJPsi~\mathrm{[nb]}$}  \\
\specialrule{\lightrulewidth}{\cmidrulesep}{0pt}
\multicolumn{11}{ c }{ $1.0 < \PtJPsi < \unit[2.0]{GeV}$} \\
\specialrule{\lightrulewidth}{0pt}{\cmidrulesep}
$0.30$ &$\div$& $0.45$ & & $0.375$  & & $ 14.9 $ & $\pm$ & $ 2.1 $ & $\pm$ & $ 1.3 $ \\
$0.45$ &$\div$& $0.60$ & & $0.525$  & & $ 28.3 $ & $\pm$ & $ 3.1 $ & $\pm$ & $ 2.5 $ \\
$0.60$ &$\div$& $0.75$ & & $0.675$  & & $ 31.8 $ & $\pm$ & $ 3.4 $ & $\pm$ & $ 2.9 $ \\
$0.75$ &$\div$& $0.90$ & & $0.825$  & & $ 33.6 $ & $\pm$ & $ 4.0 $ & $\pm$ & $ 3.0 $ \\
\specialrule{\lightrulewidth}{\cmidrulesep}{0pt}
\multicolumn{11}{ c }{ $2.0 < \PtJPsi < \unit[3.0]{GeV}$} \\
\specialrule{\lightrulewidth}{0pt}{\cmidrulesep}
$0.30$ &$\div$& $0.45$ & & $0.375$  & & $~~5.1 $ & $\pm$ & $ 0.8 $ & $\pm$ & $ 0.5 $ \\
$0.45$ &$\div$& $0.60$ & & $0.525$  & & $ 11.6 $ & $\pm$ & $ 1.4 $ & $\pm$ & $ 1.0 $ \\
$0.60$ &$\div$& $0.75$ & & $0.675$  & & $ 14.1 $ & $\pm$ & $ 1.6 $ & $\pm$ & $ 1.3 $ \\
$0.75$ &$\div$& $0.90$ & & $0.825$  & & $ 13.1 $ & $\pm$ & $ 1.8 $ & $\pm$ & $ 1.2 $ \\
\specialrule{\lightrulewidth}{\cmidrulesep}{0pt}
\multicolumn{11}{ c }{ $3.0 < \PtJPsi < \unit[4.5]{GeV}$} \\
\specialrule{\lightrulewidth}{0pt}{\cmidrulesep}
$0.30$ &$\div$& $0.45$ & & $0.375$  & & $~~2.60 $ & $\pm$ & $ 0.42 $ & $\pm$ & $ 0.23 $ \\
$0.45$ &$\div$& $0.60$ & & $0.525$  & & $~~6.05 $ & $\pm$ & $ 0.73 $ & $\pm$ & $ 0.54 $ \\
$0.60$ &$\div$& $0.75$ & & $0.675$  & & $~~5.71 $ & $\pm$ & $ 0.71 $ & $\pm$ & $ 0.51 $ \\
$0.75$ &$\div$& $0.90$ & & $0.825$  & & $~~5.32 $ & $\pm$ & $ 0.80 $ & $\pm$ & $ 0.48 $ \\
\specialrule{\lightrulewidth}{\cmidrulesep}{0pt}
\multicolumn{11}{ c }{ $\PtJPsi > \unit[4.5]{GeV}$} \\
\specialrule{\lightrulewidth}{0pt}{\cmidrulesep}
$0.30$ &$\div$& $0.45$ & & $0.375$  & & $~~1.10 $ & $\pm$ & $ 0.20 $ & $\pm$ & $ 0.1 $ \\
$0.45$ &$\div$& $0.60$ & & $0.525$  & & $~~1.30 $ & $\pm$ & $ 0.20 $ & $\pm$ & $ 0.1 $ \\
$0.60$ &$\div$& $0.75$ & & $0.675$  & & $~~1.11 $ & $\pm$ & $ 0.17 $ & $\pm$ & $ 0.1 $ \\
$0.75$ &$\div$& $0.90$ & & $0.825$  & & $~~1.30 $ & $\pm$ & $ 0.24 $ & $\pm$ & $ 0.1 $ \\
\bottomrule
\end{tabular}
\end{center}
\caption{Measured differential photoproduction cross sections in the kinematic
range $\PtJPsi > \unit[1]{GeV}$ and $60 < \Wgp < \unit[240]{GeV}$ as a function of the elasticity \ZJPsi in bins of the transverse momentum of the \JPsi meson.}
\label{tab:xsecs:gammap:Z_Pt}
\end{table}



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


\begin{table}[ptb]
\begin{center}
\begin{tabular}{r c r p{0.5cm} c p{0.5cm}  l p{0.1cm} l p{0.1cm} l}
\multicolumn{11}{ c }{Inelastic \JPsi Electroproduction} \\
\toprule
\multicolumn{3}{ c }{$\Qsquared~\mathrm{[GeV^2]}$} & \multicolumn{3}{ c }{$\left< \Qsquared \right>~\mathrm{[GeV^2]}$} & \multicolumn{5}{ c }{$ \mathrm{d}\sigma_\mathrm{ep}/\mathrm{d}\Qsquared~\mathrm{[pb/GeV^2]}$}  \\
\midrule
$~~3.6 $&$\div$& $~~6.5$ & & $~~4.9$  & & $  14.98 $ & $\pm$ & $ 1.97  $ & $\pm$ & $ 1.27   $ \\
$~~6.5 $&$\div$& $12.0 $ & & $~~8.6$  & & $ ~~6.33 $ & $\pm$ & $ 0.75  $ & $\pm$ & $ 0.54  $ \\
$12.0$  &$\div$& $20.0 $ & & $15.0$   & & $ ~~2.11 $ & $\pm$ & $ 0.33  $ & $\pm$ & $ 0.18  $ \\
$20.0$  &$\div$& $40.0 $ & & $26.7$   & & $ ~~0.74 $ & $\pm$ & $ 0.12  $ & $\pm$ & $ 0.06 $ \\
$40.0$  &$\div$& $100.0$ & & $53.0$   & & $ ~~0.141$ & $\pm$ & $ 0.029 $ & $\pm$ & $ 0.012 $ \\
\toprule
\multicolumn{3}{ c }{$\PtStarJPsiSquare~\mathrm{[GeV^2]}$} &  \multicolumn{3}{ c }{$\left< \PtStarJPsiSquare \right>~\mathrm{[GeV^2]}$} & \multicolumn{5}{ c }{$ \mathrm{d}\sigma_\mathrm{ep}/\mathrm{d}\PtStarJPsiSquare~\mathrm{[pb/GeV^2]}$}  \\
\midrule
$1.0 $ &$\div$& $  2.2$ & & $~~1.6$   & & $ 15.5   $ & $\pm$ & $ 2.7    $ & $\pm$ & $ 1.3 $ \\
$2.2 $ &$\div$& $  3.7$ & & $~~2.9$   & & $ 11.0   $ & $\pm$ & $ 2.1    $ & $\pm$ & $ 0.9 $ \\
$3.7 $ &$\div$& $  6.4$ & & $~~4.9$   & & $ ~~8.7  $ & $\pm$ & $ 1.4    $ & $\pm$ & $ 0.7 $ \\
$6.4 $ &$\div$& $  9.6$ & & $~~7.8$   & & $ ~~5.90 $ & $\pm$ & $ 0.92   $ & $\pm$ & $ 0.50 $ \\
$9.6 $ &$\div$& $ 13.5$ & & $11.2$  & & $ ~~3.23   $ & $\pm$ & $ 0.53   $ & $\pm$ & $ 0.27 $ \\
$13.5$ &$\div$& $ 20.0$ & & $16.0$  & & $ ~~1.69   $ & $\pm$ & $ 0.27   $ & $\pm$ & $ 0.14 $ \\
$20.0$ &$\div$& $ 40.0$ & & $25.7$  & & $ ~~0.576  $ & $\pm$ & $ 0.083  $ & $\pm$ & $ 0.049 $ \\
$40.0$ &$\div$& $100.0$ & & $51.0$  & & $ ~~0.055  $ & $\pm$ & $ 0.012  $ & $\pm$ & $ 0.005 $ \\
\toprule
\multicolumn{3}{ c }{\ZJPsi} & \multicolumn{3}{ c }{$\left< \ZJPsi \right>$}  & \multicolumn{5}{ c }{$ \mathrm{d}\sigma_\mathrm{ep}/\mathrm{d}\ZJPsi~\mathrm{[pb]}$}  \\
\midrule
$0.30$ &$\div$& $0.45$ & & $0.375$ & & $ 150  $ & $\pm$ & $ 26  $ & $\pm$ & $ 13 $ \\
$0.45$ &$\div$& $0.60$ & & $0.525$ & & $ 158  $ & $\pm$ & $ 22  $ & $\pm$ & $ 14 $ \\
$0.60$ &$\div$& $0.75$ & & $0.675$ & & $ 280  $ & $\pm$ & $ 31  $ & $\pm$ & $ 24 $ \\
$0.75$ &$\div$& $0.90$ & & $0.825$ & & $ 239  $ & $\pm$ & $ 29  $ & $\pm$ & $ 20 $ \\
\toprule
\multicolumn{3}{ c }{$\Wgp~\mathrm{[GeV]}$} & \multicolumn{3}{ c }{$\left< \Wgp \right>~\mathrm{[GeV]}$} & \multicolumn{5}{ c }{$ \mathrm{d}\sigma_\mathrm{ep}/\mathrm{d}\Wgp~\mathrm{[pb/GeV]}$ } \\
\midrule
$~~60$&$\div$& $~~80$ & & $~~69$     & & $ 0.89  $ & $\pm$ & $ 0.16  $ & $\pm$ & $ 0.08 $ \\
$~~80$&$\div$& $100$  & & $~~89$    & & $ 1.03  $ & $\pm$ & $ 0.15  $ & $\pm$ & $ 0.09 $ \\
$100$ &$\div$&$ 120$  & & $110$  & & $ 0.77  $ & $\pm$ & $ 0.12  $ & $\pm$ & $ 0.007 $ \\
$120$ &$\div$&$ 140$  & & $130$  & & $ 0.75  $ & $\pm$ & $ 0.11  $ & $\pm$ & $ 0.06 $ \\
$140$ &$\div$&$ 160$  & & $150$  & & $ 0.71  $ & $\pm$ & $ 0.11  $ & $\pm$ & $ 0.06 $ \\
$160$ &$\div$&$ 180$  & & $170$  & & $ 0.55  $ & $\pm$ & $ 0.10  $ & $\pm$ & $ 0.05 $ \\
$180$ &$\div$&$ 210$  & & $194$  & & $ 0.42  $ & $\pm$ & $ 0.09  $ & $\pm$ & $ 0.04 $ \\
$210$ &$\div$&$ 240$  & & $224$  & & $ 0.30  $ & $\pm$ & $ 0.10  $ & $\pm$ & $ 0.03 $ \\
\bottomrule
\end{tabular}
\end{center}
\caption{Measured differential electroproduction cross sections in the kinematic
range $3.6 < \Qsquared < \unit[100]{GeV^2}, \PtStarJPsi > \unit[1]{GeV}$ and $0.3 < \ZJPsi < 0.9$ as function of the four momentum transfer \Qsquared, the squared transverse momentum of the \JPsi meson in the photon proton rest frame \PtStarJPsiSquare, the elasticity \ZJPsi and the photon proton centre-of-mass energy \Wgp.}
\label{tab:xsecs:dis:Q2e}
\label{tab:xsecs:dis:PtStar2}
\label{tab:xsecs:dis:z}
\label{tab:xsecs:dis:Wgp}
\end{table}


\begin{table}[ptb]
\begin{center}
\begin{tabular}{r c r p{0.5cm} c p{0.5cm} l p{0.1cm} l p{0.1cm} l}
\multicolumn{11}{ c }{Inelastic \JPsi Electroproduction} \\
\toprule
\multicolumn{3}{ c }{$\PtStarJPsiSquare~\mathrm{[GeV^2]}$} & \multicolumn{3}{ c }{$\left< \PtStarJPsiSquare \right>\mathrm{[GeV^2]}$} & \multicolumn{5}{ c }{$ \mathrm{d}\sigma_\mathrm{ep}/\mathrm{d}\PtStarJPsiSquare~\mathrm{[nb/GeV^2]}$}  \\
\midrule
\multicolumn{11}{ c }{ $0.30 < \ZJPsi < 0.60$} \\
\midrule
$ 1.0$ &$\div$& $ 4.0$ & & $ 2.2$  & & $ 5.5   $ & $\pm$ & $ 1.1   $ & $\pm$ & $ 0.5 $ \\
$ 4.0$ &$\div$& $ 9.0$ & & $ 5.6$  & & $ 3.0   $ & $\pm$ & $ 0.6   $ & $\pm$ & $ 0.3 $ \\
$ 9.0$ &$\div$& $20.0$ & & $11.3$  & & $ 0.89  $ & $\pm$ & $ 0.17  $ & $\pm$ & $ 0.08 $ \\
$20.0$ &$\div$& $60.0$ & & $27.0$  & & $ 0.11  $ & $\pm$ & $ 0.02  $ & $\pm$ & $ 0.01 $ \\
\midrule
\multicolumn{11}{ c }{ $0.60 < \ZJPsi < 0.75$} \\
\midrule
$ 1.0$ &$\div$& $ 4.0$ & & $ 2.3$  & & $ 3.7   $ & $\pm$ & $ 0.7   $ & $\pm$ & $ 0.3 $ \\
$ 4.0$ &$\div$& $ 9.0$ & & $ 5.7$  & & $ 2.7   $ & $\pm$ & $ 0.4   $ & $\pm$ & $ 0.2 $ \\
$ 9.0$ &$\div$& $20.0$ & & $11.3$  & & $ 0.92  $ & $\pm$ & $ 0.15  $ & $\pm$ & $ 0.08 $ \\
$20.0$ &$\div$& $60.0$ & & $27.0$  & & $ 0.13  $ & $\pm$ & $ 0.03  $ & $\pm$ & $ 0.01 $ \\
\midrule
\multicolumn{11}{ c }{ $0.75 < \ZJPsi < 0.90$} \\
\midrule
$ 1.0$ &$\div$& $ 4.0$ & & $ 2.3$  & & $ 3.3    $ & $\pm$ & $ 0.7  $ & $\pm$ & $ 0.3 $ \\
$ 4.0$ &$\div$& $ 9.0$ & & $ 5.7$  & & $ 2.6    $ & $\pm$ & $ 0.5  $ & $\pm$ & $ 0.2 $ \\
$ 9.0$ &$\div$& $20.0$ & & $11.5$  & & $ 0.67   $ & $\pm$ & $ 0.13  $ & $\pm$ & $ 0.06 $ \\
$20.0$ &$\div$& $60.0$ & & $27.0$  & & $ 0.026  $ & $\pm$ & $ 0.013  $ & $\pm$ & $ 0.002 $ \\
\toprule
\multicolumn{3}{ c }{$\ZJPsi$} & \multicolumn{3}{ c }{$\left< \ZJPsi \right>$} & \multicolumn{5}{ c }{$\mathrm{d}\sigma_\mathrm{ep}/\mathrm{d}\ZJPsi~\mathrm{[nb]}$}  \\
\midrule
\multicolumn{11}{ c }{ $1.0 < \PtStarJPsi < \unit[2.0]{GeV}$} \\
\midrule
$0.30$ &$\div$& $0.45$ & & $0.375$  & & $ 60.0  $ & $\pm$ & $ 17.0  $ & $\pm$ & $ 5.1 $ \\
$0.45$ &$\div$& $0.60$ & & $0.525$  & & $ 48.0  $ & $\pm$ & $ 11.4  $ & $\pm$ & $ 4.1 $ \\
$0.60$ &$\div$& $0.75$ & & $0.675$  & & $ 74.6  $ & $\pm$ & $ 12.8  $ & $\pm$ & $ 6.3 $ \\
$0.75$ &$\div$& $0.90$ & & $0.825$  & & $ 66.8  $ & $\pm$ & $ 12.9  $ & $\pm$ & $ 5.7 $ \\
\midrule
\multicolumn{11}{ c }{ $2.0 < \PtStarJPsi < \unit[3.5]{GeV}$} \\
\midrule
$0.30$ &$\div$& $0.45$ & & $0.375$  & & $ 62.4  $ & $\pm$ & $ 15.1  $ & $\pm$ & $ 5.3 $ \\
$0.45$ &$\div$& $0.60$ & & $0.525$  & & $ 67.1  $ & $\pm$ & $ 13.1  $ & $\pm$ & $ 5.7 $ \\
$0.60$ &$\div$& $0.75$ & & $0.675$  & & $ 115.3  $ & $\pm$ & $ 16.2  $ & $\pm$ & $ 9.8 $ \\
$0.75$ &$\div$& $0.90$ & & $0.825$  & & $ 105.0  $ & $\pm$ & $ 16.7  $ & $\pm$ & $ 8.9 $ \\
\midrule
\multicolumn{11}{ c }{ $3.5 < \PtStarJPsi < \unit[10.]{GeV}$} \\
\midrule
$0.30$ &$\div$& $0.45$ & & $0.375$  & & $ 28.4  $ & $\pm$ & $ 6.7  $ & $\pm$ & $ 2.4 $ \\
$0.45$ &$\div$& $0.60$ & & $0.525$  & & $ 41.9  $ & $\pm$ & $ 7.5  $ & $\pm$ & $ 3.6 $ \\
$0.60$ &$\div$& $0.75$ & & $0.675$  & & $ 79.6  $ & $\pm$ & $ 10.6  $ & $\pm$ & $ 6.8 $ \\
$0.75$ &$\div$& $0.90$ & & $0.825$  & & $ 58.9  $ & $\pm$ & $ 9.5  $ & $\pm$ & $ 5.0 $ \\
\bottomrule
\end{tabular}
\end{center}
\caption{Measured differential electroproduction cross sections in the kinematic
range $3.6 < \Qsquared < \unit[100]{GeV^2}, \PtStarJPsi > \unit[1]{GeV}$ and $60 < \Wgp < \unit[240]{GeV}$ as a function of the squared transverse momentum in the photon proton rest frame \PtStarJPsiSquare in bins of the elasticity \ZJPsi and the elasticity \ZJPsi in bins of the transverse momentum in the photon proton rest frame \PtStarJPsi.}
\label{tab:xsecs:dis:ptstar2_Z}
\label{tab:xsecs:dis:Z_PtStar}
\end{table}



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

\clearpage



\begin{table}[htbp]
\vspace{3.5cm}
\begin{center}
\begin{tabular}{r c l c l l}
\multicolumn{6}{c}{Inelastic \JPsi Photoproduction} \\
\multicolumn{6}{c}{{Helicity Frame}} \\
\toprule 
\multicolumn{3}{c}{\PtJPsi [GeV]} & $\left< \PtJPsi \right>$ [GeV] & \multicolumn{1}{c}{$\alpha$} & \multicolumn{1}{c}{$\nu$}\\
\midrule
$1.0$ &$\div$& $2.0 $ & $1.45$ & $ +0.54~^{+ 0.27}_{- 0.24}$ & $ +0.25~^{+ 0.20}_{- 0.20}$\\
$2.0$ &$\div$& $3.0 $ & $2.46$ & $ -0.15~^{+ 0.24}_{- 0.21}$ & $ -0.74~^{+ 0.40}_{- 0.16}$\\
$3.0$ &$\div$& $4.5 $ & $3.65$ & $ -0.18~^{+ 0.26}_{- 0.23}$ & $ -0.04~^{+ 0.32}_{- 0.34}$\\
$4.5$ &$\div$& $10.0$ & $6.21$ & $ -0.28~^{+ 0.32}_{- 0.26}$ & $ +0.59~^{+ 0.31}_{- 0.36}$\\
\midrule
\multicolumn{3}{c}{\ZJPsi} & $\left< \ZJPsi \right>$ & \multicolumn{1}{c}{$\alpha$} & \multicolumn{1}{c}{$\nu$} \\
\hline
$0.30$ &$\div$& $0.45$ & $0.375$  & $ -0.65~^{+ 0.24}_{- 0.21}$ & $ -0.28~^{+ 0.34}_{- 0.35}$\\
$0.45$ &$\div$& $0.60$ & $0.525$  & $ +0.35~^{+ 0.25}_{- 0.22}$ & $ +0.40~^{+ 0.23}_{- 0.24}$\\
$0.60$ &$\div$& $0.75$ & $0.675$  & $ -0.18~^{+ 0.23}_{- 0.21}$ & $ +0.01~^{+ 0.24}_{- 0.25}$\\
$0.75$ &$\div$& $0.90$ & $0.825$  & $ +0.71~^{+ 0.19}_{- 0.40}$ & $ -0.10~^{+ 0.31}_{- 0.32}$\\
\toprule
\multicolumn{6}{c}{{Collins-Soper Frame}} \\
\midrule 
\multicolumn{3}{c}{\PtJPsi [GeV]} & $\left< \PtJPsi \right>$ [GeV] & \multicolumn{1}{c}{$\alpha$} & \multicolumn{1}{c}{$\nu$} \\
\midrule
$1.0$ &$\div$& $2.0 $ & $1.45$  & $ +0.25~^{+ 0.18}_{- 0.17}$ & $ +0.41~^{+ 0.15}_{- 0.16}$ \\
$2.0$ &$\div$& $3.0 $ & $2.46$  & $ -0.26~^{+ 0.17}_{- 0.15}$ & $ -0.42~^{+ 0.29}_{- 0.31}$\\
$3.0$ &$\div$& $4.5 $ & $3.65$  & $ -0.02~^{+ 0.23}_{- 0.20}$ & $ -0.31~^{+ 0.30}_{- 0.31}$\\
$4.5$ &$\div$& $10.0$ & $6.21$  & $ +0.19~^{+ 0.39}_{- 0.32}$ & $ +0.09~^{+ 0.33}_{- 0.34}$\\

\midrule
\multicolumn{3}{c}{\ZJPsi} & $\left< \ZJPsi \right>$ & \multicolumn{1}{c}{$\alpha$} & \multicolumn{1}{c}{$\nu$}\\
\midrule
%

$0.30 $&$\div$& $0.45$ & $0.375$  & $ +0.47~^{+ 0.34}_{- 0.28}$ & $ -0.18~^{+ 0.26}_{- 0.26}$\\
$0.45$ &$\div$& $0.60$ & $0.525$  & $ -0.02~^{+ 0.18}_{- 0.16}$ & $ +0.24~^{+ 0.18}_{- 0.19}$\\
$0.60 $&$\div$& $0.75$ & $0.675$  & $ -0.00~^{+ 0.18}_{- 0.16}$ & $ -0.16~^{+ 0.23}_{- 0.23}$\\
$0.75$ &$\div$& $0.90$ & $0.825$  & $ -0.02~^{+ 0.23}_{- 0.19}$ & $ +0.50~^{+ 0.26}_{- 0.28}$\\

\bottomrule
\end{tabular}
\end{center}
\caption{Measured polarisation parameters in the helicity and the Collins-Soper
frame as function of \PtJPsi and \ZJPsi in the kinematic range $\PtJPsi >
\unit[1]{GeV}, 60 < \Wgp < \unit[240]{GeV}$ and $0.3 < \ZJPsi < 0.9$.}
\label{tab:pol:heli}
\label{tab:pol:cs}
\end{table}
\renewcommand{\arraystretch}{1.0}
	

\end{document}

