%================================================================
% LaTeX file with prefered layout for H1 paper drafts
% use: dvips -D600 file-name
%================================================================
%%%%%%%%%%%%% LATEX HEADER
%%%%%%%%%%%%% DO NOT DELET NOT CHANGE ANYTHING IN THE HEADER
%%%%%%%%%%%%% UNLESS CLEARLY RECOMMENDED
\documentclass[12pt]{article}
\usepackage{epsfig}
\usepackage{amsmath}
\usepackage{hhline}
\usepackage{amssymb}
\usepackage{times}
\usepackage{cite}
%%%%%%%%%%%% Comment the next two lines to remove the line numbering
%\usepackage[]{lineno}
%\linenumbers
%%%%%%%%%%%%
\renewcommand{\topfraction}{1.0}
\renewcommand{\bottomfraction}{1.0}
\renewcommand{\textfraction}{0.0}
\newlength{\dinwidth}
\newlength{\dinmargin}
\setlength{\dinwidth}{21.0cm}
\textheight23.5cm \textwidth16.0cm
\setlength{\dinmargin}{\dinwidth}
\setlength{\unitlength}{1mm}
\addtolength{\dinmargin}{-\textwidth}
\setlength{\dinmargin}{0.5\dinmargin}
\oddsidemargin -1.0in
\addtolength{\oddsidemargin}{\dinmargin}
\setlength{\evensidemargin}{\oddsidemargin}
\setlength{\marginparwidth}{0.9\dinmargin}
\marginparsep 8pt \marginparpush 5pt
\topmargin -42pt
\headheight 12pt
\headsep 30pt \footskip 24pt
\parskip 3mm plus 2mm minus 2mm
%%%%%%%%%%%%%%%% END OF LATEX HEADER
%===============================title page=============================
\begin{document}  
%%%%%%%%%%%%%%%% Pre-defined commands, you can use for the most obvious notations
\newcommand{\pom}{{I\!\!P}}
\newcommand{\reg}{{I\!\!R}}
\newcommand{\slowpi}{\pi_{\mathit{slow}}}
%\newcommand{\gevsq}{\mathrm{GeV}^2}
\newcommand{\fiidiii}{F_2^{D(3)}}
\newcommand{\fiidiiiarg}{\fiidiii\,(\beta,\,Q^2,\,x)}
\newcommand{\n}{1.19\pm 0.06 (stat.) \pm0.07 (syst.)}
\newcommand{\nz}{1.30\pm 0.08 (stat.)^{+0.08}_{-0.14} (syst.)}
\newcommand{\fiidiiiful}{F_2^{D(4)}\,(\beta,\,Q^2,\,x,\,t)}
\newcommand{\fiipom}{\tilde F_2^D}
\newcommand{\ALPHA}{1.10\pm0.03 (stat.) \pm0.04 (syst.)}
\newcommand{\ALPHAZ}{1.15\pm0.04 (stat.)^{+0.04}_{-0.07} (syst.)}
\newcommand{\fiipomarg}{\fiipom\,(\beta,\,Q^2)}
\newcommand{\pomflux}{f_{\pom / p}}
\newcommand{\nxpom}{1.19\pm 0.06 (stat.) \pm0.07 (syst.)}
\newcommand {\gapprox}
   {\raisebox{-0.7ex}{$\stackrel {\textstyle>}{\sim}$}}
\newcommand {\lapprox}
   {\raisebox{-0.7ex}{$\stackrel {\textstyle<}{\sim}$}}
\def\gsim{\,\lower.25ex\hbox{$\scriptstyle\sim$}\kern-1.30ex%
\raise 0.55ex\hbox{$\scriptstyle >$}\,}
\def\lsim{\,\lower.25ex\hbox{$\scriptstyle\sim$}\kern-1.30ex%
\raise 0.55ex\hbox{$\scriptstyle <$}\,}
\newcommand{\pomfluxarg}{f_{\pom / p}\,(x_\pom)}
\newcommand{\dsf}{\mbox{$F_2^{D(3)}$}}
\newcommand{\dsfva}{\mbox{$F_2^{D(3)}(\beta,Q^2,x_{I\!\!P})$}}
\newcommand{\dsfvb}{\mbox{$F_2^{D(3)}(\beta,Q^2,x)$}}
\newcommand{\dsfpom}{$F_2^{I\!\!P}$}
\newcommand{\gap}{\stackrel{>}{\sim}}
\newcommand{\lap}{\stackrel{<}{\sim}}
\newcommand{\fem}{$F_2^{em}$}
\newcommand{\tsnmp}{$\tilde{\sigma}_{NC}(e^{\mp})$}
\newcommand{\tsnm}{$\tilde{\sigma}_{NC}(e^-)$}
\newcommand{\tsnp}{$\tilde{\sigma}_{NC}(e^+)$}
\newcommand{\st}{$\star$}
\newcommand{\sst}{$\star \star$}
\newcommand{\ssst}{$\star \star \star$}
\newcommand{\sssst}{$\star \star \star \star$}
\newcommand{\tw}{\theta_W}
\newcommand{\sw}{\sin{\theta_W}}
\newcommand{\cw}{\cos{\theta_W}}
\newcommand{\sww}{\sin^2{\theta_W}}
\newcommand{\cww}{\cos^2{\theta_W}}
\newcommand{\trm}{m_{\perp}}
\newcommand{\trp}{p_{\perp}}
\newcommand{\trmm}{m_{\perp}^2}
\newcommand{\trpp}{p_{\perp}^2}
\newcommand{\alp}{\alpha_s}

\newcommand{\alps}{\alpha_s}
\newcommand{\sqrts}{$\sqrt{s}$}
\newcommand{\LO}{$O(\alpha_s^0)$}
\newcommand{\Oa}{$O(\alpha_s)$}
\newcommand{\Oaa}{$O(\alpha_s^2)$}
\newcommand{\PT}{p_{\perp}}
\newcommand{\JPSI}{J/\psi}
\newcommand{\sh}{\hat{s}}
%\newcommand{\th}{\hat{t}}
\newcommand{\uh}{\hat{u}}
\newcommand{\MP}{m_{J/\psi}}
%\newcommand{\PO}{\mbox{l}\!\mbox{P}}
\newcommand{\PO}{I\!\!P}
\newcommand{\xbj}{x}
\newcommand{\xpom}{x_{\PO}}
\newcommand{\ttbs}{\char'134}
\newcommand{\xpomlo}{3\times10^{-4}}  
\newcommand{\xpomup}{0.05}  
\newcommand{\dgr}{^\circ}
\newcommand{\pbarnt}{\,\mbox{{\rm pb$^{-1}$}}}
\newcommand{\gev}{\,\mbox{GeV}}
\newcommand{\WBoson}{\mbox{$W$}}
\newcommand{\fbarn}{\,\mbox{{\rm fb}}}
\newcommand{\fbarnt}{\,\mbox{{\rm fb$^{-1}$}}}
\newcommand{\dsdx}[1]{$d\sigma\!/\!d #1\,$}
\newcommand{\eV}{\mbox{e\hspace{-0.08em}V}}
%
% Some useful tex commands
%
\newcommand{\qsq}{\ensuremath{Q^2} }
\newcommand{\gevsq}{\ensuremath{\mathrm{GeV}^2} }
\newcommand{\et}{\ensuremath{E_t^*} }
\newcommand{\rap}{\ensuremath{\eta^*} }
\newcommand{\gp}{\ensuremath{\gamma^*}p }
\newcommand{\dsiget}{\ensuremath{{\rm d}\sigma_{ep}/{\rm d}E_t^*} }
\newcommand{\dsigrap}{\ensuremath{{\rm d}\sigma_{ep}/{\rm d}\eta^*} }

%%% Dstar stuff
\newcommand{\dstar}{\ensuremath{D^*}}
\newcommand{\dstarp}{\ensuremath{D^{*+}}}
\newcommand{\dstarm}{\ensuremath{D^{*-}}}
\newcommand{\dstarpm}{\ensuremath{D^{*\pm}}}
\newcommand{\zDs}{\ensuremath{z(\dstar )}}
\newcommand{\Wgp}{\ensuremath{W_{\gamma p}}}
\newcommand{\ptds}{\ensuremath{p_t(\dstar )}}
\newcommand{\etads}{\ensuremath{\eta(\dstar )}}
\newcommand{\ptj}{\ensuremath{p_t(\mbox{jet})}}
\newcommand{\ptjn}[1]{\ensuremath{p_t(\mbox{jet$_{#1}$})}}
\newcommand{\etaj}{\ensuremath{\eta(\mbox{jet})}}
\newcommand{\detadsj}{\ensuremath{\eta(\dstar )\, \mbox{-}\, \etaj}}

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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Our macros
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\newcommand{\m}{\ensuremath{\mathrm{m}}}
\newcommand{\cm}{\ensuremath{\mathrm{cm}}}
\newcommand{\um}{\ensuremath{\mathrm{\mu m}}}
\newcommand{\us}{\ensuremath{\mathrm{\mu s}}}
\newcommand{\pb}{\ensuremath{\mathrm{pb}}}
\newcommand{\GeV}{\ensuremath{\mathrm{GeV}}}
\newcommand{\MeV}{\ensuremath{\mathrm{MeV}}}
\newcommand{\der}{\ensuremath{\mathrm{d}}}
\newcommand{\ns}{\ensuremath{\mathrm{ns}}}
\newcommand{\xxx}{\ensuremath{\mathrm{xxx}}}
\newcommand{\rem}[1]{{\bfseries \itshape #1}}
\newcommand{\sub}[1]{\ensuremath{_\mathrm{#1}}}
\newcommand{\super}[1]{\ensuremath{^\mathrm{#1}}}
\newcommand{\pt}{p\sub T}
\newcommand{\ptrel}{p\sub T \super{rel}}
\newcommand{\pz}{p\sub z}
\newcommand{\ptmu}{p\sub T ^\mu}
\newcommand{\ptjet}{p\sub T \super{jet}}
\newcommand{\ptjets}{p\sub T \super{jet 1 (2)}}
\newcommand{\ptjetone}{p\sub T \super{jet 1}}
\newcommand{\ptjettwo}{p\sub T \super{jet 2}}
\newcommand{\etamu}{\eta^\mu}
\newcommand{\etajet}{\eta\super{jet}}
\newcommand{\etajets}{\eta\super{jet 1 (2)}}
\newcommand{\dEdx}{\der E / \der x}
\newcommand{\mK}{m\sub {K}}
\newcommand{\mKK}{m\sub {KK}}
\newcommand{\mphi}{m\sub {\phi}}
\newcommand{\nS}{n\sub {S}}
\newcommand{\Nbin}{N\sub {bin}}
\newcommand{\Nb}{N\sub {b}}
\newcommand{\fb}{f\sub {b}}
\newcommand{\effrec}{\epsilon\sub {rec}}
\newcommand{\efftrig}{\epsilon\sub {trig}}
\newcommand{\sigmaep}{\sigma\sub{ep}}
\newcommand{\sigmagammap}{\sigma\sub{\gamma p}}
\newcommand{\sigmagammapT}{\sigma\sub{\gamma p}\super{T}}
\newcommand{\sigmagammapL}{\sigma\sub{\gamma p}\super{L}}
\newcommand{\phiL}{\phi\sub{\gamma/e}\super{L}}
\newcommand{\phiT}{\phi\sub{\gamma/e}\super{T}}
\newcommand{\Phigammae}{\Phi\sub{\gamma/e}}
\newcommand{\Lumi}{\ensuremath{\mathcal{L}}}
\newcommand{\Lint}{\ensuremath{\int {\mathcal{L} \der t}}}
\newcommand{\MY}{\ensuremath{M\sub{Y}}}
\newcommand{\Ntag}{\ensuremath{N\sub{tag}}}
\newcommand{\Nuntag}{\ensuremath{N\sub{untag}}}
\newcommand{\ftag}{\ensuremath{f\sub{tag}}}
\newcommand{\epsel}{\ensuremath{\epsilon\sub{el}}}
\newcommand{\epspd}{\ensuremath{\epsilon\sub{pd}}}
\newcommand{\Nel}{\ensuremath{N\sub{el}}}
\newcommand{\Npd}{\ensuremath{N\sub{pd}}}
%\newcommand{\artitle}[1]{{\itshape``#1''}}
\newcommand{\artitle}[1]{}
\newcommand{\ysigma}{\ensuremath{y \sub \Sigma}}
\newcommand{\xgamma}{\ensuremath{x \sub \gamma \super{obs}}}
\newcommand{\deltaphi}{\ensuremath{\delta\phi\sub{jets}}}
\newcommand{\mur}{\ensuremath{\mu\sub{r}}}
\newcommand{\muf}{\ensuremath{\mu\sub{f}}}

\begin{titlepage}

% \noindent
% Date:          \today      \\
% Version:       0.95      \\
% Editors:       Mira Kr\"amer $<$Mira.Kraemer@desy.de$>$, Benno List
% $<$Benno.List@desy.de$>$    \\
% Referees:  Daniel Pitzl $<$Daniel.Pitzl@desy.de$>$, Wulfrin Bartel  $<$Wulfrin.Bartel@desy.de$>$       \\
% Comments by: Monday, July 21, 2008.
% 
\noindent
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%% For conference papers  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%% coment the header and fill the right conference
\begin{center} %%% you may want to use this line for working versions
\begin{small}
\begin{tabular}{llrr}
Submitted to & & &
\epsfig{file=H1logo_bw_small.epsi
,width=2.cm} \\[.2em] \hline
\multicolumn{3}{l}{{\bf
                ICHEP08},
                Philadelphia, PA} 
                 & Abstract:   {\bf 852}    \\
%                 & Parallel Session & {\bf x}   &\\ 
                 \hline
 & \multicolumn{3}{r}{\footnotesize {\it Electronic Access:
    www-h1.desy.de/h1/www/publications/conf/conf\_list.html}} \\[.2em]
\end{tabular}
\end{small}
\end{center}


\vspace{2cm}
\begin{center}
\begin{Large}

{\bf A Measurement of Beauty Photoproduction
Through Decays to Muons and Jets at HERA-II }

\vspace{2cm}

H1 Collaboration

\end{Large}
\end{center}

\vspace{2cm}

\begin{abstract}
The photoproduction of beauty quarks in ep collisions 
has been measured using a data sample of $170\,\pb^{-1}$ collected
with the H1 detector at HERA-II in the years 2006 and 2007.
% Beauty dijet events were identified by the presence of 
% muons in the final state with large transverse momentum to
% a jet or a large impact parameter.
This measurement follows closely a measurement performed at HERA-I,
 where beauty photoproduction events were investigated with
two jets and a muon in the final state, and
beauty events were identified using the muons relative transverse momentum to 
a jet and its impact parameter.
Visible cross sections were measured differentially in the transverse
momenta of the highest energy jet ($\ptjetone$) and the muon ($\ptmu$), 
the pseudorapidity of the muon ($\etamu$) and of the photon's momentum fraction
$\xgamma$ entering the hard interaction.
% Differential cross sections have been measured as function of 
% the transverse momentum of the leading jet and the muon, the muon's
% pseudorapidity, the photon's momentum fraction $\xgamma$ that enters the hard
% interaction, and the azimuthal angle difference between the two leading jets.
The measurements are found to be well described by QCD calculations at NLO.
\end{abstract}

\vspace{1.5cm}
\vspace{1.5cm}

% \begin{center}
% To be submitted to \EJC \;\; or \PLB
% \end{center}

\end{titlepage}

%          THE PAPER DRAFTS HAVE NO AUTHORLIST
%
%          FOR PAPER ISSUED FOR THE FINAL READING 
%          COPY THE AUTHOR AND INSTITUTE LISTS 
%          INTO YOUR AREA
%
% from /h1/iww/ipublications/h1auts.tex 
%          AND UNCOMMENT THE NEXT THREE LINES 
%
%\begin{flushleft}
%  \input{h1auts}
%\end{flushleft}
%
% Please not that the author list may need re-formatting.

\newpage

\section{Introduction}

The production of beauty quarks in $ep$ collisions has been investigated
in considerable detail at HERA-I. In several analyses it was observed that
measured cross sections where significantly above the predictions from 
 perturbative QCD calculations in
next-to-leading order (NLO). The data from HERA-II with its larger
statistics makes it possible to repeat these measurements with increased
accuracy. 

This measurement follows closely a measurement performed with data from HERA-I
\cite{Aktas:2005zc}, where beauty photoproduction events were investigated with
two jets and a muon in the final state.
Visible cross sections were measured differentially in the transverse
momenta of the highest energy jet ($\ptjetone$) and the muon ($\ptmu$), 
the pseudorapidity of the muon ($\etamu$) and of the photon's momentum fraction
$\xgamma$ entering the hard interaction.

The result of the HERA-I measurement was that NLO calculations 
describe the data reasonably well, except
for the lowest bin of the muon and jet transverse momentum, $\ptmu$ and
$\ptjetone$, where the data were significantly above the predictions.
A similar measurement has also been 
made by the ZEUS collaboration \cite{Chekanov:2003si}, in a similar but
not identical phase space. Here, good agreement was found with QCD
predictions, also at low jet and muon $\pt$.

The measurement presented here uses the HERA-II data set to measure the
same cross sections as in the previous H1 publication with increased
statistics and correspondingly smaller errors.
In addition, the differential cross section as a function
of the azimuthal angle difference $\deltaphi$ between the two leading jets
has also been measured.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Data Sample}

The data was collected with the H1 detector in the years 2006 and 2007,
when HERA collided electron and positron beams
of an energy of $E\sub e=27.55\,\GeV$ with protons of 
$E\sub p=920\,\GeV$.

A detailed description of the H1 detector can be found elsewhere
\cite{H1detector}. The main components utilized in this analysis were
the Central Jet Chamber (CJC), the Liquid Argon Calorimter (LAr), the rear Spaghetti
Calorimeter (SpaCal) \cite{Appuhn:1996na}, the Central Muon Detector (CMD),
the Central Silicon Tracker (CST) \cite{bib:cst} and the Fast Track Trigger
(FTT) \cite{bib:ftt}. 

The total integrated luminosity corresponds to $170\,\pb^{-1}$,
of which $49\,\pb^{-1}$ were collected in the year 2006 running with electrons,
and $121\,\pb^{-1}$ were collected in the years 2006/07 running with positrons.

The events were triggered on the first trigger Level by demanding
a track segment in the muon system and track activity in the central jet chamber.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{QCD Models}

The Monte Carlo generators PYTHIA 6.2 \cite{bib:pythia} 
and Cascade 2.0 \cite{bib:pythia}
were used for the simulation of signal and background distributions.
Event samples from both generators were simulated with a detailed 
detector simulation based on GEANT \cite{bib:geant} and subjected to the same
reconstruction program as real data.

The PYTHIA event samples were generated in an inclusive production mode
({\tt MSTP(14)= 30}), where direct and resolved processes are calculated
using massless matrix elements. The CTEQ6L \cite{Pumplin:2002vw} and
SAS-1D \cite{Schuler:1996fc}
parton density sets were used for the proton and the photon,
respectively.
The fragmentation of heavy quarks to hadrons was simulated using the
fragmentation function by Peterson {\it et al.} \cite{Peterson:1982ak}
with a parameter $\epsilon\sub{b (c)} = 0.0069 (0.058)$ for beauty
(charm) quarks. More details on the parameter settings can be found in
our previous publication \cite{Aktas:2005zc}.

For cross checks, additional  Monte Carlo samples were generated using
the Cascade program, which in contrast to PYTHIA is based on $k\sub T$
factorization and the CCFM evolution rather than collinear factorization
and DGLAP evolution. The proton parton density set A0 \cite{Jung:2004gs} is used for the
unintegrated gluon density in the proton. 

Both programs, PYTHIA and Cascade, use Leading Order (LO), i.e.
$O (\alpha\sub s)$, QCD matrix elements for the hard scattering,
augmented by parton showers that approximate the effect of additional 
multiple gluon emission. 
 
To compare the data to Next-to-Leading Order (NLO) QCD, i.e.
$O (\alpha\sub s^2)$, the FMNR program \cite{Frixione:1995qc} has been used,
which is based on the NLO calculation by Nason, Dawson, and Ellis
\cite{bib:nde}.
This program provides weighted events with two or three partons that
result from the hard interaction, but no parton showers. The
implementation of beauty hadron formation and decay to leptons is described in  
our earlier publication \cite{Aktas:2005zc}.
The parton states were subjected to the same jet algorithm as the hadrons
in the PYTHIA and Cascade event samples, to form observables (in
particular, the $\pt$ of the most energetic jet $\ptjetone$, $\xgamma$,
and $\deltaphi$)
on parton level.
We have used the CTEQ5F4 \cite{Lai:1999wy} (GRV-G HO \cite{Gluck:1991jc}) 
parton density set for the proton
(photon).


Hadronisation corrections were calculated using the PYTHIA program;
for events generated with PYTHIA, the jet related observables were
calculated once using the final state partons, and a second time from
the resulting hadrons. From this, for each observable 
a migration matrix was determined that describes how the events in 
bins of the parton level observables migrate to the hadron level bins.
This is a refinement of the procedure used in our previous publication
\cite{Aktas:2005zc}, where only the bin-by-bin ratio of events on parton
and hadron level was determined and applied to the parton level
distributions. The resulting differences in the predicted hadron level
cross sections are, however, minor compared to the theoretical
uncertainties.

The theory uncertainties were determined by varying the input beauty
mass\footnote{
The QCD calculation uses the pole mass. The most recent evaluation by
the Particle Data Group \cite{Yao:2006px} 
yields $\overline m\sub b = 4.20\pm0.07\,\GeV$
for the $\overline{\mathrm{MS}}$ running mass, corresponding to 
$m\sub b = 4.80\pm0.07\,\GeV$ for the pole mass. 
Here, we continue to use the
mass range $m\sub b = 4.75\pm0.25\,\GeV$ for better comparison with
other publications.
}
$m\sub{b}$ from the nominal value $m\sub{b}=4.75\,\GeV$ up and down by
$0.25\,\GeV$. In addition, the renormalization and factorization scales
$\mur$ and $\muf$ were varied independently in the range
$\mu_0 / 2 \le \mur, \muf \le 2\,\mu_0$, with the constraint 
$1/2 \le \mur/\muf \le 2$.
This is the procedure adopted by the HERA-LHC workshop \cite[p. 406]{Alekhin:2005dy}.
This is again a refinement compared to our previous publication,
where $\mur$ and $\muf$ were kept equal and were varied 
at the same time as the beauty mass.
For each bin on hadron level, the deviations due to the beauty mass 
variation and the largest deviation due to the scale variation in the
upward and downward direction were determined and added in quadrature for the
total model uncertainty, following the prescription in \cite{Alekhin:2005dy}. 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Analysis Method}

Events were selected with the following experimental cuts:
\begin{itemize}
\item No electron candidate with an energy above $E>6\,\GeV$ 
  is allowed to be found in the detector.
  This removes most deep inelastic scattering (DIS) events with a momentum transfer 
  $Q^2$ of $Q^2 \gapprox 2.5\,\GeV^2$. 
  Monte Carlo simulations were used to correct for the remaining DIS background 
  with $Q^2>1\,\GeV^2$. 
\item The inelasticity $y$ as measured from the hadronic activity in the 
  detector using $y\sub h~=~\sum (E\sub i-p\sub{z,i})/(2 E_e)$
  must lie in the range $0.2 < y\sub h < 0.8$.
\item Two jets with transverse momenta $\ptjet > 7 (6)\,\GeV$
  for the highest (second highest) $\pt$ jet
  have to be found within the pseudorapitity range $-2.5 < \etajet< 2.5$.
  Jets were identified by the inclusive $k\sub t$ jet algorithm
  \cite{bib:ktalgorithm} in the $\pt$ recombination scheme, 
  with a distance parameter 
  $R = 1.0$.
  The jet algorithm was applied in the laboratory frame, using
  hadronic final state objects combining track and calorimeter information.
\item A muon, identified in the central muon system, with transverse momentum
  $\ptmu > 2.5\,\GeV$ has to be found in a pseudorapidity range
  $-0.55 < \etamu < 1.1$. The muon has to be associated to one of the two
  highest energy jets.
\end{itemize}

To ensure a good muon reconstruction and suppress events with cosmic
muons the following quality cuts were applied to the muon track:
\begin{itemize}
  \item The muon
  must be identified in the Central Muon Detector CMD, i.e. it must have
  at least three hits in the ten layers of limited streamer tubes of the
  CMD. In events with several muon candidates, 
  the one with  the highest $\ptmu$
  was selected.
  \item The $\chi^2$ probability for the match between the CJC and the CMD tracks must
  be larger than $5\,\%$.
  \item The signed impact parameter of the muon\footnote{The sign of the impact parameter is
  determined in relation to the jet axis and defined such that it is positive for
  muons originating from a secondary vertex displaced along the jet direction.
  } must lie in the range $-0.05 < \delta < 0.1\,\cm$ 
  (this cut
  reduces background from cosmic muons and inflight decays of pions and
  kaons).
  The impact parameter is
   calculated with respect to the primary vertex of the event,
   which is determined from all tracks except the muon track.
  \item To reduce background from cosmic muons, 
  the timing information of the CJC track must be consistent
  within $3\,\sigma$ (corresponding to $2.8\,\ns$) with the 
  average event time\footnote{For each event, a global $T_0$ is calculated from 
  the CJC tracks; this $T_0$ is averaged over many events to give the
  average event time.}.
  \item At least one hit in the CST must be linked to the CJC muon track.
  \item Events are rejected as cosmics
   where a second track in the opposite hemisphere 
   fitting the muon candidate in $\phi$ and $\pt$ 
   is found. 
\end{itemize}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Cross Section Definition and Measurement}

We have measured the cross section of the process $ep \to e b \bar b X \to ejj\mu X'$,
i.e. beauty production with the formation of two jets and 
the subsequent decay of a beauty hadron to a muon.
The muon may be produced by a direct semileptonic decay of a beauty hadron,
from a cascade decay, where the charm hadron decays semileptonically,
or from a $J/\psi$ or $\psi'$ decay.
The muon must be associated to either of the two highest $\pt$ jets.

The visible range is defined by

\begin{eqnarray*}
        &Q^2&< 1\,\GeV^2\\
  0.2 < &y& < 0.8 \\
        &\ptmu& > 2.5\,\GeV\\
  -0.55 < & \etamu& < 1.1\\
        &\ptjets& > 7 (6) \,\GeV  \\   
   -2.5 < &\etajets& < 2.5.   
\end{eqnarray*}

Events that pass the reconstruction cuts were binned differentially in one
of the following quantities:
\begin{itemize}
  \item The transverse momentum of the muon $\ptmu$,
  \item the transverse momentum $\ptjetone$ of the highest $\pt$ jet,
  \item the pseudorapidity $\etamu$ of the muon,
  \item \xgamma, the momentum fraction of the photon
  entering the hard interaction,
  \item and $\deltaphi$, the difference in azimuthal angle between the two jets.
\end{itemize}

To extract the beauty fraction, two quantities were used:
\begin{itemize}
\item $\ptrel$ is the transverse momentum of the muon with respect to the 
  axis of the most energetic jet, 
 reconstructed without the muon four vector.
\item $\delta$ is the impact parameter of the muon with respect to the primary
  vertex of the event.
\end{itemize}

For each bin, the two-dimensional distribution of these quantities
were fitted with three template distributions derived
from Monte Carlo simulations; these templates were generated separately
for events containing only light ($u$, $d$, and $s$) quarks, charm
quarks, and beauty quarks, respectively.
The fit takes into account the statistical uncertainties from
the data sample and the Monte Carlo templates
\cite{Barlow:1993dm, bib:root}.

The result of the fit is the relative amount of beauty induced events
$\fb$ in each analysis bin, from which the observed number of beauty events
$\Nb$ in the bin is calculated as
$
  \Nb = \fb \cdot \Nbin,
$
where $\Nbin$ is the total number of data events observed in the respective bin;
the data were then corrected for effects of detector resolution by a
matrix unfolding procedure, with a migration matrix determined from Monte Carlo
simulations.

\subsection{Control Distributions}

Fig.~\ref{fig:1} shows control distributions of several important kinematic
quantities for the selected sample: the muon transverse
momentum $\ptmu$ and its pseudorapidity $\etamu$, 
the transverse momenta $\ptjets$ of the highest and the
second highest $\pt$ jets,  and the observable $\xgamma$.
The data are compared to the PYTHIA Monte Carlo samples that were used
to correct the data. The PYTHIA sample was reweighted as a function of the
inelasticity $y$ to improve the description of the data, and the
relative fractions of light (uds), c and b quarks were adjusted to the
values observed in the fit to the data.
Overall, the description is very satisfactory.

Fig.~\ref{fig:2} shows the distributions of the quantities used to
extract the beauty fraction from the data:
the impact parameter
$\delta$ of the muon track and the transverse muon momentum
$\ptrel$ relative to the axis of the associated jet.
Both quantities are described quite well by the Monte Carlo simulation.
In particular, the impact parameter distribution, which is very
sensitive to the detector resolution, is very well described, in the
region $\delta<0$ that is dominated by resolution effects as well as in
the region $\delta>0$, which shows the tails due to long-lived particles
from charm and beauty decays.
This description has been achieved by careful tuning of the simulation
software with regard to signal heights, noise levels and dead strips
in the CST \cite{bib:mira}, the inclusion of effects from alignment
imperfections, and the description of the dead material in front of the
CST and the CJC. Further smearing of measured track parameters in the
simulation is not necessary.

As can be seen from Fig.~\ref{fig:2}, events with beauty quarks, 
which already contribute at the level of $\approx 30\,\%$ to the sample, are
enriched by cuts on the impact parameter or the relative $\pt$.
Fig.~\ref{fig:3} illustrates the effect of the cut on one variable on
the other one: The impact parameter distribution is shown for a beauty
enriched sample with $\ptrel>1.2\,\GeV$, and vice versa the $\ptrel$
distribution is shown for events with $0.01 < \delta < 0.1\,\cm$.
The Monte Carlo simulation describes these subsamples quite
satisfactorily, which demonstrates that the $\delta$ and $\ptrel$
distributions of the beauty and charm quark samples are very well
modeled.


\subsection{Systematic Uncertainties}

A number of sources of systematic errors were considered.
The following uncertainties contribute to the overall normalization uncertainty:
\begin{itemize}
\item The trigger efficiency,
\item the identification of the muon in the instrumented iron,
\item the efficiency to reconstruct the muon track in the 
  tracking system and to link a CST hit,
\item the luminosity.
\end{itemize}
Additional uncertainties affect the data differently in
various bins:
\begin{itemize}
\item The impact parameter resolution,
\item the reconstruction of the jet axis, which has an impact on the measurement of
  $\ptrel$,
\item the energy scale for hadrons of the calorimeter,
\item the model uncertainties, estimated by using the Monte Carlo generator CASCADE
  \cite{bib:cascade} instead of PYTHIA 6.2 \cite{bib:pythia},
\item the uncertainty from the fragmentation process, estimated by using the Lund \cite{Andersson:1983jt}
 instead of the Peterson \cite{Peterson:1982ak} fragmentation function,
\item the uncertainty from the fragmentation fractions of $c$ and $b$ quarks into hadrons,
their branching ratios and lifetimes,
\item and the uncertainty on the modelling of $\pi$ and $K$ inflight
decays.
\end{itemize}

The resulting systematic uncertainty is $12\,\%$.


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

The total cross section for the process $ep \to e b \bar b X \to ejj\mu X'$
in the visible range given by
$        Q^2< 1\,\GeV^2,
  0.2 < y < 0.8,
        \ptmu > 2.5\,\GeV,
  -0.55 < \etamu < 1.1,
        \ptjets > 7 (6) \,\GeV, $ and  
$   -2.5 < \etajets < 2.5 $ 
has been measured to be
$$
  \sigma\sub{vis}\,(ep \to e b \bar b X \to ejj\mu X')  
  = 31.4 \pm 1.3 (stat.) \pm 3.8 (syst.) \,\pb^{-1}.
$$
This result is somewhat lower than the published result \cite{Aktas:2005zc} from HERA-I,
which is
$$
  \sigma\sub{vis}\,(ep \to e b \bar b X \to ejj\mu X') 
  = 38.4 \pm 3.4 (stat.) \pm 5.4 (syst.) \,\pb^{-1},
$$
but compatible within errors.

In comparison, the FMNR calculation yields
$$
  \sigma\sub{vis}\super{FMNR}\,(ep \to e b \bar b X \to ejj\mu X') 
  = 25.3^{+6.4}_{-4.7} \,\pb^{-1},
$$
in agreement with the data,
and the PYTHIA prediction is
$$
  \sigma\sub{vis}\super{FMNR}\,(ep \to e b \bar b X \to ejj\mu X') 
  = 21.7 \,\pb^{-1}.
$$


The differential cross sections also tend to be lower, but compatible, with the HERA-I
results. 
The largest discrepancies are observed for the differential measurements
in the lowest bins of $\ptmu$ and $\ptjetone$, respectively.
For these bins, the ratio between the HERA-II and the HERA-I measurements is 
about $2.5\,\sigma$ below unity if systematic uncertainties that do not cancel
between both measurements are taken into account. 
It has been checked that this discrepancy is not caused by differences in the analysis
method between the HERA-I and HERA-II analyses, such as different fitting procedures and
different binnings in the two-dimensional fits used to extract the beauty fraction in
each bin, or the fact the the new analysis uses an unfolding procedure to extract the
differential cross section in   $\ptjetone$; we therefore attribute the difference to a
statistical fluctuation.

Fig.~\ref{fig:4}a) shows the differential cross section as a function of the muon's
pseudorapidity, $\etamu$. The measured values show little dependence on $\etamu$; 
in particular, there is no drop of the cross section in the forward direction,
in contrast to PYTHIA and FMNR expectations.

In Fig.~\ref{fig:4}b), the dependence on the muon's transverse momentum, $\ptmu$, is
compared to the predictions. PYTHIA gives a steeper fall off than the data, while
the NLO calculation describes the shape rather well.

The cross section also falls with the transverse momentum of the leading jet, $\ptjetone$,
albeit less steeply, as can be seen from 
Fig.~\ref{fig:4}c). The PYTHIA and FMNR calculations predict a rather
similar behaviour and are both in agreement with the data.

Fig.~\ref{fig:4}d) shows the cross section as function of $\xgamma$, which approximates
the photon's momentum entering the hard interaction. This observable is sensitive to the
relative amounts of beauty produced in direct ($\xgamma \gapprox 0.75$) and resolved 
($\xgamma \lapprox 0.75$) photon processes. Here, the NLO calculation gives a good
description of the data,
while PYTHIA underestimates the relative amount of events at high $\xgamma$.

Finally, Fig.~\ref{fig:4}e) shows the difference in azimuthal angle, $\deltaphi$,
between the two jets. In leading order QCD the two outgoing quarks must be exactly
opposite in azimuthal angle, corresponding to $\deltaphi = 180^\circ$. 
Thus, values of $\deltaphi$ substantially lower than $180^\circ$ are indicative of the
presence of further final state gluons, and therefore this quantity is very sensitive to
the description of gluon emission. PYTHIA, which employs parton showers to simulate 
the effect of multiple gluon emission, describes the shape of this observable
reasonably 
well, as does the fixed order calculation of FMNR, which allows at most one
hard gluon in the final state.

Overall, the data are reasonably well described  in shape by the predictions from PYTHIA,
but lie approximatly a factor $1.4$ above the PYTHIA prediction.

The NLO predictions, derived with the FMNR program, also lie systematically below the
data, but also describe the differential distributions well in shape.
In particular, the deficiency in the lowest bins 
of $\ptmu$ and $\ptjetone$ is not substantially larger than in the other bins,
in contrast to the findings of the HERA-I analysis. 

This observation is in agreement with the result from an analysis \cite{Chekanov:2003si}
of HERA-I data by ZEUS, where good agreement with the FMNR predictions over the whole range
in $\ptjetone$ and $\ptmu$ is found.
In the ZEUS analysis,
 the same process is measured, albeit in a slightly different phase
space (no cut on $\etamu$ is made in the ZEUS analysis). 
  


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Conclusions}

We have performed a measurement of the photoproduction of beauty quarks,
using events where at least one beauty hadron decays with a muon in the final state, and
two jets are visible in the detector.
Beauty events are separated from background by means of the transverse momentum
$\ptrel$ of the muon with respect to a jet, and by its impact parameter $\delta$,
utilizing the large lifetime of beauty hadrons. 
The measurement extends over the phase space region defined by
$        Q^2< 1\,\GeV^2,
  0.2 < y < 0.8,
        \ptmu > 2.5\,\GeV,
  -0.55 < \etamu < 1.1,
        \ptjets > 7 (6) \,\GeV, $ and  
$   -2.5 < \etajets < 2.5.  $ 

The visible cross section has been measured to be
$$
  \sigma\sub{vis}\,(ep \to e b \bar b X \to ejj\mu X')  
  = 31.4 \pm 1.3 (stat.) \pm 3.8 (syst.) \,\pb^{-1}.
$$
A NLO QCD calculation is in agreement with this measurement within the theoretical uncertainties.

Differential cross sections have been measured as function of the observables
$\etamu$, $\ptmu$, $\ptjetone$, $\xgamma$, and $\deltaphi$. 
The shape of these distributions is reasonably well described by 
the NLO QCD calculation as well as the PYTHIA LO Monte Carlo program.

At low values of $\ptmu$ and $\ptjetone$, the new measurement lies lower than 
the previous HERA-I measurement published by H1, and is thus better described by the NLO
predictions than the previous measurement. 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section*{Acknowledgements}

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

This analysis was initiated by Beate Naroska,
who passed away much too early.
We dedicate this work to her memory.

%\\(You can of course thank people external to H1 for discussions and help when they contributed to clarifying the physics analysis, provided software programs, theory calculations etc.)



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{References}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{thebibliography}{99}

\bibitem{Aktas:2005zc}
  A.~Aktas {\it et al.}  [H1 Collaboration],
  \artitle{Measurement of beauty production at HERA using events with muons and jets,}
  Eur.\ Phys.\ J.\  C {\bf 41} (2005) 453
  [hep-ex/0502010].
  %%CITATION = EPHJA,C41,453;%%
  
\bibitem{Chekanov:2003si}
  S.~Chekanov {\it et al.}  [ZEUS Collaboration],
  \artitle{Bottom photoproduction measured using decays into muons in dijet events  in $ep$ collisions at $\sqrt s = 318\,{\mathrm{GeV}}$,}
  Phys.\ Rev.\  D {\bf 70} (2004) 012008
  [Erratum-ibid.\  D {\bf 74} (2006) 059906]
  [hep-ex/0312057].
  %%CITATION = PHRVA,D70,012008;%%


\bibitem{H1detector}
I.~Abt {\it et al.}  [H1 Collaboration],
\artitle{The H1 Detector at HERA,}
Nucl.\ Instrum.\ Meth.\ A {\bf 386} (1997) 310 and 348.
%%CITATION = NUIMA,A386,310;%%
%%CITATION = NUIMA,A386,348;%%

\bibitem{Appuhn:1996na}
R.~D.~Appuhn {\it et al.}  [H1 SpaCal Group],
\artitle{The H1 lead/scintillating-fibre calorimeter,}
Nucl.\ Instrum.\ Meth.\ A {\bf 386} (1997) 397.
%%CITATION = NUIMA,A386,397;%%
    
\bibitem{bib:cst}
%\bibitem{Pitzl:2000wz}
  D.~Pitzl {\it et al.},
  \artitle{The H1 silicon vertex detector,}
  Nucl.\ Instrum.\ Meth.\  A {\bf 454} (2000) 334
  [hep-ex/0002044];
  %%CITATION = NUIMA,A454,334;%%
  \\
%\bibitem{List:2002bd}
  B.~List,
  \artitle{The H1 central silicon tracker,}
  Nucl.\ Instrum.\ Meth.\  A {\bf 501} (2001) 49.
  %%CITATION = NUIMA,A501,49;%% 
     
\bibitem{bib:ftt}
%\bibitem{Baird:2001xc}
  A.~Baird {\it et al.},
  \artitle{A fast high resolution track trigger for the H1 experiment,}
  IEEE Trans.\ Nucl.\ Sci.\  {\bf 48} (2001) 1276
  [hep-ex/0104010];
  %%CITATION = HEP-EX 0104010;%%
  \\
%\bibitem{Baird:2001sz}
  A.~Baird {\it et al.}  [H1 Collaboration],
  \artitle{A fast track trigger for the H1 Collaboration,}
  Nucl.\ Instrum.\ Meth.\ A {\bf 461} (2001) 461;
  %%CITATION = NUIMA,A461,461;%%
\\
%\bibitem{Schoning:2004jp}
  A.~Sch\"oning  [H1 Collaboration],
  \artitle{A fast track trigger for the H1 collaboration,}
  Nucl.\ Instrum.\ Meth.\ A {\bf 518} (2004) 542;
  %%CITATION = NUIMA,A518,542;%%  
  \\
  N.~Berger {\it et al.},
  \artitle{First Results from the First Level of the H1 Fast Track Trigger,}
  Nucl.\ Sc.\ Symp.\ Conf.\ Rec., 2004 IEEE,
  {\bf 3} (2004) 1976.
  % 16-22 Oct. 2004


%  \bibitem{RAPGAP} 
% H.~Jung,
% %``Hard diffractive scattering in high-energy e p collisions and the Monte Carlo generation RAPGAP,''
% Comput.\ Phys.\ Commun.\  {\bf 86} (1995) 147; \\
% \; RAPGAP 3.2 program manual unpublished \\
% \;http://www-h1.desy.de/$\sim$jung/RAPGAP.html.

\bibitem{bib:pythia} 
%\bibitem{Sjostrand:2000wi}
  T.~Sj\"ostrand {\it et al.}, % P.~Eden, C.~Friberg, L.~L\"onnblad, G.~Miu, S.~Mrenna and E.~Norrbin,
  \artitle{High-energy-physics event generation with PYTHIA 6.1,}
  Comput.\ Phys.\ Commun.\  {\bf 135} (2001) 238
  [hep-ph/0010017].
  %%CITATION = CPHCB,135,238;%%
  
\bibitem{bib:geant}
  R.~Brun {\it et al.},
  CERN-DD/EE--84--1 (1987).
  
%\cite{Pumplin:2002vw}
\bibitem{Pumplin:2002vw}
  J.~Pumplin {\it et al.}, %, D.~R.~Stump, J.~Huston, H.~L.~Lai, P.~Nadolsky and W.~K.~Tung,
  \artitle{New generation of parton distributions with uncertainties from global  QCD
  analysis,}
  JHEP {\bf 0207} (2002) 012
  [hep-ph/0201195].
  %%CITATION = JHEPA,0207,012;%%

\bibitem{Schuler:1996fc}
  G.~A.~Schuler and T.~Sj\"ostrand,
  \artitle{Parton distributions of the virtual photon,}
  Phys.\ Lett.\  B {\bf 376} (1996) 193
  [hep-ph/9601282].
  %%CITATION = PHLTA,B376,193;%%  

\bibitem{Jung:2004gs}
  H.~Jung,
  \artitle{Un-integrated PDFs in CCFM,}
  hep-ph/0411287 (2004).
  %%CITATION = HEP-PH/0411287;%%

\bibitem{Peterson:1982ak}
  C.~Peterson {\it et al.}, %D.~Schlatter, I.~Schmitt and P.~M.~Zerwas,
  \artitle{Scaling violations in inclusive $e^+ e^-$ annihilation spectra,}
  Phys.\ Rev.\  D {\bf 27} (1983) 105.
  %%CITATION = PHRVA,D27,105;%%  
  
\bibitem{Frixione:1995qc}
  S.~Frixione, P.~Nason and G.~Ridolfi,
  \artitle{Differential distributions for heavy flavor production at HERA,}
  Nucl.\ Phys.\  B {\bf 454} (1995) 3
  [hep-ph/9506226].
  %%CITATION = NUPHA,B454,3;%%  

  
\bibitem{bib:nde}
%\bibitem{Nason:1987xz}
  P.~Nason, S.~Dawson and R.~K.~Ellis,
  \artitle{The total cross-section for the production of heavy quarks in hadronic
  collisions,}
  Nucl.\ Phys.\  B {\bf 303} (1988) 607;\\
  %%CITATION = NUPHA,B303,607;%%
%
%\bibitem{Nason:1989zy}
  {\it idem,} %P.~Nason, S.~Dawson and R.~K.~Ellis,
  \artitle{The one particle inclusive differential cross-section for heavy quark
  production in hadronic collisions,}
  Nucl.\ Phys.\  B {\bf 327} (1989) 49
  [Erratum-ibid.\  B {\bf 335} (1990) 260];\\
  %%CITATION = NUPHA,B327,49;%%  
%
%\bibitem{Ellis:1988sb}
  R.~K.~Ellis and P.~Nason,
  \artitle{QCD radiative corrections to the photoproduction of heavy quarks,}
  Nucl.\ Phys.\  B {\bf 312} (1989) 551.
  %%CITATION = NUPHA,B312,551;%%

\bibitem{Yao:2006px}
  W.~M.~Yao {\it et al.}  [Particle Data Group],
  \artitle{Review of particle physics,}
  J.\ Phys.\ G {\bf 33} (2006) 1.
  %%CITATION = JPHGB,G33,1;%%
  
\bibitem{Lai:1999wy}
  H.~L.~Lai {\it et al.}  [CTEQ Collaboration],
  \artitle{Global {QCD} analysis of parton structure of the nucleon: CTEQ5 parton
  distributions,}
  Eur.\ Phys.\ J.\  C {\bf 12} (2000) 375
  [hep-ph/9903282].
  %%CITATION = EPHJA,C12,375;%%

\bibitem{Gluck:1991jc}
  M.~Gl\"uck, E.~Reya and A.~Vogt,
  \artitle{Photonic parton distributions,}
  Phys.\ Rev.\  D {\bf 46} (1992) 1973.
  %%CITATION = PHRVA,D46,1973;%%  
    
\bibitem{Alekhin:2005dy}
  S.~Alekhin {\it et al.},
  \artitle{HERA and the LHC - A workshop on the implications of HERA for LHC  physics:
  Proceedings Part B,}
  hep-ph/0601013 (2006).
  %%CITATION = HEP-PH/0601013;%%  

\bibitem{bib:ktalgorithm}
%\bibitem{Ellis:1993tq}
  S.~D.~Ellis and D.~E.~Soper,
  \artitle{Successive combination jet algorithm for hadron collisions,}
  Phys.\ Rev.\  D {\bf 48} (1993) 3160
  [hep-ph/9305266];
  %%CITATION = PHRVA,D48,3160;%%
  \\
%\bibitem{Cacciari:2005hq}
  M.~Cacciari and G.~P.~Salam,
  \artitle{Dispelling the $N^3$ myth for the $k\sub t$ jet-finder,}
  Phys.\ Lett.\  B {\bf 641} (2006) 57 
  [hep-ph/0512210].
  %%CITATION = PHLTA,B641,57;%%  
  

%%\bibitem{Bassler:1994uq}
%%  U.~Bassler and G.~Bernardi,
%%  \artitle{On the kinematic reconstruction of deep inelastic scattering at HERA: The Sigma method,}
%%  Nucl.\ Instrum.\ Meth.\  A {\bf 361} (1995) 197
%%  [hep-ex/9412004].
%%  xxCITATION = NUIMA,A361,197;xx


\bibitem{Barlow:1993dm}
  R.~J.~Barlow and C.~Beeston,
  \artitle{Fitting using finite Monte Carlo samples,}
  Comput.\ Phys.\ Commun.\  {\bf 77} (1993) 219.
  Implemented in class TFractionFitter in the Root program.  
  %%CITATION = CPHCB,77,219;%%  

\bibitem{bib:root}
  R.~Brun {\it et al.,}
  {\it ``ROOT, an object oriented data analysis framework, user's guide 5.14,''}
  available at {\tt ftp://root.cern.ch/root/doc/Users\_Guide\_5\_14.pdf}.

\bibitem{bib:mira}
  M.~Kr\"amer,
  {\it ``Tuning the vertex detector simulation of H1,''}
  Proceedings of the International School of Subnuclear Physics:
   $45^{th}$ course: Search for the `totally unexpected' in the LHC era.
   Erice, Sicily, 29.8.-7.9.2007.

\bibitem{bib:cascade} 
%\bibitem{Jung:2000hk}
  H.~Jung and G.~P.~Salam,
  \artitle{Hadronic final state predictions from CCFM: The hadron-level Monte  Carlo
  generator CASCADE,}
  Eur.\ Phys.\ J.\  C {\bf 19} (2001) 351
  [hep-ph/0012143];\\
  %%CITATION = EPHJA,C19,351;%%
%\bibitem{Jung:2001hx}
  H.~Jung,
  \artitle{The CCFM Monte Carlo generator CASCADE,}
  Comput.\ Phys.\ Commun.\  {\bf 143} (2002) 100
  [hep-ph/0109102].
  %%CITATION = CPHCB,143,100;%%
 
\bibitem{Andersson:1983jt}
  B.~Andersson, G.~Gustafson and B.~S\"oderberg,
  \artitle{A general model for jet fragmentation,}
  Z.\ Phys.\  C {\bf 20} (1983) 317.
  %%CITATION = ZEPYA,C20,317;%%  

\end{thebibliography}

\clearpage 

\begin{figure}[tbp]
\begin{center}
\setlength{\unitlength}{1cm}
\begin{picture}(15.3,18)
 \put(0,0){\epsfig{file=H1prelim-08-071.fig1f.eps,width=0.48\textwidth}}
 \put(7.8,0){\epsfig{file=H1prelim-08-071.fig1e.eps,width=0.48\textwidth}}
 \put(0,6){\epsfig{file=H1prelim-08-071.fig1c.eps,width=0.48\textwidth}}
 \put(7.8,6){\epsfig{file=H1prelim-08-071.fig1d.eps,width=0.48\textwidth}}
 \put(0,12){\epsfig{file=H1prelim-08-071.fig1a.eps,width=0.48\textwidth}}
 \put(7.8,12){\epsfig{file=H1prelim-08-071.fig1b.eps,width=0.48\textwidth}}
 \put(6.5,17.2){a)}
 \put(14.3,16.7){b)}
 \put(6.5,11.2){c)}
 \put(14.3,11.2){d)}
 \put(6.5,5.2){e)}
 \put(13.8,5.2){f)}

\end{picture}
\end{center}
    \caption{Distributions of a) the muon transverse
momentum $\ptmu$, 
b) the pseudorapidity of the muon $\etamu$, c)
and d) the transverse momenta $\ptjets$ of the highest and the
second-highest $\pt$ jets, respectively, e) the observable $\xgamma$, and
f) the azimuthal angle between the jets, $\deltaphi$.
Included in the figure are the estimated contributions of events
arising from beauty quarks (dark grey line), charm quarks (light grey line) and
light quarks (dotted line). The shapes of the distributions
from the different sources are taken from the PYTHIA Monte
Carlo simulation and their relative fractions are determined
from a fit to the two-dimensional data distribution of $\ptrel$
and the impact parameter $\delta$ (see text).}
 \label{fig:1}
\end{figure}

\begin{figure}[tbp]
\begin{center}
\setlength{\unitlength}{1cm}
\begin{picture}(15,6)
 \put(0,0){\epsfig{file=H1prelim-08-071.fig2a.eps,width=0.48\textwidth}}
 \put(7.5,0){\epsfig{file=H1prelim-08-071.fig2b.eps,width=0.48\textwidth}}
 \put(6.4,5.3){a)}
 \put(14,5.2){b)}
\end{picture}
\end{center}
    \caption{Distributions of a) the impact parameter
$\delta$ of the muon track and b) the transverse muon momentum
$\ptrel$ relative to the axis of the associated jet.
    Included in the figure are the estimated contributions of events
arising from beauty quarks (dark grey line), charm quarks (light gerey line) and
light quarks (dotted line). The shapes of the distributions
from the different sources are taken from the PYTHIA Monte
Carlo simulation and their relative fractions are determined
from a fit to the two-dimensional data distribution of $\ptrel$
t and the impact parameter $\delta$ (see text).}
 \label{fig:2}
\end{figure}

\begin{figure}[tbp]
\begin{center}
\setlength{\unitlength}{1cm}
\begin{picture}(15,6)
 \put(0,0){\epsfig{file=H1prelim-08-071.fig3a.eps,width=0.48\textwidth}}
 \put(7.5,0){\epsfig{file=H1prelim-08-071.fig3b.eps,width=0.48\textwidth}}
 \put(6.5,5.3){a)}
 \put(14,5.2){b)}
\end{picture}
\end{center}
    \caption{
    Distributions in the restricted sample
of a) the impact parameter $\delta$ for events with 
$\ptrel > 1.2\,\GeV$
and b) the transverse muon momentum $\ptrel$
relative to the jet
axis for tracks with impact parameter $\delta > 0.01\,\cm$. The predictions
for the contributions to the restricted sample
    from beauty quarks (dark grey line), charm quarks (light grey line) and
light quarks (dotted line). The shapes of the distributions
from the different sources are taken from the PYTHIA Monte
Carlo simulation and their relative fractions are determined
from a fit to the two-dimensional data distribution of $\ptrel$
and the impact parameter $\delta$ (see text).}
 \label{fig:3}
\end{figure}


\begin{figure}[tbp]
\setlength{\unitlength}{1cm}
\begin{center}
\begin{picture}(15,17)
 \put(0,0){\epsfig{file=H1prelim-08-071.fig4e.eps,width=0.48\textwidth}}
 \put(0,6){\epsfig{file=H1prelim-08-071.fig4c.eps,width=0.48\textwidth}}
 \put(7.5,6){\epsfig{file=H1prelim-08-071.fig4d.eps,width=0.48\textwidth}}
 \put(0,12){\epsfig{file=H1prelim-08-071.fig4a.eps,width=0.48\textwidth}}
 \put(7.5,12){\epsfig{file=H1prelim-08-071.fig4b.eps,width=0.48\textwidth}}
 \put(6.4,17.2){a)}
 \put(13.9,17.2){b)}
 \put(6.4,11.2){c)}
 \put(13.9,11.2){d)}
 \put(6.4,5.2){e)}

\end{picture}
\end{center}
    \caption{
Differential cross sections for the photoproduction process
$ep \to eb\bar b X \to ejj\mu X$ in the kinematic range 
$Q^2 < 1\,\GeV^2$, $0.2 < y < 0.8$, $\ptmu > 2.5\,\GeV$, 
$0.55 < \etamu < 1.1$, $\ptjets > 7(6)\,\GeV$ and $|\etajets| < 2.5$. The cross sections are
shown as functions of a) the muon pseudorapidity $\etamu$, b) the
muon transverse momentum $\ptmu$, c) the jet transverse momentum
$\ptjetone$ of the highest transverse momentum jet, d) the photon's momentum
fraction $\xgamma$ entering the hard interaction, 
and e) the azimuthal angle difference $\deltaphi$ between the jets.
 The inner error bars show the statistical error,
the outer error bars represent the statistical and systematic
uncertainties added in quadrature. The NLO QCD predictions
are corrected to the hadron
level (solid line) using the PYTHIA generator. The shaded
band around the hadron level prediction indicates the systematic
uncertainties as estimated from scale variations (see text).
Predictions from the Monte Carlo generator PYTHIA (dotted line) are also
shown.}
 \label{fig:4}
\end{figure}

\end{document}




