%
% LaTeX file with instructions for contributions to the proceedings of 
% the workshop on 'Monte Carlo Generators for HERA Physics' 1998/99
% ===========================================================================
% Contains instructions on formats, submission (with deadline) and plans 
% for proceedings on www and in print. Run LaTeX on this file to get 
% nicely formatted version. 
% Should be used as a template (delete present text and insert your text) 
% to get the right format. 
% 
% NOTE: 1. Format: DIN A4, with DESY standard margins
%       2. Fontsize: 12pt 
%       3. Reduction 2 A4 pages --> 1 A4 page to be made in printed proceedings
%
% Alternatives:
% 1. LaTeX2e (prefered); usage at DESY with alternatives for including figures:
%    \documentclass[12pt]{article} \usepackage{dina4p} 
%    \usepackage{psfig}  \usepackage{epsfig}  \usepackage{graphics}
% 2. Simple LaTeX with DINA4 page definition at DESY:
%    \documentstyle[12pt,dina4p]{article} 
% 3. Simple LaTex with explicit DINA4 page definition from DESY style files:
     \documentclass[12pt]{article}   
     \usepackage{epsfig}
     \newlength{\dinwidth}                       
     \newlength{\dinmargin}                      
     \setlength{\dinwidth}{21.0cm}               
     \textheight24.2cm \textwidth17.0cm          
     \setlength{\dinmargin}{\dinwidth}           
     \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 \footheight 12pt \footskip 24pt % <--- Use this with LaTeX2e
    \headsep 30pt \footskip 24pt                  % <--- Use this with LaTeX
%
% Use always these additional commands:
% ------------------------------------
\parskip 2mm plus 2mm minus 1mm     % Add space between paragraphs
\renewcommand{\topfraction}{1.0}    % These three commands assure that floats
\renewcommand{\bottomfraction}{1.0} %  (figures, tables) can cover a whole page
\renewcommand{\textfraction}{0.0}   %  and no text is required
%
% Optional useful macros for Journals:
%-- refs according to JHEP
\newcommand\hepph[1]{hep-ph/#1}
\newcommand\hepth[1]{hep-th/#1}
\newcommand\hepex[1]{hep-ex/#1}
\newcommand\adp[3]{{{\it Adv. Phys. }{\bf #1} (#2) #3}}
\newcommand\ap[3]{{{\it Ann. Phys. (NY) }{\bf #1} (#2) #3}}
\newcommand\arnps[3]{{{\it Ann. Rev. Nucl. Part. Sci. }{\bf #1} (#2) #3}}
\newcommand\cpc[3]{{{\it Comput. Phys. Commun.}{\bf #1} (#2) #3}}
\newcommand\cmp[3]{{{\it Comm. Math. Phys. }{\bf #1} (#2) #3}}
\newcommand\jmp[3]{{{\it J. Math. Phys. }{\bf #1} (#2) #3}}
\newcommand\jgp[3]{{{\it J. Geom. Phys. }{\bf #1} (#2) #3}}
\newcommand\cqg[3]{{{\it Class. and Quant. Grav. }{\bf #1} (#2) #3}}
\newcommand\hpa[3]{{{\it Helv. Phys. Acta }{\bf #1} (#2) #3}}
\newcommand\jhep[3]{{{\it J. High Energy Phys. }{\bf #1} (#2) #3}}
\newcommand\npa[3]{{{\it Nucl. Phys. }{\bf A #1} (#2) #3}}
\newcommand\npb[3]{{{\it Nucl. Phys. }{\bf B #1} (#2) #3}}
\newcommand\npps[3]{{{\it Nucl. Phys. }{\bf #1} {\it(Proc. Suppl.)} (#2) #3}}
\newcommand\pla[3]{{{\it Phys. Lett. }{\bf A #1} (#2) #3}}
\newcommand\plb[3]{{{\it Phys. Lett. }{\bf B #1} (#2) #3}}
\newcommand\pr[3]{{{\it Phys. Rev. }{\bf #1} (#2) #3}}
\newcommand\pra[3]{{{\it Phys. Rev. }{\bf A #1} (#2) #3}}
\newcommand\prb[3]{{{\it Phys. Rev. }{\bf B #1} (#2) #3}}
\newcommand\prc[3]{{{\it Phys. Rev. }{\bf C #1} (#2) #3}}
\newcommand\prd[3]{{{\it Phys. Rev. }{\bf D #1} (#2) #3}}
\newcommand\pre[3]{{{\it Phys. Rev. }{\bf E #1} (#2) #3}}
\newcommand\prep[3]{{{\it Phys. Rep. }{\bf #1} (#2) #3}}
\newcommand\prl[3]{{{\it Phys. Rev. Lett. }{\bf #1} (#2) #3}}
\newcommand\ptp[3]{{{\it Prog. Theor. Phys. }{\bf #1} (#2) #3}}
\newcommand\rmp[3]{{{\it Rev. Mod. Phys. }{\bf #1} (#2) #3}}
\newcommand\zpc[3]{{{\it Z. Physik }{\bf C #1} (#2) #3}}
\newcommand\mpla[3]{{{\it Mod. Phys. Lett. }{\bf A #1} (#2) #3}}
\newcommand\mplb[3]{{{\it Mod. Phys. Lett. }{\bf B #1} (#2) #3}}
\newcommand\sjnp[3]{{{\it Sov. J. Nucl. Phys. }{\bf #1} (#2) #3}}
\newcommand\jetp[3]{{{\it Sov. Phys. JETP }{\bf #1} (#2) #3}}
\newcommand\zetf[3]{{{\it Zh. Eksp. Teor. Fiz. }{\bf #1} (#2) #3}}
\newcommand\yf[3]{{{\it Yad. Fiz. }{\bf #1} (#2) #3}}
\newcommand\jetpl[3]{{{\it JETP Lett. }{\bf #1} (#2) #3}}
\newcommand\joth[3]{{{\it J. Operator Theory }{\bf #1} (#2) #3}}
\newcommand\ibid[3]{{\it ibid. }{\bf #1} (#2) #3}
\newcommand\ijmpa[3]{{{\it Int. J. Mod. Phys. }{\bf A #1} (#2) #3}}
\newcommand\ijmpb[3]{{{\it Int. J. Mod. Phys. }{\bf B #1} (#2) #3}}
\newcommand\epjc[3]{{{\it Eur. Phys. J. }{\bf C #1} (#2) #3}}

%\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} A}
%\def\NPB{{\em Nucl. Phys.} B}
%\def\PLB{{\em Phys. Lett.}  B}
%\def\PRL{\em Phys. Rev. Lett.}
%\def\PRD{{\em Phys. Rev.} D}
%\def\ZPC{{\em Z. Phys.} C}
%
% User commands may be inserted here:
% ----------------------------------
\def\ga{\gamma}

\newcommand{\Pom}{I$\!$P}                % gives pomeron symbol
\def\lsim{\mathrel{\rlap{\lower4pt\hbox{\hskip1pt$\sim$}}
    \raise1pt\hbox{$<$}}}                % less than or approx. symbol
\def\gsim{\mathrel{\rlap{\lower4pt\hbox{\hskip1pt$\sim$}}
    \raise1pt\hbox{$>$}}}                % greater than or approx. symbol
%
\renewcommand{\thefootnote}{\fnsymbol{footnote}}
% Start of document
% -----------------
\begin{document}
\begin{flushright}
  Bicocca--FT--99--14\\
  \\
  May 1998
\end{flushright}

\vspace*{10mm}
\begin{center}  \begin{Large} \begin{bf}
Small-$x$ final states and the CCFM equation\footnote{Work supported in
  part by the E.U. QCDNET 
 contract FMRX-CT98-0194}\footnote{Contribution to workshop on Monte
 Carlo generators for HERA physics.}\\ 
  \end{bf}  \end{Large}
  \vspace*{5mm}
  \begin{large}
G.P. Salam$^a$\\
  \end{large}
\end{center}
$^a$ INFN --- Sezione di Milano, and Universit\`a di Milano--Bicocca,
Milano, Italy\\
%
\begin{quotation}
\noindent
{\bf Abstract:} The status of the Milan group's work on CCFM-based
phenomenology at small $x$, and possible directions of future
investigation, are discussed.
\end{quotation}
%

In order to study final states in QCD it is generally necessary to
take into account the coherence of QCD radiation \cite{coherence}. At
small $x$, coherence is embodied in the CCFM equation \cite{CCFM}.  In
\cite{BMSS} a first analysis of the CCFM equation was carried out with
regard to small-$x$ phenomenology at HERA. The information gained was
to have served as a basis for a CCFM-based Monte Carlo event
generator.  The main results are summarised here, together with some
discussion of possible future work.


As a first step, we studied the freedom available in implementing the
CCFM equation, and noted that the non-Sudakov form factor (which
resums virtual corrections at small $x$) is not uniquely defined.
There are two possible choices for it. One leads to small NLL
corrections compared to BFKL \cite{BFKL,NLL}, the other to large NLL
corrections (much bigger that the true NLL corrections to BFKL).

Soft emissions (those from the $1/(1-z)$ part of the splitting
function) were neglected throughout the analysis.  Formally this
should have been acceptable for the observables we considered. Their
inclusion would have been another source of NLL corrections.

As a phenomenological constraint on the free parameters of our
evolution we fitted the small-$x$ $F_2$ structure function \cite{F2}.
It was found that of the two possible implementations of the
non-Sudakov form factor, only the one with the large corrections was
able to reproduce the $F_2$ data. We then examined other, more
exclusive, observables such as the transverse energy flow \cite{Et},
charged-particle transverse-momentum spectra \cite{pt} and the
forward-jet cross section \cite{fjZS,fjH1}. No hadronisation was
implemented. Nor did we include the contribution to the final state
from the quark box.

All final-state observables were systematically low (especially the
$E_t$ flow and the forward-jet cross section, which were off by a
factor of $2$). In the case of the $E_t$ flow this might be
explainable by the absence of the contributions from soft emissions
and hadronisation corrections.  Less so for the other observables. As
an example, figure~\ref{fig:fj} shows our results for the forward-jet
cross section, together with an indication of the uncertainties, as
compared to the data.

\begin{figure}[htbp]
  \begin{center}
    \epsfig{file=fj-ccfm.eps,width=0.8\textwidth}
    \caption{Forward-jet cross section from \cite{BMSS} as compared to
      ZEUS data \cite{fjZS}. The differences between the three
      histograms are indicative of the uncertainty on the
      theoretical result.\vspace{-0.5cm}}
    \label{fig:fj}
  \end{center}
\end{figure}

What is our understanding as to why we failed to reproduce the data?
We suspect that the choice of implementation of the non-Sudakov form
factor is largely to blame. We needed a `source of subleading
corrections' in order to be able to fit $F_2$. This was given to us by
our modification of the non-Sudakov form factor. But we subsequently
discovered that this modification leads to strange behaviours, such as
small-$x$ cross sections (e.g.\ for forward jets) which first decrease
as $x$ decreases, and only at very small $x$ start to increase: our
source of subleading corrections was too large and had the wrong
characteristics.

Does this mean that the CCFM equation is inconsistent with the data?
Not necessarily.  It should be remembered that we did not include the
soft emissions. It is quite possible that they alone would have been a
sufficient source of subleading corrections, in order to allow a fit
to the structure function data, without requiring a modification of
the non-Sudakov form factor (the soft emissions lead to somewhat
smaller, but still sizeable NLL corrections). They would also have
increased activity in the final state. Such an approach might
therefore have led to much better agreement with the data.  (This
would have been close in physics content to the SMALLX \cite{SMALLX}
program, which is currently being studied by Hannes Jung \cite{Jung},
who finds that it does actually reproduce a range of data relatively
well).

We have not however taken the seemingly logical step of implementing
the soft emissions within our approach. Why? Two practical issues are
involved. One is undoubtedly the difficulty involved --- the inclusion
of soft emissions complicates enormously the numerical solution of the
CCFM equation.\footnote{It should be pointed out that the recent
  results on the LL equivalence of CCFM (including soft emissions) and
  BFKL final states \cite{gps99} is not sufficient to allow one to use
  the BFKL equation as the basis of a Monte Carlo --- indeed as
  discussed in \cite{dis99}, at subleading orders, in DIS, the BFKL
  equation has serious pathologies in its final states.}
Additionally, as mentioned above, there is already a program in
existence, SMALLX, which implements CCFM including soft emissions.
There, the difficulty of implementing soft emissions is traded for
difficulties of use (it is a forward-evolution Monte Carlo generating
weighted events, with considerable fluctuations in the weights), and
restricted flexibility in certain areas such as the inclusion of the
appropriate scale for the running coupling.


But there are also physics issues of relevance. In the past year the
NLL BFKL corrections have become available \cite{NLL} and there has
been significant progress in understanding the physical origins of the
numerically dominant parts of the corrections \cite{NLLresum}. The
CCFM equation (even with soft emissions and running coupling) is
missing two of the main physical parts of the NLL corrections. There
is a contribution related to the finite part of the gluon splitting
function (while the $z\to0$ and $z\to1$ singular parts are included by
the CCFM equation, the part which is finite at both $z=0$ and $z=1$ is
not). This might be quite straightforwardly included in the CCFM
equation. But there is also a component of the NLL corrections which
ensures the symmetry of the evolution --- the requirement that
evolution from a small (transverse) scale to a large scale give the
same answer as evolution from a large scale to a small scale.  The
CCFM equation does not satisfy this symmetry, corresponding
technically to the absence of the component of the NLL corrections
which goes as $1/(1-\ga)^3$ ($\ga$ is the Mellin transform variable
conjugate to the transverse momentum). Nor is it immediately clear how
to introduce this symmetry without modifying the whole structure of
the CCFM equation. (Such a modification has been proposed, in the form
of the LDC model \cite{LDC}, but it changes the small-$x$ leading
logarithms).


Given that without these contributions, in terms of the physics
included, it would be difficult to improve significantly over the
already-existing SMALLX program, we have decided to place on hold our
plans for the construction of a small-$x$ Monte Carlo event generator.
Such a project is of considerable importance. But an important (and
ongoing) part of it is to understand how to include the various
physical components of the NLL corrections outlined above, into an
equation such as the CCFM equation which is suitable for final-state
studies.

\section*{Acknowledgements}
The work on implementing a CCFM description of final states was done
in collaboration with Giulio Bottazzi, Giuseppe Marchesini and Massimo
Scorletti. I would like also to thank Yuri Dokshitzer, Marcello
Ciafaloni and Dimitri Colferai for discussions.

\begin{thebibliography}{99}
\bibitem{coherence}
       A.H. Mueller, \plb{104}{1981}{161}; 
       B.I. Ermolaev and V.S. Fadin, \jetpl{33}{1981}{285}; 
       Yu.L.\ Dokshitzer, V.S. Fadin and V.A. Khoze, \zpc{15}{1982}{325}; 
       A.\ Bassetto, M.\ Ciafaloni and G.\ Marchesini and
       A.H. Mueller, \npb{207}{1982}{189}; 
       A.\ Bassetto, M.\ Ciafaloni and G.\ Marchesini,
       \prep{100}{1983}{201}; 
       Yu.L.\ Dokshitzer, V.A.\ Khoze, S.I.\ Troyan and A.H.\ Mueller,
       {\it Basics of Perturbative QCD}, Editions Frontieres, Paris,
       France 1991; 
       Yu.L. Dokshitzer, L.V. Gribov, V.A. Khoze and
       S.I. Troyan, \plb{202}{1988}{276};  
       L.V. Gribov, Yu.L. Dokshitzer, V.A. Khoze and S.I. Troyan,
       \jetp{67}{1988}{1303}; 
       M. Ciafaloni, A. Bassetto, G. Marchesini, \prep{100}{1983}{201}.
\bibitem{CCFM} M. Ciafaloni, \npb{296}{1987}{249};
  S. Catani, F. Fiorani and G. Marchesini, \plb{234}{1990}{339}; S.
  Catani, F. Fiorani and G. Marchesini, \npb{336}{1990}{18}.
  
\bibitem{BMSS} G. Bottazzi, G. Marchesini, G.P. Salam and M.
  Scorletti, \npb{505}{1997}{366} [\hepph{9702418}]; G. Bottazzi, G.
  Marchesini, G.P. Salam, M. Scorletti, \jhep{9812}{1998}{011},
  [\hepph{9810546}].
\bibitem{BFKL} L.N. Lipatov, Sov. J. Phys. 23 (1976) 338; E.A. Kuraev,
  L.N. Lipatov and V.S. Fadin, Sov. Phys. JETP 45 (1977) 199; Ya.
  Balitskii and L.N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822.
\bibitem{NLL} M. Ciafaloni and G. Camici, \plb{430}{1998}{127}
  [\hepph{9803389}];
  V.S. Fadin and L.N. Lipatov, \plb{429}{1998}{127} [\hepph{9802290}]; 
  and references therein.

\bibitem{F2} H1 Collaboration (S. Aid et al.), \npb{470}{1996}{3}
  [\hepex{9603004}]; 
 ZEUS Collaboration (M. Derrick et al.),
  \zpc{72}{1996}{399} [\hepex{9607002}].
\bibitem{Et} H1 collaboration, contrib.\ paper pa02-073 to ICHEP'96,
  Warsaw 1996.
\bibitem{pt} H1 collaboration (C. Adloff et al.),
  \npb{485}{1997}{3} [\hepex{9610006}].
\bibitem{fjZS} ZEUS Collaboration (J. Breitweg  et al.), \epjc{6}{239}{1999}
  [\hepex{9805016}].
\bibitem{fjH1}  H1 Collaboration (C. Adloff et al.),
  \npb{538}{3}{1999} [\hepex{9809028}]. 

\bibitem{SMALLX}  G. Marchesini and B.R. Webber, \npb{386}{1992}{215}.

\bibitem{Jung} H. Jung, these proceedings.

\bibitem{gps99} G.P. Salam,
  \jhep{9903}{1999}{009} [\hepph{9902324}]. 
  
\bibitem{dis99} G.P. Salam, talk given at the 7th International Workshop
  on Deep Inelastic Scattering and QCD (DIS 99), Zeuthen, Germany,
  1999, \hepph{9905308}.


\bibitem{NLLresum}  G.P. Salam, \jhep{9807}{1998}{019}
  [\hepph{9806482}]; M. Ciafaloni and D. Colferai,
  \plb{452}{1999}{372} [\hepph{9812366}]; M. Ciafaloni, D. Colferai
  and G.P. Salam, \hepph{9905566}.

\bibitem{LDC} B. Andersson, G. Gustafson, J. Samuelsson,
  \npb{467}{1996}{443}; H. Kharraziha and L. L\"onnblad,
  \jhep{9803}{1998}{006} [\hepph{9709424}].

\end{thebibliography}
\end{document}

