\documentclass[12pt]{iopart}
% Uncomment next line if AMS fonts required
\usepackage{iopams} 
\usepackage{my} 
\usepackage{epsfig} \usepackage{graphics}
 
%
% 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
%
% User commands may be inserted here:
% ----------------------------------
\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
%
\newcommand{\cA}{{\cal A}}
\newcommand{\as}{\alpha_\mathrm{s}}
\newcommand{\asb}{{\bar \alpha}_\mathrm{s}}
\newcommand{\ga}{\gamma}
\newcommand{\om}{\omega}
\newcommand{\order}[1]{{\cal O}\left(#1\right)}
\newcommand{\tq}{{q}_t}
\newcommand{\tqn}{{q}_{t\;n}}
\newcommand{\Pmax}{\bar{q}}
\newcommand{\kt}{k_{t}}
\newcommand{\ktp}{k_{t}^{\prime}}
\newcommand{\SMALLXC}{SMALLXa,SMALLXb}
\newcommand{\CCFM}{CCFMa,CCFMb,CCFMc,CCFMd}
\newcommand{\BFKL}{BFKLa,BFKLb,BFKLc}
\newcommand{\LDC}{LDCa,LDCb,LDCc,LDCd}
\newcommand{\alphasb}{\bar{\alpha}_s}
%\newcommand{\JETSET}{Jetseta,Jetsetb,Jetsetc}
\newcommand{\JETSET}{Jetsetnew}
\newcommand{\LEPTO}{Ingelman_LEPTO65}
\newcommand{\PYTHIAMC}{Jetsetc}
\newcommand{\RAPGAPMC}{RAPGAP206}
\newcommand{\DGLAP}{DGLAPa,DGLAPb,DGLAPc,DGLAPd}
%\newcommand{\SMMOD}{$SMMOD$}
%\newcommand{\CASCADE} {CASCADE}
\def\CASCADE{{\sc Cascade}}
\def\SMALLX{{\sc Smallx}}
\def\RAPGAP{{\sc Rapgap}}
\def\PYTHIA{{\sc Pythia}}
\def\HERWIG{{\sc Herwig}}
% Start of document
% -----------------
\begin{document}
\title[Unintegrated parton densities applied to  
heavy quark production in CCFM ]
{Unintegrated parton densities  
applied to heavy quark production  
in the CCFM approach}
\author{H. Jung 
\footnote[3]{Hannes.Jung@desy.de}
}
\address{   Physics Department, Lund University, 
Box 118, S-221~00 Lund, Sweden }
%
\begin{abstract}
 The application of $k_t$ - factorization supplemented with the
CCFM small-$x$ evolution equation to heavy quark production  
 is discussed.
The $b\bar{b}$
production cross sections at the TEVATRON
can be consistently described using the $k_t$ -
factorization formalism together with the unintegrated gluon density obtained
within the CCFM evolution approach from a fit to HERA $F_2$ data. 
Special attention is drawn to the comparison with measured 
visible cross sections. The visible measured cross sections at 
HERA are compared to the hadron level Monte Carlo generator \CASCADE.
\end{abstract}
%----------------------------------------------------------------------
\section{Introduction}
The calculation of inclusive quantities, like the structure function
$F_2(x,Q^2)$ at HERA, performed in NLO QCD is in perfect agreement with the
measurements.  
However, Catani argues, that the NLO approach,
although phenomenologically successful for $F_2(x,Q^2)$, is
not fully satisfactory from a theoretical viewpoint, because 
{\it ``the truncation of the
splitting functions at a fixed perturbative order is equivalent to assuming that
the dominant dynamical mechanism leading to scaling violations is the evolution
of parton cascades with strongly ordered transverse 
momenta"}~\cite{catani-feb2000}.
As soon as exclusive quantities like jet or heavy quark
production  are investigated, the agreement between NLO coefficient functions
convoluted with NLO DGLAP~\cite{\DGLAP}
 parton densities and the data is not at all
satisfactory: large so-called $K$-factors (normalization factors)
\cite{CDF_bbar,D0_bbar,H1_bbar,ZEUS_bbar} are needed
to bring the NLO calculations close to the data ($K \sim 2-4$ for bottom
production at the TEVATRON), indicating that in the calculations a significant
part of the cross section is still missing.
\par 
At small $x$ the structure function $F_2(x,Q^2)$ is proportional
to the sea quark density, and the sea-quarks are driven via the DGLAP evolution
equations by the gluon density. The standard QCD fits determine the parameters
of the initial parton distributions at a starting scale $Q_0$. With help of the
DGLAP evolution equations these parton distributions are then evolved to any
other scale $Q^2$, with the splitting functions still truncated at fixed
$O(\alpha_s)$ (LO) or $O(\alpha_s ^2)$ (NLO).
Any physics process in the fixed order scheme is then calculated via 
collinear factorization into the coefficient functions
$C^a(\frac{x}{z})$ and collinear (independent of $k_t$) 
parton density functions:
$f_a(z,Q^2)$:
\begin{equation}
\sigma = \sigma_0 \int \frac{dz}{z} C^a(\frac{x}{z}) f_a(z,Q^2)
\label{collinear-factorisation}
\end{equation}
At large energies (small $x$) the evolution of parton densities proceeds over a large
region in rapidity $\Delta y \sim \log(1/x)$ and effects of finite transverse
momenta of the partons may become increasingly important.
Cross sections can then be $k_t$ - factorized~\cite{CCH}
into an off-shell ($k_t$ dependent) partonic cross section
$\hat{\sigma}(\frac{x}{z},k_t) $
and a $k_t$ - unintegrated parton density function 
${\cal F}(z,k_t)$:
\begin{equation}
 \sigma  = \int 
\frac{dz}{z} d^2k_t \hat{\sigma}(\frac{x}{z},k_t) {\cal F}(z,k_t)
\label{kt-factorisation}
\end{equation}
The unintegrated gluon density ${\cal F}(z,k_t)$ is 
described by the BFKL~\cite{\BFKL}
 evolution equation in the region of asymptotically large energies (small $x$). 
 An appropriate description valid for
both small and large $x$ is given by the CCFM evolution
equation~\cite{\CCFM}, resulting in an unintegrated gluon density 
${\cal A} (x,\kt,\Pmax ) $, which is a function also of the 
additional evolution scale $\Pmax $ described below.
\par
In \cite{catani-dis96} Catani argues that
by explicitly carrying out the $k_t$ integration in eq.(\ref{kt-factorisation})
one can obtain a form fully consistent with collinear factorization: the
coefficient functions and also the DGLAP splitting functions leading to 
$f_a(z,Q^2)$ are no longer evaluated in fixed order perturbation theory but
supplemented with the all-order resummation of the $\as \log 1/x$ contribution  
at small $x$. 
\par
In this paper  heavy quark production at the TEVATRON and at HERA is  
investigated using the $k_t$ - factorization approach.
The unintegrated gluon density has been obtained 
previously in \cite{jung_salam_2000} from a
CCFM fit to the HERA structure function $F_2(x,Q^2)$. All free
parameters are thus fixed and absolute predictions for bottom production can be
made. 
As both TEVATRON and HERA data can be described, this shows for the first time
evidence for the universality of the unintegrated CCFM gluon distribution.
\par
First the basic features of the CCFM evolution
equation are recalled and 
the unintegrated gluon density is investigated. Then 
the calculations for $b \bar{b}$ production at the TEVATRON 
is presented as well
as calculations of the visible cross section for $b\bar{b}$
production at HERA.
%======================================================================
\section{The CCFM evolution equation}
\label{sec:CCFMEquation}

\begin{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
% begin picture kinematic variables 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
\begin{center} 
\epsfig{figure=ladder-ccfm_ep.eps,width=6cm,height=6cm}
\end{center} 
\caption{{\it Kinematic variables for multi-gluon emission. The $t$-channel gluon
four-vectors are given by $k_i$ and the gluons emitted in the initial state
cascade have four-vectors $p_i$. The upper angle for any emission is obtained
from the quark box, as indicated with $\Xi$.  
\label{CCFM_variables} }} 
\end{figure}
A solution of the CCFM evolution equation, which properly describes the
inclusive structure function $F_2(x,Q^2)$ and also typical small $x$ final
state processes at HERA has been presented in detail in ~\cite{jung_salam_2000}. 
Figure~\ref{CCFM_variables} shows the pattern of QCD initial-state
radiation in a small-$x$ lepto-production process,  together with labels for the
kinematics.  According to the CCFM evolution equation, the emission of
partons during the initial cascade is only allowed in an
angular-ordered region of phase space. The maximum allowed angle $\Xi$
is defined by the hard scattering quark box, producing the heavy quark pair.
 In terms
of Sudakov variables the quark pair momentum is written as:
\begin{equation}
p_q + p_{\bar{q}} = \Upsilon (p_p + \Xi p_e) + Q_t
\end{equation}
where $p_e$ ($p_p$) are the incoming electron (proton) momenta,
respectively and $Q_t$ is the transverse momentum of the quark pair
in the laboratory frame.
Similarly, the momenta $p_i$ of the gluons emitted during the initial
state cascade are given by (here treated massless):
\begin{equation}
p_i = \upsilon_i (p_p + \xi_i p_e) + p_{ti} \;  , \;\; 
\xi_i=\frac{p_{ti}^2}{s \upsilon_i^2},
\end{equation}
with $\upsilon_i = (1 - z_i) x_{i-1}$, $x_i = z_i x_{i-1}$ and
$s=(p_e+p_p)^2$ being the squared center of mass energy.
The variable $\xi_i$ is connected to the angle of the emitted gluon
with respect to the incoming proton and $x_i$ and $\upsilon_i$ are the
momentum fractions of the exchanged and emitted gluons, while $z_i$ is
the momentum fraction in the branching $(i-1) \to i$ and $p_{ti}$ is
the transverse momentum of the emitted gluon $i$.
\par
The angular-ordered region is then specified by (Fig.~\ref{CCFM_variables}):
\begin{equation}
\xi_0 < \xi_1< \cdots < \xi_n < \Xi
\end{equation}
which becomes:
\begin{equation}
z_{i-1} q_{i-1} < q_{i} 
\end{equation}
where the rescaled transverse momenta $q_{i}$ of the emitted
gluons is defined by:
\begin{equation}
 q_{i} = x_{i-1}\sqrt{s \xi_i} = \frac{p_{ti}}{1-z_i}
\end{equation}
\par
The CCFM equation for the unintegrated gluon density can be
written~\cite{CCFMd,jung_salam_2000,Salam,Martin_Sutton}
as an integral equation:
\begin{eqnarray}
{\cal A} (x,\kt,\Pmax ) = & {\cal A}_0 (x,\kt,\Pmax ) + 
\nonumber \\ & 
 \int \frac{dz }{z} 
\int \frac{d^2 q}{\pi q^{2}} \Theta(\Pmax - zq) \Delta_s(\Pmax ,z q) 
\tilde{P}(z,q,\kt) {\cal A}\left(\frac{x}{z},\ktp,q\right) 
\label{CCFM_integral} 
\end{eqnarray}  
with $\vec{k}'_{t} = | \vec{k}_{t} + (1-z) \vec{q}|$ and $\Pmax$ being
the upper scale for the last angle of the emission: $\Pmax > z_n q_n$,
$q_n > z_{n-1} q_{n-1}$, ..., $q_{1} > Q_0$. Here $q$
is used as a shorthand notation for the 2-dimensional vector of the
rescaled transverse momentum $\vec{q}\equiv\vec{q}_t=\vec{p}_t/(1-z)$.  The
splitting function $\tilde{P}(z,q,\kt)$ 
and the Sudakov form factor $\Delta_s(\Pmax ,z q)$ are given explicitly
in~\cite{jung_salam_2000}.


%----------------------------------------------------------------------
\subsection{The unintegrated gluon density}
In ~\cite{jung_salam_2000} 
the unintegrated gluon density $x {\cal A}(x,k_{t}^2,\Pmax)$ has been
obtained from a fit to the structure function 
$F_2(x,Q^2)$~\footnote[2]{A Fortran program for the unintegrated gluon
density  $x {\cal A}(x,k_{t}^2,\Pmax)$ can be obtained from 
\protect\cite{CASCADEMC}}. 
\par
In Fig.~\ref{gluon} the CCFM unintegrated gluon density distribution 
as a function of  $x$ and  $k_t^2$ is shown and compared to
\begin{equation}
{\cal F}(x,k_t^2) \simeq \left.
\frac{d xG(x,\mu^2)}{d\mu^2}\right|_{\mu^2 = k_t^2}
\end{equation}
with $xG(x,\mu)$ being the collinear gluon density of 
GRV~98~\protect\cite{GRV98} in LO and NLO.  
\begin{figure}[htb] 
  \vspace*{1mm} 
  \begin{center}
\epsfig{figure=gluon_bbarpaper.eps,width=13cm,height=12cm} 
\caption{{\small 
The  {\it CCFM} $k_t$ dependent (unintegrated) gluon density
~\protect\cite{jung_salam_2000,CASCADEMC} at $\Pmax =10$ GeV 
 as a function of
$x$ for different values of $k_t^2$ (upper)
and as a function of 
$k_t^2$ for different values of $x$ (lower) compared to
$\frac{dxG(x,\mu^2)}{d\mu^2}$  with $xG(x,\mu)$ being the
collinear gluon density of 
GRV~98~\protect\cite{GRV98} in LO and NLO. 
 }}\label{gluon} 
 \end{center}
\end{figure} 
\par
The unintegrated gluon density can be related to the integrated one by:
\begin{equation}
\left. xG(x,\mu) \right|_{\mu = \Pmax}
 \simeq \int_0 ^{\Pmax^2} d k_t^2 x {\cal A}(x,k_{t}^2,\Pmax)
\label{ccfm_intglu}
\end{equation}
Here the dependence on the scale of the maximum angle
$\Pmax$  is made explicit: the
evolution proceeds up to the maximum angle $\Pmax$, which plays the role of the
evolution scale in the collinear parton densities. This becomes
 obvious since 
\begin{equation}
\Pmax^2 = x_g^2 \Xi s=
y x_g s 
 = \hat{s} + Q_t^2
\end{equation}
The last expression is derived by using 
$p_Q + p_{\bar{Q}} \simeq x_g p_p + y p_e + Q_t$,
$\Xi  \simeq y/x_g$ and $\hat{s} = y x_g s - Q_t^2$.
This can be compared to a possible choice of 
the renormalization and factorization scale $\mu^2$  
in the collinear approach with $\mu^2 = Q_t^2 + 4 \cdot m_Q^2$ and the
similarity between $\mu$ and $\Pmax$ becomes obvious.
\par
In Fig.~\ref{intglu} the CCFM gluon density integrated over $k_t$
according to eq.(\ref{ccfm_intglu}) is compared to
the gluon densities of GRV~98~\cite{GRV98} in LO and NLO.
\begin{figure}[htb] 
  \vspace*{1mm} 
  \begin{center}
\epsfig{figure=gluonint.eps,width=13cm,height=12cm} 
\caption{{\small 
The  {\it CCFM} gluon density (solid line) integrated over $k_t$
as a function of $x$ for different values of $Q$. For comparison the
GRV~98~\protect\cite{GRV98} gluon density in LO 
(dashed line) and NLO (dotted line) is also shown.
 }}\label{intglu} 
 \end{center}
\end{figure} 
It is interesting to note that the CCFM gluon density is flat for $x\to 0$ at
the input scale $Q=1$ GeV. 
Even at larger scales the collinear gluon densities rise faster with decreasing
$x$ than the CCFM gluon density. However, 
 after evolution and convolution with the off-shell matrix element 
 the scaling violations of $F_2(x,Q^2)$ and the rise of $F_2$
towards small $x$ is reproduced,
 as shown in \cite[Fig.~4 therein]{jung_salam_2000}. 
A similar trend
is observed in the collinear fixed order calculations, when going from LO to
NLO: at NLO the gluon density is less steep at small $x$, because part of the
$x$ dependence is already included in the NLO $P_{qg}$ splitting function, as
argued in \cite{catani-feb2000}.
\par
From Fig.~\ref{intglu} one can see, that the integrated gluon density from CCFM
is larger in the medium $x$ range, than the ones from the collinear approach.
Due to an additional $1/k_t^2$ suppression in the off-shell matrix elements, the
gluon density obviously 
has to be larger to still reproduce the same cross section. In 
addition  only gluon
ladders are considered in the $k_t$ - factorization approach used here,
which means that the sea quark contribution to the
structure function $F_2(x,Q^2)$ comes entirely from boson gluon fusion, without
any contribution from the intrinsic quark sea.
One also should remember, that the relation in eq.(\ref{ccfm_intglu}) is only 
approximately true, since the gluon density itself is not a 
physical observable.
%----------------------------------------------------------------------
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{{\boldmath $b \bar{b}$} production at the TEVATRON}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
The cross section for $b \bar{b}$ production in $p\bar{p}$ collision at
$\sqrt{s}=1800$~GeV is calculated with 
\CASCADE\ ~\cite{jung_salam_2000,CASCADEMC}, which is a Monte Carlo
implementation of the CCFM approach described above.
The off-shell matrix element as given
in~\cite{CCH} for heavy quarks is used with $m_b=4.75$ GeV. The scale
$\mu$ used in
$\alpha_s(\mu^2)$ is set to $\mu^2= m_T^2 = m_b^2 + p_T^2$ (as in
\cite{jung_salam_2000}), with $p_T$ being
the transverse momentum of the heavy quarks in the $p\bar{p}$ 
center-of-mass frame.
In Fig.~\ref{d0_bbar} the prediction for the cross section for $b\bar{b}$
production with pseudo-rapidity
$|y^b| < 1 $ is shown as a function of $p_T^{min}$ and
compared to the measurement of D0~\cite{D0_bbar}. Also shown is the NLO
prediction from \cite[taken from~\cite{D0_bbar}]{Frixione-Mangano}. 
\begin{figure}[htb]
\begin{minipage}{0.45\textwidth}
  \begin{center}
\epsfig{figure=d0bbar_add.epsi,width=7.5cm,height=7.5cm}
\caption{{\it 
Cross section for $b\bar{b}$ production with $|y^b| < 1 $ as a function of
$p_T^{min}$. Shown are the D0~\protect\cite{D0_bbar} 
data points, the fixed
order NLO prediction, and the prediction of \CASCADE.
    }}\label{d0_bbar}
    \end{center}
\end{minipage} 
\hspace*{0.8cm}
\begin{minipage}{0.45\textwidth}
\vskip 0.4cm
\epsfig{figure=cdf_add.epsi,width=7.5cm,height=7.5cm}
\caption{{\it 
Cross section for $b\bar{b}$ production with $|y_b| < 1 $
and $p_T^{min}(b) > 6.5 $~GeV as a function of
$p_T^{min}(\bar{b})$. Shown are the CDF~\protect\cite{CDF_bbar} 
data points, the fixed
order NLO prediction, and the prediction of \CASCADE.
    }}\label{cdf_bbar}
\end{minipage}   
\end{figure}
In Fig.~\ref{cdf_bbar} the measured
cross section of CDF ~\cite{CDF_bbar} is shown 
for 3 values of $p_T^{min}(\bar{b}) $
with the kinematic constraint of  $|y^b|,|y^{\bar{b}}| < 1 $ and 
$p_T^{min}(b) > 6.5 $~GeV  together
with the prediction from \CASCADE\ and the NLO calculation
from~\cite[taken from~\cite{CDF_bbar}]{Mangano-Nason-Ridolfi}.
In all cases the NLO calculation used $m_b=4.75$ GeV and the factorization and
renormalization scales were set to $\mu^2= m_T^2 = m_b^2 + p_T^2$. Both, D0 and
CDF measurements are above the NLO predictions by a factor of $\sim 2$. The
\CASCADE\  predictions are in reasonable agreement with the measurements.
A similarly good description of the D0 and CDF data 
has been obtained in \cite{Hagler_bbar} using also $k_t$ - factorization but
supplemented with a BFKL type unintegrated gluon density.
 It is
interesting to note, that the CCFM unintegrated gluon density has been obtained
from inclusive $F_2(x,Q^2)$ at HERA. In this sense, the prediction of \CASCADE\ 
is a parameter free prediction of the $b\bar{b}$ cross section in $p\bar{p}$
collisions. This also shows for the first time 
evidence for the universality of unintegrated CCFM gluon distribution.
\begin{figure}[htb]
  \vspace*{2mm}
  \begin{center}
\epsfig{figure=kin_bbar_sel.eps,width=14cm,height=8cm}
\caption{{\it 
Comparison of $x_n$ and $k_{tn}$ distributions of the gluons entering the
process $g +g \to b + \bar{b}$. Shown are the predictions from \CASCADE\ 
representing  $k_t$-factorization with the CCFM unintegrated gluon density
and also from \PYTHIA\ representing the collinear approach supplemented with
initial and final state DGLAP parton showers to cover the phase space of $p_t$
ordered QCD cascades. 
    }}\label{kin_bbar}
    \end{center}
\end{figure}
\par
In Fig.~\ref{kin_bbar} the 
$x_n$ and $k_{tn}$ distributions  of the gluons
entering the hard scattering process 
$g +g \to b + \bar{b}$ (see Fig.~\ref{CCFM_variables}) are shown
and the predictions from the $k_t$ - factorization approach (\CASCADE\ ) 
are compared to the
standard collinear approach (here \PYTHIA\ ~\cite{\PYTHIAMC} with LO
  $gg \to  b  \bar{b}$ matrix elements 
supplemented with DGLAP parton showers to simulate higher order 
effects). Whereas the $x_n$ distributions agree reasonably well, a significant
difference is observed in the  $k_{tn}$ distribution. However this is not
surprising: 
the $k_t$ factorization approach includes
a large part of the fixed order NLO corrections (in collinear factorization).
Such corrections are $gg \to Q \bar{Q} g$, where the final state gluon 
can have any kinematically allowed transverse momentum, 
which could be regarded as a first step
toward a non-$p_T$ ordered QCD cascade.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{{\boldmath $b \bar{b}$} production at HERA}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
In \cite{jung_salam_2000} the prediction of \CASCADE\
for the total $b\bar{b}$ cross section was compared to
the extrapolated measurements of the
H1~\cite{H1_bbar} and ZEUS~\cite{ZEUS_bbar} experiments at HERA. 
Since \CASCADE\ generates full hadron level events, 
a direct comparison with measurements can be done,
before extrapolating the measurement over the full phase space 
to the total $b\bar{b}$ 
cross section. ZEUS~\cite{ZEUS_bbar} has measured the dijet 
cross section which can be attributed to bottom production 
by demanding an electron 
inside one of the jets. In the kinematic range of
$Q^2 < 1$~GeV$^2$, $0.2 < y< 0.8$, at least two jets with
$E_T^{jet1(2)}>7(6)$~GeV and $|\eta^{jet}|<2.4$ and a prompt electron with
$p_T^{e^-} > 1.6$~GeV and $|\eta^{e^-}|<1.1$, 
ZEUS~\cite{ZEUS_bbar} quotes the cross section as: 
\begin{equation}
\sigma^{b \to e^-}_{e^+p \to e^+ + dijet +e^- +X } =24.9 \pm 6.4^{+4.2}_{-7.3}
\mbox{ pb }\hspace{1cm}\mbox{(ZEUS~\cite{ZEUS_bbar})}
\end{equation}
with the statistical (first) and systematic (second) error given.
Within the same kinematic region and applying the same jet algorithm 
\CASCADE\  predicts:
\begin{equation}
\sigma^{b \to e^-}_{e^+p \to e^+ + dijet +e^- +X } = 20.3 ^{+1.6}_{-1.9}
\mbox{ pb,}
\end{equation}
where the error reflects the variation of $m_b=4.75\mp 0.25$~GeV. 
This value agrees with the measurement within the
statistical error. The  $b\bar{b}$ cross section, extrapolated from the measured
jet cross section to the region 
 $p^b_t>5$~GeV and $|\eta^b| <2$ using  \CASCADE\, is given by:
\begin{equation}
\sigma(ep \to e' b \bar{b}X)= 1.07 \pm 0.27\; ^{+0.18}_{-0.3}\; \mathrm{nb}
\hspace{1cm}\mbox{(ZEUS using \CASCADE\ ),}
\end{equation}
which can be compared with the 
\CASCADE\  prediction 
$\sigma(ep \to e' b \bar{b}X) = 0.87 \pm 0.08 \;\mathrm{nb}$
using $m_b=4.75 \mp 0.25$~GeV.
The NLO prediction is $\sigma(ep \to e' b \bar{b}X) = 0.64 \;\mathrm{nb}$.
When using \HERWIG\ ~\cite{Herwig} for the extrapolation, ZEUS quotes an
extrapolated cross section of~\cite{ZEUS_bbar}:
$$\sigma
= 1.6\pm0.4(stat.)^{+0.3}_{-0.5}(syst.)^{+0.2}_{-0.4}(ext.) \; \mathrm{nb}$$
The quoted extrapolation uncertainty includes
the extrapolation with \PYTHIA , which is similar to the one obtained with 
\CASCADE .
Thus the extrapolated cross section includes large model
uncertainties and the comparison of extrapolated cross sections 
with predictions from NLO calculations are questionable. This also shows the
advantage of a full hadron level simulation in form of a Monte Carlo program 
over a fixed order parton level prediction. 
\begin{figure}[htb]
\vspace*{2mm}
\begin{center}
\epsfig{figure=zeusbbar_xgamma.eps,width=10cm,height=10cm}
\end{center}
\vskip -1cm
\caption{{\it The differential cross section $d\sigma/dx^{obs}_{\gamma}$
for heavy quark decays as measured by ZEUS~\protect\cite{ZEUS_bbar} and compared
to \CASCADE. Shown are the charm, bottom and the sum of both contributions to
the cross section (Note there is no additional $K$ factor applied).
}}\label{zeus_xgamma}
\end{figure}
\par
It is interesting to note, that the \CASCADE\
results agree well with the \PYTHIA\ results
for the di-jet plus electron cross section, 
if heavy quark excitation is included,
 as well as with the extrapolation factor.
As in the case of $b\bar{b}$ production at the TEVATRON higher order 
QCD effects are important, which are already included in the 
$k_t$ - factorization approach.
\par
In Fig.~\ref{zeus_xgamma} within the same kinematic range the differential 
cross section for heavy quark decays (charm and bottom)
as a function of $x_{\gamma}$ predicted by \CASCADE\  
is compared with the measurement of ZEUS~\cite{ZEUS_bbar}.
A similar trend as in the case of charm photo-production is observed,
namely a significant fraction of the cross section with $x_{\gamma}<1$.
In LO in the collinear factorization
approach this is attributed to resolved photon processes. However, in $k_t$ -
factorization, the $x_{\gamma}$ distribution is explained naturally because 
gluons in the initial state need not to be radiated 
in a $p_t$ ordered region
and therefore can give rise to a high $p_t$ jet with transverse momentum larger
than that of the heavy quarks.
\par
The prediction of \CASCADE\  has also been compared to the measurement of 
H1~\cite{H1_bbar} for electro-production cross section in 
$ Q^2 < 1$ GeV$^2$, $0.1< y < 0.8$, $p_{\perp}^{\mu} < 2$ GeV and $35^o <
\theta^{\mu} < 130^o$:
$$\sigma(ep \to e' b \bar{b}X \to \mu X')
= 0.176\pm0.016(stat.)^{+0.026}_{-0.017}(syst.)\; \mathrm{nb}$$
This cross section already includes the extrapolation from measured jets to
the muon.
In the same kinematic range \CASCADE\ predicts:
\begin{equation*}
\sigma(ep \to e' b \bar{b}X \to \mu X') = 
 0.066^{+0.009}_{-0.007} \;\mathrm{nb}.
\end{equation*} 
which is a factor $\sim 2.6 $ below the measurement. It is interesting to note,
that the ratio 
\begin{equation*}
R(\mbox{H1})=\frac{\sigma(\mbox{\CASCADE })}{\sigma(\mbox{NLO})}= 
\frac{0.066}{0.054} = 1.2 
\end{equation*} 
is similar to the one obtained for the ZEUS extrapolated measurements: 
$R(\mbox{ZEUS}) = 1.4 $. 
The measured cross sections of H1 and ZEUS cannot be
compared directly because different kinematic ranges and decay channels were
used. To
compare both experiments, a ratio $R_{MC}$ is defined:
\begin{eqnarray}
R_{MC}= & \frac{\sigma_{measured}}{\sigma_{MC}}
\end{eqnarray} 
Using \CASCADE , the ratios are:
\begin{eqnarray}
 R_{MC}(\mbox{\small H1}) & = \frac{\sigma_{measured}}
 {\sigma_{ MC}} = 2.7 \pm 0.25^{+0.4}_{-0.26}& \\
 R_{MC}(\mbox{\small ZEUS})& =\frac{\sigma_{measured}}
 {\sigma_{MC}} = 1.2 \pm 0.32 ^{+0.21}_{-0.37}&
\end{eqnarray} 
where the error comes only from $\sigma_{measured}$. This exercise shows
that both experiments differ by a factor 
of more than $2$ in the published measurements, when compared to \CASCADE .

\section{Conclusion}
%   
Bottom production at the TEVATRON can be reasonably well
described using the $k_t$-factorization approach with off-shell matrix elements
for the hard scattering process. One essential ingredient for the satisfactory
description is the unintegrated gluon density, which was obtained by a CCFM
evolution fitted to structure function data at HERA, showing
evidence for the universality of the unintegrated gluon density.
The comparison with the data was performed with the \CASCADE\ Monte
Carlo event generator, which implements $k_t$-factorization together with the
CCFM unintegrated gluon density.
\par
Measurements of bottom production at HERA are also
compared to predictions from \CASCADE. 
The visible dijet plus electron cross section attributed to
$b$-production as measured by ZEUS could be reproduced within the statistical
error. It was pointed out, that the extrapolation from the measured to the total
$b\bar{b}$ cross section contains large model dependencies. If \CASCADE\ or
\PYTHIA\ was used for extrapolation, the cross section was found to agree with
\CASCADE\ and even with NLO calculations within the combined statistical and
systematic error. 
\par
However, the situation is different with the H1 measurement: the visible
muon cross section is already a factor of 2.7 
above the prediction from \CASCADE. 
Comparing both HERA measurements with Monte Carlo predictions of \CASCADE, it
could be shown that both experiments, H1 and ZEUS, differ in their measurements
by a factor of more than two. Further measurements, also differential, are     
desirable to clarify the situation at HERA.
\par
In general the $k_t$-factorization approach has now proven to be successful even
in a kinematic region, where typical small-$x$ effects are expected to be small.
It is the advantage of that approach that important parts of NLO and even NNLO
contributions are consistently included due to the off-shell gluons, which
enter into the hard scattering process.  
%\par
  
\section{Acknowledgments}
Many thanks go  to the organizers of this very nice workshop: G.~Grindhammer,
B.~Kniehl, G.~Kramer and W.~Ochs.
I am grateful to G.~Salam for all his patience and his help in all
different kinds of discussions concerning CCFM and a backward
evolution approach. 
I am grateful to E.~Elsen, L.~Gladilin, G.~Ingelman, L.~J\"onsson,
P.~Schleper for careful reading and comments on the manuscript. 
No words are adequate for everything I have with Antje.
I also want to thank the DESY directorate for
hospitality and support.

\begin{thebibliography}{10}

\bibitem{catani-feb2000}
S.~Catani, Aspects of QCD, from the Tevatron to LHC,  in {\em Proceedings of
  the International Workshop {\em Physics at TeV Colliders}} (Les Houches,
  France, 8-18 June, 1999), \mbox.

\bibitem{DGLAPa}
V.~Gribov, L.~Lipatov, {\em Sov. J. Nucl. Phys.} {\bf 15}\nolinebreak
  [2]\,(1972)\nolinebreak [2]\,438 and 675.

\bibitem{DGLAPb}
L.~Lipatov, {\em Sov. J. Nucl. Phys.} {\bf 20}\nolinebreak
  [2]\,(1975)\nolinebreak [2]\,94.

\bibitem{DGLAPc}
G.~Altarelli, G.~Parisi, {\em Nucl. Phys. {\bf B}} {\bf 126}\nolinebreak
  [2]\,(1977)\nolinebreak [2]\,298.

\bibitem{DGLAPd}
Y.~Dokshitser, {\em Sov. Phys. JETP} {\bf 46}\nolinebreak
  [2]\,(1977)\nolinebreak [2]\,641.

\bibitem{CDF_bbar}
\mbox{CDF} Collaboration; F. Abe~et al., {\em Phys. Rev. {\bf D}} {\bf
  55}\nolinebreak [2]\,(1997)\nolinebreak [2]\,2546.

\bibitem{D0_bbar}
\mbox{D0} Collaboration; B. Abbott~et al., {\em Phys. Lett. {\bf B}} {\bf
  487}\nolinebreak [2]\,(2000)\nolinebreak [2]\,264.

\bibitem{H1_bbar}
\mbox{H1} Collaboration; C. Adloff~et al., {\em Phys. Lett. {\bf B}} {\bf
  467}\nolinebreak [2]\,(1999)\nolinebreak [2]\,156, and erratum ibid.

\bibitem{ZEUS_bbar}
\mbox{ZEUS} Collaboration; J. Breitweg~et al., {\em Eur. Phys. J. {\bf C}} {\bf
  18}\nolinebreak [2]\,(2001)\nolinebreak [2]\,625.

\bibitem{CCH}
S.~Catani, M.~Ciafaloni, F.~Hautmann, {\em Nucl. Phys. {\bf B}} {\bf
  366}\nolinebreak [2]\,(1991)\nolinebreak [2]\,135.

\bibitem{BFKLa}
E.~Kuraev, L.~Lipatov, V.~Fadin, {\em Sov. Phys. JETP} {\bf 44}\nolinebreak
  [2]\,(1976)\nolinebreak [2]\,443.

\bibitem{BFKLb}
E.~Kuraev, L.~Lipatov, V.~Fadin, {\em Sov. Phys. JETP} {\bf 45}\nolinebreak
  [2]\,(1977)\nolinebreak [2]\,199.

\bibitem{BFKLc}
Y.~Balitskii, L.~Lipatov, {\em Sov. J. Nucl. Phys.} {\bf 28}\nolinebreak
  [2]\,(1978)\nolinebreak [2]\,822.

\bibitem{CCFMa}
M.~Ciafaloni, {\em Nucl. Phys. {\bf B}} {\bf 296}\nolinebreak
  [2]\,(1988)\nolinebreak [2]\,49.

\bibitem{CCFMb}
S.~Catani, F.~Fiorani, G.~Marchesini, {\em Phys. Lett. {\bf B}} {\bf
  234}\nolinebreak [2]\,(1990)\nolinebreak [2]\,339.

\bibitem{CCFMc}
S.~Catani, F.~Fiorani, G.~Marchesini, {\em Nucl. Phys. {\bf B}} {\bf
  336}\nolinebreak [2]\,(1990)\nolinebreak [2]\,18.

\bibitem{CCFMd}
G.~Marchesini, {\em Nucl. Phys. {\bf B}} {\bf 445}\nolinebreak
  [2]\,(1995)\nolinebreak [2]\,49.

\bibitem{catani-dis96}
S.~Catani, $k_t$-factorisation and perturbative invariants at small $x$,  in
  {\em Proceedings of the International Workshop on Deep Inelastic Scattering,
  {\em DIS 96}} (Rome, Italy, 15-19 April, 1996), \mbox.

\bibitem{jung_salam_2000}
H.~Jung, G.~Salam, {\em Eur. Phys. J. {\bf C}} {\bf 19}\nolinebreak
  [2]\,(2001)\nolinebreak [2]\,351, \mbox.

\bibitem{Salam}
G.~Bottazzi, G.~Marchesini, G.~Salam, M.~Scorletti, {\em JHEP} {\bf
  12}\nolinebreak [2]\,(1998)\nolinebreak [2]\,011, .

\bibitem{Martin_Sutton}
J.~Kwiecinski, A.~Martin, P.~Sutton, {\em Phys. Rev. {\bf D}} {\bf
  52}\nolinebreak [2]\,(1995)\nolinebreak [2]\,1445.

\bibitem{CASCADEMC}
H.~Jung, {\em The CCFM Monte Carlo generator {\sc Cascade} for lepton - proton
  and proton - proton collisions}, Lund University, 2001,
  \mbox, \mbox{DESY 01-114},
  \verb+http://www.quark.lu.se/~hannes/cascade/+.

\bibitem{GRV98}
M.~Gl\"uck, E.~Reya, A.~Vogt, {\em Eur. Phys. J. {\bf C}} {\bf 5}\nolinebreak
  [2]\,(1998)\nolinebreak [2]\,461, \mbox.

\bibitem{Frixione-Mangano}
S.~Frixione, M.~Mangano, {\em Nucl. Phys. {\bf B}} {\bf 483}\nolinebreak
  [2]\,(1997)\nolinebreak [2]\,321.

\bibitem{Mangano-Nason-Ridolfi}
M.~Mangano, P.~Nason, G.~Ridolfi, {\em Nucl. Phys. {\bf B}} {\bf
  373}\nolinebreak [2]\,(1992)\nolinebreak [2]\,295.

\bibitem{Hagler_bbar}
P.~Hagler et~al., {\em Phys. Rev.} {\bf D62}\nolinebreak
  [2]\,(2000)\nolinebreak [2]\,071502.

\bibitem{Jetsetc}
T.~Sj\"ostrand, {\em Comp. Phys. Comm.} {\bf 82}\nolinebreak
  [2]\,(1994)\nolinebreak [2]\,74.

\bibitem{Herwig}
G.~Marchesini et~al., {\em Comp. Phys. Comm.} {\bf 76}\nolinebreak
  [2]\,(1992)\nolinebreak [2]\,465, .

\end{thebibliography}
 
\end{document} 




