%================================================================
% 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{G\eV}^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{G\eV}}
\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{G\eV}^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^*} }

% BtoMuX stuff
\newcommand{\ptjetone}{\ensuremath{P_{T}^{\text{jet1}}}}
\newcommand{\ptjettwo}{\ensuremath{P_{T}^{\text{jet2}}}}
\newcommand{\ptmuon}{\ensuremath{P_{T}^{\mu}}}
\newcommand{\ptmuonjet}{\ensuremath{P_{T}^{\mu\text{jet}}}}
\newcommand{\ptotherjet}{\ensuremath{P_{T}^{\text{other jet}}}}
\newcommand{\etajet}{\ensuremath{\eta^{\text{jet}}}}
\newcommand{\etamuon}{\ensuremath{\eta^{\mu}}}
\newcommand{\phimuon}{\ensuremath{\phi^{\mu}}}
\newcommand{\ptrel}{\ensuremath{P_{T}^{rel}}}
\newcommand{\muonjet}{\ensuremath{\mu\text{-jet}}}
\newcommand{\xgamma}{\ensuremath{x_{\gamma}^{\text{obs}}}}
\newcommand{\dphijets}{\ensuremath{\Delta\phi^{\text{jj}}}}
\newcommand{\impactpar}{\ensuremath{\delta}}





%%% 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.}

\begin{titlepage}

\noindent
\begin{flushleft}
{\tt DESY 12-059    \hfill    ISSN 0418-9833} \\
{\tt April 2012}                  \\
\end{flushleft}


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

{\bf Measurement of Beauty and Charm Photoproduction using Semi-muonic Decays in Dijet Events at HERA }

\vspace{1cm}%2cm

H1 Collaboration

\end{Large}
\end{center}

\vspace{1cm} %2cm

\begin{abstract}
\noindent Measurements of cross sections for beauty and charm events with dijets and a muon in the photoproduction regime at HERA are presented. The data were collected with the H1 detector and correspond to an integrated luminosity of $179$ pb$^{-1}$. Events with dijets of transverse momentum \ptjetone $> 7$ G\eV \ and \ptjettwo $> 6$ G\eV \ in the pseudorapidity range $-1.5<\etajet < 2.5$ in the laboratory frame are selected in the kinematic region of photon virtuality \qsq $< 2.5$ \gevsq and inelasticity $0.2 < y < 0.8$. One of the two selected jets must be associated to a muon with \ptmuon $> 2.5$ G\eV\  in the pseudorapidity range $-1.3 <$ \etamuon $< 1.5$.  The fractions of beauty and charm events are determined using  the impact parameters of the muon tracks with respect to the primary vertex and their transverse momentum relative to the axis of the associated jet. 
Both variables are reconstructed using the H1 vertex detector. The measurements are in agreement with QCD predictions at leading and next-to-leading order. 
\end{abstract}

\vspace{0.5cm}%1.5cm

\begin{center}
Submitted to \EJC \;\; 
\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}
%-- H1AUTS Author list by names 
%-- Status: Tue Apr  3 13:47:10 CEST 2012  Number of authors = 197 

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

%\newpage

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

\bigskip
 $ ^{42}$ Also at Physics Department, National Technical University,
          Zografou Campus, GR-15773 Athens, Greece \\
 $ ^{43}$ Also at Rechenzentrum, Universit\"at Wuppertal,
          Wuppertal, Germany \\
 $ ^{44}$ Also at IPNL, Universit\'e Claude Bernard Lyon 1, CNRS/IN2P3,
          Villeurbanne, France \\
 $ ^{45}$ Also at CERN, Geneva, Switzerland \\
 $ ^{46}$ Also at Max-Planck-Institut f\"ur Physik, M\"unchen, Germany \\
 $ ^{47}$ Also at Comenius University, Bratislava, Slovak Republic \\
 $ ^{48}$ Also at Faculty of Physics, University of Bucharest,
          Bucharest, Romania \\
 $ ^{49}$ Also at Ulaanbaatar University, Ulaanbaatar, Mongolia \\
 $ ^{50}$ Supported by the Initiative and Networking Fund of the
          Helmholtz Association (HGF) under the contract VH-NG-401. \\
 $ ^{51}$ Absent on leave from NIPNE-HH, Bucharest, Romania \\
 $ ^{52}$ Also at  Department of Physics, University of Toronto,
          Toronto, Ontario, Canada M5S 1A7 \\
 $ ^{53}$ Also at LAPP, Universit\'e de Savoie, CNRS/IN2P3,
          Annecy-le-Vieux, France \\


\bigskip
 $ ^a$ Supported by the Bundesministerium f\"ur Bildung und Forschung, FRG,
      under contract numbers 05H09GUF, 05H09VHC, 05H09VHF,  05H16PEA \\
 $ ^b$ Supported by the UK Science and Technology Facilities Council,
      and formerly by the UK Particle Physics and
      Astronomy Research Council \\
 $ ^c$ Supported by FNRS-FWO-Vlaanderen, IISN-IIKW and IWT
      and  by Interuniversity
Attraction Poles Programme,
      Belgian Science Policy \\
 $ ^d$ Partially Supported by Polish Ministry of Science and Higher
      Education, grant  DPN/N168/DESY/2009 \\
 $ ^e$ Supported by the Deutsche Forschungsgemeinschaft \\
 $ ^f$ Supported by VEGA SR grant no. 2/7062/ 27 \\
 $ ^g$ Supported by the Swedish Natural Science Research Council \\
 $ ^h$ Supported by the Ministry of Education of the Czech Republic
      under the projects  LC527, INGO-LA09042 and
      MSM0021620859 \\
 $ ^i$ Supported by the Swiss National Science Foundation \\
 $ ^j$ Supported by  CONACYT,
      M\'exico, grant 48778-F \\
 $ ^k$ Russian Foundation for Basic Research (RFBR), grant no 1329.2008.2
      and Rosatom \\
 $ ^l$ This project is co-funded by the European Social Fund  (75\%) and
      National Resources (25\%) - (EPEAEK II) - PYTHAGORAS II \\
 $ ^m$ Supported by the Romanian National Authority for Scientific Research
      under the contract PN 09370101 \\
 $ ^n$ Partially Supported by Ministry of Science of Montenegro,
      no. 05-1/3-3352 \\
}
\end{flushleft}
%
% Please not that the author list may need re-formatting.

%\newpage

%\tableofcontents

\newpage

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Introduction}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The photoproduction of heavy quarks at the HERA electron-proton collider is of
particular interest for testing calculations in the framework of perturbative
quantum chromodynamics (pQCD). The masses $m_b$ and $m_c$ of the beauty and
charm quarks, as well as the transverse momentum of a jet, provide a hard scale,
which is necessary for the calculation of pQCD predictions \cite{mcatnlo}. The
dominant process for beauty and charm production at HERA is boson-gluon fusion
(BGF): $\gamma g \rightarrow Q\bar{Q}X$, with $Q=b, c$. In the kinematic regime
of photoproduction, where the photon virtuality \qsq is small, two classes of processes contribute to BGF. In direct photon processes, the photon emitted from the electron\footnote{In this paper the term 'electron' is used generically to refer to both electrons and positrons.} enters the hard BGF process directly. In resolved photon processes, the photon fluctuates into a hadronic state before the hard interaction and acts as a source of partons, one of which takes part in the hard interaction. At leading order (LO) pQCD resolved photon processes are expected to contribute significantly in the photoproduction region.
\par
Heavy quark photoproduction has been investigated at HERA using different analysis techniques in different regions of phase space. The lifetime or mass of heavy flavoured hadrons \cite{ZeusLifetime, H1Lifetime2010, H1Lifetime2011, H1Hera1desy05004, ZEUSHera1desy03212,ZEUSdesy08210}, semileptonic decays \cite{ZeusSemileptonic, H1Hera1desy05004, ZEUSHera1desy03212} or the full reconstruction of a $D$ or $D^*$ meson \cite{H1D, ZeusD} are exploited to perform the measurements. In general, the measured cross sections agree reasonably well with the theoretical predictions. In the H1 measurement performed with HERA I data \cite{H1Hera1desy05004}, where beauty photoproduction was investigated using two jets and a muon in the final state, the data were found to be described reasonably well by the pQCD calculations at next-to-leading order (NLO), except in the region of low transverse momentum of the muon $2.5<\ptmuon\ <3.3$ \gev \ and of low leading jet transverse momentum $7.0< \ptjetone <10.0$ \gev, where the NLO prediction was lower than the data. Comparable measurements were performed by the ZEUS collaboration \cite{ZEUSHera1desy03212,ZEUSdesy08210}, covering a similar phase space. Here, a good description by the NLO QCD predictions is found, also at low transverse momenta of the leading jet and the muon. 
\par
This paper reports on beauty and charm measurements of cross sections for
photoproduction of events with two jets and a muon, where the muon is associated
with one of the jets. The analysis exploits the lifetime and the mass of heavy
flavoured hadrons as in the former HERA I beauty production analysis
\cite{H1Hera1desy05004}. The measurement is additionally performed for charm
photoproduction. Increased statistics due to increased luminosity and extended
phase space, a better understanding and description of the H1 vertex detector
\cite{CSTTuning}, which is crucial for this ana\-lysis, as well as reduced
systematic uncertainties compared to the previous H1 ana\-lysis make a more
detailed test of pQCD predictions possible. Going beyond the earlier HERA I
analysis, the cross sections as a function of the azimuthal angular difference \dphijets\  between the two leading jets, which are sensitive to higher order corrections, are also measured. Furthermore, cross sections are measured in two different regions of phase space, which are either enriched by resolved or direct photon processes. 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{QCD Calculations}\label{MonteCarloSimulation}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The data presented here are compared with LO calculations supplemented by parton showers as well as with NLO calculations. The calculations are performed using either collinear factorisation, which is based on the DGLAP evolution equations \cite{DGLAP}, or the $k_t$-factorisation approach, which employs the CCFM \cite{CCFM} evolution equations. In the collinear approach transverse momenta obtained through initial state QCD evolution are neglected and all the transverse momenta are generated in the hard scattering process, i.e. the partons entering the hard interaction are collinear with the proton. Effects from non-zero transverse momenta of the gluons enter only at NLO. In the $k_t$-factorisation ansatz the transverse momenta of gluons $k_t$ entering the hard interaction are already included at leading order both in the off-shell matrix element and the $k_t$-dependent unintegrated gluon density \cite{uPDFs}. Therefore, corrections appearing only at higher order in collinear factorisation are partially included at LO in the $k_t$-factorisation approach.
\par
For beauty and charm photoproduction two classes of processes occur, the direct photon processes and the resolved photon processes.
The distinction between these two classes depends on the factorisation scheme and the order in which the calculation is performed. 
\par
The production of heavy quarks is calculated either in the massive scheme, where heavy quarks are produced only perturbatively via boson gluon fusion, or in the massless scheme, where heavy quarks are treated as massless partons. These two schemes are expected to be appropriate in different regions of phase space \cite{schemes}: the massive scheme is expected to be reliable when the transverse momentum $P_T$ of the heavy quarks is of similar size compared to the heavy quark mass $m_Q$, whereas the massless scheme is expected to be valid for $P_T\gg m_Q$.
\par
An overview of the parameters used in the Monte Carlo (MC) generators is given in table \ref{table_qcdsettings}. The following MC generators are used:

\begin{description}
\item{\textsc{Pythia}: The MC program \textsc{Pythia} 6.4 \cite{Pythia} is based on LO QCD matrix elements with leading-log parton showers in the collinear factorisation approach. \textsc{Pythia} includes both direct photon gluon fusion and resolved photon processes. In the resolved photon processes either a beauty or a charm quark or a gluon from the photon enters the hard scattering. In the inclusive mode of \textsc{Pythia} used here beauty and charm quarks are treated as massless partons in all steps of the calculation in both types of processes. Three samples are generated containing photoproduction events for the processes $ep\rightarrow b\mu X$, $ep\rightarrow c\mu X$ and $ep\rightarrow qjX$ where $q$ is a light quark of flavour $u$, $d$ or $s$ and $j$ denotes a jet. The latter sample is generated without specifically requiring a muon in order to use it for studying the background arising from muon candidates which originate from sources other than beauty or charm processes.
The hadronisation process is simulated according to the Lund string model  \cite{LundFrag}. For the systematic uncertainty arising from the fragmentation model, additional samples are generated using the Peterson fragmentation function \cite{PeteFrag} for heavy quarks.
}
%\end{description}

%\begin{description}
\item{\textsc{Cascade}: The MC program \textsc{Cascade} 2.0 \cite{Cascade} is used for simulating events based on LO QCD calculations in the $k_t$-factorisation approach. Only the direct boson gluon fusion process is implemented using off-shell matrix elements. Higher order QCD corrections are simulated with initial state parton showers applying the CCFM evolution \cite{CCFM}. Here, two samples containing the processes $ep\rightarrow b\mu X$ and $ep\rightarrow c\mu X$ are generated. The unintegrated PDFs of the proton from set A0 \cite{CCFMSetA0} are used. The hadronisation of partons is performed with the Lund string model as implemented in \textsc{Pythia}. 
}
%\end{description}

%\begin{description}
\item{\textsc{Herwig}: The MC program \textsc{Herwig} 6.510 \cite{herwig} is used to simulate events based on collinear factorisation and massless LO QCD calculations. \textsc{Herwig} includes both direct photon gluon fusion and resolved photon processes. The hadronisation of partons is performed using the cluster fragmentation~\cite{cluster}. 
}
%\end{description}
%\begin{description}
\item{MC@NLO: The MC@NLO program for HERA \cite{mcatnlo} combines a collinear NLO calculation of heavy flavour production in the massive approach \cite{fmnr} with parton showers and hadronisation as described in \cite{mcatnlohad}. The direct and resolved parts of the cross section are calculated separately. MC@NLO uses parton showers applying the DGLAP evolution to simulate higher order contributions and the cluster fragmentation~\cite{cluster} as implemented in \textsc{Herwig} \cite{herwig}. The quark masses are set to $m_c=1.5$ G\eV \ and $m_b=4.75$ G\eV. The central value of the renormalisation scale $\mu_{R}$ is set to $\sqrt{m_Q^2+P_T^2}$, while the factorisation scale $\mu_F$ is $2\mu_R$. As an estimate of the theoretical uncertainties on the NLO QCD predictions the scales $\mu_R$ and $\mu_F$ are varied separately by factors of $0.5$ and $2$, while $m_b$ is changed by \mbox{$\pm0.25$ G\eV} and $m_c$ is changed by $\pm0.2$ G\eV. The resulting variations are added in quadrature to obtain the resulting theoretical uncertainty. 
}
\end{description}

\renewcommand{\arraystretch}{1.18} % this increase the table in vertical direction, very useful!

\begin{table}

\centering
\begin{tabular}{|c|c|c|c|c|}

\hline
               &\textsc{Pythia} & \textsc{Cascade} & \textsc{Herwig} & MC@NLO \\ \hline \hline
 Version & 6.4 & 2.0 & 6.510 & \textsc{Herwig} 6.510     \\  
 Evolution scheme & DGLAP & CCFM &  DGLAP  &  DGLAP    \\ \hline                          
 $m_b$ [G\eV] & $4.75$ & $4.75$ & $4.75$ & $4.75$  \\
 $m_c$ [G\eV] & $1.5$   & $1.5$ & $1.5$  & $1.5$   \\ \hline
 Proton PDF & CTEQ6L1\cite{CTEQ6}   & A0\cite{CCFMSetA0} & HERAPDF 1.0\cite{HERAPDF} & HERAPDF 1.0   \\
 Photon PDF & SaS2D\cite{Sas2d}   &  -  & GRV-G LO\cite{GRVGLO} & GRV-G LO   \\ \hline
& & & & \\[-0.5cm]
Renorm. scale $\mu_R$ & $\sqrt{P_T^2+\frac{1}{2}(m_1^2+m_2^2)}$  & $\sqrt{\hat{s}+Q_{\bot}^2}$  & $\sqrt{m_Q^2+P_T^2}$  & $\sqrt{m_Q^2+P_T^2}$ \\
Factor. scale $\mu_F$ & $\mu_R$ & $\mu_R$ &$\mu_R$ & $2\mu_R$\\ \hline
Fragmentation & Lund  & Lund & cluster & cluster \\    \hline
                           & $a=0.437$  &  $a=0.437$ & - & - \\
                            &  $b=0.850$ & $ b=0.850$ & - & -\\
                           & Peterson & - & - & - \\
                           & $\epsilon_b=0.006$ & - & - & -\\ 
                          & $\epsilon_c=0.06$ & - & - & -\\ \hline
 \end{tabular}
\caption{Parameters used in the QCD calculations \textsc{Pythia}, \textsc{Cascade}, \textsc{Herwig}, and MC@NLO. Here, $P_T$ stands for the transverse momentum, $m_1, \ m_2$ are the masses of the two outgoing partons from the hard process, and $m_Q$ denotes the heavy quark mass. The variable $\hat{s}$ represents the invariant mass of the $Q\bar{Q}$ system and $Q_{\bot}$ stands for its transverse momentum.}\label{table_qcdsettings}
\end{table}

\textsc{Pythia} and \textsc{Cascade} are used to simulate detector effects in order to determine the acceptance and the efficiency and to estimate the systematic uncertainties associated with the measurement. The generated events are passed through a detailed simulation of the detector response based on the GEANT simulation program \cite{Brun:1987ma} and are processed using the same reconstruction and analysis chain as is used for the data.  

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{H1 Detector}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Only a short description of the H1 detector is given here including the most relevant detector components for this analysis. A more complete description may be found elsewhere \cite{H1detector, Appuhn:1996na}. A right-handed coordinate system is employed at H1, with its origin at the nominal interaction vertex, its $z$-axis pointing in the proton beam direction and its $x(y)$ axis pointing in the horizontal (vertical) direction. Polar ($\theta$) and azimuthal ($\phi$) angles are measured with respect to this reference system. The pseudorapidity $\eta$ is related to the polar angle $\theta$ by $\eta=-\ln\tan(\theta/2)$.
\par
Charged particles are measured in the central tracking detector (CTD) with a transverse momentum resolution of $\sigma(P_T)/P_T\approx0.5\%P_T/\text{GeV} \oplus 1.5\%$ \cite{CTDResolution}. This device consists of two cylindrical drift chambers (CJC) interspersed with a drift chamber designed to improve the $z$-coordinate reconstruction. A multiwire proportional chamber mainly used for triggering is located in front of the inner CJC. The CTD is operated in a uniform solenoidal $1.16 \ \text{T}$ magnetic field, enabling the momentum measurement of charged particles over the polar angular range $20^{\circ}<\theta<160^{\circ}$. The efficiency for finding tracks in the CTD is greater than $99\%$.
\par
The CTD tracks are linked to hits in the vertex detector, the central silicon tracker (CST) \cite{CST}, to provide precise spatial track reconstruction. The CST consists of two layers of double-sided silicon strip detectors surrounding the beam pipe, covering an angular range of $30^{\circ}<\theta<150^{\circ}$ for tracks passing through both layers. The information of the $z$-coordinate of the CST hits is not used in the analysis presented in this paper. For CTD tracks with CST hits in both layers the transverse distance of closest approach (DCA) to the nominal vertex in $x-y$, averaged over the azimuthal angle, is measured to have a resolution of $43\ \mu$m $\oplus$ $51\ \mu$m$/(P_T/\text{GeV})$, where the first term represents the intrinsic resolution (including alignment uncertainty) and the second term is the contribution from multiple scattering in the beam pipe and the CST. The efficiency for linking hits in both layers of the CST to a CTD track is around $84\%$. 
\par
The track detectors are surrounded in the forward and central directions ($4^{\circ}<\theta<154^{\circ}$) by a finely grained liquid argon calorimeter (LAr) and in the backward region ($153^{\circ}<\theta<178^{\circ}$) by a lead-scintillating fibre calorimeter (SpaCal) both with electromagnetic and hadronic sections. These calorimeters provide energy and angular reconstruction for final state particles from the hadronic system. In the LAr electromagnetic shower energies are measured with a precision of $\sigma(E)/E=11\%/\sqrt{E/\text{GeV}}\oplus 1\%$ and hadronic energies with $\sigma(E)/E=50\%/\sqrt{E/\text{GeV}} \oplus 1\%$, as determined in test beam measurements. The energy resolution for electromagnetic showers in the SpaCal is $\sigma(E)/E=7.1\%/\sqrt{E/\text{GeV}}\oplus1\%$, as determined in test beam measurements \cite{SpaCalTestBeam}.
\par
The calorimeters are surrounded by the muon system. The central muon detector (CMD) is integrated in the iron return yoke of the superconducting coil and consists of $64$ modules, which are grouped in the forward endcap ($5^{\circ}\leq\theta\leq35^{\circ}$), the forward and backward barrel ($35^{\circ}\leq\theta\leq130^{\circ}$) and the backward endcap ($130^{\circ}\leq\theta\leq175^{\circ}$). Muon candidates are identified by requiring a geometrical matching of a CMD track segment with a CTD track.
\par
The luminosity determination is based on the measurement of the Bethe-Heitler process $ep\rightarrow ep\gamma$, where the photon is detected in a calorimeter located downstream of the interaction point in the electron beam direction at $z=-103\ \text{m}$.


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Experimental Method}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The data were collected with the H1 detector at the HERA collider during the years $2006$ and $2007$ and correspond to an integrated luminosity of $\mathcal{L}=179 \ \text{pb}^{-1}$. 
The beam energies were $E_e=27.6 \ \text{G\eV}$ and $E_p=920 \ \text{G\eV}$ for electrons and protons, respectively, resulting in a centre-of-mass energy of $\sqrt{s}\approx 320 \ \text{G\eV}$.  The  trigger requires a track segment in the muon system and track activity in the central jet chamber. A detailed account of the present analysis can be found in \cite{MiraKraemerThesis}. A summary of the kinematic range and the definition of the measurement is given in table \ref{CutSummaryTable}. 

\subsection{Photoproduction Event Selection}
Events in the photoproduction regime are selected by requiring that no isolated high energy electromagnetic cluster, consistent with a signal from a scattered electron, is detected in the LAr and SpaCal calorimeters. This limits the photon virtuality to values of \qsq$< 2.5 \ \gevsq$. 
The inelasticity $y$ is reconstructed using the relation $y=\sum_h(E-P_z)/2E_e$ \cite{JacquetBlondel}. Here, the sum includes all particles of the hadronic final state (HFS), while $E$ denotes their energies and $P_z$ stands for the $z$-components of their momenta. The HFS particles are reconstructed using a combination of tracks and calorimeter deposits in an energy flow algorithm that avoids double counting \cite{HFS}. The inelasticity in this analysis is restricted to $0.2 < y < 0.8$. 

\subsection{Muon Reconstruction and Selection}\label{subsec:MuonRec}
Muon candidates are identified as track segments in the barrel and endcap parts of the instrumented iron. The iron track segments must be well matched to a track reconstructed in the CTD. At least two CST hits in the $r-\phi$ plane have to be associated with the muon track. The combined CTD-CST track in $r-\phi$ is required to have a fit probability of at least $10\%$. The muon momentum is reconstructed using the CTD-CST track information. The CST hit requirements for the muon track restrict the allowed range of $ep$ interaction vertices along the $z$-axis to $|z_{vtx}|\leqslant 20 \ \text{cm}$. Events are selected with at least one muon candidate reconstructed in the instrumented iron having a pseudorapidity within $-1.3<\etamuon<1.5$ and a transverse momentum of $\ptmuon>2.5 \ \text{G\eV}$. If more than one muon candidate is found, the one with the highest transverse momentum is selected and other candidates are ignored. In $5.4\%$ of the events after the full selection more than one muon is found. 

\subsection{Jet Reconstruction and Selection}
Jets are reconstructed using the inclusive longitudinally invariant $k_T$ algorithm in the massless $P_T$ recombination scheme and with the distance parameter $R_0=1$ in the $\eta-\phi$ plane \cite{JetAlgo}. The algorithm is applied in the laboratory frame using all reconstructed HFS particles including the muon candidate. A jet is defined as a \muonjet\  if the selected muon candidate lies within a cone of radius 1 about the jet axis in the $\eta-\phi$ plane. The efficiency for this matching amounts to about $90\%$ and is consistent for data and MCs.  
The jet with the highest $P_T$ is referred to as jet1, while the second highest is called jet2. Events with at least two jets are selected, where the leading two jets are required to be in the angular range $-1.5<\etajet<2.5$ and to have a transverse momentum of $\ptjetone>7\ \text{G\eV}$ and $\ptjettwo>6\ \text{G\eV}$. One of the two selected jets must be classified as a \muonjet. The Monte Carlo simulation is used to define hadron level jets, which consist of stable particles including neutrinos, but excluding the scattered electron, before they are passed through the simulation of the detector response. 


\subsection{Separation of Direct and Resolved Processes}
The fraction of the photon energy entering the hard interaction is estimated using the observable $\xgamma$:
\begin{equation}\nonumber
\xgamma=\frac{\sum_{jet1}(E-P_z)+\sum_{jet2}(E-P_z)}{\sum_{HFS}(E-P_z)},
\end{equation}
where the sums in the numerator run over the particles associated with the two jets and the one in the denominator over all detected hadronic final state particles. For direct processes \xgamma\  approaches unity, as the hadronic final state consists only of the two hard jets selected in the present analysis and the proton remnant in the forward region only has a minor contribution to $\sum_{HFS}(E-P_z)$. In resolved processes \xgamma\  can have smaller values. 

\begin{table}
\centering
\begin{tabular}{|c|c|}
\hline 
& Photoproduction of $b(c)\rightarrow \mu jj X$  \\
\hline
\hline 
Kinematic range &  \qsq$< 2.5 \ \gevsq$  \\
 & $0.2 < y < 0.8$ \\
\hline
Event selection & \ptmuon $>2.5 \ \text{G\eV}$ \\
& $-1.3<\etamuon<1.5$\\
& \ptjetone $> 7 \ \text{G\eV}$ \\
& \ptjettwo $>6 \ \text{G\eV}$\\
& $-1.5<\etajet<2.5$ \\
\hline
Event sample & $N_{events}=6807$ \\
& $\mathcal{L}=179 \ \text{pb}^{-1}$ \\
\hline
\end{tabular}
\caption{Definition of the kinematic range of the measurement and event yield for the data sample collected in the years $2006$ and $2007$. The variables are measured in the laboratory frame.}\label{CutSummaryTable}
\end{table}

\subsection{Flavour Separation}
The flavour of an event is defined as the hadron flavour of the $\muonjet$. The measured cross sections are proportional to the rate of events with a muon and a dijet system rather than the rate of muons or jets. 
\par
The separation of $b$, $c$ and light quark ($uds$) events is only briefly described here. The procedure closely follows that described in \cite{H1Hera1desy05004}. The separation is performed using the properties of the muon track associated to the \muonjet. The impact parameter $\delta$ of a muon track is the transverse distance of closest approach (DCA) of the muon track to the beam spot point, which is the position of the beam interaction region in $x$ and $y$. The beam spot is derived from tracks with CST hits averaging over many events and is updated regularly to account for drifts during beam storage. Muon tracks with $\delta > 0.1 \ \text{cm}$ are rejected to suppress contributions from the decays of long-lived strange particles. If the angle $\alpha$ between the azimuthal angle of the \muonjet \ and the line joining the primary vertex to the point of  muon closest approach is less than $90^{\circ}$, the impact parameter is defined as positive. It is defined as negative otherwise. The transverse momentum \ptrel \ of the muon relative to the \muonjet \ axis is also sensitive to the quark content of the event sample and used together with the impact parameter for the flavour separation.
\par
The fractions of events with beauty, charm and light quarks, $f_{b}$, $f_{c}$ and $f_{l}$, are obtained by a binned likelihood fit \cite{BarlowFit} in the $\impactpar - \ptrel$ plane. Following \cite{BarlowFit}, a likelihood ratio is calculated based on Poisson statistics. The fit is performed separately for each individual measurement bin, while the total cross sections are determined using the fractions obtained from a fit to the complete event sample. The $uds$ (light), $c$ and $b$ \textsc{Pythia} Monte Carlo simulation samples are used as templates. Only the statistical errors of the data and the Monte Carlo simulations are considered in the fit. As a cross check, all fits are also performed using one-dimensional distributions of \ptrel\ and \impactpar\ separately. These two one-dimensional fits give a compatible $b$ fraction in all measurement bins within the statistical error. The one-dimensional \ptrel\ fit does not allow a determination of the charm fraction.
\par
The two distributions that are used in the flavour separation are shown in figure \ref{fig:plotsflavourseparation}. The distribution of the impact parameter \impactpar \ shown in figure~\ref{fig:plotsflavourseparation} (a) is symmetric for $uds$ events while $b$ and $c$ events contribute more at large positive values of \impactpar. Therefore, the fit of this variable can distinguish the three different quark contributions. As shown in figure \ref{fig:plotsflavourseparation} (a), the sum of all three fitted contributions in the Monte Carlo simulation is able to describe the data quite well. This description is achieved by a better understanding of the detector and an improved detector simulation with regard to signal heights, noise levels and dead strips in the CST \cite{CSTTuning}, the inclusion of effects from alignment imperfections, and the description of the dead material in front of the CST and CJC. Therefore, a further smearing of measured track parameters in the simulation is not necessary, as it was done in the former H1 analysis using HERA I data \cite{H1Hera1desy05004}.
The \ptrel \ distribution is shown in figure \ref{fig:plotsflavourseparation} (b). The $uds$ and $c$ distributions are very similar and peak at low values of \ptrel, while the $b$ events contribute more at higher values of \ptrel. The sum of all quark contributions in the Monte Carlo simulation is able to describe the data reasonably well.
\par
The fitted parameters $f_{b}$, $f_{c}$ and $f_{l}$ for the whole kinematic range are:
\begin{align}\nonumber
f_b&=(26.0\pm1.2)\%,\\ \nonumber
f_c&=(48.6\pm2.5)\%,\\ \nonumber
f_l&=(25.3\pm2.6)\%. \nonumber
\end{align}
The $\chi^2/ndf$ is found to be $0.76$ for the total sample. 
\par
Control distributions for the data sample in comparison to the Monte Carlo simulations are shown in figure \ref{fig:controlplots}. All selection cuts are applied. The data are compared with the MC contributions from beauty, charm and light quark events and their sum with the relative fractions taken from the fit as discussed above. The number of events in the simulation is normalised to the one of the data. It is observed that the shapes of the MC contributions are rather similar for beauty, charm and light quark events for the distributions shown here. For the determination of the detector corrections a reweighting of \dphijets\  and the transverse momentum of the leading selected track in the event is performed on hadron level in the Monte Carlo simulation to provide a better description of the data. Only small deviations between data and MC are observed, such as in the forward \etamuon \ region in figure \ref{fig:controlplots} (b) and in the very small \dphijets \ region in \ref{fig:controlplots} (d). 

\subsection{Cross Section Determination}
Total and differential visible beauty and charm $ep$ cross sections are measured in the photoproduction regime. The fitted fractions $f_b$ and $f_c$ are converted to cross sections $\sigma_{b(c)}$ in each bin using

\begin{equation}
\sigma_{b(c)}=\frac{f_{b(c)}N_{Data}N_{b(c)}^{MCgen}}{\mathcal{L}N_{b(c)}^{MCrec}}.
\end{equation}
Here, $N_{Data}$ and $N_{b(c)}^{MCrec}$ represent the number of data or Monte Carlo simulation events passing all selection cuts on reconstruction level. The variable $N_{b(c)}^{MCgen}$ stands for the events selected in the Monte Carlo simulation on hadron level and $\mathcal{L}$ denotes the luminosity of the data. The differential cross sections are obtained by dividing by the bin width.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Fake Muon Rate}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
All backgrounds to semi-muonic $b$ and $c$ decays are called {\it fake muons} here. 
These contributions are modelled using MC simulation and originate mainly
from $uds$ events, with a small fraction from $b$ and $c$ events.  
Three sources of fake muons are considered:
\begin{itemize}
\item{Hadrons which reach the muon detector and are misidentified as muons. According to the fully inclusive \textsc{Pythia} MC, $1.6\%$ of the selected $b$ events and $2.3\%$ of the selected $c$ events originate from hadron misidentification. } 
\item{Muons which do not originate from $b$ or $c$ hadron decays, but from other hadrons such as kaons and pions. This background source is denoted as {\it inflight decay} in the following. 
According to the full inclusive \textsc{Pythia} MC, $0.9\%$ of the $b$ and $0.7\%$ of the selected $c$ events originate from inflight decays.}
\item{Cosmic ray muons which coincide in time with real $ep$ events. About $1\%$ of the selected muon candidates are 
rejected as cosmic ray muons based on timing information from the CTD \cite{MiraKraemerThesis}. The remaining background from cosmic ray muons is negligible.}
\end{itemize}
The probability to fake a muon depends on the particle species.  The kinematic distributions of different particle species differ, but the $uds$ MC is used as one single template. Therefore the fake probabilities for the most important particle species are studied in data and MC.
The fake muon probabilities from misidentification and inflight decays $P^{\mu}_{h}$ for hadron $h$ are defined as: 
\begin{equation}\nonumber
P^{\mu}_h=\frac{\# \ \text{fake muons}}{\# \ \text{all hadrons}}.
\end{equation}
The fake muon probability is investigated for pions originating from $K^0_S\rightarrow \pi^+\pi^-$ decays and amounts to $P^{\mu}_{\pi}\approx0.05$ in the data. For protons the decay channels $\Lambda\rightarrow p\pi^-$ and $\Lambda\rightarrow \bar{p}\pi^+$ are used and the fake muon probability is found to be $P^{\mu}_{p}\approx0.04$ in the data. For $K^{\pm}$ mesons from the decay $D^{*\pm}\rightarrow D^0\pi_{slow}^{\pm}\rightarrow (K^{\mp}\pi^{\pm})\pi^{\pm}_{slow}$ a fake muon probability of $P^{\mu}_{K{\pm}}\approx0.01$ is measured in the data. It is observed that the pion and proton fake muon probabilities in the data are not described by the Monte Carlo simulation. They are reweighted in the Monte Carlo simulation by factors of $2.0$ and $1.9$, respectively, to match the data. The $K^{\pm}$ fake muon probability in the data and the simulation are in agreement. The misidentified muon events remain in the event sample.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Systematic Uncertainties}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The following uncertainties are taken into account in order to evaluate the systematic errors.  
\begin{itemize}
\item{The trigger efficiencies are determined using independent trigger channels in the DIS regime since no independent triggers exist in photoproduction. The uncertainty is estimated by the difference between the efficiency found in the data and the simulation. The Monte Carlo simulation is reweighted to match the data trigger efficiency.}
\item{The efficiency for the identification of the muons is determined using a high statistics sample of events of elastically produced $J/\psi$ mesons \cite{MartinGoettlichThesis}. The efficiency is known to a precision of $4\%$.}
\item{The track efficiency of the CTD is known to $\pm1\%$ and that of the CST to $\pm2\%$. The uncertainty due to the track efficiencies is estimated by varying the efficiencies of the CTD and CST correspondingly.}
\item{The integrated luminosity is known to a precision of $4\%$.}
\item{The uncertainty due to the resolution of the impact parameter $\delta$ of the muon tracks is estimated by varying the resolution by an amount that encompasses any difference between the data and the simulation. This is achieved by applying an additional Gaussian smearing in the Monte Carlo simulation of $200 \ \mu\text{m}$ to $5\%$ of randomly selected tracks and $12\ \mu\text{m}$ to the rest.}
\item{The uncertainty on the cross section arising from the uncertainty on the reconstruction of $\phi_{jet}$ is estimated by shifting its value by $\pm2^{\circ}$. }
\item{The uncertainty arising from the hadronic energy scale is estimated by changing it by $\pm1\%$ for the complete hadronic final state.}
\item{The dependence of the measurement on the physics model used for the templates representing different QCD evolution schemes is estimated by replacing the \textsc{Pythia} $b$ and $c$ Monte Carlo templates with \textsc{Cascade}.}
\item{The uncertainty on the cross section arising from the uncertainty of the parton fragmentation model is estimated by replacing the \textsc{Pythia} Monte Carlo samples using the Lund fragmentation function with samples based on the Peterson fragmentation function. }
\item{The uncertainty arising from fake muon background is estimated by not applying the weights that have been found in the fake rate probabilities for $K$ and $\Lambda$ decays.}
\item{The impact of the reweighting on the cross sections are investigated and found to be negligible.}

\end{itemize}
\par
The effect of the listed experimental uncertainties are estimated by varying the relevant variables in the Monte Carlo simulation or by modifying the corresponding efficiencies in the cross section calculation. The difference between the obtained cross sections with and without the change result in the measurement systematic uncertainties, which are summarised in table \ref{table_syserr}.
The individual effects of the above experimental uncertainties are combined in quadrature, yielding a total uncertainty of $10.5\%$ and $10.4\%$ on the measured $b$ and $c$ cross sections, respectively. The systematic uncertainties as derived from the integrated sample are applied to each analysis bin in order to avoid statistical fluctuations. The largest contribution to this uncertainty for the $c$ measurement arises from systematics attributed to the hadronic energy scale ($5\%$). The systematic errors of the $b$ analysis are not dominated by a single source. 
\begin{table}
\centering
\begin{tabular}{|l|c|c|}
 \hline
 Systematic error source & Beauty $\Delta\sigma/\sigma[\%]$& Charm $\Delta\sigma/\sigma[\%]$\\ \hline 
 Trigger efficiency &  $4$ & $4$ \\ 
 Muon identification & $4$ & $4$ \\
 Track finding efficiency & $3$ & $3$\\
 Luminosity & $4$ & $4$ \\ 
 $\delta$ Resolution & $3$ & $2$  \\ 
Jet axis & $4$  & $2$ \\
 Hadronic energy scale & $3$ & $5$ \\  
 Physics model &  $3$ &  $1$ \\
 Fragmentation &  $3$ & $4$ \\
 Fake muon background & $1$ & $1$ \\ \hline \hline
Total & $10.5$  & $10.4$ \\  \hline
\end{tabular}
\caption{Summary of the systematic uncertainties of the beauty and charm cross sections.}\label{table_syserr}
\end{table}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Results}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The cross sections for $b$ and $c$ in photoproduction using semi-muonic decays in dijet events are measured. The cross sections are determined for the phase space defined by the kinematic range and the event selection cuts presented in table \ref{CutSummaryTable}. The measured cross sections are compared to the expectations of the MC programs {\sc Pythia}, {\sc Cascade}, {\sc Herwig}, and MC@NLO. The total measured and predicted cross sections are listed in table \ref{table_totalresults}. 
{\sc Pythia} shows the highest normalisation of the three LO MCs, while the normalisation of {\sc Cascade} is below the one of {\sc Pythia} and {\sc Herwig} has the lowest normalisation. For the beauty measurement, the {\sc Pythia} prediction is closest to the data and gives the best description of the three LO MC predictions.
The beauty and charm data cross sections tend to be underestimated by the MC@NLO predictions but are in agreement within the errors. The precision of the measured cross sections are much higher than the ones of the theory predictions shown here.
\par
The beauty and charm cross sections are measured differentially as a function of the transverse momentum of the leading jet \ptjetone\ and of the muon \ptmuon, the pseudorapidity of the muon \etamuon, the momentum fraction \xgamma\  carried by the photon entering the hard interaction and the azimuthal angular difference \dphijets\  between the two leading jets. The measurements are performed for the full sample, as well as for direct and resolved enriched processes separately. The distinction is performed by the variable $\xgamma$, which leads to enriched resolved processes in the region $\xgamma\leq0.75$ and direct photon enriched processes for $\xgamma>0.75$. The cross sections for the beauty measurements are given in tables \ref{tab:beautyfull}-\ref{tab:beautyresolved} and shown in figures \ref{fig:beautyxsecfullsample}-\ref{fig:beautyxsecresolvedsample}. 
\begin{table}
\centering
\begin{tabular}{|c|c|c|}
\hline
                   & $\sigma_{\text{vis}}(ep\rightarrow eb\bar{b}X\rightarrow ejj\mu X')$ [pb]  & $\sigma_{\text{vis}}(ep\rightarrow ec\bar{c}X\rightarrow ejj\mu X')$ [pb]  \\ \hline  \hline
H1 Data  & $43.3\pm 2.1\ (\text{stat.})\pm 4.5\ (\text{sys.})$ & $81.3\pm 4.3\ (\text{stat.})\pm 8.5\ (\text{sys.})$ \\ \hline
\textsc{Pythia} & $35.3$  & $94.3$ \\
\textsc{Cascade} &  $29.0$ & $76.8$ \\
\textsc{Herwig}  & $20.6$ & $58.5$ \\
MC@NLO &  $33.4^{+7.1}_{-9.2}$   &  $58.6^{+29.5}_{-11.2}$   \\ \hline
\end{tabular}
\caption{Total visible measured beauty and charm cross sections along with their statistical and systematic errors. The total predictions from \textsc{Pythia}, \textsc{Cascade}, \textsc{Herwig}, and MC@NLO are also shown. The MC@NLO predictions are given with their theoretical uncertainties. }\label{table_totalresults}
\end{table}
\par
In the case of beauty production the models provide a good description of the measured cross sections in terms of shape in all distributions.  For the LO MCs this is true for the full sample, as well as for the direct and resolved enriched regions. 
The cross sections for direct enriched processes are well described in shape, but tend to be underestimated by MC@NLO, while for resolved enriched processes a reasonable agreement is observed both in shape and normalisation.
\par
In the analysis of semi-muonic $b$ decays in dijet events with HERA I data \cite{H1Hera1desy05004} an excess of data compared to the NLO predictions of the FMNR program was observed in the first bin of \ptmuon\ and \ptjetone. In this analysis the NLO predictions are provided by MC@NLO which is based on the FMNR parton level calculations. 
%The predictions from FMNR are found to agree with those of MC@NLO within the theoretical uncertainties. The excess in the first bin of  \ptmuon\ and \ptjetone, which was seen in the HERA I analysis, is not confirmed. Especially the first bin of \ptjetone\ does not show such an excess in the full sample. 
Also in this analysis, the NLO predictions lie below the data in the
first bin of \ptmuon\ and \ptjetone, but they are consistent with the data
within $2 \sigma$ of the experimental and theoretical uncertainty.
Whereas in enriched direct processes, the data tend to be underestimated by the MC@NLO prediction in the first bin of \ptmuon \ and the first bin of \ptjetone, in enriched resolved processes, no such effect is visible. 
\par
The cross sections as a function \dphijets \ show a significant contribution away from the back-to-back configuration at \dphijets$\simeq 180^{\circ}$. Such a configuration can be described by models which include significant contributions from higher order QCD radiation or a transverse momentum of the gluon in the initial state. This distribution is reasonably well described by all models. In direct and resolved enriched processes, this observation also holds.  
\par
The measured charm cross sections are presented in tables \ref{tab:charmfull}-\ref{tab:charmresolved} and figures \ref{fig:charmxsecfullsample}-\ref{fig:charmxsecresolvedsample}. The distributions are reasonably well described by all models.  
Similar to recent observations in H1 measurements of the photoproduction of $D^*$ mesons \cite{2011DStar}, the central value of the MC@NLO calculations tend to be lower than the measured charm cross sections. 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Conclusions}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Beauty and charm photoproduction cross sections for events with dijets and a muon are measured using the data collected by the H1 detector at HERA. Compared to the previous H1 beauty measurement \cite{H1Hera1desy05004}, the analysis profits from a three times larger luminosity of the data sample, an extended phase space as well as improved understanding of the H1 vertex detector.  The flavour composition of the event sample is determined by the transverse momentum of the muon relative to the jet axis of its associated jet \ptrel\ and by its impact parameter \impactpar. Total visible  and differential cross sections are measured, and the results are compared to leading order QCD models provided by {\sc Pythia}, {\sc Cascade}, and {\sc Herwig}, as well as to the next-to-leading order  calculations provided by MC@NLO. At low values of \ptmuon\ and \ptjetone, the present beauty measurement does not show a significant excess as observed by the previous H1 measurement with respect to the NLO calculation. In general the predictions are in reasonable agreement with the beauty and charm measurements.  


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

We are grateful to the HERA machine group whose outstanding
efforts have made this experiment possible. 
We thank the engineers and technicians for their work in constructing and
maintaining the H1 detector, our funding agencies for 
financial support, the
DESY technical staff for continual assistance
and the DESY directorate for support and for the
hospitality which they extend to the non-DESY 
members of the collaboration. We would like to thank T. Toll for fruitful discussions and for his help with the MC@NLO predictions.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{thebibliography}{99}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\bibitem{mcatnlo}
T.~Toll and S.~Frixione, Phys. Lett. {\bf B703} (2011) 452 [arXiv:1106.1614].

\bibitem{ZeusLifetime}
S.~Chekanov {\it et al.} [ZEUS Collaboration], Eur. Phys. J. {\bf C63} (2009) 171 [arXiv:0812.3775].

\bibitem{H1Lifetime2010}
F.D.~Aaron {\it et al.} [H1 Collaboration], Eur. Phys. J. {\bf C65} (2010) 89 [arXiv:0907.2643]. 

\bibitem{H1Lifetime2011}
F.D.~Aaron {\it et al.} [H1 Collaboration], Eur. Phys. J. {\bf C71} (2011) 1509 [arXiv:1008.1731].

\bibitem{H1Hera1desy05004}
A.~Aktas {\it et al.}  [H1 Collaboration],
Eur. Phys. J. {\bf C41} (2005) 453 [hep-ex/052010].

\bibitem{ZEUSHera1desy03212}
S.~Chekanov {\it et al.} [ZEUS Collaboration],
Phys. Rev. {\bf D70} (2004) 012008 [hep-ex/0312057]; \\
Erratum.ibid. {\bf D74} (2006) 059906.

\bibitem{ZEUSdesy08210}
S.~Chekanov {\it et al.} [ZEUS Collaboration], JHEP {\bf 0904} (2009) 133 [arXiv:0901.2226].

\bibitem{ZeusSemileptonic}
S.~Chekanov {\it et al.} [ZEUS Collaboration], Eur. Phys. J. {\bf C65} (2010) 65 [arXiv:0904.3487].

\bibitem{H1D}
F.D.~Aaron {\it et al.} [H1 Collaboration], Phys. Lett. {\bf B686} (2010) 91 [arXiv:0911.3989]. 

\bibitem{ZeusD}
H.~Abramowicz {\it et al.} [ZEUS Collaboration], JHEP {\bf 1011} (2010) 009 [arXiv:1007.1945].


\bibitem{CSTTuning}
M.~Kr\"amer, Proceedings of 'International School of Subnuclear Physics, 45th course: Search for the Ôtotally unexpected
in the LHC era', Erice, Sicily (2007).


\bibitem{DGLAP} 
V.N.~Gribov and L.N.~Lipatov, 
Yad. Fiz. {\bf 15} (1972) 781 [Sov.J.Nucl.Phys. {\bf 15} (1972) 438]; \\
G.~Altarelli and G.~Parisi, 
Nucl. \ Phys. {\bf B126} (1977) 298;\\
Y.L.~Dokshitzer, 
Sov. Phys. JETP {\bf 46} (1977) 641 [Zh. Eksp. Teor. Fiz. {\bf 73} 1216].

\bibitem{CCFM} 
M.~Ciafaloni,
Nucl. Phys. {\bf B296} (1988) 49; \\
S.~Catani, F.~Fiorani, and G.~Marchesini,
Phys. Lett. {\bf B234} (1990) 339; \\
G.~Marchesini, 
Nucl. Phys. {\bf B445} (1995) 49. 

\bibitem{uPDFs}
J.~Collins and H.~Jung, Proceedings of 'HERA and the LHC', eds. H.~Jung and A.~De Roeck, CERN-2005-014, DESY-PROC-2005-001 (2005) 261 [hep-ph/0508280].

\bibitem{schemes}
W.K.~Tung, S.~Kretzer, and C.~Schmidt, J. Phys. {\bf G28} (2002) 983 [hep-ph/0110247].



\bibitem{Pythia} 
T.~Sj\"ostrand, S.~Mrenna, and P.Z.~Skands,
%\textsc{Pythia} 6.4 Physics and Manual, 
JHEP {\bf 0605} (2006) 026 [hep-ph/0603175].

\bibitem{LundFrag}
B.~Andersson, G.~Gustafson, G.~Ingelmann, and T.~Sj\"ostrand,
Phys. Rept. {\bf 97} (1983) 31.

\bibitem{PeteFrag}
C.~Peterson, D.~Schlatter, I.~Schmitt, and P.M.~Zerwas,
Phys. Rev. {\bf D27} (1983) 105.       

\bibitem{Cascade} 
H.~Jung {\it et al}, \textsc{Cascade} 2.2, Eur. Phys. J. {\bf C70} (2010) 1237 [arXiv:1008.0152].

\bibitem{CCFMSetA0}
H.~Jung, Proceedings of 'XII International Workshop in Deep-Inelastic Scattering (DIS 2004)', eds. D.~Bruncko, J.~Ferencei, and P.~Str\'{\i}\v{z}enec, (2004) 299 [hep-ph/0411287].

\bibitem{herwig}
G.~Corcella {\it et al.}, JHEP {\bf 0101} (2001) 010 [hep-ph/0011363].

\bibitem{cluster}
B.R.~Webber, Nucl. Phys. {\bf B238} (1984) 492; \\
G.~Marchesini and B.R.~Webber, Nucl. Phys. {\bf B310} (1988) 461.

\bibitem{fmnr}
S.~Frixione, P.~Nason, and G.~Ridolfi, Nucl. Phys. {\bf B454} (1995) 3 [hep-ph/9506226];\\
S.~Frixione, M.L.~Mangano, P.~Nason, and G.~Ridolfi, Phys. Lett. {\bf B348} (1995) 633 [hep-ph/9412348].


\bibitem{mcatnlohad}
S.~Frixione and B.R.~Webber, JHEP {\bf 0206} (2002) 029 [hep-ph/0204244]; \\
S.~Frixione, P.~Nason, and B.R.~Webber, JHEP {\bf 0308} (2003) 007 [hep-ph/0305252].

\bibitem{Brun:1987ma}
\textsc{Geant} 3, R.~Brun  {\it et al.}, 
%, F.~Bruyant, M.~Maire, A.~C.~McPherson and P.~Zanarini,
%``Geant3,''
CERN-DD/EE/84-1 (1987).
%\href{http://www.slac.stanford.edu/spires/find/hep/www?r=cern-dd\%2Fee\%2F84-1}{SPIRES entry}   



\bibitem{CTEQ6}
J.~Pumplin {\it et al},
JHEP {\bf 0207} (2002) 012 [hep-ph/0201195].  

\bibitem{Sas2d}
G.A.~Schuler and T.~Sj\"ostrand, Z. Phys. {\bf C68} (1995) 607.     


\bibitem{HERAPDF}
F.D.~Aaron {\it et al.} [H1 and ZEUS Collaboration], JHEP {\bf 1001} (2010) 109 [arXiv:0911.0884].

\bibitem{GRVGLO}
M.~Gl\"uck, E.~Reya, and A.~Vogt, Phys. Rev. {\bf D45} (1992) 3986; \\
M.~Gl\"uck, E.~Reya, and A.~Vogt, Phys. Rev. {\bf D46} (1992) 1973.

\bibitem{H1detector}
I.~Abt {\it et al.}  [H1 Collaboration],
%``The H1 Detector At Hera,''
Nucl.\ Instrum.\ Meth.\ {\bf A386} (1997) 310 and 348.

\bibitem{Appuhn:1996na}
R.D.~Appuhn {\it et al.} [H1 SpaCal Group],
%``The H1 lead/scintillating-fibre calorimeter,''
Nucl.\ Instrum.\ Meth.\ {\bf A386} (1997) 397.
%%CITATION = NUIMA,A386,397;%%

\bibitem{CTDResolution}
C.~Kleinwort, Proceedings of 'First LHC Detector
Alignment Workshop', eds. S.~Blusk {\it et al.}, CERN-2007-004 (2007) 41.

\bibitem{CST}
D.~Pitzl {\it et al.},
Nucl.\ Instrum.\ Meth.\ {\bf A454} (2000) 334 [hep-ex/0002044];\\
B.~List,
Nucl.\ Instrum.\ Meth.\ {\bf A501} (2001) 49.

\bibitem{SpaCalTestBeam}
T.~Nicholls {\it et al.} [H1 SpaCal Group], Nucl. Instrum. Methods {\bf A374} (1996) 149.

\bibitem{MiraKraemerThesis}
M.~Kr\"amer, 'Measurement of Charm and Beauty Cross Sections in Photoproduction Using Events with Muons and Dijets at HERA', PhD thesis, Univ. Hamburg (2009) DESY-THESIS-2009-030, \\ available at 
http://www-h1.desy.de/publications/theses$\_$list.html.

\bibitem{JacquetBlondel}
A.~Blondel and F.~Jacquet, Proceedings of 'Study of an ep Facility for Europe', 
ed. U.~Amaldi, DESY 79/48 (1979) 391.

\bibitem{HFS}
M.~Peez, 'Recherche de d\'eviations au Mod\`ele Standard dans les processus de grande \'energie transverse sur le collisionneur \'electron-proton HERA', PhD thesis (in French), Universit\'e de Lyon (2003), DESY-THESIS-2003-023 \\
available at http://www-h1.desy.de/publications/theses$\_$list.html;\\
S.~Hellwig, 'Untersuchung der $D^*$-$\pi_{slow}$ Double Tagging Methode in Charmanalysen', Dipl. thesis (in German), Univ. Hamburg (2004) \\
available at http://www-h1.desy.de/publications/theses$\_$list.html.

\bibitem{JetAlgo}
S.D.~Ellis and D.E.~Soper, Phys. \ Rev. {\bf D48} (1993) 3160 [hep-ph/9305266];\\
S.~Catani, Y.L.~Dokshitzer, M.H.~Seymour, and B.R.~Webber, Nucl. \ Phys. \ {\bf B406} (1993) 187.

\bibitem{BarlowFit}
R.J.~Barlow and C.~Beeston,
%``Fitting using finite Monte Carlo samples,''
Comput.\ Phys.\ Commun.\  {\bf 77} (1993) 219.
%%CITATION = CPHCB,77,219;%%

\bibitem{MartinGoettlichThesis}
M.~G\"ottlich, 'Study of Charm and Beauty Production at HERA/H1 using Dilepton Events', PhD thesis, Univ. Hamburg (2007), DESY-THESIS-2007-012, \\ available at 
http://www-h1.desy.de/publications/theses$\_$list.html.

\bibitem{2011DStar}
F.D.~Aaron {\it et al.} [H1 Collaboration], Eur. Phys. J. {\bf C72} (2012) 1995 [arXiv:1203.1170]. 

    
\end{thebibliography}

\clearpage


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

\renewcommand{\arraystretch}{1.19} % this increase the table in vertical direction, very useful!
\begin{table}[tb]
  \begin{center}
    \begin{tabular}{|rr|rrr||c|}
      \hline
      \multicolumn{6}{|c|}{\bf \boldmath H1 Beauty Dijet Muon Cross Sections}   \\
      \multicolumn{6}{|c|}{$\sigma_{vis}(ep\rightarrow eb\bar{b}X\rightarrow ejj\mu X)$}   \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\ptmuon\ range} & \dsdx{\ptmuon} & stat. & sys. & $f_{b}\pm$stat.  \\  
      \multicolumn{2}{|c|}{[G\eV]} & \multicolumn{3}{c||}{[pb/G\eV]} &  \\ 
      \hline
      $2.5$ & $3.3$ & $25.3$ & $2.2$ & $2.8$ & $0.21\pm0.02$ \\ 
     $3.3$ & $4.7$ & $12.0$ & $0.9$ & $1.3$ & $0.28\pm0.02$\\ 
     $4.7$ & $15.0$ & $0.772$ & $0.078$ & $0.085$ & $0.31\pm0.03$ \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\etamuon\ range} & \dsdx{\etamuon} & stat. & sys. & $f_{b}\pm$stat.   \\  
      \multicolumn{2}{|c|}{} & \multicolumn{3}{c||}{[pb]} &  \\ 
      \hline
      $-1.3$ & $-0.3$ & $9.9$ & $1.2$ & $1.1$ & $0.23\pm0.03$ \\ 
     $-0.3$ & $0.0$ & $19.6$ & $2.2$ & $2.2$ & $0.24\pm0.03$  \\ 
     $0.0$ & $0.3$ & $21.7$ & $2.3$ & $2.4$ & $0.25\pm0.03$ \\ 
     $0.3$ & $0.6$ & $23.4$ & $2.4$ &$2.6$ & $0.30\pm0.03$ \\ 
     $0.6$ & $1.5$ & $14.3$ & $1.4$ & $1.6$ & $0.26\pm0.03$ \\ 
      \hline
      \hline
      \multicolumn{2}{|c|}{\ptjetone\ range} & \dsdx{\ptjetone} & stat. & sys. & $f_{b}\pm$stat.  \\  
      \multicolumn{2}{|c|}{[G\eV]} & \multicolumn{3}{c||}{[pb/G\eV]} &  \\ 
      \hline
     $7$ & $11$ & $4.53$ & $0.41$ & $0.50$ & $0.24\pm0.02$  \\ 
    $11$ & $15$& $3.40$ & $0.27$ & $0.37$ & $0.25\pm0.02$ \\ 
    $15$ & $38$ & $0.469$ & $0.038$ & $0.052$ & $0.30\pm0.02$  \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\dphijets\ range} & \dsdx{\dphijets} & stat. & sys. & $f_{b}\pm$stat.  \\  
      \multicolumn{2}{|c|}{[deg]} & \multicolumn{3}{c||}{[pb/deg]} &  \\ 
      \hline
   $0$ & $155$ & $0.0576$ & $0.0063$ & $0.0061$ & $0.23\pm0.03$ \\ 
$155$ & $173$ & $1.01$ & $0.07$ & $0.11$ & $0.26\pm0.02$  \\ 
$173$ & $180$ & $2.17$ & $0.17$ & $0.24$ & $0.26\pm0.02$ \\ 
      \hline
      \hline
      \multicolumn{2}{|c|}{\xgamma\ range} & \dsdx{\xgamma} & stat. & sys. & $f_{b}\pm$stat.  \\  
      \multicolumn{2}{|c|}{} & \multicolumn{3}{c||}{[pb]} &  \\ 
      \hline
$0.0$ & $0.4$ & $10.6$ & $2.7$ & $1.2$ & $0.17\pm0.04$  \\ 
$0.4 $ & $ 0.75 $ & $ 35.6 $ & $ 3.4 $ & $ 3.9 $ & $ 0.23\pm0.02$   \\ 
$0.75 $ & $ 1.0 $ & $ 103.5 $ & $ 5.9 $ & $ 11.4 $ & $ 0.29\pm0.02$ \\ 
      \hline
    \end{tabular}
    \caption{Bin averaged differential cross sections for beauty photoproduction of dijet events using semi-muonic decays in bins of \ptmuon, \etamuon, \ptjetone, \dphijets, and \xgamma\  with their statistical and systematic uncertainties. The fit parameter $f_b$ is given including its statistical error. }
    \label{tab:beautyfull}
  \end{center}
\end{table}


\begin{table}[tb]
  \begin{center}
    \begin{tabular}{|rr|rrr||c|}
      \hline
      \multicolumn{6}{|c|}{\bf \boldmath H1 Beauty Dijet Muon Cross Sections, \bf \boldmath \xgamma$>0.75$}   \\
      \multicolumn{6}{|c|}{$\sigma_{vis}(ep\rightarrow eb\bar{b}X\rightarrow ejj\mu X)$}   \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\ptmuon\ range} & \dsdx{\ptmuon} & stat. & sys. & $f_{b}\pm$stat.  \\  
      \multicolumn{2}{|c|}{[G\eV]} & \multicolumn{3}{c||}{[pb/G\eV]} &  \\ 
      \hline
     $ 2.5 $ & $ 3.3 $ & $ 14.8 $ & $ 1.5 $ & $ 1.6 $ & $ 0.24\pm0.02 $  \\ 
$3.3 $ & $ 4.7 $ & $ 6.7 $ & $ 0.6 $ & $ 0.7 $ & $ 0.28\pm0.03 $  \\ 
$4.7 $ & $ 15.0 $ & $ 0.487 $ & $ 0.068 $ & $ 0.054 $ & $ 0.31\pm0.04 $ \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\etamuon\ range} & \dsdx{\etamuon} & stat. & sys. & $f_{b}\pm$stat.  \\  
      \multicolumn{2}{|c|}{} & \multicolumn{3}{c||}{[pb]} &  \\ 
      \hline
      $-1.3 $ & $ -0.3 $ & $ 6.5 $ & $ 0.9 $ & $ 0.7 $ & $0.24\pm0.03$  \\ 
$-0.3 $ & $ 0.0 $ & $ 11.3 $ & $ 1.6 $ & $ 1.2 $ & $0.24\pm0.03$ \\ 
$0.0 $ & $ 0.3 $ & $ 14.4 $ & $ 1.8 $ & $ 1.6 $ & $0.29\pm0.03$ \\ 
$0.3 $ & $ 0.6 $ & $ 13.9 $ & $ 1.7 $ & $ 1.5 $ & $0.32\pm0.04$  \\ 
$0.6 $ & $ 1.5 $ & $ 6.1 $ & $ 1.0 $ & $ 0.7 $ & $0.25\pm0.04$ \\ 
      \hline
      \hline
      \multicolumn{2}{|c|}{\ptjetone\ range} & \dsdx{\ptjetone} & stat. & sys. & $f_{b}\pm$stat.   \\  
      \multicolumn{2}{|c|}{[G\eV]} & \multicolumn{3}{c||}{[pb/G\eV]} &  \\ 
      \hline
     $ 7 $ & $ 11 $ & $ 1.81 $ & $ 0.24 $ & $ 0.20 $ & $0.21\pm0.03$ \\ 
$11 $ & $ 15 $ & $ 2.05 $ & $ 0.21 $ & $ 0.23 $ & $0.26\pm0.03$ \\ 
$15 $ & $ 38 $ & $ 0.361 $ & $ 0.032 $ & $ 0.040 $ & $0.35\pm0.03$ \\ 
      \hline
      \hline
      \multicolumn{2}{|c|}{\dphijets\ range} & \dsdx{\dphijets} & stat. & sys. & $f_{b}\pm$stat.  \\  
      \multicolumn{2}{|c|}{[deg]} & \multicolumn{3}{c||}{[pb/deg]} & \\ 
      \hline
 $  0 $ & $ 155 $ & $ 0.0084 $ & $ 0.0030 $ & $ 0.0016 $ & $0.18\pm0.04$ \\ 
$155 $ & $ 173 $ & $ 0.617 $ & $ 0.052 $ & $ 0.068 $ & $0.28\pm0.02$  \\ 
$173 $ & $ 180 $ & $ 1.640 $ & $ 0.144 $ & $ 0.180 $ & $0.30\pm0.03$  \\ 
      \hline
    \end{tabular}
    \caption{Bin averaged differential cross sections for beauty photoproduction of dijet events using semi-muonic decays for $\xgamma>0.75$ in bins of \ptmuon, \etamuon, \ptjetone, and \dphijets \ with their statistical and systematic uncertainties. The fit parameter $f_b$ is given including its statistical error. }
    \label{tab:beautydirect}
  \end{center}
\end{table}

\begin{table}[tb]
  \begin{center}
    \begin{tabular}{|rr|rrr||c|}
      \hline
     \multicolumn{6}{|c|}{\bf \boldmath H1 Beauty Dijet Muon Cross Sections, \bf \boldmath \xgamma$\leq0.75$}   \\
      \multicolumn{6}{|c|}{$\sigma_{vis}(ep\rightarrow eb\bar{b}X\rightarrow ejj\mu X)$}   \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\ptmuon\ range} & \dsdx{\ptmuon} & stat. & sys. & $f_{b}\pm$stat. \\  
      \multicolumn{2}{|c|}{[G\eV]} & \multicolumn{3}{c||}{[pb/G\eV]} &  \\ 
      \hline
      $2.5 $ & $ 3.3 $ & $ 11.3 $ & $ 1.7 $ & $ 1.2 $ & $0.18\pm0.03$ \\ 
$3.3 $ & $ 4.7 $ & $ 5.47 $ & $ 0.67 $ & $ 0.60 $ & $0.29\pm0.03$ \\ 
$4.7 $ & $ 15.0 $ & $ 0.243 $ & $ 0.053 $ & $ 0.027 $ & $0.27\pm0.06$ \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\etamuon\ range} & \dsdx{\etamuon} & stat. & sys. & $f_{b}\pm$stat.   \\  
      \multicolumn{2}{|c|}{} & \multicolumn{3}{c||}{[pb]} &  \\ 
      \hline
      $-1.3 $ & $ -0.3 $ & $ 2.58 $ & $ 1.38 $ & $ 0.29 $ & $ 0.15\pm0.08$ \\ 
$-0.3 $ & $ 0.0 $ & $ 7.98 $ & $ 1.55 $ & $ 0.88 $ & $ 0.23\pm0.04$  \\ 
$0.0 $ & $ 0.3 $ & $ 5.51 $ & $ 1.64 $ & $ 0.61 $ & $ 0.14\pm0.04$ \\ 
$0.3 $ & $ 0.6 $ & $ 7.68 $ & $ 1.77 $ & $ 0.85 $ & $ 0.22\pm0.05$  \\ 
$0.6 $ & $ 1.5 $ & $ 6.38 $ & $ 1.08 $ & $ 0.70 $ & $ 0.21\pm0.04$ \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\ptjetone\ range} & \dsdx{\ptjetone} & stat. & sys. & $f_{b}\pm$stat.  \\  
      \multicolumn{2}{|c|}{[G\eV]} & \multicolumn{3}{c||}{[pb/G\eV]} &  \\ 
      \hline
      $7 $ & $ 11 $ & $ 2.82 $ & $ 0.33 $ & $ 0.31 $ & $0.26\pm0.03$ \\ 
$11 $ & $ 15 $ & $ 1.116 $ & $ 0.188 $ & $ 0.123 $ & $0.20\pm0.03$ \\ 
$15 $ & $ 38 $ & $ 0.095 $ & $ 0.018 $ & $ 0.010 $ & $0.19\pm0.04$  \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\dphijets\ range} & \dsdx{\dphijets} & stat. & sys. & $f_{b}\pm$stat.   \\  
      \multicolumn{2}{|c|}{[deg]} & \multicolumn{3}{c||}{[pb/deg]} &  \\ 
      \hline
     $ 0 $ & $ 155 $ & $ 0.0365 $ & $ 0.0063 $ & $ 0.0040 $ & $0.23\pm0.04$ \\ 
$155 $ & $ 173 $ & $ 0.462 $ & $ 0.056 $ & $ 0.051 $ & $0.26\pm0.03$ \\ 
$173 $ & $ 180 $ & $ 0.320 $ & $ 0.012 $ & $ 0.035 $ & $0.11\pm0.04$  \\
      \hline
    \end{tabular}
    \caption{Bin averaged differential cross sections for beauty photoproduction of dijet events using semi-muonic decays for $\xgamma\leq0.75$ in bins of \ptmuon, \etamuon, \ptjetone, and \dphijets \ with their statistical and systematic uncertainties. The fit parameter $f_b$ is given including its statistical error. }
    \label{tab:beautyresolved}
  \end{center}
\end{table}


\begin{table}[tb]
  \begin{center}
    \begin{tabular}{|rr|rrr||c|}
      \hline
      \multicolumn{6}{|c|}{\bf \boldmath H1 Charm Dijet Muon Cross Sections}   \\
      \multicolumn{6}{|c|}{$\sigma_{vis}(ep\rightarrow ec\bar{c}X\rightarrow ejj\mu X)$}   \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\ptmuon\ range} & \dsdx{\ptmuon} & stat. & sys. & $f_{c}\pm$stat.  \\  
      \multicolumn{2}{|c|}{[G\eV]} & \multicolumn{3}{c||}{[pb/G\eV]} & \\ 
      \hline
      $2.5 $ & $ 3.3 $ & $ 49.1 $ & $ 4.0 $ & $ 5.1 $ & $0.50\pm0.04$\\ 
$3.3 $ & $ 4.7 $ & $ 18.3 $ & $ 1.7 $ & $ 1.9 $ & $0.47\pm0.04$ \\ 
$4.7 $ & $ 15.0 $ & $ 0.854 $ & $ 0.126 $ & $ 0.088 $ & $0.31\pm0.04$ \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\etamuon\ range} & \dsdx{\etamuon} & stat. & sys. & $f_{c}\pm$stat.  \\  
      \multicolumn{2}{|c|}{} & \multicolumn{3}{c||}{[pb]} &  \\ 
      \hline
      $-1.3 $ & $ -0.3 $ & $ 20.3 $ & $ 2.5 $ & $ 2.1 $ & $0.44\pm0.05$\\ 
$-0.3 $ & $ 0.0 $ & $ 37.0 $ & $ 4.6 $ & $ 3.9 $ & $0.45\pm0.05$ \\ 
$0.0 $ & $ 0.3 $ & $ 42.3 $ & $ 4.6 $ & $ 4.4 $ & $0.48\pm0.05$ \\ 
$0.3 $ & $ 0.6 $ & $ 38.0 $ & $ 4.6 $ & $ 4.0 $ & $0.49\pm0.06$ \\ 
$0.6 $ & $ 1.5 $ & $ 23.6 $ & $ 3.1 $ & $ 2.5 $ & $0.46\pm0.05$ \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\ptjetone\ range} & \dsdx{\ptjetone} & stat. & sys. & $f_{c}\pm$stat.  \\  
      \multicolumn{2}{|c|}{[G\eV]} & \multicolumn{3}{c||}{[pb/G\eV]} & \\ 
      \hline
      $7 $ & $ 11 $ & $ 11.8 $ & $ 0.9 $ & $ 1.2 $ & $0.55\pm0.04$\\ 
$11 $ & $ 15 $ & $ 5.22 $ & $ 0.59 $ & $ 0.54 $ & $0.39\pm0.04$ \\ 
$15 $ & $ 38 $ & $ 0.657 $ & $ 0.066 $ & $ 0.068 $ & $0.51\pm0.05$ \\ 

      \hline
      \hline
      \multicolumn{2}{|c|}{\dphijets\ range} & \dsdx{\dphijets} & stat. & sys. & $f_{c}\pm$stat.   \\  
      \multicolumn{2}{|c|}{[deg]} & \multicolumn{3}{c||}{[pb/deg]} & \\ 
      \hline
     $0 $ & $ 155 $ & $ 0.0820 $ & $ 0.0119 $ & $ 0.0085 $ & $0.37\pm0.05$ \\ 
$155 $ & $ 173 $ & $ 2.10 $ & $ 0.15 $ & $ 0.22 $ & $0.52\pm0.04$ \\ 
$173 $ & $ 180 $ & $ 4.06 $ & $ 0.35 $ & $ 0.43 $ & $0.50\pm0.04$ \\ 

      \hline
      \hline
      \multicolumn{2}{|c|}{\xgamma\ range} & \dsdx{\xgamma} & stat. & sys. & $f_{c}\pm$stat. \\  
      \multicolumn{2}{|c|}{} & \multicolumn{3}{c||}{[pb]} & \\ 
      \hline
     $ 0.0 $ & $ 0.4 $ & $ 12.3 $ & $ 5.0 $ & $ 1.3 $ & $0.22\pm0.09$ \\ 
$0.4 $ & $ 0.75 $ & $ 63.5 $ & $ 5.8 $ & $ 6.6 $ & $0.50\pm0.04$ \\ 
$0.75 $ & $ 1.0 $ & $ 206.7 $ & $ 13.8 $ & $ 21.5 $ & $0.51\pm0.03$ \\
      \hline
    \end{tabular}
    \caption{Bin averaged differential cross sections for charm photoproduction of dijet events using semi-muonic decays in bins of \ptmuon, \etamuon, \ptjetone, \dphijets, and \xgamma \ with their statistical and systematic uncertainties. The fit parameter $f_c$ is given including its statistical error. }
    \label{tab:charmfull}
  \end{center}
\end{table}

\begin{table}[tb]
  \begin{center}
    \begin{tabular}{|rr|rrr||c|}
      \hline
      \multicolumn{6}{|c|}{\bf \boldmath H1 Charm Dijet Muon Cross Sections, \bf \boldmath \xgamma$>0.75$}   \\
      \multicolumn{6}{|c|}{$\sigma_{vis}(ep\rightarrow ec\bar{c}X\rightarrow ejj\mu X)$}   \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\ptmuon\ range} & \dsdx{\ptmuon} & stat. & sys. & $f_{c}\pm$stat.  \\  
      \multicolumn{2}{|c|}{[G\eV]} & \multicolumn{3}{c||}{[pb/G\eV]} & \\ 
      \hline
     $2.5 $ & $ 3.3 $ & $ 31.2 $ & $ 3.2 $ & $ 3.2 $ & $0.55\pm0.05$  \\ 
$3.3 $ & $ 4.7 $ & $ 11.6 $ & $ 1.4 $ & $ 1.2 $ & $0.49\pm0.06$ \\ 
$4.7 $ & $ 15.0 $ & $ 0.368 $ & $ 0.119 $ & $ 0.038 $ & $0.20\pm0.07$ \\ 

      \hline
      \hline
      \multicolumn{2}{|c|}{\etamuon\ range} & \dsdx{\etamuon} & stat. & sys. & $f_{c}\pm$stat.  \\  
      \multicolumn{2}{|c|}{} & \multicolumn{3}{c||}{[pb]} & \\ 
      \hline
     $ -1.3 $ & $ -0.3 $ & $ 14.6 $ & $ 2.2 $ & $ 1.5 $ & $0.45\pm0.06$ \\ 
$-0.3 $ & $ 0.0 $ & $ 21.9 $ & $ 4.0 $ & $ 2.3 $ & $0.41\pm0.07$ \\ 
$0.0 $ & $ 0.3 $ & $ 27.0 $ & $ 3.4 $ & $ 2.8 $ & $0.50\pm0.06$ \\ 
$0.3 $ & $ 0.6 $ & $ 25.8 $ & $ 3.0 $ & $ 2.7 $ & $0.56\pm0.06$ \\ 
$0.6 $ & $ 1.5 $ & $ 9.4 $ & $ 1.9 $ & $1.0 $ & $0.39\pm0.08$\\ 
      \hline
      \hline
      \multicolumn{2}{|c|}{\ptjetone\ range} & \dsdx{\ptjetone} & stat. & sys. & $f_{c}\pm$stat.  \\  
      \multicolumn{2}{|c|}{[G\eV]} & \multicolumn{3}{c||}{[pb/G\eV]} & \\ 
      \hline
     $ 7 $ & $ 11 $ & $ 6.6 $ & $ 0.6 $ & $ 0.7 $ & $0.57\pm0.05$ \\ 
$11 $ & $ 15 $ & $ 3.45 $ & $ 0.42 $ & $ 0.36 $ & $0.42\pm0.05$ \\ 
$15 $ & $ 38 $ & $ 0.377 $ & $ 0.048 $ & $ 0.039 $ & $0.42\pm0.05$ \\ 
      \hline
      \hline
      \multicolumn{2}{|c|}{\dphijets\ range} & \dsdx{\dphijets} & stat. & sys. & $f_{c}\pm$stat.  \\  
      \multicolumn{2}{|c|}{[deg]} & \multicolumn{3}{c||}{[pb/deg]} & \\ 
      \hline
     $0 $ & $ 155 $ & $ 0.0424 $ & $ 0.0075 $ & $ 0.0045 $ & $0.46\pm0.08$ \\ 
$155 $ & $ 173 $ & $ 1.202 $ & $ 0.122 $ & $ 0.125 $ & $0.46\pm0.05$ \\ 
$173 $ & $ 180 $ & $ 3.00 $ & $ 0.29 $ & $ 0.31 $ & $0.53\pm0.05$ \\ 

      \hline
    \end{tabular}
    \caption{Bin averaged differential cross sections for charm photoproduction of dijet events using semi-muonic decays for $\xgamma>0.75$ in bins of \ptmuon, \etamuon, \ptjetone, and \dphijets \ with their statistical and systematic uncertainties. The fit parameter $f_c$ is given including its statistical error. }
    \label{tab:charmdirect}
  \end{center}
\end{table}

\begin{table}[tb]
  \begin{center}
    \begin{tabular}{|rr|rrr||c|}
      \hline
       \multicolumn{6}{|c|}{\bf \boldmath H1 Charm Dijet Muon Cross Sections, \bf \boldmath \xgamma$\leq0.75$}   \\
      \multicolumn{6}{|c|}{$\sigma_{vis}(ep\rightarrow ec\bar{c}X\rightarrow ejj\mu X)$}   \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\ptmuon\ range} & \dsdx{\ptmuon} & stat. & sys. & $f_{c}\pm$stat.  \\  
      \multicolumn{2}{|c|}{[G\eV]} & \multicolumn{3}{c||}{[pb/G\eV]} & \\ 
      \hline
      $2.5 $ & $ 3.3 $ & $ 17.6 $ & $ 2.7 $ & $ 1.8 $ & $0.43\pm0.06$ \\ 
$3.3 $ & $ 4.7 $ & $ 6.7 $ & $ 1.0 $ & $ 0.7 $ & $0.45\pm0.06$\\ 
$4.7 $ & $ 15.0 $ & $ 0.424 $ & $ 0.068 $ & $ 0.044 $ & $0.45\pm0.07$ \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\etamuon\ range} & \dsdx{\etamuon} & stat. & sys. & $f_{c}\pm$stat.  \\  
      \multicolumn{2}{|c|}{} & \multicolumn{3}{c||}{[pb]} & \\ 
      \hline
     $ -1.3 $ & $ -0.3 $ & $ 4.16 $ & $ 1.25 $ & $ 0.43 $ & $0.32\pm0.10$\\ 
$-0.3 $ & $ 0.0 $ & $ 12.7 $ & $ 2.5 $ & $ 1.3 $ & $0.45\pm0.08$ \\ 
$0.0 $ & $ 0.3 $ & $ 12.3 $ & $ 2.6 $ & $ 1.3 $ & $0.37\pm0.08$ \\
$0.3 $ & $ 0.6 $ & $ 9.9 $ & $ 2.7 $ & $ 1.0 $ & $0.32\pm0.09$ \\ 
$0.6 $ & $ 1.5 $ & $ 11.7 $ & $ 1.8 $ & $ 1.2 $ & $0.43\pm0.06$ \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\ptjetone\ range} & \dsdx{\ptjetone} & stat. & sys. & $f_{c}\pm$stat.  \\  
      \multicolumn{2}{|c|}{[G\eV]} & \multicolumn{3}{c||}{[pb/G\eV]} &\\ 
      \hline
      $7 $ & $ 11 $ & $ 4.58 $ & $ 0.51 $ & $ 0.48 $ & $0.50\pm0.05$ \\ 
$11 $ & $ 15 $ & $ 1.216 $ & $ 0.375 $ & $ 0.126 $ & $0.24\pm0.07$ \\ 
$15 $ & $ 38 $ & $ 0.265 $ & $ 0.037 $ & $ 0.028 $ & $0.65\pm0.08$  \\
      \hline
      \hline
      \multicolumn{2}{|c|}{\dphijets\ range} & \dsdx{\dphijets} & stat. & sys. & $f_{c}\pm$stat.  \\  
      \multicolumn{2}{|c|}{[deg]} & \multicolumn{3}{c||}{[pb/deg]} & \\ 
      \hline
      $0 $ & $ 155 $ & $ 0.0320 $ & $ 0.0092 $ & $ 0.0047 $ & $0.24\pm0.07$ \\ 
$155 $ & $ 173 $ & $ 0.947 $ & $ 0.102 $ & $ 0.098 $ & $0.64\pm0.07$\\ 
$173 $ & $ 180 $ & $ 0.834 $ & $ 0.185 $ & $ 0.087 $ & $0.35\pm0.07$ \\ 

      \hline
    \end{tabular}
    \caption{Bin averaged differential cross sections for charm photoproduction of dijet events using semi-muonic decays for $\xgamma\leq0.75$ in bins of \ptmuon, \etamuon, \ptjetone, and \dphijets \ with their statistical and systematic uncertainties. The fit parameter $f_c$ is given including its statistical error. }
    \label{tab:charmresolved}
  \end{center}
\end{table}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}
\setlength{\unitlength}{1cm}
\centering
\begin{picture}(15.,8.)(0.,0.)
\put(-1.5,0.){\epsfig{file=d12-059f1a.eps,width = 9cm}}
\put(7.7,0.){\epsfig{file=d12-059f1b.eps,width = 9cm}}
\put(0.,7.2){(a)}
\put(9.2,7.2){(b)}
\end{picture}
\caption{The impact parameter \impactpar \ and \ptrel \ distributions for the total event sample. The data are compared to the \textsc{Pythia} predictions of the different quark contributions and their sum. The MC predictions are obtained from the fit to the total data sample.}
\label{fig:plotsflavourseparation}
\end{figure}

\begin{figure}
\setlength{\unitlength}{1cm}
\centering
\begin{picture}(15.,21.)(0.,0.)
\put(-1.5,15.){\epsfig{file=d12-059f2a.eps,width = 9cm}}
\put(7.7,15.){\epsfig{file=d12-059f2b.eps,width=9cm}}
\put(-1.5,7.5){\epsfig{file=d12-059f2c.eps,width = 9cm}}
\put(7.7,7.5){\epsfig{file=d12-059f2d.eps,width = 9cm}}
\put(-1.5,0.){\epsfig{file=d12-059f2e.eps,width = 9cm}}
\put(7.7,0.){\epsfig{file=d12-059f2f.eps,width = 9cm}}
\put(0.,22.1){(a)}
\put(9.2,22.1){(b)}
\put(0.,14.6){(c)}
\put(9.2,14.6){(d)}
\put(0.,7.2){(e)}
\put(9.2,7.2){(f)}
\end{picture}
\caption{Distributions of kinematic variables for the selected event sample. $P_T^{leading track}$ is the transverse momentum of the track with the highest transverse momentum in the event. The other notations are described in the text. The data are compared to the \textsc{Pythia} predictions of the different quark contributions and their sum. The MC predictions are obtained from the fit to the total data sample.}
\label{fig:controlplots}
\end{figure}


\begin{figure}
\setlength{\unitlength}{1cm}
\centering
\begin{picture}(15.,21.)(0.,0.)
\put(-1.5,15.){\epsfig{file=d12-059f3a.eps,width = 8.5cm}}
\put(7.7,15.){\epsfig{file=d12-059f3b.eps,width=8.5cm}}
\put(-1.5,7.5){\epsfig{file=d12-059f3c.eps,width = 8.5cm}}
\put(7.7,7.5){\epsfig{file=d12-059f3d.eps,width = 8.5cm}}
\put(-1.5,0.){\epsfig{file=d12-059f3e.eps,width = 8.5cm}}
\put(0.,22.1){(a)}
\put(9.2,22.1){(b)}
\put(0.,14.6){(c)}
\put(9.2,14.6){(d)}
\put(0.,7.2){(e)}
\end{picture}
\caption{The differential cross sections for beauty photoproduction of dijet events with a muon as a function of \ptmuon, \etamuon, \ptjetone, \dphijets, and \xgamma. The inner error bars show the statistical errors, the outer error bars represent the statistical and systematic errors added in quadrature. The data are compared with the predictions from the LO models \textsc{Pythia},  \textsc{Cascade} and \textsc{Herwig} as well as to the NLO predictions of MC@NLO. The theoretical uncertainties of MC@NLO are given as shaded band.}
\label{fig:beautyxsecfullsample}

\end{figure}


\clearpage

\begin{figure}
\setlength{\unitlength}{1cm}
\centering
\begin{picture}(15.,15.)(0.,0.)
\put(-1.5,7.5){\epsfig{file=d12-059f4a.eps,width = 8.5cm}}
\put(7.7,7.5){\epsfig{file=d12-059f4b.eps,width=8.5cm}}
\put(-1.5,0.){\epsfig{file=d12-059f4c.eps,width = 8.5cm}}
\put(7.7,0.){\epsfig{file=d12-059f4d.eps,width = 8.5cm}}
\put(0.,14.6){(a)}
\put(9.2,14.6){(b)}
\put(0.,7.2){(c)}
\put(9.2,7.2){(d)}
\end{picture}
\caption{The differential cross sections for beauty photoproduction of dijet events using semi-muonic decays for $\xgamma>0.75$ as a function of \ptmuon, \etamuon, \ptjetone, and \dphijets. For details see caption of figure~\ref{fig:beautyxsecfullsample}.}
\label{fig:beautyxsecdirectsample}
\end{figure}

\clearpage

\begin{figure}
\setlength{\unitlength}{1cm}
\centering
\begin{picture}(15.,15.)(0.,0.)
\put(-1.5,7.5){\epsfig{file=d12-059f5a.eps,width = 8.5cm}}
\put(7.7,7.5){\epsfig{file=d12-059f5b.eps,width=8.5cm}}
\put(-1.5,0.){\epsfig{file=d12-059f5c.eps,width = 8.5cm}}
\put(7.7,0.){\epsfig{file=d12-059f5d.eps,width = 8.5cm}}
\put(0.,14.6){(a)}
\put(9.2,14.6){(b)}
\put(0.,7.2){(c)}
\put(9.2,7.2){(d)}
\end{picture}
\caption{The differential cross sections for beauty photoproduction of dijet events using semi-muonic decays for $\xgamma\leq0.75$ as a function of \ptmuon, \etamuon, \ptjetone, and \dphijets. For details see caption of figure~\ref{fig:beautyxsecfullsample}.}
\label{fig:beautyxsecresolvedsample}
\end{figure}

\clearpage

\begin{figure}
\setlength{\unitlength}{1cm}
\centering
\begin{picture}(15.,21.)(0.,0.)
\put(-1.5,15.){\epsfig{file=d12-059f6a.eps,width = 8.5cm}}
\put(7.7,15.){\epsfig{file=d12-059f6b.eps,width=8.5cm}}
\put(-1.5,7.5){\epsfig{file=d12-059f6c.eps,width = 8.5cm}}
\put(7.7,7.5){\epsfig{file=d12-059f6d.eps,width = 8.5cm}}
\put(-1.5,0.){\epsfig{file=d12-059f6e.eps,width = 8.5cm}}
\put(0.,22.1){(a)}
\put(9.2,22.1){(b)}
\put(0.,14.6){(c)}
\put(9.2,14.6){(d)}
\put(0.,7.2){(e)}
\end{picture}
\caption{The differential cross sections for charm photoproduction of dijet events using semi-muonic decays as a function of \ptmuon, \etamuon, \ptjetone, \dphijets, and \xgamma. For details see caption of figure~\ref{fig:beautyxsecfullsample}.}\label{fig:charmxsecfullsample}
\end{figure}


\clearpage

\begin{figure}
\setlength{\unitlength}{1cm}
\centering
\begin{picture}(15.,15.)(0.,0.)
\put(-1.5,7.5){\epsfig{file=d12-059f7a.eps,width = 8.5cm}}
\put(7.7,7.5){\epsfig{file=d12-059f7b.eps,width=8.5cm}}
\put(-1.5,0.){\epsfig{file=d12-059f7c.eps,width = 8.5cm}}
\put(7.7,0.){\epsfig{file=d12-059f7d.eps,width = 8.5cm}}
\put(0.,14.6){(a)}
\put(9.2,14.6){(b)}
\put(0.,7.2){(c)}
\put(9.2,7.2){(d)}
\end{picture}
\caption{The differential cross sections for charm photoproduction of dijet events using semi-muonic decays for $\xgamma>0.75$ as a function of \ptmuon, \etamuon, \ptjetone, and \dphijets. For details see caption of figure~\ref{fig:beautyxsecfullsample}.}
\label{fig:charmxsecdirectsample}
\end{figure}

\clearpage

\begin{figure}
\setlength{\unitlength}{1cm}
\centering
\begin{picture}(15.,15.)(0.,0.)
\put(-1.5,7.5){\epsfig{file=d12-059f8a.eps,width = 8.5cm}}
\put(7.7,7.5){\epsfig{file=d12-059f8b.eps,width=8.5cm}}
\put(-1.5,0.){\epsfig{file=d12-059f8c.eps,width = 8.5cm}}
\put(7.7,0.){\epsfig{file=d12-059f8d.eps,width = 8.5cm}}
\put(0.,14.6){(a)}
\put(9.2,14.6){(b)}
\put(0.,7.2){(c)}
\put(9.2,7.2){(d)}
\end{picture}
\caption{The differential cross sections for charm photoproduction of dijet events using semi-muonic decays for $\xgamma\leq0.75$ as a function of \ptmuon, \etamuon, \ptjetone, and \dphijets.  For details see caption of figure~\ref{fig:beautyxsecfullsample}.}
\label{fig:charmxsecresolvedsample}
\end{figure}






\end{document}




