%================================================================
% LaTeX file with prefered layout for H1 paper drafts
% use: dvips -D600 file-name
%================================================================
\documentclass[12pt]{article}
%\usepackage{epsfig}
\usepackage{lineno}
\usepackage{graphicx}
\usepackage{amsmath}
\usepackage{hhline}
\usepackage{amssymb}
\usepackage{times}
\usepackage{cite}
\usepackage{array}
\usepackage{multirow}

\usepackage{color}
\definecolor{rltred}{rgb}{0.75,0,0}
\definecolor{rltgreen}{rgb}{0,0.5,0}
\definecolor{rltblue}{rgb}{0,0,0.75}

\newif\ifpdf
\ifx\pdfoutput\undefined
    \pdffalse          % we are not running PDFLaTeX
\else
    \pdfoutput=1       % we are running PDFLaTeX
    \pdftrue
\fi

\ifpdf
\linenumbers           % Remove for final publication!!!!!     
\usepackage{thumbpdf}
\usepackage[pdftex,
        colorlinks=true,
        urlcolor=rltblue,       % \href{...}{...} external (URL)
        filecolor=rltgreen,     % \href{...} local file
        linkcolor=rltred,       % \ref{...} and \pageref{...}
        pdftitle={Example File for pdflatex},
        pdfauthor={H1},
        pdfsubject={Template file},
        pdfkeywords={High-Energy Physics, Particle Physics},
%        pagebackref,
        pdfpagemode=None,
        bookmarksopen=true]{hyperref}
 
\fi

\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

\newcommand{\myref}[2]{%
\ifpdf \href{#1}{#2}\else #2 \fi%
}

%===============================title page=============================
\begin{document}  


% The rest
\newcommand{\pom}{{I\!\!P}}
\newcommand{\reg}{{I\!\!R}}
\newcommand{\slowpi}{\pi_{\mathit{slow}}}
%\newcommand{\gevsq}{\mathrm{GeV}^2}
\newcommand{\fiidiii}{F_2^{D(3)}}
\newcommand{\fiidiiiarg}{\fiidiii\,(\beta,\,Q^2,\,x)}
\newcommand{\n}{1.19\pm 0.06 (stat.) \pm0.07 (syst.)}
\newcommand{\nz}{1.30\pm 0.08 (stat.)^{+0.08}_{-0.14} (syst.)}
\newcommand{\fiidiiiful}{F_2^{D(4)}\,(\beta,\,Q^2,\,x,\,t)}
\newcommand{\fiipom}{\tilde F_2^D}
\newcommand{\ALPHA}{1.10\pm0.03 (stat.) \pm0.04 (syst.)}
\newcommand{\ALPHAZ}{1.15\pm0.04 (stat.)^{+0.04}_{-0.07} (syst.)}
\newcommand{\fiipomarg}{\fiipom\,(\beta,\,Q^2)}
\newcommand{\pomflux}{f_{\pom / p}}
\newcommand{\nxpom}{1.19\pm 0.06 (stat.) \pm0.07 (syst.)}
\newcommand {\gapprox}
   {\raisebox{-0.7ex}{$\stackrel {\textstyle>}{\sim}$}}
\newcommand {\lapprox}
   {\raisebox{-0.7ex}{$\stackrel {\textstyle<}{\sim}$}}
\def\gsim{\,\lower.25ex\hbox{$\scriptstyle\sim$}\kern-1.30ex%
\raise 0.55ex\hbox{$\scriptstyle >$}\,}
\def\lsim{\,\lower.25ex\hbox{$\scriptstyle\sim$}\kern-1.30ex%
\raise 0.55ex\hbox{$\scriptstyle <$}\,}
\newcommand{\pomfluxarg}{f_{\pom / p}\,(x_\pom)}
\newcommand{\dsf}{\mbox{$F_2^{D(3)}$}}
\newcommand{\dsfva}{\mbox{$F_2^{D(3)}(\beta,Q^2,x_{I\!\!P})$}}
\newcommand{\dsfvb}{\mbox{$F_2^{D(3)}(\beta,Q^2,x)$}}
\newcommand{\dsfpom}{$F_2^{I\!\!P}$}
\newcommand{\gap}{\stackrel{>}{\sim}}
\newcommand{\lap}{\stackrel{<}{\sim}}
\newcommand{\fem}{$F_2^{em}$}
\newcommand{\tsnmp}{$\tilde{\sigma}_{NC}(e^{\mp})$}
\newcommand{\tsnm}{$\tilde{\sigma}_{NC}(e^-)$}
\newcommand{\tsnp}{$\tilde{\sigma}_{NC}(e^+)$}
\newcommand{\st}{$\star$}
\newcommand{\sst}{$\star \star$}
\newcommand{\ssst}{$\star \star \star$}
\newcommand{\sssst}{$\star \star \star \star$}
\newcommand{\tw}{\theta_W}
\newcommand{\sw}{\sin{\theta_W}}
\newcommand{\cw}{\cos{\theta_W}}
\newcommand{\sww}{\sin^2{\theta_W}}
\newcommand{\cww}{\cos^2{\theta_W}}
\newcommand{\trm}{m_{\perp}}
\newcommand{\trp}{p_{\perp}}
\newcommand{\trmm}{m_{\perp}^2}
\newcommand{\trpp}{p_{\perp}^2}
\newcommand{\alp}{\alpha_s}

\newcommand{\alps}{\alpha_s}
\newcommand{\sqrts}{$\sqrt{s}$}
\newcommand{\LO}{$O(\alpha_s^0)$}
\newcommand{\Oa}{$O(\alpha_s)$}
\newcommand{\Oaa}{$O(\alpha_s^2)$}
\newcommand{\PT}{p_{\perp}}
\newcommand{\JPSI}{J/\psi}
\newcommand{\sh}{\hat{s}}
%\newcommand{\th}{\hat{t}}
\newcommand{\uh}{\hat{u}}
\newcommand{\MP}{m_{J/\psi}}
%\newcommand{\PO}{\mbox{l}\!\mbox{P}}
\newcommand{\PO}{I\!\!P}
\newcommand{\xbj}{x}
\newcommand{\xpom}{x_{\PO}}
\newcommand{\ttbs}{\char'134}
\newcommand{\xpomlo}{3\times10^{-4}}  
\newcommand{\xpomup}{0.05}  
\newcommand{\dgr}{^\circ}
\newcommand{\pbarnt}{\,\mbox{{\rm pb$^{-1}$}}}
\newcommand{\gev}{\,\mbox{GeV}}
\newcommand{\WBoson}{\mbox{$W$}}
\newcommand{\fbarn}{\,\mbox{{\rm fb}}}
\newcommand{\fbarnt}{\,\mbox{{\rm fb$^{-1}$}}}
%
% Some useful tex commands
%
\newcommand{\qsq}{\ensuremath{Q^2} }
\newcommand{\gevsq}{\ensuremath{\mathrm{GeV}^2} }
\newcommand{\et}{\ensuremath{E_t^*} }
\newcommand{\rap}{\ensuremath{\eta^*} }
\newcommand{\gp}{\ensuremath{\gamma^*}p }
\newcommand{\dsiget}{\ensuremath{{\rm d}\sigma_{ep}/{\rm d}E_t^*} }
\newcommand{\dsigrap}{\ensuremath{{\rm d}\sigma_{ep}/{\rm d}\eta^*} }
% 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
%%%Date:        \today       \\
%%%%Version:     0.0          \\
%%%Editors:     \myref{mailto:mehta@mail.desy.de}{A.~Mehta}, \myref{mailto:thompspd@mail.desy.de}{P.~Thompson}       \\
%%%Referees:    \myref{mailto:cmk@mppmu.mpg.de}{C.~Kiesling}, \myref{mailto:sefkow@mail.desy.de}{F.~Sefkow}     \\
%%%%Comments by  Deadline       

%%%\vspace{2cm}

\vspace{1cm}
\begin{center}
{\large \bf H1prelim-04-072 } 
\end{center}
\vspace{1cm}



\begin{center}
\begin{Large}
  
  {\boldmath {\bf Measurement of $F_2^{c\bar{c}}$ and $F_2^{b\bar{b}}$
      at High $Q^2$}}

\vspace{2cm}

H1 Collaboration

\end{Large}
\end{center}

\vspace{2cm}

\begin{abstract}
  \noindent
  A measurement is presented of the charm and beauty cross sections at
  for $Q^2 > 150$~GeV and $0.1<y<0.7$, using a method based on the
  distance of closest approach of tracks, reconstructed using the
  central silicon tracker, to the production vertex.  The data are
  divided into 3 $x$-$Q^2$ bins and values for $F_2^{c\bar{c}}$ and
  $F_2^{b\bar{b}}$ are evaluated.
\end{abstract}

\vspace{1.5cm}

\begin{center}
Prepared for DIS 2004, Strbske Pleso, Slovakia 
\end{center}


\begin{center}
%To be submitted to 
\end{center}

\end{titlepage}

%
%          COPY THE AUTHOR AND INSTITUTE LISTS 
%       AT THE TIME OF THE T0-TALK INTO YOUR AREA
%
% from /h1/iww/ipublications/h1auts.tex 

%\begin{flushleft}
%  \input{h1auts}
%\end{flushleft}

\newpage

Cross sections for  $Q^2>150$~${\rm GeV}^2$, $0.1<y<0.7$:
\begin{equation*}
 \sigma_{c\bar{c}} =  431 \pm 59 \pm 69 \ {\rm pb}
\end{equation*}
\begin{equation*}
 \sigma_{b\bar{b}} = 45 \pm 11 \pm 11 \ {\rm pb}
\end{equation*}


Compared with results of a next to leading order QCD fit
 which was fitted to inclusive cross sections from H1~\cite{H1NCCC}:
\begin{equation*}
 \sigma_{c\bar{c}} =  455  \ {\rm pb}
\end{equation*}
\begin{equation*}
 \sigma_{b\bar{b}} =  51.8  \ {\rm pb}
\end{equation*}

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

\bibitem{H1NCCC}
C.~Adloff {\it et al.}  
\ifpdf \href{http://www-h1.desy.de/publications/H1publication.short_list.html}{[H1 Collaboration],}
\else [H1 Collaboration],
\fi
Eur.\ Phys.\ J {\bf C 30} (2003) 1-32
\ifpdf
\href{http://arXiv.org/ps/hep-ex/0304003}{[hep-ex/0304003.}
\else
[hep-ex/0304003].
\fi
                             

\end{thebibliography}


\newpage

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[hhh]
  \begin{center} \includegraphics[width=\textwidth]{H1prelim-04-072.fig1.eps}
  \caption{The number of reconstructed central silicon tracker (CST)
  tracks per jet (or per event if there are no jets with
  $P_T>5$~GeV). Each CST track is required to have at least 2 CST hits
  and $P_T>0.5$ GeV. Also included is the expectation from the RAPGAP
  Monte Carlo, showing the contribution from the various quark
  flavours after applying the scale factors obtained from the fit of the
  significance distributions of the data (see figure~\ref{fig:s1negsub}).}
\label{fig:csttracks} 
  \end{center}
\end{figure}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[hhh]
  \begin{center}
    \includegraphics[width=\textwidth]{H1prelim-04-072.fig2.eps}
    \caption{The transverse momentum distribution of all CST tracks.
      Also included is the expectation from the RAPGAP Monte
      Carlo, showing the contribution from the various quark flavours after
      applying the scale factors obtained  from the fit to the
  significance distributions of the data.}
    \label{fig:pttracks} 
  \end{center}
\end{figure}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[hhh]
  \begin{center}
    \includegraphics[width=\textwidth]{H1prelim-04-072.fig3.eps}
    \caption{The theta distribution  of all CST tracks.
      Also included is the expectation from the RAPGAP Monte
      Carlo, showing the contribution from the various quark flavours after
      applying the scale factors obtained  from the fit to the
  significance distributions of the data.}
    \label{fig:thetatracks} 
  \end{center}
\end{figure}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[hhh]
  \begin{center}
    \includegraphics[width=\textwidth]{H1prelim-04-072.fig4.eps}
    \caption{The theta distribution  of all jets which contain
      at least 1 reconstructed CST track within a cone of radius 1.0.
      Also included is the expectation from the RAPGAP Monte Carlo after
      applying the scale factors obtained  from the fit to the
  significance distributions of the data.}
      \label{fig:jettheta} 
  \end{center}
\end{figure}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[hhh]
  \begin{center}
    \includegraphics[width=\textwidth]{H1prelim-04-072.fig5.eps}
    \caption{The transverse momentum distribution  of all jets which contain
      at least 1 reconstructed CST track within a cone of radius 1.0.
      If there are no reconstructed jets the complete hadronic final
      state $P_T$ is used.  Also included is the expectation from the
      RAPGAP Monte Carlo, showing the contribution from the various
      quark flavours after
      applying the scale factors obtained  from the fit to the
  significance distributions of the data.}
\label{fig:jetpt}
  \end{center}
\end{figure}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[hhh]
  \begin{center}
    \includegraphics[width=\textwidth]{H1prelim-04-072.fig6.eps}
    \caption{The distance of closest approach of a track 
      to the vertex (DCA) in the $XY$ plane for all CST tracks.  Also
      included is the expectation from the RAPGAP Monte Carlo, showing
      the contribution from the various quark flavours after
      applying the scale factors obtained  from the fit to the
  significance distributions of the data.}
    \label{fig:dca} 
  \end{center}
\end{figure}




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[hhh]
  \begin{center}
    \includegraphics[width=\textwidth]{H1prelim-04-072.fig7.eps}
    \caption{The signed significance ($DCA/ \sigma (DCA)$ 
%%%, where $\sigma (DCA)$ is the error on the DCA,
      distribution for all CST tracks, 
where $\sigma (DCA)$ is the error on the DCA.
The sign is defined as positive
      (negative) if the azimuthal difference between the jet-axis and the
      vector between the vertex and the point of closest approach of the track is less (greater) than $90^\circ$. If there are no reconstructed
jets in the event the azimuthal difference is taken with respect to
the complete hadronic final state. Also included is the
      expectation from the RAPGAP Monte Carlo, showing the contribution
      from the various quark flavours after
      applying the scale factors obtained  from the fit to the
  significance distributions of the data.}
    \label{fig:sigall} 
  \end{center}
\end{figure}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[hhh]
  \begin{center}
    \includegraphics[width=\textwidth]{H1prelim-04-072.fig8.eps}
    \caption{The signed significance $S_1=  DCA /\sigma (DCA) $ 
      distribution for all jets that contain  1 reconstructed CST
      track. Also included is the
      expectation from the RAPGAP Monte Carlo, showing the contribution
      from the various quark flavours after
      applying the scale factors obtained  from the fit to the
  significance distributions of the data.}
    \label{fig:s1} 
  \end{center}
\end{figure}



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[hhh]
  \begin{center} \includegraphics[width=\textwidth]{H1prelim-04-072.fig9.eps}
  \caption{The signed significance $S_2= DCA / \sigma (DCA) $
  distribution of the track with the second largest absolute value of
  significance for all jets that contain 2 or more reconstructed CST
  tracks. Also included is the expectation from the RAPGAP Monte
  Carlo, showing the contribution from the various quark flavours
  after applying the scale factors obtained from the fit to the
  significance distributions of the data.}  \label{fig:s2}
  \end{center}
\end{figure}




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\newpage
\begin{figure}[hhh]
\vspace{-2.2cm}
  \begin{center}
  \includegraphics[width=\textwidth]{H1prelim-04-072.fig10a.eps}
  \includegraphics[width=\textwidth]{H1prelim-04-072.fig10b.eps}
  \caption{The distibutions (a) $S_1$ and (b) $S_2$, after subtracting
  the bins with equal magnitude but negative sign from the postive.
  Also included is the results from the fit to the data of the Monte Carlo 
 distributions of events arising from 
  $c$ quarks or $b$ quarks. The light quark constribution is also shown.}
  \label{fig:s1negsub} \end{center}
  \vspace{-21.5cm} \hspace{2.5cm} (a) 

 \vspace{9.5cm} \hspace{2.5cm} (b) 

\end{figure}




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[hhh]
  \begin{center}
    \includegraphics[width=\textwidth]{H1prelim-04-072.fig11a.eps}
    \includegraphics[width=\textwidth]{H1prelim-04-072.fig11b.eps}
    \caption{ The invariant mass distribution of all
      track in jets for which (a) $5<S_2<10$ and (b) $-10<S_2<-5$ . 
        Also included is the
      expectation from the RAPGAP Monte Carlo, showing the
      contribution from the various quark flavours after
      applying the scale factors obtained  from the fit to the
  significance distributions of the data.}
    \label{fig:trackmass} 
  \end{center}
  \vspace{-21.0cm} \hspace{2.5cm} (a) 

 \vspace{10cm} \hspace{2.5cm} (b) 

  \vspace{11cm}
\end{figure}



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[hhh]
  \begin{center}
    \includegraphics[width=\textwidth]{H1prelim-04-072.fig12.eps}
    \caption{A display of an example data event showing the
      reconstructed tracks around the primary vertex position in
      the $x-y$ plane . The  primary vertex position is at $x=y=0$.
 The jet axis is shown by the dash-dotted
      line. The tracks are shown by dashed lines, with the
      bands indicating the errors. This event shows a clear interection of
tracks at about 4~mm from the primary vertex, giving good evidence
for a secondary vertex formed by long lived heavy quark hadrons.
This event has $S_2=7.2$.}
    \label{fig:event} 
  \end{center}
  
  \vspace{11cm}
\end{figure}


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[hhh]
  \begin{center}
    \includegraphics[width=\textwidth]{H1prelim-04-072.fig13.eps}
    \caption{The ratio $f^{c\bar{c}}= F_2^{c\bar{c}}/F_2$ and $f^{b\bar{b}}=F_2^{b\bar{b}}/F_2$ shown as a function
of $x$ for 2 different $Q^2$ values. Also shown is a prediction of the 
next to leading order QCD fit which was obtained using inclusive cross sections
from H1~\cite{H1NCCC}.}
    \label{fig:fraccb} 
  \end{center}
\end{figure}




%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[hhh]
  \begin{center}
    \includegraphics[width=\textwidth]{H1prelim-04-072.fig14.eps}
    \caption{The measured $F_2^{c\bar{c}}$ shown as a function
of $x$ for 2 different $Q^2$ values. Also shown is a prediction of the 
next to leading order QCD fit which was obtained using inclusive cross sections
from H1~\cite{H1NCCC}.}
    \label{fig:f2c} 
  \end{center}
\end{figure}



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}[hhh]
  \begin{center}
    \includegraphics[width=\textwidth]{H1prelim-04-072.fig15.eps}
    \caption{The measured $F_2^{b\bar{b}}$ shown as a function
of $x$ for 2 different $Q^2$ values. Also shown is a prediction of the 
next to leading order QCD fit which was obtained using inclusive cross sections
from H1~\cite{H1NCCC}.}
    \label{fig:f2b} 
  \end{center}
\end{figure}


\end{document}