\section{The H1 Detector}\label{detector}
A detailed description of the H1 apparatus can be found elsewhere~\cite{h1nim}.
The following section briefly describes the components of the detector relevant
to this analysis.

The H1 liquid
argon~(LAr) calorimeter~\cite{larc} was used to measure positrons scattered into
the central and forward (proton direction) parts of the H1
detector and also to determine the hadronic energy flow. The
 calorimeter extends over
the polar angle range $4^\circ < \theta <  154^\circ$ with full azimuthal
coverage, where $\theta$ is defined with respect to the incoming proton
direction. It consists of
an electromagnetic section with lead absorbers and a hadronic section with
steel absorbers. Both sections are highly segmented in the transverse and
longitudinal directions.
Energy resolutions for electrons~\cite{h1po} and charged pions~\cite{h1pi}
of
$\sigma_{E_e}/E_e\approx 0.12/\sqrt{E_e\ [\,\mbox{GeV}]} \oplus 0.01$ and
$\sigma_{E_h}/E_h\approx 0.50/\sqrt{E_h\ [\,\mbox{GeV}]} \oplus 0.02$,
 respectively, have been
established in test beams. The uncertainties in the absolute
electromagnetic and hadronic energy scales are $3$\% and
$4$\%, respectively, for the present data sample.

The backward electromagnetic lead-scintillator calorimeter
(BEMC) was used to
measure the properties of the scattered positron for polar angles in the
range $155^\circ<\theta<176^\circ$. An energy resolution of
$0.10/\sqrt{E_e [\,\mbox{GeV}]}\oplus 0.42/E_e$[GeV]$\oplus0.03$
  has been achieved~\cite{h1bemc}. The absolute
electromagnetic scale has been determined to a precision of
$1$\%~\cite{f2pap}. Since it consists of only one interaction length
of material, the hadronic response of the BEMC is poor and approximately
30\% of incident hadrons
leave no significant energy deposition.
Consequently, a large
scale uncertainty of 20\% exists for hadronic measurements made with this
device. For results presented here based on 1994 data, the BEMC was
used both to measure properties of the scattered positron and
the hadronic energy flow.

The BEMC was replaced in 1995 by the SPACAL  as the
main rear calorimeter in the H1 detector. This contains  electromagnetic
and hadronic sections, achieving energy resolutions of $0.075/\sqrt{E_e\
[\,\mbox{GeV}]}\oplus 0.025$
and $30\%$
 for positrons and hadrons,
respectively~\cite{shad}. The absolute electromagnetic energy scale
is known to a precision of 2\% and the hadronic energy scale to
7\%~\cite{spaca}. For the analysis presented here,
the SPACAL is used to measure hadronic energy flow for events in which the
scattered positron is reconstructed in the LAr.

A series of tracking chambers are in place which provide the measurement
of the polar angle of the scattered positron in this analysis. Forward and
central tracking chambers
cover polar angle ranges of $5^\circ<\theta<25^\circ$ and $25^\circ<\theta<155^\circ$,
respectively, and provide a measurement of the primary vertex position and
the polar
angle of the scattered positron. In the backward region
($155^\circ<\theta<176^\circ$) the scattered positron is measured by the
Backward
Proportional Chamber (BPC) which lies in front of the BEMC.
These detectors allow the polar angle of the scattered positron
to be measured with a precision of $1$ mrad~\cite{h1f297}.

The transverse energy~ flow measurement is extended in the proton direction by a small
calorimeter (PLUG) with copper absorber and silicon pad readout, covering the
region between the beam-pipe and the LAr cryostat ($0.7^\circ < \theta <
3.3^\circ$).
Owing to the large amount of
passive material  between the PLUG and the interaction point, which varies
between $0.8$ and $5$ interaction lengths, less than $40$\% of the energy
measured
in
the PLUG  originates from the primary vertex. Using a full simulation of the
material distribution in and around the H1 detector, energy loss
corrections have been
determined. The precision of these corrections has been studied using
 $ep$ data~\cite{enrico}.    An energy scale uncertainty of 26\% and
an energy resolution of $ 1.5/\sqrt{E_h\ [\,\mbox{GeV}]}$
have been achieved.

\section{Event Selection}\label{select}
In this analysis three independent event samples are considered: two from
the 1994 and one from the 1996 data taking period.
The 1994 low {$Q^2$}  samples consist of events in which a scattered
positron is found in the BEMC and correspond to
 running periods in which
 the $ep$ interaction point was at the nominal position and in which
it was
shifted by about 67 cm in the proton direction, $z$. This shift allowed
lower {$Q^2$}  values, down to 2.5 GeV$^2$, to be accessed.
The 1994 data samples together permit the range
$2.5$~GeV$^2 < Q^2 <  100$
GeV$^2$ to be studied.  The integrated luminosities of the nominal and shifted
vertex data samples  are $2.7$ pb$^{-1}$ and $0.058$ pb$^{-1}$, respectively.

 The high {$Q^2$}  sample consists of events
collected in 1996 and corresponds to an integrated luminosity
of $8.2$ pb$^{-1}$.  For this
sample, the scattered positron is detected in the
LAr calorimeter and  {$Q^2$}  is required to be larger  than $100$ GeV$^2$.

DIS events are selected by demanding a well-reconstructed scattered positron
with an energy larger than $12$ GeV.  The event kinematics are determined using
the scattered positron energy, $E_e'$, and polar angle, $\theta_e$
 ({electron
method}):
  $ Q^{2} = 4 E_{e} E_{e}^{'}\; {\rm cos}^{2}(\theta_{e} /2)$
and $ y_e = 1 - (E_{e}^{'} / E_{e})\; {\rm sin^2}(\theta_{e} /2)$
 where $E_e$ is the incident positron beam energy.
 The
scaling variable Bjorken-$x$ is related to these quantities via the square
of the centre of mass energy $s$: $x = Q^{2} / (ys)$, and the
hadronic invariant mass squared is $W^2=sy-Q^2$.

 In addition, the following selection criteria are
applied to suppress non-$ep$ and photoproduction background and to
maintain
optimal  resolution in the kinematic variables:

\begin{itemize}
\item A reconstructed event vertex must be found within $\pm 30$ cm in
 $z$ of the nominal interaction point.
\item The longitudinal momentum balance must be within   $30 \mbox{ GeV} <
      \sum_i E_i - P_{zi} < 70$ GeV, where the sum runs over all energy
      deposits in the LAr and the backward calorimeter.

\item $W_e^2$, the invariant mass squared of the hadronic final state
determined using the scattered positron
 is required to be larger
than $4\,400\,$GeV$^2$.

\item The invariant mass of the hadronic final state determined using
the hadronic energy deposited in the calorimeters, $W_h$, is
also required to satisfy $W_h^2 > 4\,400\,$GeV$^2$.

\item To reduce background from photoproduction interactions, $y_e < 0.6$
is required.
\end{itemize}

For the low {$Q^2$}  sample the positron candidate is required to be found within
the angular acceptance of the BEMC ($157^\circ<\theta_e<173^\circ$ for the
nominal vertex data and
$164^\circ<\theta_e<176^\circ$  for the shifted vertex data). In addition,
positron identification cuts to suppress photoproduction background are made
using the information obtained from the BEMC cluster radius and the matching of
the cluster position with that of a charged particle measured in the BPC.

For the high {$Q^2$}  sample the scattered positron candidate is required to be
found within the LAr calorimeter ($\theta_e<150^\circ$).  Since the high
{$Q^2$}
analysis is especially sensitive to QED radiation collinear with the $e^+$ beam,
further selection criteria are included:
\begin{itemize}
\item $0.5 < y_h/y_e < 1.3$. The variable $y_h$ is given by the Jacquet
       Blondel method~\cite{jablo}: $y_h=\frac{1}{2}\Sigma/E_e$ where
       $\Sigma=\sum_i(E_i-p_{zi})$ and $i$ refers to all  energy clusters
       detected in the calorimeters except that due to the scattered positron.
\item In addition to the {electron method}, the $\Sigma$ {method
}~\cite{sigma} is
        used to calculate $x$ and {$Q^2$}  and a  requirement is made that both
       methods should prescribe the same kinematic bin for a given event.
       According to the $\Sigma$ {method} {$Q^2$} , $y$ and $x$ are given by
       \\
       $Q^2_{\Sigma}={E{_e'}^2} {\sin}^2 \theta_e/(1-y_{\Sigma})$,
       $y_{\Sigma}=\Sigma /(\Sigma+E_e'(1-\cos \theta_e))$ and
       $x_{\Sigma}=Q_{\Sigma}^2/(s y_{\Sigma})$.
\end{itemize}

After the application of these cuts, the
efficiency for selecting events for every interval in $x$ and $Q^2$ used
in this study is greater than 75\%.

The high {$Q^2$}  data from 1994 have also been studied.
The transverse energy~ flow measurements made with the two high {$Q^2$}  data sets are
consistent within their statistical errors.
Because of the different correction
procedures necessary for the 1994 and the 1996 data samples and the
significantly
larger amount of data available to
 the 1996  high {$Q^2$}  analysis, the two
measurements have not been combined and for the high {$Q^2$}  region
only results based on 1996 data are presented here.

\section{Experimental Method}

The data are corrected bin-by-bin for QED radiation and detector
effects using the event generator {\sc Django}~\cite{DJANGO}
together with a full simulation of the detector response.
To demonstrate the accuracy of the simulation and the understanding of the
response of the H1 calorimeters the distributions of transverse momentum
balance, $p_T^h/p_T^e$,  are shown in
Fig.~\ref{ptbal}
 for the different data samples used in this analysis. Here,
$p_T^h$ is defined as the sum of
the azimuthal four-momentum vector components $x$ and $y$ of each energy
deposition $i$ in the LAr and SPACAL calorimeters according to
$(p_T^{h})^{2}={(\Sigma}_i
 p_{x,i})^2 +{(\Sigma}_i p_{y,i})^2$. Similarly, $p_T^e$
is the transverse momentum of the scattered positron.
 Predictions made with the {\sc Django} model are also
shown and they describe the data well.

\begin{figure}
\unitlength1cm
  \begin{picture}(6,5)  \put(0,-4.25){
     \put(0.2,4){\psfig{file=nvd_bal.eps,width=5cm,clip=}}
     \put(5.5,4){\psfig{file=svd_bal.eps,width=5cm,clip=}}
     \put(10.7,4.15){\psfig{file=highb.eps,width=5.08cm,clip=}}
                                }
  \end{picture}
  \caption{\label{ptbal} The distribution of $p_T^h/p_T^e$
   measured  with the 1994 nominal vertex data-set
 (left), the 1994 shifted vertex data-set (middle)
and the 1996 data-set (right). The
data are compared with predictions from {\sc Django}.}
\end{figure}

The work presented in this paper comprises studies of the
production of transverse energy $E_T^*$ as a function of pseudorapidity $\eta^*$
and the kinematics of the $ep$ scattering  process. Measurements are made of
 the distribution  $1/Nd{\mbox{$E_T^*$}} /d\eta^*$ in which, for a specified kinematic range,
$N$ is the total number of DIS events and $d{\mbox{$E_T^*$}} /d\eta^*$ is the
sum of the transverse energies of each particle $i$ per unit of
pseudorapidity. The transverse energy,  $E_{Ti}^*$, of a particle $i$
 with energy $E_i^*$ and polar angle $\theta_i^*$ is defined as
$E_{Ti}^*=E_i^*\sin\theta_i^*$. Results are given in the
hadronic centre of  mass
 system (hCMS) for which the incoming photon direction defines the $+z^*$
direction.\footnote{All quantities presented in the hCMS frame are denoted by the superscript *}

Whereas in the laboratory frame (LAB) a
sizeable contribution to the transverse energy is produced by the
kinematic recoil from the scattered lepton, in the hCMS transverse
energy
is due largely to perturbative QCD and fragmentation effects. The Lorentz
transformation from the LAB into the hCMS therefore  serves to isolate the
physically interesting part of the $E_T$ flow~\hspace{-.85ex}.

At high {$Q^2$}  the Lorentz transformation from the LAB to the hCMS requires great
care. Any deviation from the true transformation produces an artificial
contribution to the transverse energy~ seen in the hCMS. This is a more severe problem at high
{$Q^2$}  than at low {$Q^2$}  since the mean particle transverse energy
increases with $Q$ in the LAB but depends
only weakly on $Q$ in  the
hCMS. For the high {$Q^2$}  data sample the maximum $\langle E_T\rangle$
per unit pseudorapidity is
typically an order of magnitude larger in the LAB than in the
hCMS (about 20~GeV compared to 2~GeV). In practice the transformation is
calculated from the energy and direction of the scattered positron. The
main sources of error are the energy resolution and
calibration of the calorimeters and QED radiation.  Although the energy
resolution and calibration of the H1 LAr
calorimeter are well understood and the influence of QED radiation  is well
described by {\sc Django}, these remain an important source of systematic
bias.
To minimise their influence on the measurement of the transverse energy
for the high {$Q^2$}  data,  the method described below  was used.
This exploits the precise information which is available on the direction
of the scattered
lepton.

A deviation in the energy measurement of the scattered positron  results
in
an artificial component in the momentum of any particle in the hCMS which
lies within the scattering plane
defined by the incoming proton and the scattered positron. To suppress
this artificial momentum, the variable used to measure transverse energy~
in the high {$Q^2$}  data set is redefined as
$E_{\perp}^*=E_T^* |\sin \varphi^*|\frac{\pi}{2}$, where $\varphi^*$
is the
 azimuthal angle w.r.t.\ the lepton scattering plane.
Integrated over
$\varphi^*$ this is again equal to $E_T^*$, assuming an isotropic
distribution of the true transverse energy~ around the proton in the hCMS in $\varphi^*$.
Therefore the
measurements of
\mean{E_T^*} and \mean{E_{\perp}^*} are equivalent for an isotropic $\varphi^*$
distribution.
Although $\varphi^*$
asymmetries are predicted in pQCD from processes such as boson-gluon fusion
and QCD-Compton scattering they  introduce biases of less than 1\%
in the transverse energy flow at high $Q^2$. This is
a significantly
smaller effect than the biases due to poor {\mbox{$E_T^*$}}  resolution which would be
introduced were the redefinition not used.

\section{Systematic Errors}

There are several sources of systematic effects which may
affect the measurements presented in this paper. The systematic
errors owing to these have been investigated~\cite{enrico,fabian} and are
discussed below.

\begin{itemize}
   \item The hadronic energy scales of the calorimeters are known to an
         accuracy of $4$\% (LAr), $26$\% (PLUG), $7$\% (SPACAL) and $20$\%
         (BEMC). This directly gives the uncertainties for the energy
         measurements presented here.
   \item For the scattered positron the energy calibration uncertainty is
        $1$\% in the BEMC region (low {$Q^2$}  sample) and  $3$\% for the
         LAr calorimeter (high {$Q^2$}  sample). Resulting errors on
         measured {{\mbox{$E_T^*$}} } are typically
         $2$\% at low {$Q^2$}   and $6$\% at high {$Q^2$} . At high {$Q^2$}  and at large
         values of  $\eta^*$ this error can become large ($\le 27\%$).
         Section \ref{tef} describes in more detail the specific
         problems related to the high {$Q^2$}
         measurements.
   \item Uncertainties owing to the reconstruction of the kinematic
         variables $x$ and
         {$Q^2$}  are typically of the order of $2$\% and $1.5\%$ for the low
         {$Q^2$}   and high {$Q^2$}  spectra, respectively. The improved precision at high {$Q^2$}
         is due to the additional
         cuts applied to this data sample (see section \ref{select}).
   \item The model dependence of the bin-by-bin correction is estimated
         using the event generators
         {\sc Ariadne}~\cite{ariadne},
         {\sc Herwig}~\cite{herw} and {\sc Lepto}~\cite{lepto}  with two
parton density  functions (GRV~\cite{GRV} and MRSH~\cite{MRSH}).
         The differences between the models are typically $3\%$ at
high {$Q^2$}   and
$4.5$\% at low
         $Q^2$. Larger differences are found for the two $\eta^*$ bins measured
         with the PLUG calorimeter ($20$\%). The larger
 errors in the  PLUG pseudorapidity region are due to the uncertainties
in the modelling of the physics in this region.
   \item The transverse energy~ measurement is strongly reliant on a correct simulation
of     the inactive material in and around the H1 detector, particularly
         in the forward region.
          The sensitivity of the measurements to the
         assumed material distribution is estimated by varying the
         amount and location of the inactive material in the simulation.
         The transverse energy flows measured with the 1994 and 1996
         data, for which different configurations of dead material
          were present, are also compared.
         This gives rise to  typical  uncertainties
         of the order of $4\%$ for the
data samples. These differences are
included as systematic errors. For the
         points measured with the PLUG
         calorimeter a different procedure is used. Using the shifted
vertex  and the nominal vertex data sets energy flow is measured
in approximately the same detector volume in the PLUG although
incident particles pass through different amounts of inactive material.
 After applying the dead material
correction,  the change in total
energy measured in the PLUG  is
compatible with Monte Carlo  expectations to within $5$\%. The magnitude
of
the
total energy shift observed in data
between the nominal and shifted vertex event samples, $11$\%, is taken as
 the systematic error on this measurement due to the influence of
 dead material.
   \item Another possible source of error is the simulation of the hadronic
         shower within the LAr and the PLUG calorimeters. This has been
         investigated by comparing the simulation programs {\sc
Gheisha}~\cite{geisha}  and {\sc Calor}~\cite{calor}. The differences are
         small for the LAr calorimeter ($3\%$) but large for the PLUG
         calorimeter  ($20$\%).
   \item Photoproduction background is only important at low $Q^2$.
         The typical uncertainty on the measured points owing to this
source
         is  $1\%$ and has been estimated with the
         {\sc Phojet}~\cite{phojet} program.
\end{itemize}

\input{qcd_final}

\section{Results}

\subsection{Transverse Energy}\label{tef}

In addition to the selection criteria outlined in section \ref{select},  a
further cut on the total energy produced in the angular region  $4.4^\circ\! <
\theta <15^\circ$ is also imposed. The energy in this angular region is
required to be larger than $0.5$ GeV  ($E_{fwd}$ cut). This condition was
introduced for the low {$Q^2$}  analysis  to suppress diffractive-like
events, which are characterised
by an absence of hadronic energy in the
forward region~\cite{h1r,zr}. It
allows a
comparison with  the predictions of QCD based models which is less
dependent on uncertainties
introduced  by the attempt to model  diffractive processes.

In Fig.~\ref{lpl} the $E_T$ flow~ in the hCMS  is presented in 17 regions
of low $Q^2$ and $x$.
The values  of average {$Q^2$}
 cover a range from $3.2$ GeV$^2$ to $70$ GeV$^2$
and the mean values of $x$ extend
 from  $0.08\cdot 10^{-3}$ to $ 7\cdot 10^{-3}$.
Here, and in all following figures
in which
two sets of error bars are displayed,
the total errors (outer error bars) are the result of
adding in quadrature
the statistical errors (inner error bars) and the
systematic errors.
\begin{figure}[p]\unitlength1cm
  \begin{picture}(16,18)
     \put(0.2,0){\epsfig{file=fig1.eps,width=14.7cm,clip=}}
  \end{picture}
  \caption{\label{lpl}
  The inclusive transverse energy flow $1/Nd{\mbox{$E_T^*$}} /d\eta^*$
  at different values of $x$ and {$Q^2$}  for the low
   {$Q^2$}  sample. Note that the errors on all of the measurements made at the
two lowest values of $\eta^*$ in each $x$ and {$Q^2$}  interval
are highly correlated and largely independent of the errors
at larger values of $\eta^*$.
The data are compared to four
  QCD based models.
 The arrows
   mark the average position of the origin of the Breit frame
($\frac{1}{2}\ln{\frac{1}{{\langle}x{\rangle}}-1} $).}
\end{figure}

The data exhibit a mean transverse energy of approximately $2$ GeV per unit of
pseudorapidity. The $E_T$ flow~ shows a plateau for {$Q^2$}  values below about
$10$ GeV$^2$ in the current region ($\eta^* > 0$), although as
{$Q^2$}  increases the distributions become peaked in this region. Using the
PLUG calorimeter it is possible  to measure transverse energy~ in the vicinity of
the proton remnant; these data can be seen as the two  points at smallest
$\eta^*$.
Here
the transverse energy~ flow tends to be about one half of
 that which is measured in the current region.

The systematic errors on the {\mbox{$E_T^*$}}  measurements shown in Fig.~\ref{lpl} are
correlated.
The measurements in the range $-1<\eta^*<-4$ were made predominantly using
the LAr calorimeter and  suffer from a
10\% normalisation uncertainty owing to the LAr hadronic energy
 scale and the model dependence of the correction procedure. These
 are uncorrelated with the errors on the measurements made using  the
PLUG and the BEMC calorimeters. For each $x$ and $Q^2$ interval
the PLUG was used exclusively for
  the data points at the two
lowest values of $\eta^*$ and the BEMC was the main calorimeter used for
measurements in three highest $\eta^*$ bins.  The errors on the measurements
provided by both of these calorimeters are  predominantly normalisation uncertainties.

In Fig.~\ref{lpl} four QCD
 models,
{\sc Ariadne}, {\sc Herwig}, {\sc Lepto} and {\sc Rapgap}
are compared to the data.
These models give an acceptable overall
description of the measured  transverse energy flow with the exceptions
that  {\sc Herwig} fails to describe the
shape of the distribution for values of $Q^2$ above about 7 GeV$^2$ and
shows a peak
at $\eta^* \approx 1$ and {\sc Lepto} fails to match the data
in the central pseudorapidity region for the lowest values of
$x$ and $Q^2$.
{\sc Ariadne} and {\sc Rapgap} give the best description of
all of the models in the lowest
$x$ and {$Q^2$}  bins although each tends to overestimate the transverse
energy flow in the vicinity of the proton remnant.

\begin{figure}[p]\unitlength1cm
   \begin{picture}(16,20)
      \put(-0.1,0){\epsfig{file=fig3.eps,width=16cm,clip=}}
   \end{picture}
   \caption{\label{gpl}
   The inclusive transverse energy flow $1/Nd{\mbox{$E_T^*$}} /d\eta^*$
   at different values of $x$ and {$Q^2$}  for the high
   {$Q^2$}  sample. The data are compared to four QCD based models. The arrows
   mark the average position of the origin of the Breit frame
($\frac{1}{2}\ln{\frac{1}{{\langle}x{\rangle}}-1} $).
      }
\end{figure}

\begin{figure}[!ht]\unitlength1cm
   \begin{picture}(16,7.5)
      \put(-0.2,-.1){\epsfig{file=gpl_one.eps,width=15cm,clip=}}
   \end{picture}
   \caption{\label{gone}
The inclusive transverse energy flow $1/Nd{\mbox{$E_T^*$}} /d\eta^*$
for two selected kinematic bins from Figs. \ref{lpl}
and ~\ref{gpl}.
The influence of GAL string reinteractions and parton showers
on the  expected  transverse energy flow are shown.
      }
\end{figure}

The measurements at high {$Q^2$}  are shown in Fig.~\ref{gpl}.
The $E_T$ flow~ is presented in 7 regions
in $x$ and {$Q^2$}  and compared again to the four QCD models.
At high {$Q^2$}  the experimental resolution
in $\eta^*$ becomes strongly dependent on $x$ and {$Q^2$}  and the $\eta^*$
bin widths are adjusted accordingly.
The  average {$Q^2$}  values cover a range from $175$ GeV$^2$
to $2\,200$
GeV$^2$ and average $x$ from 0.0043 to 0.11. For
consistency with the low {$Q^2$}  results the $E_{fwd}$ cut is applied,
although its
influence is small at high $Q^2$.
  As is the case for the low {$Q^2$}  spectra,  the transverse energy~ flow is seen to  peak
in the
current region. Here, the measured  transverse energy~  is significantly higher  than
at low
{$Q^2$} .
In this kinematic range all QCD models describe the data well
 with the
exception of {\sc Herwig} which
produces insufficient transverse energy~ over the {$Q^2$}  range shown
and does not
describe the shape of the energy flows.

At high {$Q^2$}  the positions of the maxima of the $E_T$ flow~ coincide
approximately
with the average origin of the Breit frame at $\eta^*_{O-BF}=
\frac{1}{2}\ln{\frac{1}{{\langle}x{\rangle}}-1}$. The
Breit frame (BF) is defined as the
frame in which
the exchanged boson carries the
momentum $(0,0,0,-Q)$.
In \cite{OBF} it was argued that at high {$Q^2$}  the maximum of the radiated
transverse energy should coincide with $\eta^*_{O-BF}$.
This behaviour can be seen in our data. In
Fig.~\ref{lpl} and Fig.~\ref{gpl} the position of $\eta^*_{O-BF}$
is marked by an arrow.

In Fig.~\ref{gone} the transverse energy flow is shown for
two {$Q^2$}  bins (\mean{Q^2}$=3.2$ GeV$^2$ and \mean{Q^2}$=253$ GeV$^2$)
compared to the predictions of {\sc Lepto} with different parameter
settings. The effects of using  GAL based string reinteractions
are shown as
the dashed
line.
The solid line represents {\sc Lepto} predictions when the GAL
model is not used. The
GAL approach
leads to a deficit of transverse energy production at low {$Q^2$}
although its effect is small at high $Q^2$. The expectations of
{\sc Lepto} without parton showers (PS) are shown as the dotted curves.
The predictions lie below the data for most of the pseudorapidity
range measured, illustrating the sensitivity of our measurements to
 pQCD processes, as  modelled here using a parton shower algorithm.

\subsection{Dependence on  Bjorken-{\boldmath $x$} }

It has so far proved impossible to distinguish between
the different QCD evolution schemes mentioned in section \ref{QCDpre}
using inclusive
 structure function measurements. It  has
been suggested that the hadronic final state  may prove to be
a more sensitive testing  ground~\cite{mueller}. Measurements of energy
flow~\cite{mk2}, jets and leading particles~\cite{dave} in
the region near
the  proton remnant have indicated that DGLAP  parton evolution fails to
produce sufficient QCD radiation at low $x$.  On the other hand,
BFKL evolution predicts a  rise of hard parton emissions at
low $x$ in the central pseudorapidity region of the hCMS.

\begin{figure}[t]\unitlength1cm
  \begin{picture}(16.,20.5)
    \put(0,0){\epsfig{file=fig7.eps,width=14.25cm,clip=}}
  \end{picture}
  \caption{\label{etx}
   Variation of mean $E_T^*$ in the central pseudorapidity region
($-0.5<\eta^*<0.5$) with $x$ in different regions of {$Q^2$}
compared to four QCD based models and an analytical BFKL calculation.}
\end{figure}

To study these QCD evolution effects,
the dependence of the mean transverse energy~ in the central
pseudorapidity region ($-0.5< \eta^* < 0.5$) with $x$ is shown in
Fig.~\ref{etx}
for different ranges of $Q^2$. The  central region has been chosen since
it
is less affected by the hard scattering
and
still lies  within the acceptance of the H1 detector. A  rise
in the measured transverse energy~ is observed as $x$ decreases.  The general behaviour
of the data can be understood as a rise of the average transverse energy
with hadronic mass, $W$, which increases as $x$ decreases at fixed $Q^2$.

The QCD based Monte Carlo models exhibit reasonable agreement with the
data over most of the kinematic range presented here.
However, the shape of the $x$ dependence as predicted by
the DGLAP based model {\sc
Lepto} does not follow closely the data at the smallest
{$Q^2$}  and $x$ values shown in Fig.~\ref{etx}. The calculations of the
CCFM based
LDCMC show the same behaviour in this region.
Analytical BFKL calculations at the parton
level~\cite{etbfkl} describe the $x$
dependence better and are closer to the data at the lowest
values of {$Q^2$}  and $x$.
However, uncertainties due to hadronisation corrections and
missing higher orders in the calculations preclude
an interpretation of the data as a
signal of the onset of BFKL dynamics.

The effects of a resolved component of the virtual photon provide another possible explanation
of the observed increase of transverse energy~ production with decreasing $x$.
Calculations  using {\sc Rapgap} including  resolved and direct virtual
photon
components  describe well the $x$ dependence and the amount of
transverse energy~ flow. The CDM approach of {\sc Ariadne} is also
in good agreement with the data.

\subsection{Dependence on  {\boldmath $Q^2$} and {\boldmath
$W$}}\label{ww}

\begin{figure}[t]\unitlength1cm
   \begin{picture}(16,9.9)
      \put(0,0){\epsfig{file=fig4.eps,width=16.cm,clip=}}
   \end{picture}
   \caption{\label{wdep}
   Variation of the average $E_T^*$ in the central pseudorapidity region
($-0.5<\eta^*<0.5$)
   with $W$. The present high {$Q^2$}  data are shown together with earlier H1
   results at lower $Q^2$. Also shown are
    hadron-hadron results
   ($p\bar{p}$ and $pp$ from UA1, NA22 and AFS).
   }
\end{figure}

\begin{figure}[ht]\unitlength1cm
   \begin{picture}(15.5,17.)
      \put(-0.15,0){\epsfig{file=fig5.eps,width=15.5cm,clip=}}
   \end{picture}
   \caption{\label{wpl}
   The inclusive transverse energy flow
   $1/Nd{\mbox{$E_T^*$}} /d\eta^*$ in different regions of {$Q^2$}
    for \mean{W} $\approx$ 185 GeV.
    The data are
   compared to three QCD based models.  The two pairs of dashed
vertical lines
denote the central ($-0.5<\eta^*<0.5$) and photon fragmentation
($2<\eta^*<3$) bins.  }
\end{figure}

As mentioned in the introduction, photon-proton scattering
is similar to a hadron-hadron process, as it can be described as
the  fluctuation of a photon into a
 hadronic system which then scatters off the proton. In a previous
publication
on transverse energy flow, H1 showed that this  picture is not only valid
in  photoproduction
but also in DIS~\cite{andre1}.
 A similar level of \mean{E_T^*} in the central pseudorapidity region and
a $W$ dependence similar to that observed in
 hadron-hadron interactions were found.  This is consistent
 with the Bjorken-Kogut picture in which hadronic
final state quantities in the central pseudorapidity region are expected
to be insensitive to the nature of the colliding particles~\cite{bjoko}
and
depend only on the centre of mass energy.

\begin{figure}[t]\unitlength1cm
   \begin{picture}(16.,12.25)
      \put(0,0){\epsfig{file=fig6.eps,width=16.cm,clip=}}
   \end{picture}
   \caption{\label{www}
   Variation of the mean $E_T^*$ in the central pseudorapidity bin
($-0.5<\eta^*<0.5$)
   and the photon fragmentation bin ($2<\eta^*<3$) with $Q^2$. The data
   are compared to three QCD based models.
}
\end{figure}

Fig.~\ref{wdep} shows the $W$ dependence of the new H1 high {$Q^2$}
data
in the central pseudorapidity region
($|\eta^*|<0.5$) together with our previously published low {$Q^2$}  results.
Data from $p\bar{p}$ and $pp$ experiments~\cite{ua1,hd}, taken with
non-single
diffractive triggers,  are also shown. This trigger requirement
excludes elastic and quasi-elastic scattering events.

A $W$ dependence compatible with that of the hadron-hadron data
is found in the H1 data although the mean transverse energy~ is larger at high
$Q^2$. For a
given value of $Q^2$, the
increase in \mean{E_T^*} seen in the H1 data  is consistent with a linear dependence on $\ln W$.

In order to investigate further the {$Q^2$}  dependence,  data
 have been selected
from a small range in $W$ ($165$ GeV $<W<$ $213$ GeV).
 To ensure consistency with the previous study~\cite{andre1},
the $E_{fwd}$ cut is removed
for this part of the analysis and  data points and model predictions
 in all subsequent plots are shown with this condition dropped. This has
an effect of, at most, 10\% for
the lowest {$Q^2$}  points and is negligible for the high {$Q^2$}  measurements.
These  fixed-$W$ transverse energy~ flow measurements are shown in
Fig.~\ref{wpl}. From these distributions, it can be seen that the maximum
of the $E_T$ flow~ not only increases with {$Q^2$}  but also that its position
moves towards the central pseudorapidity region.

To study the ${\mbox{$E_T^*$}} $ flow from the fragmenting photon
and in the central pseudorapidity region, two slices in $\eta^*$
are studied, namely the
 central bin ($-0.5<\eta^*<0.5$)  and  the so-called
photon fragmentation  bin ($2<\eta^*<3$). These are delineated in Fig.~\ref{wpl}.
The definition of the  photon fragmentation bin is somewhat arbitrary.
The chosen
bin in $\eta^*$ is expected to be dominated by the fragmentation of the
hadronic fluctuation of the photon.  Since
the rapidity range within which particles of mass $m$ can be produced is
limited to $\pm \ln W/m$, and the width of the photon fragmentation
peak is $\sim \ln (Q^2/m^2)$~\cite{bjoko},  the
peak position of the transverse energy~
flow associated with
the fragmenting photon moves towards the central pseudorapidity bin
 as {$Q^2$}  increases for fixed $W$. This can be seen in Fig.~\ref{wpl}.

Fig.~\ref{www} shows that the measured transverse energy~ in the photon fragmentation bin
rises with $Q^2$. Qualitatively,
this is consistent with the Generalised Vector Dominance
Model~\cite{gvdm,sak,gor}, in
which mesons with higher masses  can
contribute to the scattering at higher $Q^2$. A more
formal approach~\cite{wave} has shown that the photon wave function allows
for $q\bar{q}$ states with higher $p_T$ at higher $Q^2$.

Fig.~\ref{www} also provides the first experimental evidence for a
rise in the transverse energy with $Q^2$ in the central pseudorapidity
bin. This result was not observed in  previous work~\cite{andre1},
which used less than one tenth of the data available for this analysis
and which showed no significant $Q^2$ dependence in this region.

  Comparisons are made with {\sc Ariadne}, {\sc Lepto} and {\sc Rapgap}.
As they are used here, {\sc Lepto} and {\sc Rapgap}
do not include diffractive events. Therefore, the predictions of
these models are expected to lie above
the data  in
Fig.~\ref{www}.
However,
as the
fraction of diffractive  events
is
less
than $10$\%,  with an average transverse energy production of about one
half
of that observed in non-diffractive interactions~\cite{difET},
 the bias is not larger than $10$\%. A comparison
of {\sc Rapgap} predictions including pomeron exchange processes
with those {\sc Rapgap} calculations shown here also supports this.
All of the QCD based models describe the rise with {$Q^2$}
of the measured average transverse energy in the central and
photon fragmentation regions although {\sc Rapgap}
predicts too much  transverse energy in the central
pseudorapidity range.

\section*{Summary}

Measurements of energy production transverse to the photon and
proton directions in the
hadronic centre of mass
system have been presented using deep-inelastic scattering data
taken in positron-proton collisions by the H1 Collaboration at HERA. The
measurements cover  more than 8 units of pseudorapidity and extend over
the kinematic range $3.2 < Q^2< 2\,200$ GeV$^2$, $x> 8 \cdot
10^{-5}$ and $66<W<233$ GeV.

For {$Q^2$}  values below about 10 GeV$^2$ the inclusive
transverse energy flow $1/Nd{\mbox{$E_T^*$}} /d\eta^*$ shows a plateau-like behaviour
in the current hemisphere with typical values of about 2 GeV.
As {$Q^2$}  increases, the distribution shows a peak of increasing
magnitude.
For values of {$Q^2$}   greater than approximately 70 GeV$^2$, the
peak position  coincides with the origin of the Breit frame. Transverse
energy flow
in the vicinity of the proton remnant is observed to be
approximately half of that in the current hemisphere, albeit
with large systematic uncertainties.

Transverse energy production in the
central pseudorapidity region rises with increasing $W$ which is
consistent with observations made in
 hadron-hadron experiments. However, for the first time
there
is evidence of an increase in the level of transverse energy in the
central rapidity region with $Q^2$.

Four QCD based models have been compared to the data.
  At low $x$ and  $Q^2$,
predictions made using approaches based on DGLAP and CCFM evolution,
implemented within  the {\sc Lepto}
and {\sc LDCMC} models, respectively,  are not able to
fully describe  the transverse
energy flow in the central
pseudorapidity region.  {\sc Rapgap}, which includes resolved
virtual photon processes, and the Colour Dipole Model, as implemented
in {\sc Ariadne}, give a reasonable description of the data
in this region. Furthermore, the {\sc Lepto} model can only
describe the measurements over the full current hemisphere at low {$Q^2$}
 if a General Area Law based string reinteractions scheme is not used.  At
the highest
{$Q^2$}
values, all models describe the data well except for {\sc Herwig} which
provides a relatively poor description over most of the kinematic range
considered in this analysis.

\section*{Acknowledgements}
We are grateful to the HERA machine group whose outstanding efforts have made
and continue to make this experiment possible. We thank the engineers and
technicians for their work constructing and maintaining the H1 detector, our
funding agencies for financial support, the DESY technical staff for continual
assistance and the DESY directorate for the hospitality which they extend
to the non-DESY members of the collaboration.

\include{bibli_final}

\end{document}

\section{Introduction}
Measurements of hadronic final state quantities
are extremely useful in investigating the different QCD processes that
occur in the wide range of phase space made accessible by the
 $ep$ collider HERA.
 One such quantity is transverse energy,
measurements of
which contain global information about charged and neutral particles
and cover a wider pseudorapidity range than  equivalent available charged
track analyses~\cite{ua1,charged,mk1}.

Electron-proton scattering for values of $Q^2$, the virtuality of the
exchanged boson,
significantly above $1\,$
GeV$^2$ is
usually considered
as a deep inelastic scattering (DIS) process, in which an exchanged
 boson directly couples to a parton in the proton.
This approach successfully describes inclusive
cross-section measurements~\cite{h1f297,zeusf297},
provided that
appropriate parton distribution functions are used to describe the partonic
content of the proton. Several approaches, based  on the
DGLAP~\cite{DGLAP},
BFKL~\cite{bfkl} and CCFM~\cite{CCFM} equations are available for
the QCD evolution of these parton densities to an appropriate scale before
the interaction with the exchanged boson.

Extending this picture to describe the hadronic final state
introduces a number of further complications. It becomes not only
necessary to understand what happens to the parton involved in the
partonic scattering process in more detail, but also to understand the
effects of the interaction on the entire proton. Thus the influence of
the evolution process leading to the parton undergoing the hard
scattering must be modelled, as must the behaviour of the proton
remnant and the fragmentation process.
Measurements of transverse energy flow~\cite{mk1,mk2} have proven useful
in discriminating between the different approaches used
in these QCD models.
For example, it was shown that early Monte Carlo models based on
DGLAP evolution tended to produce insufficient transverse energy in the
region near the proton remnant for values of the Bjorken
scaling variable $x$ of less than about $10^{-3}$.

Recent measurements of jet and leading particle production in
DIS~\cite{tania,zeusfj,dave} suggest that
the description of the data provided by DGLAP based models
can be improved by allowing the virtual photon to have structure. That
is, in addition to ``direct photon'' events, in which the entire momentum
of the virtual photon enters the hard scattering process, ``resolved
photon'' events are allowed. In these, the virtual photon is considered
to have a developed partonic structure and one of the partons in the
photon participates in the hard interaction. The significance of this
idea for the description of the transverse energy flow, in the framework
of the above models of DIS, is also investigated in this paper.

The deep inelastic $ep$ scattering may also be viewed in the rest frame of
the proton. In this frame
  the photon can be considered to fluctuate into a hadronic object which
subsequently interacts with the proton,
even for $Q^2$ values of up to $1\,000$ GeV$^2$~\cite{bjoko,duca,shaw}.
In the proton rest frame and for values of $Q^2$ above several GeV$^2$
the fluctuation time~\cite{ioffe,duca}
is given  by $\tau \approx 1/(xM_p)$, where $M_p$ is the mass of the hadronic target.
For values of  $x$ of less than $10^{-2}$ the virtual photon can fluctuate into
 and exist as a hadronic object over a distance of $10$ to $1\,000$ fm,
which
is far larger than the size of the proton. Under the na\"{\i}ve assumption
that this fluctuation behaves like a
single hadron then
$ep$ interactions should resemble hadron-hadron
scattering.
This approach was tested in an earlier publication~\cite{andre1},
in which it was demonstrated that the average transverse energy flow in
the central rapidity region in the hadronic centre of mass system showed
no significant
$Q^2$ dependence. This is consistent with observations  made in
hadron-hadron scattering, in which particle production in the central
region is largely independent
of the nature of the hadronic object and depends only
on the total centre of mass energy~\cite{na22}. Conversely,
transverse energy production in the  expected photon fragmentation region
was found to
be  strongly dependent on $Q^2$.
The
data collected by H1 in recent years now allow
this picture to be studied with far greater precision than was possible in the
previous work~\cite{andre1}.

  It is the aim of this analysis to study transverse energy production
within both the traditional DIS framework and in the picture given
in the proton rest frame. Distributions showing the
dependence of transverse energy on pseudorapidity, $x$, $Q^2$ and
the total hadronic mass, $W$,
are
presented and  quantitative comparisons of QCD based models are made
with the data. A
more qualitative approach is adopted in examining the
interpretation of our data as
a hadron-hadron scattering process.  The $Q^2$ and $W$
dependence of transverse energy production in the central pseudorapidity
region and a region associated with the fragmenting photon  is
investigated and
compared to hadron-hadron data and earlier H1 results in
both photoproduction and DIS.

The analysis presented here is based on data taken at
the HERA collider for which $820$ GeV protons were collided with
$27.5$ GeV positrons at a centre of mass energy of $300$ GeV.
Transverse energy production in the hadronic centre of mass system
is studied in  the kinematic range
$3.2$ GeV$^2$ $< Q^2<$ $2\,200$ GeV$^2$,
$x> 8 \cdot 10^{-5}$ and $66<W<233$ GeV. An order of magnitude more data
are
used than
 in the previous H1 measurements~\cite{mk2,andre1}. Furthermore, the phase
space region has been extended into the proton remnant fragmentation
region and now covers more than 8 units of pseudorapidity.

\section{QCD-based Models}\label{QCDpre}

Although progress has been made on several other
fronts~\cite{evs,mlla,lattice,neweth} it is still the case that most QCD
based
predictions of the hadronic final state are produced with Monte Carlo
methods which use  the following prescription. Phenomenologically derived
parton distribution functions, evolved to the relevant scale, are used to
determine the properties of the partons  emerging from the initial state
hadron, or photon. Some of these partons
undergo a hard scatter, the cross-section for which
is calculated using leading order (LO) QCD. The resulting partons
undergo fragmentation to produce the observed hadrons.

The proton parton distribution functions are reasonably well constrained
by inclusive measurements of DIS (in particular $F_2$), and data from
hadron-hadron scattering experiments. The effects
on the hadronic final state of changing the
input parton distribution functions within the limits allowed by
the inclusive DIS data are small compared to the effects arising from
varying other aspects of the calculations. The proton parton
distribution functions used here are the CTEQ4L~\cite{CTEQ,pdflib}
set. For calculations involving resolved virtual photon processes,
 the SAS-G virtual photon parton parameterisation~\cite{sasgam} is
used.

Different QCD evolution equations are known. They
have been derived in the
Leading-Log-Approximation and  are expected to be valid in certain regions
of Bjorken $x$. The DGLAP evolution
equation effectively resums the leading
$\log Q^2$ terms which corresponds to the
strong ordering in transverse momentum of successive
parton emissions and is applicable at large $Q^2$.
The BFKL approach sums leading $\log(1/x)$
terms and is expected to become significant at small $x$.
Successive parton emissions in this approach have strongly decreasing
longitudinal momenta, but are not ordered in transverse momentum.
This latter feature is emulated in parton emissions produced
by the Colour Dipole Model (CDM)~\cite{dipole} in which
partons are radiated from colour dipoles produced in
the hard interaction~\cite{cdmbfkl}.
 The CCFM evolution equation forms a bridge between the DGLAP and
BFKL approaches and resums the leading $\log 1/x$ and $\log Q^2$
terms both for inclusive and non-inclusive quantities.
To obtain CCFM based hadronic final state predictions,  the Linked Dipole
Chain~\cite{ldc} model is used here. This is a reformulation of
the CCFM equation and redefines the separation of initial and final
state  QCD emission using the CDM.

 The fragmentation of the produced partons typically involves a
showering process followed by a hadronisation phase.
 The models to which the data are compared here use the
string~\cite{string} or cluster~\cite{cluster}
approach to hadronisation.

Although the evolution equations discussed above can be derived within the
picture of a parton cascade, neither the DGLAP nor BFKL
approximations describe the details of the hadronic final state. It is
therefore necessary to reformulate these inclusive equations as an
iterative process, the shower algorithm, to produce
parton emissions in the Monte Carlo models.\footnote{The CCFM
approach is an iterative process by definition.}
The amount of transverse energy depends on these shower algorithms
and
hence
the
expected levels
of $E_T$
are
sensitive to
the
different evolution
schemes.

The Monte Carlo models used here are based on various combinations
of the above ideas. They are described briefly in the following text as
are any changes to the default settings and modifications of
the Monte Carlo programs.

{\sc Ariadne}~\cite{ariadne} version 4.10. This is an implementation
of the CDM. In order to generate an increase of the transverse energy with
$Q^2$ the
program has been modified; an additional switch has been
introduced~\cite{arpr}
and is set to MHAR(151)=2. This changes the phase space
restriction for
radiation from an extended source. In addition, two control
parameters were changed from their default values. The parameter settings used
are PARA(10)=1.5 and PARA(15)=0.5~\cite{tune}. These alter the suppression
of
radiation from the proton remnant and the struck quark, respectively.
The expectations of this model are marked ``{\sc Ariadne} 4.10 mod'' on
all figures.

The LDCMC~\cite{ldcmc}, version 1.0. This is based on the Linked Dipole
Chain implementation of CCFM evolution and uses the
CDM to simulate the effect of higher order emissions.
It is incorporated in the framework of the {\sc Ariadne} package. All
parameters used by both programs are set to their default values,
with the exception of those listed above. The expectations of this model
are marked ``LDCMC'' on all figures.

{\sc Lepto} 6.51~\cite{lepto}. {\sc Lepto} matches exact first order QCD
matrix elements to DGLAP based leading log parton showers.
It also allows two different methods of non-perturbative rearrangement of
the event colour topology. One way is via ``soft colour interactions''
(SCI)~\cite{SCI}, in which low momentum gluons are exchanged between
partons in the proton remnant and those which are perturbatively produced
and this is the default option for this version of {\sc Lepto}.
However,  this scheme leads to an excess of soft particle production at
high $Q^2$~\cite{mcworkshop}. In a different approach, a string
reinteraction scheme based on a  Generalised Area Law (GAL) is now
available~\cite{GAL}. This allows interactions between the colour
strings connecting the final state partons, leading to a reduction
in the total string area.
To use this option,
retuned values of parameters relating to the hadronisation and
 the parton shower schemes have to be
used.\footnote{The following parameter settings are used with
the GAL string re-interaction option: LST(34)=3, PARJ(42)=0.45,
PARJ(82)=2.0, PYPAR(22)=4.0, and PARL(7)=0.10}
However, as is demonstrated later in this report, the implementation
of this scheme leads to a worsening of the description of
our data.  Because of this,  the GAL option is not used in
comparisons with most of the spectra presented here and model expectations
which use neither GAL nor SCI are marked ``{\sc Lepto} 6.51 mod'' on
these figures.

{\sc Rapgap} 2.06/48~\cite{rapgap}. {\sc Rapgap} also matches exact first
order QCD matrix
elements to DGLAP based leading log parton showers. In addition
to the direct photon processes simulated by {\sc Lepto}, {\sc Rapgap}
simulates resolved photon interactions in which the virtual photon is
assumed to have structure, parameterised using an implementation of the
SAS-G virtual photon parton distributions. The {\sc Rapgap}
package allows a choice of renormalisation and
factorisation scales. For predictions presented here,  these are set to
$p_T^2 + Q^2$, where $p_T$ is the transverse momentum of the partons in
the hard scattering.
The default value for the cut-off necessary to regulate the
matrix elements is $p_T^2 > {\rm PT2CUT} = 1.0 \,$GeV$^2$.
This leads to a predicted DIS cross-section which is larger than
measurements
for $Q^2 > 200\,$GeV$^2$, so here ${\rm PT2CUT} = 4.0\,$GeV$^2$ is used.
The expectations of this model are marked ``{\sc Rapgap} 2.06/48 dir + res'' on
all figures. If only direct photon interactions are
allowed, the results of {\sc Rapgap} are very similar to those of {\sc Lepto}.

To perform the hadronisation step,  all of the above models use string
fragmentation as implemented in JETSET 7.4~\cite{jetset74}.

{\sc Herwig} 5.9~\cite{herw}.
 This model is based on leading log parton
showers with
matrix element corrections~\cite{mecor} and implements a
 cluster hadronisation scheme. The expectations of this model are marked
``{\sc Herwig} 5.9'' where shown.

The {\sc Phojet}~\cite{phojet} Monte Carlo program is
used to calculate the amount of background due to photoproduction
processes. This generator contains LO QCD matrix elements
for hard
subprocesses, a  parton shower model and a
phenomenological description of  soft processes.

As an alternative to the Monte Carlo based approach,
        analytical BFKL calculations~\cite{etbfkl}
 are available to
 predict energy flow in the forward region.
  These calculations are
based on asymptotic expressions derived from the BFKL equation
at LO and do not include hadronisation
effects. The predictions of the calculations
are marked ``BFKL Partons''
where shown.


\section{Introduction}

The study of double spin asymmetries in appropriate  final states is
a very
promising way to measure the polarized gluon distribution function in
the nucleon. The  favored processes are those
for which the involved production cross sections and subprocess asymmetries
can be calculated in the framework of perturbative QCD
and be predicted with small uncertainties.
Recently, the measurement of $J/\psi$ photoproduction asymmetries
was suggested as a tool for direct access to the  polarized gluon density
\cite{HERMES,COMPASS}; the  analysis of asymmetries at the
 subprocess level  was based on results obtained in the framework
 of the color singlet model (CSM) \cite{CSM,Guillet}.
However, the CSM can not  be considered as a  sufficiently
reliable model to describe heavy quarkonium production and decay processes.
The measured cross section of prompt $J/\psi$ and $\psi'$ production at
large $p_T$ at the Tevatron exceeds the CSM predictions by  more than one order magnitude \cite{Tevatron}.
In photon-proton collisions, even after taking into account
next-to-leading order corrections \cite{Kraemer}, the CSM does not
describe satisfactorily different distributions in $J/\psi$
production \cite{H1}.

The factorization approach (FA) based on
non-relativistic QCD (NRQCD) represents a more  reliable framework to
study heavy quarkonium production and decay processes \cite{NRQCD}.
In an extension of the  CSM \cite{CSM}, the NRQCD FA
implies that  the heavy quarkonium can be produced through color octet
quark-antiquark  states, as well. The production of a heavy quark-antiquark
 pair can be
calculated perturbatively because the hard scale is
determined by the mass of heavy quark,
whereas the evolution of the $(Q\bar Q)$ pair to the final hadronic state
is parameterized by nonperturbative  matrix elements,
$\langle{\cal O}^{\psi}_n\rangle$.
The relative importance of the different
intermediate states, both  color octet and color singlet,
is   defined by  velocity scaling rules \cite{LMNMH}.
The color singlet long distance parameters,
being related to the quarkonium wave function at the origin, can be calculated
in the framework of  nonrelativistic models,
 while the color octet parameters have to be extracted
from  experimental data.
Once this is done,  predictions for other processes can be made.
 Unfortunately, the present-day theoretical
uncertainties are still too large, preventing a determination of  these
 parameters with  enough accuracy
 to check the universality of the NRQCD FA  \cite{BK}-\cite{KK}.

In the framework of the color octet model the main contribution to $J/\psi$
photoproduction is coming from the color octet states
$^1S_0^{(8)}$ and $^3P_J^{(8)}$ \cite{CK,KLS}. Throughout this paper we
exploit the fact that the corresponding long distance matrix elements can
not be determined  separately from the
prompt $J/\psi$ production data at the Tevatron;
it is only possible to extract their combination.

Spin effects in heavy quarkonium production, such as
$J/\psi$ \cite{BK,BKV} and $\Upsilon$ polarization \cite{KLST}-\cite{T}
in unpolarized experiments and spin asymmetries in hadroproduction
\cite{TT,GM} can provide a possibility to check the NRQCD factorization
approach. As ratios of  cross sections, spin-dependent
 observables do not strongly depend neither on the mass of heavy quarks,
nor on the  factorization and renormalization scale.
These parameters are crucial in the prediction of cross sections
and they usually give rise to large uncertainties in the extraction of
color octet parameters \cite{BK}.

In the present paper we consider
double spin asymmetries in $J/\psi$ photoproduction in the framework of the
NRQCD FA. We investigate the uncertainties depending on the
long distance color octet parameters and the possibility to use $J/\psi$
production asymmetries for a determination of the  polarized gluon density.
We will argue that a measurement of the polarized gluon density via double spin asymmetries in
charmonium photoproduction can potentially provide an additional test of
the NRQCD FA, or alternatively, the polarized gluon distribution can be
measured once the FA is confirmed.

The deviation from the CSM prediction can also
be used as a test of an  alternative model for heavy quarkonium production
proposed by
Hoyer and Peigne \cite{HP} which reduces to the CSM in the case of
$J/\psi$ photoproduction.

In the next section we consider $J/\psi$ photoproduction and electroproduction
asymmetries at  the subprocess level.
 The uncertainties in the color octet long
distance parameters are discussed in Sec. 3. The
double spin asymmetries for various sets of long distance matrix elements
and different parameterizations of polarized parton densities are calculated
at   $\sqrt{s}=10$ GeV which can be considered as representative for
both the HERMES and the COMPASS experiment \cite{HERMES,COMPASS}.
We also estimate the expected double spin asymmetries in $J/\psi$
electroproduction for a possible polarized electron-proton collider with
energy $\sqrt{s}=25$ GeV \cite{EPIC}.

\section{Double Spin Asymmetries for $J/\psi$  Subprocesses}

\subsection{$J/\psi$ Photoproduction}
The double-spin asymmetry $A_{LL}$ for  inclusive $J/\psi$ photoproduction
 is defined as
\begin{eqnarray}
A_{LL}^{J/\psi}(\gamma p) = \frac{
 d\sigma(\overrightarrow{\gamma}+\overrightarrow{p}\to J/\psi+X)
-d\sigma(\overrightarrow{\gamma}+\overleftarrow{p}\to J/\psi+X)}
{d\sigma(\overrightarrow{\gamma}+\overrightarrow{p}+\to J/\psi+X)
+d\sigma(\overrightarrow{\gamma}+\overleftarrow{p}\to J/\psi+X)}=
\frac{d\Delta\sigma/d^3p}{d\sigma/d^3p}~,
\end{eqnarray}
where the arrows over $\gamma$ and $p$ denote the helicity projection on
the direction of the corresponding momenta of the initial particles.

Throughout this paper we use the standard notation of spectroscopy
 to define the  quantum numbers of intermediate quark-antiquark states,
 $^{(2S+1)}L_J^{(1,8)}$. Here $S$ denotes the spin
of quark-antiquark state, $J$ is the total angular momentum, the superscripts
$(1)$ and $(8)$ correspond to color singlet and color octet states
respectively.
The  long distance parameter
$\langle{\cal O}_{1,8}^{H}(^{(2S+1)}L_J)\rangle$ describes the transition
from the  $^{(2S+1)}L_J^{(1,8)}$-state to a hadron $H$.

In  lowest order of $\alpha_s$ only  color octet states contribute to
$J/\psi$ production, namely through the following subprocesses:
\begin{eqnarray}
\gamma+g &\to& ^1S_0^{(8)}~,\\
\gamma+g &\to& ^3P_{0,2}^{(8)}
\end{eqnarray}
These heavy quarkonium states  are produced at the kinematical
endpoint in $2\to1$ subprocesses,  i.e. at  $z\simeq1$, where
$z= P_{\psi}\cdot P_{p}/P_{\gamma}P_{p}$.
The double-spin asymmetries for these subprocesses were calculated
in Ref. \cite{GM}. Uncontrolled uncertainties at $z\simeq1$ arise
from higher order terms in the NRQCD velocity expansion;
in this particular
case the higher order terms are enhanced by the kinematical factor $1/(1-z)$
which makes all predictions unreliable close to the kinematical boundary
 \cite{BRW,FM}.
In the same kinematics, the $J/\psi$ is produced via the diffractive process,
 which is an additional source of uncertainty.
Another kinematical characteristic of the $2\to1$ subprocesses, Eqs. (2)-(3),
is the small transverse momentum of the produced $J/\psi$; the transverse
momentum of the final hadronic state is determined mainly by internal
motion of the colliding partons and by the
recoil of emitted soft gluons in the hadronization phase.

All these uncertainties are smaller in $J/\psi$ production at larger transverse
momenta. The leading order contribution in this case
 comes from  $2\to2$ subprocesses which are of the order of
 $\alpha\alpha_s^2$:
\begin{eqnarray}
\gamma+g &\to& J/\psi +g~,\\
\gamma+g(q) &\to& ^1S_0^{(8)}+g(q),\\
\gamma+g(q) &\to& ^3S_1^{(8)}+g(q),\\
\gamma+g(q) &\to& ^3P_J^{(8)}+g(q).
\end{eqnarray}
The double spin asymmetries for  $J/\psi$ production in the
framework of the CSM, whose corresponding subprocess is described by Eq. (4),
  were studied in Ref. \cite{Guillet}.
We calculated the cross sections for the production of a color octet
heavy quark-antiquark pair, Eqs. (5)-(7),
for different helicity states of initial photon and gluon (quark).
he expressions for $\Delta\sigma$ for photon-gluon subprocesses which
give the dominant contribution to the $J/\psi$ production
are listed in the Appendix\footnote{ The FORTRAN
codes can be provided on request}.
Our results after averaging the spins of the initial particles  are in
agreement with those  used  for the calculation of the $J/\psi$ photoproduction
cross section in the NRQCD FA \cite{CK}\footnote{
 We are grateful to Michael Kr\"amer for providing us with the
analytical results for subprocess cross sections.}.

\begin{figure}[th]
\vspace{-4mm}
\centering
\begin{minipage}[c]{7.5cm}
\centering
\epsfig{file=fig1a.eps,width=7.5cm}
\end{minipage}
\hspace*{0.5cm}
\begin{minipage}[c]{7.5cm}
\centering
\epsfig{file=fig1b.eps,width=7.5cm}
\end{minipage}
\caption{\small
Different contributions to double spin asymmetries in quarkonium
photoproduction,
 a) vs. $z$  and  b) vs. $p_T$, all for the GS LO
 parameterization (set A). The solid line
corresponds to the color singlet $^3S_1$ state,  the dashed line to the
color octet
$^1S_0^{(8)}$ state, the dash-dotted line to the color octet $^3P_J^{(8)}$
states and the dotted line to the color octet $^3S_1^{(8)}$ state.
}
\end{figure}

The resulting double spin asymmetries for the photoproduction of various
heavy quark-antiquark  states at $\sqrt{s}=10$ GeV
are shown in Figs.  1(a) and 1(b). They
are  calculated for the Gehrmann-Stirling (GS) leading order (LO)
parameterization (set A) of the
polarized parton distribution functions (PDFs) \cite{GS}.
For the unpolarized parton distribution functions the Gl\"uck-Reya-Vogt (GRV)
 LO parameterization
evaluated at $Q^2=p_T^2+4m_c^2$ was used \cite{GRV}.
As can be seen, the asymmetries for the $^1S_0^{(8)}$ (dashed lines) and
$^3P_J^{(8)}$ (dash-dotted lines)
octet states are
 rather different from the CSM prediction (solid line).
Therefore, the predicted asymmetry for $J/\psi$ production  in the NRQCD FA
is expected to be sensitive to long distance color octet parameters.

\subsection{$J/\psi$ Electroproduction}

We also calculated the double spin asymmetries in heavy quarkonium
electroproduction.
The virtuality of the colliding photon can be used as an additional large scale in
the process.
Requiring  $Q^2>>4 m_c^2$,  some theoretical uncertainties connected with
the enhancement of higher order corrections in the NRQCD
expansion can be avoided.
\cite{BRW,FM}. At the same time it is possible to reduce the contribution
from higher twist effects and from $J/\psi$ diffractive production ( for a
detailed discussion see Ref. \cite{FM}).

In the  framework of the CSM, $J/\psi$ can be produced
  in electron-proton collisions
 only in order $\alpha^2\alpha_s^2$
subprocesses, $e+g\to e+(c\bar c)+g$. The asymmetries for these subprocesses
were calculated in Ref. \cite{Guillet}.
However,
the main contribution to the total $J/\psi$ electroproduction cross section
comes not from CSM, but from the leading order
$\alpha^2\alpha_s$ subprocesses,
similar to (2),(3), when
color octet states can be produced only with small transverse momenta
in the center-of-mass system of virtual photon and proton.
In this paper we consider the double spin asymmetries only for these
leading order subprocesses.
The next-to-leading order color octet part, $O(\alpha^2\alpha_s^2)$,
is even smaller than the color singlet model prediction \cite{FM} and does
 not significantly change
the asymmetry calculated in the framework of the CSM in Ref. \cite{Guillet}.

In the NRQCD factorization approach in  leading order  QCD
$J/\psi$ can be produced through  two color octet states,
 $^1S_0$ and $^3P_J$.
The expression for the corresponding cross section was obtained
in Ref. \cite{FM}:
\begin{eqnarray}
&&\sigma (l+p\to l+J/\psi+X) =
\int{\frac{ dQ^2}{Q^2}} \int{d y \frac{G(x)}{y S} }
\; \frac{2\pi^2 \alpha_s(\mu^2) \alpha^2 e_c^2}{m_c(Q^2+4m_c^2)}
~~~~~~~~~~~~~~\nonumber \\
&&~~~~~~~~~~~~\times \Biggl\{\frac{1+(1-y)^2}{y}
 \times  \left[\langle {\cal O}_8^{J/\psi}(^1S_0)\rangle +
\frac{3Q^2+28 m_c^2}{Q^2+4m_c^2}
\frac{\langle{\cal O}_8^{J/\psi}(^3P_0)\rangle}{m_c^2}\right] \\
&&~~~~~~~~~~~~ -y\frac{32 m_c^2 Q^2}{(Q^2+4 m_c^2)^2}
\frac{\langle{\cal O}_8^{J/\psi}(^3P_0)\rangle}{m_c^2}\Biggr\},
\nonumber
\end{eqnarray}
where $G(x)$ is the gluon distribution in the proton, $S$ is the
invariant mass  of the electron-proton system, and the momentum fraction
of the struck
parton is given by $x=(Q^2+4 m_c^2)/y S$.

The corresponding polarized cross section has the form
\begin{eqnarray}
&&\Delta\sigma(l + p\to l +J/\psi+X) =
\frac {1}{2}(\sigma(\overrightarrow{l}\overrightarrow{p})
-\sigma(\overrightarrow{l}\overleftarrow{p})) = \nonumber \\
&&~~~~~~~~~~~~~~~\int{\frac{dQ^2}{Q^2} } \int{d y \frac{\Delta G(x)}{y S} }
\; \frac{2\pi^2 \alpha_s(\mu^2) \alpha^2 e_c^2}{m_c(Q^2+4m_c^2)}\\
&&\times \frac{1-(1-y)^2}{y}
 \left[\langle {\cal O}_8^{J/\psi}(^1S_0)\rangle +
\frac{3 Q^4+8m_c^2 Q^2-16 m_c^4}{(Q^2+4m_c^2)^2}
\frac{\langle{\cal O}_8^{J/\psi}(^3P_0)\rangle}{m_c^2}\right],
\nonumber
\end{eqnarray}
where $\Delta G(x)$ is the polarized gluon distribution  in the
proton. The double spin asymmetry is defined by the ratio
$A_{LL} = \Delta\sigma(lp)/\sigma(lp)$.

In the limit  $Q^2\to0$ the expressions for the cross sections
reduce to the convoluted photoproduction cross sections \cite{CK,GM}
with polarized and unpolarized splitting functions correspondingly:
\begin{eqnarray}
\lim_{Q^2\to0} \sigma(l+p\to J/\psi+X) &\to& \frac{\alpha}{2\pi}
\int{\frac{dQ^2}{Q^2}}\int{dy \frac{1+(1-y)^2}{y}\hat\sigma(\gamma+P\to J/\psi+X)},
\\
\lim_{Q^2\to0}\Delta\sigma(l+p\to J/\psi+X) &\to& \frac{\alpha}{2\pi}
\int{\frac{dQ^2}{Q^2}}\int{dy \frac{1-(1-y)^2}{y}\Delta\hat\sigma(\gamma+P\to J/\psi+X)}.
\end{eqnarray}

\section{Matrix Elements and Results}

In $J/\psi$ photoproduction the dominant contribution comes from the
color singlet state $^3S_1$ and the color octet states
 $^1S_0^{(8)}$ and $^3P_J^{(8)}$.
The production cross section
for the $^3S_1^{(8)}$ octet state is much smaller than the cross sections
for these two octet states \cite{CK}.
 The color singlet matrix element is connected to the nonrelativistic wave
 function at the origin and can be evaluated in the framework
of potential models. It can also be extracted from the leptonic decay width
 of the $J/\psi$.
The corresponding color octet long distance matrix elements
 can not be calculated
in  perturbative QCD (pQCD) and are usually extracted from  experimental data.
For the color octet $P$-wave states they
can be reduced to one parameter using the NRQCD spin symmetry relation:
\begin{eqnarray}
\langle{\cal O}_8^H(^3P_J)\rangle & = & (2J+1)
\langle{\cal O}_8^H(^3P_0)\rangle.
\end{eqnarray}

After using these relations there are only two matrix elements left
which give the main octet contribution to $J/\psi$ photoproduction,
$\langle{\cal O}_8^{\psi}(^1S_0)\rangle$ and
$\langle{\cal O}_8^{\psi}(^3P_0)\rangle$.
From the available  data of prompt $J/\psi$ production at the Tevatron
it is only possible to extract  the combination
of these matrix elements with rather large uncertainties \cite{CL,BK,KK}.
Apparently, at Tevatron energies  the states
 $^1S_0^{(8)}$ and $^3P_J^{(8)}$ contribute
significantly
at small transverse momenta of $J/\psi$, $p_T<5$ GeV.
Consequently, the fitted value of the combination of matrix elements
$\langle{\cal O}_8^{\psi}(^1S_0)\rangle$ and
$\langle{\cal O}_8^{\psi}(^3P_0)\rangle$, defined as \cite{KK}
\begin{eqnarray}
M_r^{J/\psi} &=& \langle{\cal O}_8^{J/\psi}(^1S_0)\rangle
+\frac{r\cdot\ \langle{\cal O}_8^{J/\psi}(^3P_0)\rangle}{m_c^2}~,
\end{eqnarray}
is very sensitive to effects
which change the $p_T$ spectrum, such as dependence on factorization and/or
renormalization scales, choice of parameterization of parton distribution
 functions
etc. \cite{BK}. The number $r$ is fitted from experimental data.

To test the sensitivity of double spin asymmetries on the color octet matrix
elements we used  two sets of parameters extracted from the
Tevatron data on prompt $J/\psi$ hadroproduction \cite{KK}, as shown in
 Table I:
\begin{center}
\begin{tabular}{c|c|c|c}
\hline \hline
          & LO   &  HO   & scaling \\
\hline
$\langle{\cal O}_1^{J/\psi}(^3S_1)\rangle$ &
$(7.63\pm0.54)\cdot10^{-1}$ GeV$^3$ & $ (1.30\pm0.09)$ GeV$^3$ &
$[m_c^3 v^3]$ \\
$\langle{\cal O}_8^{J/\psi}(^3S_1)\rangle$ &
$(3.94\pm0.63)\cdot10^{-3}$ GeV$^3$ & $ (2.73\pm0.45)\cdot10^{-3}$ GeV$^3$ &
$[m_c^3 v^7]$ \\
$M_r^{J/\psi}$  &
$(6.52\pm0.67)\cdot10^{-2}$ GeV$^3$ & $ (5.72\pm1.84)\cdot10^{-3}$ GeV$^3$ &
$[m_c^3 v^7]$ \\
$r$ & 3.47 & 3.54 &  \\
\hline
\end{tabular}
\end{center}

\noindent
 Table I. {\small Values of the long distance matrix elements from LO- and
 HO-improved fits to the CDF data.}
\vspace{0.3cm}

\noindent
The LO  set of the matrix elements corresponds to the leading-order  fit
of the prompt $J/\psi$ hadroproduction data at the Tevatron and the obtained
value for $M_r^{J/\psi}$ is the largest among existing fits in the
literature. These values of matrix elements lead to an increase
of the photoproduction  cross section as $z\to1$ and dramatically
overestimate the H1 and ZEUS data on $J/\psi$  photoproduction at HERA.

\begin{figure}[ht]
\vspace{-4mm}
\centering
\begin{minipage}[c]{7.5cm}
\centering
\epsfig{file=fig2a.eps,width=7.5cm}
\end{minipage}
\hspace*{0.5cm}
\begin{minipage}[c]{7.5cm}
\centering
\epsfig{file=fig2b.eps,width=7.5cm}
\end{minipage}
\caption{\small
  (a) $z$ and  (b) $p_T$ dependence of double spin asymmetries of
$J/\psi$ photoproduction
for LO and HO sets of color octet matrix elements. The solid curves correspond
to the CSM prediction. The dashed curves correspond to the case when
$\langle{\cal O}_8^{J/\psi}(^1P_0)\rangle = 0$, the dash-dotted ones
to  $\langle{\cal O}_8^{J/\psi}(^1S_0)\rangle = 0$.
}
\end{figure}

When  fitting the  HO set, the higher-order  corrections were taken into
account only
approximately by considering the initial and final state gluon radiation
in the framework of a Monte Carlo simulation \cite{CS,KK}.
Including effectively the primordial transverse momentum of partons
increases the production cross sections of $^1S_0^{(8)}$ and $^3P_J^{(8)}$
color octet states at small values of $p_T$. The impact of this effect is a
significant decrease of the $M_r^{J/\psi}$ combination by about one order of
magnitude, as can be seen from table I.
It is worth  noting that this
 type of HO corrections changes the photoproduction cross
sections for the color octet states only by about $20\%-30\%$
at  H1 and ZEUS energies \cite{KK}.
The HO parameters, presented in Table I, are the smallest ones from the
existing fits and are compatible with the H1 and ZEUS data on the
$z$-distribution of the $J/\psi$ {\it photo}production cross section.
We consider this set as an example of small color octet parameters.
The combination of color octet parameters $M_r^{J/\psi}$ was
 also extracted from
fixed target {\it hadro}production data \cite{BR,GuS} for $r=7$,
 $M_7^{J/\psi}\simeq 2\div3 10^{-2}$ GeV$^3$. We note that these values lie in
between the LO and HO sets of parameters, as given Table I,
We used this range  to check the sensitivity
of the expected asymmetries on the model of $J/\psi$ production.

In  Fig. 2(a) and 2(b) the expected asymmetries in $J/\psi$ photoproduction,
calculated at  $\sqrt{s}=10$ GeV
for the   GS LO  set A parameterization of polarized PDF's,
are shown in dependence on $z$ and $p_T^2$.
As can be seen  from  these figures, the uncertainty
in the asymmetry of $J/\psi$ production is large for the LO set of color octet
matrix elements because  the  values of
$\langle{\cal O}_8^{J/\psi}(^1S_0)\rangle$ and
$\langle{\cal O}_8^{J/\psi}(^3P_0)\rangle$ are not fixed separately but only
 as the combination (13).
For the HO set
the contribution to the cross section
of the $^1S_0^{(8)}$ and $^3P_J^{(8)}$ color octet states is smaller, which
reduces the deviation from the  asymmetry predicted by the CSM
especially for lower $z$-values $(z<0.7)$.
To make a choice between the different
 sets of long distance
parameters,   more accurate measurements at energies lower  than that for
H1 and ZEUS  are required in the future.
The COMPASS experiment at CERN \cite{COMPASS} and later possibly a
 presently discussed polarized electron-proton  collider (EPIC)
at  $\sqrt{s}=25$ GeV \cite{EPIC} will work at higher
luminosities and may hence  provide such  measurements.
The gluon polarization itself at that time has  probably been  measured in
semi-inclusive reactions with less uncertainties, for example in
open charm production at COMPASS  or in dijet and direct photon
 production at RHIC \cite{RHIC}.
If, in this situation,   the HO set of parameters should  fit the  data of
the lower energy experiments, the
$J/\psi$ photoproduction asymmetries can be used to test the NRQCD FA.
Alternatively, based upon the HO set of parameters and relying on the NRQCD FA
approach, the polarized gluon distribution can be measured with
acceptable theoretical uncertainties, if only $J/\psi$ events with
$z<0.7$ will be used.
\begin{wrapfigure}{l}{7.5cm}
\vspace*{-4mm}
\centering
\epsfig{file=fig3.eps,width=7.5cm}
\vspace{-7mm}
\caption{\small
The double spin asymmetries for $J/\psi$ photoproduction at
$\sqrt{s}=10$ GeV using HO long distance parameters
versus $p_T$.
The definition of lines is the same as in the previous figure.
}
\end{wrapfigure}

Using the HO set of color octet parameters in  Fig.~3 the double spin
 asymmetries for $J/\psi$ photoproduction
are presented  at $\sqrt{s}=10$ GeV
for the  different sets A, B, C of
GS NLO parameterizations for polarized PDF's \cite{GS}.
We used here different  sets of NLO parameterization to show the sensitivity of
the asymmetries to the polarized gluon distribution function, which is rather
 different for these three sets.
We note that, as expected, the    result for the GS LO  set A
parameterization  is practically the same as for the GS NLO set B.
 We show only the $p_T^2$ dependence of the asymmetries.
As  was mentioned above,   uncontrolled
corrections are expected at large $z$  from higher order terms in the velocity
expansion  because of ignoring the mass difference between the
$(c\bar c)$-pair and the final state hadron.
This type of uncertainties can be neglected in the $d\sigma/dp_T$
distribution because the integration over $z$ smears the singular region out
\cite{BRW,FM}.
 When using the HO set of  color octet parameters,
the $p_T$ behaviour of the
 double
spin asymmetry for $J/\psi$ production
is rather well described by the CSM prediction as long as color octet
$P$ states don't contribute, as can be seen from Fig. 3.
 If, in contrast, the color octet $S$-state does
not contribute, the asymmetry decreases by about $30\%$. This behaviour
appears for $p_T>1$ GeV and seems to be rather independent of the actual NLO
parameterization set used for polarized PDF's.

In $J/\psi$ photoproduction experiments at low energy
higher twist effects are expected to be
suppressed by the  heavy quark mass, which defines the hard scale of the
corresponding subprocess, ~$\Lambda_{QCD}/m_c$.
Moreover, it is not excluded that higher twist effects are suppressed even
stronger by the transverse momentum of the produced quarkonium statee,
 or by factor $(m_{Q\bar Q}^2+p_T^2)^{-1/2}$.
A study of these effects may be performed
at high luminosity experiments such as eN at DESY TeV Energy Superconductiong
Linear Accelerator TESLA,\footnote{A project of
a high luminosity experiment  scattering of  TESLA electron beam off
polarized nucleons is under consideration \cite{WD}.} where enough statistics
can be accumulated at larger values of transverse momentum.

In  Fig. 4  electroproduction asymmetries are shown at
$\sqrt{s}=25$ versus photon virtuality $Q^2$, calculated using the 'largest'
 parameterization, namely GS NLO set A.
 As can be seen from this figure,
the difference between the asymmetries for the color octet  states
 $^1S_0^{(8)}$ and $^3P_J^{(8)}$ is large for a large gluon polarization
offering a powerful tool to discriminate between the
hitherto 'coupled' contributions from color octet $S$- and $P$-wave states.
It is important to note  that
in  leading order electroproduction
large contributions from higher twist effects and diffractive $J/\psi$
production are expected. They can be  suppressed
by requiring
large enough $Q^2>>4 m_c^2$  which leads, however,
to  a rapidly falling  cross section  and
large statistical errors are expected \cite{FM}.
It is a question of attainable statistics at  EPIC
whether it may be  possible to achieve  reasonable restrictions on the values
of the color octet parameters from  $J/\psi$ electroproduction
asymmetries.
At  high energies, $\sqrt{s}>100$ GeV, which correspond to the polarized
HERA option \cite{HERA}, the total cross section is large
but the asymmetries calculated with the
commonly used parameterizations GS or Gl\"uck-Reya-Stratmann-Volgelsang
(GRSV) \cite{GS, GRSV} are expected to be very small, less than $2\%$.

\section{Conclusions}

In the present paper  double spin asymmetries in $J/\psi$
photoproduction
 and electroproduction were studied in the NRQCD factorization approach.
 The asymmetries
for the color octet states $^1S_0^{(8)}$ and $^3P_J^{(8)}$ which give
the dominant octet
 contribution to the $J/\psi$ {\it photo}production cross section, are
significantly different from the CSM prediction.
The values of the corresponding color octet long distance matrix elements
 can presently not be extracted
separately from  existing data on $J/\psi$ hadroproduction.
Two sets of long distance parameters were used to show the sensitivity
\begin{wrapfigure}{r}{7.5cm}
\vspace*{-10mm}
\centering
\epsfig{file=fig4.eps,width=7.5cm}
\vspace{-7mm}
\caption{\small
The Double spin asymmetries for $J/\psi$ electroproduction at $\sqrt{s}=25$ GeV
for GS NLO (set a) parameterization.
The dashed  line corresponds to the asymmetry of $^1S_0^{(8)}$ state and
dash-dotted line to the asymmetries of $^3P_J^{(8)}$ states.
}
\end{wrapfigure}
of the expected asymmetries to this uncertainty.
For the 'large' set of parameters extracted from CDF data \cite{KK}
the asymmetries contain very large uncertainties.
For the 'small'  set of
parameters,
when higher order QCD corrections are partly  taken into account, the
predicted asymmetries for the two extreme  choices of parameters are
for $z<0.7$ rather close
to the CSM prediction. This 'small' set of long distance parameters is
preferable to explain  H1 and ZEUS data on inelastic $J/\psi$
photoproduction. The planned high luminosity measurements at
COMPASS   and
possibly at EPIC ($\sqrt{s}=25$ GeV)
may provide the  possibility to fit the color octet parameters
in {\it electro}production at
lower energies with better accuracy.
Here even a discrimination between $S$-wave and $P$-wave color octet
states appears feasible because the CSM contribution is subleading in
electroproduction.
If the 'small' set of parameters should remain  the prefered one,
the $J/\psi$ {\it photo}production asymmetries
can be eventually used to test the NRQCD FA.
Here it is  anticipated  that the gluon polarization
will have  been  measured already with small uncertainties
through  other semi-inclusive modes in polarized
electron-proton or proton-proton collisions.
On the other hand, once the situation with the NRQCD FA has been clarified
by other experiments, such as $J/\psi$ and bottomonium
polarization measurements,
the measurement of double spin asymmetries in
$J/\psi$  {\it photo}production at HERMES and COMPASS
can be used to extract the gluon
polarization in the nucleon in  LO  pQCD using events in the range
 $0.4<z<0.7$, because for $z<0.4$
the contribution of resolved photon subprocesses becomes essential.

\vspace{0.3cm}

{\bf Acknowledgment} We are grateful to Michael Kr\"amer for stimulating
and helpful discussions. This work is supported partly by the Alexander von
Humboldt foundation and by the National Science Foundation.

\newpage
\setcounter{equation}{0}
\defA.\arabic{equation}{A.\arabic{equation}}

\noindent
{\Large\bf
Appendix}

\vspace{0.3cm}
\noindent
Here we present the polarized cross sections for
the production of $(c\bar c)$ states in photon-gluon fusion subprocesses.
The polarized cross section is defined as follows
\begin{eqnarray}
\Delta\sigma & =& \frac{\sigma(++)-\sigma(+-)}{2},
\end{eqnarray}
where the first and secod signs denote
 the helicity states of colliding photon and gluon, respectively.
The expressions for unpolarized cross sections are listed in Ref. \cite{KLS}.
All cross sections were calculated for different helicity states of
colliding particles. The obtained results were checked by comparing
the unpolarized cross sections of this work, $(\sigma(++)+\sigma(+-))/2$,
 with those
obtained in \cite{CK,KLS}.
As an additional check, we calculated the asymmetries for the
production of scalar,
 $^1S_0^{(8)}$ and $^3P_0^{(8)}$, and tensor $^3P_0^{(8)}$
octet states in the limit
$s\to M_{J/\psi}$
 (i.e., at the threshold of  quarkonium pair production).
The subprocess level asymmetries for scalar and tensor states
approach 1 and -1,
respectively. The limit $4 m^2_c\to s$ implies that the emitted gluon in the
$2\to2$ subprocess $\gamma g\to(c\bar c)g$
becomes soft; the helicity properties of the amplitude then become the same
as those of the amplitude of the $2\to1$ process $\gamma g\to(c\bar
c)$.
It is easy to show that the asymmetries for  the process $2\to1$
for scalar and tensor states are 1 and -1, respectively;
 the scalar(tensor) states
are produced only in the case when the helicities of initial
 photon and gluon are
$[++]([+-])$,  which means that they  have a total momentum
projection $J_z=0(\pm2)$.
The existence of such limits serves as an additional test for
the  analytical calculations of the polarized cross sections $\Delta\sigma$.

The spin dependent differential cross section for $J/\psi$ production  through a
color singlet state has the form
\begin{eqnarray}
\frac{\Delta\sigma(\gamma g\to J/\psi\;g)}{dt} =
\frac{(ee_cg_s)^2\langle{\cal O}_1^{J/\psi}(^3S_1)\rangle}
{16 \pi s^2} \frac{32}{27 Q^2} \big\{(st+tu+su)^2-M^2stu\bigr\},
\end{eqnarray}
with $Q=(s+t)(s+u)(t+u)$.

 For the case of an intermediate $^1S_0^{(8)}$ color octet state the
corresponding cross section reads as

\begin{eqnarray}
&&\frac{\Delta\sigma(\gamma g\to ^1S_0^{(8)} g\to J/\psi X)}{dt} =
\frac{(ee_cg_s^2)^2\langle{\cal O}_8^{J/\psi}(^1S_0)\rangle}
{16 \pi M t Q^2} 3s u(s^4-t^4-u^4+M^8).
\end{eqnarray}
For  $J/\psi$ production through the color octet state $^3S_1^{(8)}$
the polarized cross section
can be obtained from (A.2) by changing the long distance
matrix element:
$\langle{\cal O}^{J/\psi}(^3S_1)\rangle\to
15/8\langle{\cal O}_8^{J/\psi}(^3S_1)\rangle$.

Spin dependent cross sections for  $^3P_J^{(8)}$
octet states have the form:
\begin{eqnarray}
&&\frac{\Delta\sigma(\gamma g\to ^3P_0^{(8)} g\to J/\psi X)}{dt} =
\frac{(ee_cg_s^2)^2\langle{\cal O}_8^{J/\psi}(^3P_0)\rangle}{\pi s^2Q^3 t}
384\biggl\{ 4 M^6 s1^{-1} t^5 t_1 t_2
\nonumber\\
&&~~~ +4 M^6 t^4 t_2 (3 t-2 M^2)- M^2 s_1 t^3  (3 t^4+2 t^3 M^2-41 t^2 M^4
\nonumber\\
&&~~~ -14  t M^6+3 M^8)
-M^2 s s_1^2 t^3  (41 t^2-118 t M^2+68 M^4)
\nonumber\\
&&~~~ -M^2 s^2 s_1^2 t  (73 t^3-191 t^2 M^2+139 t M^4-45 M^6)
\nonumber\\
&&~~~ -4 M^2 s^3 s_1^2 t  (19 t^2-42 t M^2+24 M^4)
-9 M^2 s^4 s_1^2  (5 t^2-9 t M^2+3 M^4)
\nonumber\\
&&~~~ -12 M^2 s^5 s_1^2  t - M^2 s_1^2 t^3
 (14 t^3-33 t^2 M^2-22 t M^4+3 M^6)
\nonumber\\
&&~~~ -s_1^2 t^2 (t^2+t s+2 s^2) u_1^3
 +9 M^{10}  s^2 s_1 t +18 M^4 s s_1^2  ( M^4(t^2+s^2)+s^4)\biggr\}
\end{eqnarray}

\begin{eqnarray}
&&\frac{\Delta\sigma(\gamma g\to ^3P_1^{(8)} g\to J/\psi X)}{dt} =
\frac{(ee_cg_s^2)^2\langle{\cal O}_8^{J/\psi}(^3P_1)\rangle}{\pi s^2Q^3 t}
384\biggl\{ -4M^6 s_1^{-1} t^5 t_1 t_2
\nonumber\\
&&~~~ +4 M^4 t^5 (t^2-3 t M^2+6 M^4)+2 M^2 s_1 t^4 t_2 (28 M^4+5 t t_1)
\nonumber\\
&&~~~ +2 M^2 s s_1^2 t^3 (13 t^2+20 t M^2+3 M^4)
+2 M^2 s^2 s_1^2 t^2  (7 t^2+11 t M^2-2 M^4)
\nonumber\\
&&~~~ +12 M^2 s^3 s_1^2 t^2 u+2 M^2 s_1^2 t^4 (11 t^2+16 t M^2+23 M^4)
\nonumber\\
&&~~~ -4 s_1^2 t^2 (t^2+t s+2 s^2) u_1^3
+8 M^4 t^3 (M^4 s_1^2+2 t^3 s_1-M^6 u)\biggr\}
\end{eqnarray}

\begin{eqnarray}
&&\frac{\Delta\sigma(\gamma g\to ^3P_2^{(8)} g\to J/\psi X)}{dt} =
\frac{(ee_cg_s^2)^2\langle{\cal O}_8^{J/\psi}(^3P_2)\rangle}{\pi s^2Q^3 t}
128\biggl\{ 2 M^6 s_1^{-1} t^5 t_1 t2
\nonumber\\
&&~~~ +2 M^4 t^3 (3 t^3 t2+10 t t2 M^4-12 M^6 t_1)
+  M^2  s_1 t^4 (3 t^3+20 t^2 M^2
\nonumber\\
&&~~~ +85 t M^4+52 M^6)+M^2 s s_1^2 t^2 (11 t^3+116 t^2 M^2-55 t M^4+48 M^6)
\nonumber\\
&&~~~ +M^2 s^2 s_1^2 t^2  (t+2 M^2) (13 t+11 M^2)
-2 M^2 s^3 s_1^2 t  (t^2+39 t M^2-48 M^4)
\nonumber\\
&&~~~ -6 s^4 s_1^2 t M^2 (t+14 M^2) -24 s_1^2 M^4 s^5
\nonumber\\
&&~~~ + M^2 s_1^2 t^3
 (5 t^3+72 t^2 M^2+53 t M^4+12 M^6)
\nonumber\\
&&~~~ -2 s_1^2 t^2 (t^2+t s+2 s^2) u_1^3
+12 M^8 s^2 s_1  (2 s u+t M^2)
\nonumber\\
&&~~~ +36 M^6 s_1 (M^4 t^3+s^5-s^4 M^2+M^4 t^2 s)\biggr\}
\end{eqnarray}

\vspace{0.2cm}
The same process has been considered
 in another paper \cite{China}. We note that our
 results for unpolarized and polarized cross sections for
$^3P_0^{(8)}$ octet state production and for the spin dependent cross section
of the $^3P_2^{(8)}$ octet state are different compared to Ref. \cite{China}.
However, our
results for the sums of cross sections for $P$ states,
both polarized and unpolarized,  were confirmed by those obtained in
\cite{China}.