% %================================================================ % LaTeX file with preferred 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{rotating} \usepackage{multirow} \usepackage{bigstrut} \usepackage{color} \definecolor{rltred}{rgb}{0.75,0,0} \definecolor{rltgreen}{rgb}{0,0.5,0} \definecolor{rltblue}{rgb}{0,0,0.75} \usepackage{dcolumn} \newcolumntype{d}{D{.}{.}{1}} \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% } %\linenumbers % Remove for final publication!!!!! %===============================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} \begin{flushleft} {\tt DESY 09-096 \hfill ISSN 0418-9833} \\ {\tt October 2009} \\ \end{flushleft} \vspace{2cm} \begin{center} \begin{Large} {\bf {\boldmath {Measurement of the Charm and Beauty \\ Structure Functions using \\ the H1 Vertex Detector at HERA}}} \vspace{2cm} H1 Collaboration \end{Large} \end{center} \vspace{2cm} \begin{abstract} \noindent Inclusive charm and beauty cross sections are measured in $e^-p$ and $e^+p$ neutral current collisions at HERA in the kinematic region of photon virtuality $5 \le Q^2 \le 2000~{\rm GeV}^2$ and Bjorken scaling variable $0.0002 \le x \le 0.05$. The data were collected with the H1 detector in the years 2006 and 2007 corresponding to an integrated luminosity of 189~${\rm pb^{-1}}$. The numbers of charm and beauty events are determined using variables reconstructed by the H1 vertex detector including the impact parameter of tracks to the primary vertex and the position of the secondary vertex. The measurements are combined with previous data and compared to QCD predictions. \end{abstract} \vspace{1.5cm} \begin{center} Accepted by Eur. Phys. J.\ C. \end{center} \end{titlepage} \begin{flushleft} %-- H1AUTS Author list by names %-- Status: Tue Jun 2 10:22:10 CEST 2009 Number of authors = 247 F.D.~Aaron$^{5,49}$, %BUCH-PD 11/06 Aaron M.~Aldaya~Martin$^{11}$, %DESY-PD 10/08 Aldayamartin C.~Alexa$^{5}$, %BUCH-PD 06/06 Alexa K.~Alimujiang$^{11}$, %DESY-PD 07/08 Alimujiang V.~Andreev$^{25}$, %LPI -PD 8/88 Andreev B.~Antunovic$^{11}$, %DESY-LEFT 12/08 Antunovic A.~Asmone$^{33}$, %ROME-ST 07/2 Asmone S.~Backovic$^{30}$, %PODG-PD 03/2 Backovic A.~Baghdasaryan$^{38}$, %YERE-PD 09/03 Baghdasaryana E.~Barrelet$^{29}$, %PARI-PD 11/99 Barrelet W.~Bartel$^{11}$, %DESY-PD 8/88 Bartel K.~Begzsuren$^{35}$, %ULBA-PD 04/06 Begzsuren A.~Belousov$^{25}$, %LPI -PD 8/88 Belousov J.C.~Bizot$^{27}$, %ORSA-PD 8/88 Bizot V.~Boudry$^{28}$, %ECPL-PD 1/93 Boudry I.~Bozovic-Jelisavcic$^{2}$, %BEOG-PD 03/06 Bozovicjelisavcic J.~Bracinik$^{3}$, %BIRM-PD 01/2 Bracinik G.~Brandt$^{11}$, %DESY-PD 01/20 Brandt M.~Brinkmann$^{12}$, %HAM2-ST 02/09 Brinkmann V.~Brisson$^{27}$, %ORSA-PD 8/88 Brisson D.~Bruncko$^{16}$, %KOSI-PD 8/88 Bruncko A.~Bunyatyan$^{13,38}$, %MPIH-PD 12/95 Bunyatyan G.~Buschhorn$^{26}$, %MPIM-PD 8/88 Buschhorn L.~Bystritskaya$^{24}$, %ITEP-PD 05/99 Bystritskaya A.J.~Campbell$^{11}$, %DESY-PD 8/88 Campbella K.B. ~Cantun~Avila$^{22}$, %MEX1-ST 04/06 Cantunavila F.~Cassol-Brunner$^{21}$, %MARS-PD 12/0 Cassolbrunner K.~Cerny$^{32}$, %PRG2-ST 09/02 Cernyk V.~Cerny$^{16,47}$, %KOSI-PD 06/04 Cernyv V.~Chekelian$^{26}$, %MPIM-PD 01/90 Chekelian A.~Cholewa$^{11}$, %DESY-ST 11/05 Cholewa J.G.~Contreras$^{22}$, %MEX1-PD 04/97 Contreras J.A.~Coughlan$^{6}$, %RAL -PD 8/88 Coughlan G.~Cozzika$^{10}$, %SACL-PD 10/07 Cozzika J.~Cvach$^{31}$, %PRAG-PD 8/88 Cvach J.B.~Dainton$^{18}$, %LIVE-PD 8/88 Dainton K.~Daum$^{37,43}$, %WUPP-PD 06/96 Daum M.~De\'{a}k$^{11}$, %DESY-ST 08/06 Deak Y.~de~Boer$^{11}$, %DESY-LEFT 08/08 Deboer B.~Delcourt$^{27}$, %ORSA-PD 8/88 Delcourt M.~Del~Degan$^{40}$, %ZUTH-LEFT 09/08 Deldegan J.~Delvax$^{4}$, %BRUX-ST 10/06 Delvax E.A.~De~Wolf$^{4}$, %ANTW-PD 3/93 Dewolf C.~Diaconu$^{21}$, %MARS-PD 01/05 Diaconu V.~Dodonov$^{13}$, %MPIH-PD 04/98 Dodonov A.~Dossanov$^{26}$, %MPIM-ST 01/07 Dossanov A.~Dubak$^{30,46}$, %PODG-PD 10/03 Dubak G.~Eckerlin$^{11}$, %DESY-PD 8/88 Eckerlin V.~Efremenko$^{24}$, %ITEP-PD 8/88 Efremenko S.~Egli$^{36}$, %PSI -PD 01/01 Egli A.~Eliseev$^{25}$, %LPI -PD 01/06 Eliseev E.~Elsen$^{11}$, %DESY-PD 8/88 Elsen A.~Falkiewicz$^{7}$, %CRAC-LEFT 03/09 Falkiewicz L.~Favart$^{4}$, %BRUX-PD 8/88 Favart A.~Fedotov$^{24}$, %ITEP-PD 8/88 Fedotov R.~Felst$^{11}$, %DESY-PD 11/0 Felst J.~Feltesse$^{10,48}$, %SACL-PD 03/05 Feltesse J.~Ferencei$^{16}$, %KOSI-PD 01/05 Ferencei D.-J.~Fischer$^{11}$, %DESY-ST 03/08 Fischer M.~Fleischer$^{11}$, %DESY-PD 07/0 Fleischer A.~Fomenko$^{25}$, %LPI -PD 8/88 Fomenko E.~Gabathuler$^{18}$, %LIVE-PD 10/89 Gabathulere J.~Gayler$^{11}$, %DESY-PD 8/88 Gayler S.~Ghazaryan$^{38}$, %YERE-PD 8/88 Ghazaryan A.~Glazov$^{11}$, %DESY-PD 01/04 Glazov I.~Glushkov$^{39}$, %ZEUT-LEFT 11/08 Glushkov L.~Goerlich$^{7}$, %CRAC-PD 8/88 Goerlich N.~Gogitidze$^{25}$, %LPI -PD 8/88 Gogitidze M.~Gouzevitch$^{11}$, %DESY-PD 12/08 Gouzevitch C.~Grab$^{40}$, %ZUTH-PD 8/88 Grab T.~Greenshaw$^{18}$, %LIVE-PD 8/88 Greenshaw B.R.~Grell$^{11}$, %DESY-ST 09/04 Grell G.~Grindhammer$^{26}$, %MPIM-PD 8/88 Grindhammer S.~Habib$^{12,50}$, %HAM2-ST 12/05 Habib D.~Haidt$^{11}$, %DESY-PD 8/88 Haidt C.~Helebrant$^{11}$, %DFLC-ST 03/06 Helebrant R.C.W.~Henderson$^{17}$, %LANC-PD 8/88 Henderson E.~Hennekemper$^{15}$, %HDB2-ST 11/07 Hennekemper H.~Henschel$^{39}$, %ZEUT-PD 06/99 Henschel M.~Herbst$^{15}$, %HDB2-ST 06/08 Herbst G.~Herrera$^{23}$, %MEX2-PD 07/98 Herrera M.~Hildebrandt$^{36}$, %PSI -PD 10/99 Hildebrandtm K.H.~Hiller$^{39}$, %ZEUT-PD 8/88 Hiller D.~Hoffmann$^{21}$, %MARS-PD 10/0 Hoffmann R.~Horisberger$^{36}$, %PSI -PD 8/88 Horisberger T.~Hreus$^{4,44}$, %BRUX-ST 10/04 Hreus M.~Jacquet$^{27}$, %ORSA-PD 09/96 Jacquet M.E.~Janssen$^{11}$, %DFLC-LEFT 07/08 Janssenm X.~Janssen$^{4}$, %BRUX-PD 02/03 Janssenx L.~J\"onsson$^{20}$, %LUND-PD 8/88 Joensson A.W.~Jung$^{15}$, %HDB2-ST 11/04 Junga H.~Jung$^{11}$, %DESY-PD 07/00 Jungh M.~Kapichine$^{9}$, %JINR-PD 3/97 Kapichine J.~Katzy$^{11}$, %DESY-PD 09/1 Katzy I.R.~Kenyon$^{3}$, %BIRM-PD 8/88 Kenyon C.~Kiesling$^{26}$, %MPIM-PD 8/88 Kiesling M.~Klein$^{18}$, %LIVE-PD 8/88 Klein C.~Kleinwort$^{11}$, %DESY-PD 8/88 Kleinwort T.~Kluge$^{18}$, %LIVE-PD 05/04 Kluge A.~Knutsson$^{11}$, %DESY-PD 04/07 Knutsson R.~Kogler$^{26}$, %MPIM-ST 01/07 Kogler P.~Kostka$^{39}$, %ZEUT-PD 8/88 Kostka M.~Kraemer$^{11}$, %DESY-ST 02/06 Kraemer K.~Krastev$^{11}$, %DESY-LEFT 12/08 Krastev J.~Kretzschmar$^{18}$, %LIVE-PD 01/08 Kretzschmar A.~Kropivnitskaya$^{24}$, %ITEP-ST 07/2 Kropivnitskaya K.~Kr\"uger$^{15}$, %HDB2-PD 01/04 Kruegerk K.~Kutak$^{11}$, %DESY-PD 01/07 Kutak M.P.J.~Landon$^{19}$, %QMWC-PD 8/88 Landon W.~Lange$^{39}$, %ZEUT-PD 8/88 Lange G.~La\v{s}tovi\v{c}ka-Medin$^{30}$, %PODG-PD 06/04 Lastovickamedin P.~Laycock$^{18}$, %LIVE-PD 11/03 Laycock A.~Lebedev$^{25}$, %LPI -PD 8/88 Lebedev G.~Leibenguth$^{40}$, %ZUTH-LEFT 09/08 Leibenguth V.~Lendermann$^{15}$, %HDB2-PD 01/2 Lendermann S.~Levonian$^{11}$, %DESY-PD 8/88 Levonian G.~Li$^{27}$, %ORSA-PD 09/06 Li K.~Lipka$^{11}$, %DESY-PD 01/03 Lipka A.~Liptaj$^{26}$, %MPIM-ST 10/04 Liptaj B.~List$^{12}$, %HAM2-PD 11/99 Listb J.~List$^{11}$, %DFLC-PD 01/05 Listj N.~Loktionova$^{25}$, %LPI -PD 03/99 Loktionova R.~Lopez-Fernandez$^{23}$, %MEX2-PD 03/2 Lopezfernandez V.~Lubimov$^{24}$, %ITEP-PD 01/95 Lubimov L.~Lytkin$^{13}$, %MPIH-LEFT 06/08 Lytkine A.~Makankine$^{9}$, %JINR-PD 11/02 Makankine E.~Malinovski$^{25}$, %LPI -PD 01/89 Malinovskie P.~Marage$^{4}$, %BRUX-PD 8/88 Marage Ll.~Marti$^{11}$, %DESY-LEFT 04/09 Marti H.-U.~Martyn$^{1}$, %AAC1-PD 8/88 Martyn S.J.~Maxfield$^{18}$, %LIVE-PD 8/88 Maxfield A.~Mehta$^{18}$, %LIVE-PD 8/88 Mehta A.B.~Meyer$^{11}$, %DESY-PD 01/00 Meyeran H.~Meyer$^{11}$, %DFLC-LEFT 11/08 Meyerhe H.~Meyer$^{37}$, %WUPP-PD 8/88 Meyerhi J.~Meyer$^{11}$, %DESY-PD 8/88 Meyerj V.~Michels$^{11}$, %DESY-LEFT 08/08 Michels S.~Mikocki$^{7}$, %CRAC-PD 8/88 Mikocki I.~Milcewicz-Mika$^{7}$, %CRAC-ST 10/02 Milcewicz F.~Moreau$^{28}$, %ECPL-PD 01/90 Moreau A.~Morozov$^{9}$, %JINR-PD 06/99 Morozova J.V.~Morris$^{6}$, %RAL -PD 8/88 Morris M.U.~Mozer$^{4}$, %BRUX-PD 06/07 Mozer M.~Mudrinic$^{2}$, %BEOG-PD 01/07 Mudrinic K.~M\"uller$^{41}$, %ZUER-PD 8/88 Muellerk P.~Mur\'\i n$^{16,44}$, %KOSI-LEFT 02/09 Murin Th.~Naumann$^{39}$, %ZEUT-PD 01/89 Naumannt P.R.~Newman$^{3}$, %BIRM-PD 10/92 Newman C.~Niebuhr$^{11}$, %DESY-PD 3/93 Niebuhr A.~Nikiforov$^{11}$, %DESY-PD 05/07 Nikiforov G.~Nowak$^{7}$, %CRAC-PD 8/88 Nowakg K.~Nowak$^{41}$, %ZUER-ST 08/05 Nowakk M.~Nozicka$^{11}$, %DESY-PD 11/06 Nozicka B.~Olivier$^{26}$, %MPIM-LEFT 09/08 Olivier J.E.~Olsson$^{11}$, %DESY-PD 8/88 Olsson S.~Osman$^{20}$, %LUND-ST 02/04 Osman D.~Ozerov$^{24}$, %ITEP-ST 08/98 Ozerov V.~Palichik$^{9}$, %JINR-PD 01/04 Palichik I.~Panagoulias$^{l,}$$^{11,42}$, %DESY-ST 08/04 Panagoulias M.~Pandurovic$^{2}$, %BEOG-ST 03/06 Pandurovic Th.~Papadopoulou$^{l,}$$^{11,42}$, %DESY-PD 06/04 Papadopoulou C.~Pascaud$^{27}$, %ORSA-PD 8/88 Pascaud G.D.~Patel$^{18}$, %LIVE-PD 8/88 Patel O.~Pejchal$^{32}$, %PRG2-LEFT 10/08 Pejchal E.~Perez$^{10,45}$, %SACL-PD 10/07 Perez A.~Petrukhin$^{24}$, %ITEP-ST 01/01 Petrukhin I.~Picuric$^{30}$, %PODG-PD 01/06 Picuric S.~Piec$^{39}$, %ZEUT-ST 01/06 Piec D.~Pitzl$^{11}$, %DESY-PD 8/88 Pitzl R.~Pla\v{c}akyt\.{e}$^{11}$, %DESY-PD 10/06 Placakyte B.~Pokorny$^{12}$, %HAM2-ST 07/08 Pokorny R.~Polifka$^{32}$, %PRG2-ST 10/06 Polifka B.~Povh$^{13}$, %MPIH-PD 8/88 Povh T.~Preda$^{5}$, %BUCH-LEFT 06/08 Preda V.~Radescu$^{11}$, %DESY-PD 10/06 Radescu A.J.~Rahmat$^{18}$, %LIVE-ST 01/05 Rahmat N.~Raicevic$^{30}$, %PODG-PD 03/2 Raicevic A.~Raspiareza$^{26}$, %MPIM-LEFT 03/09 Raspiareza T.~Ravdandorj$^{35}$, %ULBA-PD 06/06 Ravdandorj P.~Reimer$^{31}$, %PRAG-PD 8/88 Reimer E.~Rizvi$^{19}$, %QMWC-PD 01/05 Rizvi P.~Robmann$^{41}$, %ZUER-PD 8/88 Robmann B.~Roland$^{4}$, %BRUX-LEFT 11/08 Roland R.~Roosen$^{4}$, %BRUX-PD 8/88 Roosen A.~Rostovtsev$^{24}$, %ITEP-PD 8/88 Rostovtsev M.~Rotaru$^{5}$, %BUCH-ST 02/07 Rotaru J.E.~Ruiz~Tabasco$^{22}$, %MEX1-ST 09/06 Ruiztabascojuliaelis Z.~Rurikova$^{11}$, %DESY-LEFT 09/08 Rurikova S.~Rusakov$^{25}$, %LPI -PD 8/88 Rusakov D.~\v S\'alek$^{32}$, %PRG2-ST 11/06 Salek D.P.C.~Sankey$^{6}$, %RAL -PD 8/88 Sankey M.~Sauter$^{40}$, %ZUTH-ST 10/05 Sauter E.~Sauvan$^{21}$, %MARS-PD 11/1 Sauvan S.~Schmitt$^{11}$, %DESY-PD 09/07 Schmittst L.~Schoeffel$^{10}$, %SACL-PD 12/98 Schoeffel A.~Sch\"oning$^{14}$, %HDB1-PD 04/09 Schoening H.-C.~Schultz-Coulon$^{15}$, %HDB2-PD 01/04 Schultzcoulon F.~Sefkow$^{11}$, %DFLC-PD 09/99 Sefkow R.N.~Shaw-West$^{3}$, %BIRM-ST 10/04 Shawwest L.N.~Shtarkov$^{25}$, %LPI -PD 8/88 Shtarkov S.~Shushkevich$^{26}$, %MPIM-ST 08/07 Shushkevich T.~Sloan$^{17}$, %LANC-PD 1/96 Sloan I.~Smiljanic$^{2}$, %BEOG-PD 03/06 Smiljanic Y.~Soloviev$^{25}$, %LPI -PD 8/88 Soloviev P.~Sopicki$^{7}$, %CRAC-ST 09/07 Sopicki D.~South$^{8}$, %DORT-PD 06/03 South V.~Spaskov$^{9}$, %JINR-PD 12/97 Spaskov A.~Specka$^{28}$, %ECPL-PD 3/95 Specka Z.~Staykova$^{11}$, %DESY-ST 08/06 Staykova M.~Steder$^{11}$, %DESY-PD 09/08 Steder B.~Stella$^{33}$, %ROME-PD 8/88 Stella G.~Stoicea$^{5}$, %BUCH-PD 02/08 Stoicea U.~Straumann$^{41}$, %ZUER-PD 8/88 Straumann D.~Sunar$^{4}$, %ANTW-ST 03/05 Sunar T.~Sykora$^{4}$, %ANTW-PD 01/06 Sykora V.~Tchoulakov$^{9}$, %JINR-PD 05/03 Tchoulakov G.~Thompson$^{19}$, %QMWC-PD 8/88 Thompsong P.D.~Thompson$^{3}$, %BIRM-PD 08/99 Thompsonp T.~Toll$^{12}$, %HAM2-ST 11/08 Toll F.~Tomasz$^{16}$, %KOSI-LEFT 12/08 Tomasz T.H.~Tran$^{27}$, %ORSA-ST 10/06 Tran D.~Traynor$^{19}$, %QMWC-PD 12/01 Traynor T.N.~Trinh$^{21}$, %MARS-LEFT 10/08 Trinh P.~Tru\"ol$^{41}$, %ZUER-PD 8/88 Truoel I.~Tsakov$^{34}$, %SOFI-PD 04/03 Tsakov B.~Tseepeldorj$^{35,51}$, %ULBA-PD 06/06 Tseepeldorj J.~Turnau$^{7}$, %CRAC-PD 8/88 Turnau K.~Urban$^{15}$, %HDB2-ST 04/05 Urbank A.~Valk\'arov\'a$^{32}$, %PRG2-PD 8/88 Valkarova C.~Vall\'ee$^{21}$, %MARS-PD 8/88 Vallee P.~Van~Mechelen$^{4}$, %ANTW-PD 12/98 Vanmechelen A.~Vargas Trevino$^{11}$, %DFLC-PD 02/07 Vargastrevino Y.~Vazdik$^{25}$, %LPI -PD 8/88 Vazdik S.~Vinokurova$^{11}$, %DESY-LEFT 10/08 Vinokurova V.~Volchinski$^{38}$, %YERE-PD 12/01 Volchinski M.~von~den~Driesch$^{11}$, %DESY-ST 06/08 Vondendriesch D.~Wegener$^{8}$, %DORT-PD 8/88 Wegener Ch.~Wissing$^{11}$, %DESY-PD 07/06 Wissing 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-PD 04/99 Zhokine T.~Zimmermann$^{40}$, %ZUTH-LEFT 01/09 Zimmermannt H.~Zohrabyan$^{38}$, %YERE-PD 11/02 Zohrabyan F.~Zomer$^{27}$, %ORSA-PD 8/88 Zomer and R.~Zus$^{5}$ %BUCH-PD 07/08 Zus %-- H1 Institutes \bigskip{\it $ ^{1}$ I. Physikalisches Institut der RWTH, Aachen, Germany$^{ a}$ \\ $ ^{2}$ Vinca Institute of Nuclear Sciences, Belgrade, Serbia \\ $ ^{3}$ School of Physics and Astronomy, University of Birmingham, Birmingham, UK$^{ b}$ \\ $ ^{4}$ Inter-University Institute for High Energies ULB-VUB, Brussels; Universiteit Antwerpen, Antwerpen; Belgium$^{ c}$ \\ $ ^{5}$ National Institute for Physics and Nuclear Engineering (NIPNE) , Bucharest, Romania \\ $ ^{6}$ Rutherford Appleton Laboratory, Chilton, Didcot, 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}$ Queen Mary and Westfield College, London, UK$^{ b}$ \\ $ ^{20}$ Physics Department, University of Lund, Lund, Sweden$^{ g}$ \\ $ ^{21}$ CPPM, CNRS/IN2P3 - Univ. Mediterranee, Marseille, France \\ $ ^{22}$ Departamento de Fisica Aplicada, CINVESTAV, M\'erida, Yucat\'an, M\'exico$^{ j}$ \\ $ ^{23}$ Departamento de Fisica, CINVESTAV, M\'exico$^{ j}$ \\ $ ^{24}$ Institute for Theoretical and Experimental Physics, Moscow, Russia$^{ k}$ \\ $ ^{25}$ Lebedev Physical Institute, Moscow, Russia$^{ e}$ \\ $ ^{26}$ Max-Planck-Institut f\"ur Physik, M\"unchen, Germany \\ $ ^{27}$ LAL, Univ Paris-Sud, CNRS/IN2P3, Orsay, France \\ $ ^{28}$ LLR, Ecole Polytechnique, IN2P3-CNRS, Palaiseau, France \\ $ ^{29}$ LPNHE, Universit\'{e}s Paris VI and VII, IN2P3-CNRS, Paris, France \\ $ ^{30}$ Faculty of Science, University of Montenegro, Podgorica, Montenegro$^{ e}$ \\ $ ^{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 University of P.J. \v{S}af\'{a}rik, Ko\v{s}ice, Slovak Republic \\ $ ^{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 DESY and University Hamburg, Helmholtz Humboldt Research Award \\ $ ^{49}$ Also at Faculty of Physics, University of Bucharest, Bucharest, Romania \\ $ ^{50}$ Supported by a scholarship of the World Laboratory Bj\"orn Wiik Research Project \\ $ ^{51}$ Also at Ulaanbaatar University, Ulaanbaatar, Mongolia \\ \bigskip $ ^a$ Supported by the Bundesministerium f\"ur Bildung und Forschung, FRG, under contract numbers 05 H1 1GUA /1, 05 H1 1PAA /1, 05 H1 1PAB /9, 05 H1 1PEA /6, 05 H1 1VHA /7 and 05 H1 1VHB /5 \\ $ ^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 PBS/DESY/70/2006 \\ $ ^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-1P05LA259 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 \\ $ ^l$ This project is co-funded by the European Social Fund (75\%) and National Resources (25\%) - (EPEAEK II) - PYTHAGORAS II \\ } \end{flushleft} \newpage \section{Introduction} The measurement of the inclusive charm ($c$) and beauty ($b$) quark cross sections and the derived structure functions $F_2^{c\bar{c}}$ and $F_2^{b\bar{b}}$ in DIS at HERA is an important test of the theory of the strong interaction, quantum chromodynamics (QCD), within the Standard Model. These measurements uniquely constrain the parton density functions (PDFs) of the proton, in particular its $b$ and $c$ content. Precise knowledge of the PDFs is for example essential at the Large Hadron Collider (LHC). The predictions of the `standard candle' QCD processes at the LHC, such as the inclusive production of $W$ and $Z$ bosons, are sensitive to the theoretical treatment of heavy quarks~\cite{collins,VFNS1,VFNS2,VFNS3,VFNS4,cteqvfns,mstw}. The $b$ quark density is important in Higgs production at the LHC in both the Standard Model and in extensions to the Standard Model such as supersymmetric models at high values of the mixing parameter $\tan\beta$\cite{bbh}. This paper reports on measurements made in neutral current deep inelastic scattering (DIS) at HERA of the charm and beauty contributions to the inclusive proton structure function $F_2$ in the range of virtuality of the exchanged photon $5 \le Q^2 \le 2000$~${\rm GeV}^2$ and Bjorken $x$ $0.0002 \le x \le 0.05$. The analysis uses the precise spatial information from the H1 vertex detector to separate events containing $c$ and $b$ flavoured hadrons from light quark events. The analysis extends to lower and higher $Q^2$ than previous H1 measurements~\cite{Aktas:2004az,Aktas:2005iw} which used a similar technique to the one used in this paper. The analysis is based on a dataset with an integrated luminosity of $189$~${\rm pb}^{-1}$, which is about three times greater than in the previous measurements. The data was recorded in the years 2006 and 2007 with $54$~${\rm pb}^{-1}$ taken in $e^-p$ mode and $135$~${\rm pb}^{-1}$ in $e^+p$ mode. The $ep$ centre of mass energy is $\sqrt{s} = 319~{\rm GeV}$, with a proton beam energy of $E_p=920~{\rm GeV}$ and electron beam energy of $E_e=27.6~{\rm GeV}$. This dataset is referred to here as HERA~II. Many details of the analysis are similar to the previous measurements~\cite{Aktas:2004az,Aktas:2005iw}, referred to here as HERA~I. The HERA~I and HERA~II measurements are combined to produce a complete HERA dataset. Measurements of the charm contribution to the proton structure function have also been made at HERA using $D$ or $D^*$ meson production~\cite{H1ZEUSDstar,H1Dstar}. There are also measurements of charm and beauty in DIS using semi-leptonic decays~\cite{H1ZEUSmuons}. Events containing heavy quarks are distinguished from those containing only light quarks using variables that are sensitive to the longer lifetimes of heavy flavour hadrons. The most important of these variables are the transverse displacement of tracks from the primary vertex and the reconstructed position of a secondary vertex in the transverse plane. For events with three or more tracks in the vertex detector the reconstructed variables are used as input to an artificial neural network. This method has better discrimination between $c$ and $b$ compared to previous methods~\cite{Aktas:2004az,Aktas:2005iw}, which used only the transverse displacement of tracks from the primary vertex. Lifetime based methods have the advantage over more exclusive methods, such as $D^*$ or muon tagging, in that a higher fraction of heavy flavour events may be used, although the background from light quark events is larger. The charm structure function $F_2^{c\bar{c}}$ and the beauty structure function $F_2^{b\bar{b}}$ are obtained from the measured $c$ and $b$ cross sections after applying small corrections for the longitudinal structure functions $F_L^{c\bar{c}}$ and $F_L^{b\bar{b}}$. \section{Monte Carlo Simulation} Monte Carlo simulations are used to correct for the effects of the finite detector resolution, acceptance and efficiency. The Monte Carlo program RAPGAP\cite{Jung:1993gf} is used to generate DIS events for the processes $ep \rightarrow eb\bar{b}X$, $ep \rightarrow ec\bar{c}X$ and $ep \rightarrow eq X$ where $q$ is a light quark of flavour $u$, $d$ or $s$. RAPGAP combines $\cal{O}$($\alpha_s$) matrix elements with higher order QCD effects modelled by parton showers. The heavy flavour event samples are generated according to the `massive' photon gluon fusion (PGF) matrix element~\cite{massive} with the mass of the $c$ and $b$ quarks set to $m_c=1.5~{\rm GeV}$ and $m_b=4.75~{\rm GeV}$, respectively. The DIS cross section is calculated using the leading order (LO) 3-flavour PDFs from MRST (MRST2004F3LO~\cite{Martin:2006qz}). The partonic system for the generated events is fragmented according to the Lund string model~\cite{Andersson:1983ia} implemented within the PYTHIA program~\cite{Sjostrand:2001yu}. The $c$ and $b$ quarks are hadronised according to the Bowler fragmentation function~\cite{bowler}. The HERACLES program\cite{Kwiatkowski:1990es} interfaced to RAPGAP calculates single photon radiative emissions off the lepton line, virtual and electroweak corrections. The Monte Carlo program PHOJET\cite{Engel:1995yd} is used to simulate the background contribution from photoproduction $\gamma p \rightarrow X$. The samples of generated events are passed through a detailed simulation of the detector response based on the GEANT3 program\cite{Brun:1978fy}, and through the same reconstruction software as is used for the data. \section{H1 Detector} \label{sec:h1detector} Only a short description of the H1 detector is given here; a more complete description may be found in\cite{Abt:1997xv}. A right handed coordinate system is employed with the origin at the position of the nominal interaction point that has its $Z$-axis pointing in the proton beam, or forward, direction and $X$ ($Y$) pointing in the horizontal (vertical) direction. The pseudorapidity is related to the polar angle $\theta$ by $\eta = - \ln \tan(\theta/2)$. % and PHI Charged particles are measured in the central tracking detector (CTD). This device consists of two cylindrical drift chambers interspersed with $Z$-chambers to improve the $Z$-coordinate reconstruction and multi-wire proportional chambers mainly used for triggering. The CTD is operated in a uniform solenoidal $1.16\,{\rm T}$ magnetic field, enabling the momentum measurement of charged particles over the polar angular range $20^\circ< \theta<160^\circ$. 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 on the $Z$-coordinate of the CST tracks 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{\rm m} \oplus 51 \ \mu{\rm m} /(P_T [{\rm GeV}])$ where $P_T$ is the transverse momentum of the particle. 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 track detectors are surrounded in the forward and central directions ($4^\circ<\theta<155^\circ$) by a fine-grained liquid argon calorimeter (LAr) and in the backward region ($153^\circ<\theta<178^\circ$) by a lead-scintillating fibre calorimeter (SPACAL)~\cite{Nicholls:1996di} with electromagnetic and hadronic sections. These calorimeters are used in this analysis to measure and identify the scattered electron\footnote{In this paper we use `electron' to also denote `positron' unless explicitly stated.} and also provide energy and angular reconstruction for final state particles from the hadronic system. Electromagnetic calorimeters situated downstream in the electron beam direction allow detection of photons and electrons scattered at very low $Q^2$. The luminosity is measured with these calorimeters from the rate of photons produced in the Bethe-Heitler process $ep\rightarrow ep\gamma$. \section{Experimental Method} \subsection{DIS Event Selection} \label{disselect} The events are triggered by requiring a compact, isolated electromagnetic cluster in either the LAr or SPACAL calorimeters with an overall trigger efficiency of almost $100$\%. The electromagnetic cluster with the highest transverse energy, which also passes stricter offline criteria is taken as the scattered electron. The $Z$-position of the interaction vertex, reconstructed by one or more charged tracks in the tracking detectors, must be within $\pm 20~{\rm cm}$ to match the acceptance of the CST. The interaction vertex approximately follows a Gaussian distribution with a standard deviation of $13~{\rm cm}$. Photoproduction events and DIS events with a hard photon radiated from the initial state electron are suppressed by requiring $\sum_{i} (E_i - p_{z,i}) >35~{\rm GeV}$. Here, $E_i$ and $p_{z,i}$ denote the energy and longitudinal momentum components of a particle and the sum is over all final state particles including the scattered electron and the hadronic final state (HFS). The HFS particles are reconstructed using a combination of tracks and calorimeter deposits in an energy flow algorithm that avoids double counting\cite{peezportheault}. The event kinematics, including the photon virtuality $Q^2$, the Bjorken scaling variable $x$ and the inelasticity variable $y$, are reconstructed with the `$e\Sigma$' method\cite{Bassler:1994uq}, which uses the scattered electron and the HFS. The variables obey the relation $x = Q^2/sy$. In order to have good acceptance in the SPACAL and to ensure that the HFS has a significant transverse momentum, events are selected for $Q^2 >4.5 \ {\rm GeV^2}$. The analysis is restricted to $y>0.07$ in order to ensure that the direction of the quark which is struck by the photon is mostly in the CST angular range. An upper $y$ cut is applied that varies from $0.5$ at low $Q^2$ to $0.85$ at high $Q^2$ in order to suppress photoproduction background. The measurement is made differentially by dividing the data into discrete $y$--$Q^2$ intervals. This binning scheme is preferable to one using $x$--$Q^2$ boundaries as it avoids event losses near the cuts on $y$. The RAPGAP Monte Carlo program is used to estimate the acceptance of $c$ ($b$) events that have a $c$ ($b$) quark with tranverse momentum greater than $0.3$~GeV and within the angular acceptance of the CST. The overall $c$ ($b$) quark acceptance is $89\%$ ($90\%$) with a minimum of $75\%$ ($75\%$) in any $y$--$Q^2$ bin. The position of the beam interaction region in $X$--$Y$ (beam spot) is calculated from information of tracks with CST hits and updated regularly to account for drifts during beam storage. The beam interaction region has an elliptical shape with a size of around $90 \ \mu {\rm m}$ in $x$ and around $22 \ \mu {\rm m}$ in $Y$. \subsection{Track Selection} \label{tracksel} The impact parameter $\delta$ of a track, which is the transverse DCA of the track to the primary vertex point (see figure~\ref{fig:alpha}), is only determined for those tracks which are measured in the CTD and have at least two CST hits linked (referred to as CST tracks). The beam spot is used as the position of the primary vertex. CST tracks are required to have a transverse momentum $P_T^{\rm track}> 0.3~{\rm GeV}$. The direction of the struck quark is used in the determination of the sign of $\delta$. The vector of the struck quark ($P_T^q$, $\theta^q$, $\phi^q$) is reconstructed as the azimuthal angle of the highest transverse momentum jet in the event. If there is no jet reconstructed in the event the vector is reconstructed from the electron~\cite{eflow} so that $P_T^q=P_T^{\rm elec}$, $\cos \theta^q =1-8E^2_ey^2/(4E^2_ey^2+Q^2(1-y))$, $\phi^q=180^\circ-\phi^{\rm elec}$, where $P_T^{\rm elec}$ and $\phi^{\rm elec}$ are the transverse momentum and the azimuthal angle in degrees of the scattered electron, respectively. Jets are reconstructed using the inclusive longitudinally invariant $k_T$ algorithm with massless $P_T$ recombination scheme and distance parameter $R_0=1$ in the $\eta-\phi$ plane\cite{KTJET}. The algorithm is run in the laboratory frame using all reconstructed HFS particles and the resultant jets are required to have transverse momentum greater than $1.5$ \ ${\rm GeV}$ and to be in the angular range $15^\circ < \theta < 155^{\rm o}$. Approximately $95\%$ ($99\%$) of $c$ ($b$) events have $\phi^q$ reconstructed from a jet, as determined from the Monte Carlo simulation. Tracks with azimuthal angle $\phi_{\rm track}$ outside $\pm 90^\circ$ of $\phi^q$ are assumed not to be associated to the struck quark and rejected. If the angle $\alpha$ between $\phi^q$ and the line joining the primary vertex to the point of DCA is less than $90^\circ$, $\delta$ is defined as positive. It is defined as negative otherwise. Figure~\ref{fig:alpha} shows an example of a track with positive $\delta$ and one with negative $\delta$. The $\delta$ distribution, shown in figure~\ref{fig:dca}, is seen to be asymmetric with positive values in excess of negative values indicating the presence of long-lived particles. It is found to be well described by the Monte Carlo simulation. CST tracks with $|\delta|>0.1~{\rm cm}$ are rejected to suppress light quark events containing long-lived strange particles. The number of reconstructed CST tracks $N_{\rm track}$ associated to the struck quark is an important quantity since for higher track multiplicities more information can be used. In the kinematic range of this measurement $22\%$ of the events have no selected track, $26\%$ of the events have $N_{\rm track}=1$, $23\%$ have $N_{\rm track}=2$ and $29\%$ have $N_{\rm track} \ge 3$. \subsection{Secondary Vertex Reconstruction} \label{secvertex} The complete set of reconstructed tracks in the event is used to simultaneously reconstruct a secondary and primary vertex in the transverse plane. Two vertices are reconstructed in each event even if they are not statistically separable. Each track is assigned a weight for each vertex, which is a measure of the probability that the track originated at that vertex\cite{svref1}. In this approach tracks are not assigned unambiguously to one vertex or the other. A simultaneous fit to the primary and secondary vertex is made, with the weights of all tracks of the event considered for the primary vertex, but only the weights of CST tracks considered for the secondary vertex. The beam spot together with its spread is used as an additional constraint to the primary vertex. The vertex configuration that minimises the global $\chi^2$ is found iteratively using deterministic annealing\cite{svref2}. The transverse distance between the secondary and primary vertex is defined as $L_{xy}$. The secondary vertex significance $S_L$ is $L_{xy} / \sigma (L_{xy})$, where $\sigma (L_{xy})$ is the uncertainty on $L_{xy}$. If the absolute difference between the azimuthal angle of the secondary vertex (taking the primary vertex as the origin) and $\phi^q$ is less than $90^\circ$, $S_L$ is signed positive and negative otherwise. A measure of the decay multiplicity $N^{\rm SV}_{\rm track}$ is made by counting the number of tracks with weight greater than $0.8$ to the secondary vertex. This method was shown in~\cite{Aktas:2004az} to yield consistent results with the default method that used track significances. \subsection{Quark Flavour Separation} \label{quarkflavourseparation} The track significance is defined as $\delta/\sigma(\delta)$, where $\sigma(\delta)$ is the uncertainty on $\delta$. The significances $S_1$, $S_2$ and $S_3$, are defined as the significance of the CST track with the highest, second highest and third highest absolute significance, respectively. The significances take the sign of $\delta$ (see section~\ref{tracksel}). Tracks with a negative sign for significance are likely not to arise from particles with a large lifetime and are used in this analysis to estimate the light quark contribution. % Tracks that do not have the same significance sign as $S_1$ are ignored. The $S_1$ distribution is used for events with one remaining CST track after this selection and the $S_2$ distribution is used if there are two remaining CST tracks. For events with three or more remaining CST tracks, where more information is available, information is combined from the significance distributions and the reconstructed secondary vertex using an artificial neural network (NN) that takes into account correlations of the input variables~\cite{jetnet}. In this way each event with at least $1$ CST track appears in exactly one distribution and the distributions are statistically independent. The $S_1$ and $S_2$ distributions are shown in figure~\ref{fig:s1s2}. The NN has inputs of $S_1$, $S_2$, $S_3$, $S_L$, $N_{\rm track}$, $N^{\rm SV}_{\rm track}$ and $P_T^{\rm track}$ of the CST tracks with the highest and second highest transverse momentum. The NN has one hidden layer with $5$ nodes. It was trained using a sample of around $5000$ Monte Carlo $b$ events as `signal' and a similar number of Monte Carlo $c$ events as `background'. No attempt is made to discriminate against the light quark events since their impact is minimized by the subtraction procedure described below. The same NN is used for all $y$--$Q^2$ bins. The distributions of the NN inputs are shown in figures~\ref{fig:nninputsa} and \ref{fig:nninputsb}. These distributions are dominated by light quark events. It can be seen that the Monte Carlo simulation gives a reasonable description of these distributions. It is also apparent that these distributions have good separation power between light, $c$ and $b$ events. The decrease in events around zero for the $S_1$ and $S_2$ distributions is due to the requirement that $|S_1|> |S_2|> |S_3|$. In order to see how well the Monte Carlo simulation describes the heavy flavour contribution of the NN input distributions, the distribution for those events with $S_2>3$ is taken and the distribution for those events with $S_2<-3$ is subtracted from it. This has the effect of greatly reducing the light quark distribution which is almost symmetric in $S_2$. This subtraction method can be used for any distribution. The subtracted distributions of the NN inputs are shown in figures~\ref{fig:nninputsnega} and~\ref{fig:nninputsnegb}. The Monte Carlo simulation gives a good description of these distributions. It can be seen that $b$ events tend to have a higher track multiplicity and more tracks at higher $P_T$. Other distributions are checked in a similar way. Distributions of variables related to the struck quark, $P_T^q$ and $\eta^q$, as well as the kinematic variables $\log Q^2$ and $\log x$ are shown in figure~\ref{fig:controlneg} for $N_{\rm track} \ge 2$. The Monte Carlo simulation gives a good description of these distributions. The output of the NN is shown in figure~\ref{fig:nnoutput}. It gives output values in the range from about $0.2$ to $0.95$. The NN output is signed according to the sign of $S_1$. It can be seen that: the light quark distribution is approximately symmetric and peaks towards low absolute values; the $c$ and $b$ distributions are asymmetric with more positive than negative entries; the $b$ events are peaked towards $1$, whereas the $c$ events are peaked towards $0$. The Monte Carlo simulation gives a good description of the distribution. Since the light quark $S_1$, $S_2$ and NN distributions are nearly symmetric around zero the sensitivity to the modelling of the light quarks can be reduced by subtracting the contents of the negative bins from the contents of the corresponding positive bins. The subtracted distributions are shown in figure~\ref{fig:s1s2nnnegsub}. The resulting distributions are dominated by $c$ quark events, with a $b$ quark fraction increasing towards the upper end of the distributions. The light quarks contribute a small fraction, although this fraction is larger than in~\cite{Aktas:2005iw}, mainly due to the lower $Q^2$ and $P_T^{\rm track}$ selections in the present analysis. The fractions of $c$, $b$ and light quarks of the data are extracted in each $y$--$Q^2$ interval using a least squares simultaneous fit to the subtracted $S_1$, $S_2$ and NN distributions (as in figure~\ref{fig:s1s2nnnegsub}) and the total number of inclusive events before any CST track selection. Only those bins in the significance distributions which have at least $25$ events before subtraction are considered in the fit, since Gaussian errors are assumed. The last fitted bin of the significance distributions, which usually has the lowest statistics, is made $3$ times as wide as the other bins (see figure~\ref{fig:s1s2nnnegsub}). If any of the bins before subtraction within the NN output range contain less than $25$ events the bin size is doubled, which ensures all bins contain at least $25$ events. The light, $c$ and $b$ Monte Carlo simulation samples are used as templates. The Monte Carlo light, $c$ and $b$ contributions in each $y$--$Q^2$ interval are scaled by factors $\rho_l$, $\rho_c$ and $\rho_b$, respectively, to give the best fit to the observed subtracted $S_1$, $S_2$ and NN distributions and the total number of inclusive events in each $y-Q^2$ interval. Only the statistical errors of the data and Monte Carlo simulation are considered in the fit. The fit to the subtracted significance and NN distributions mainly constrains $\rho_c$ and $\rho_b$, whereas the overall normalisation constrains $\rho_l$. Since the error on $\rho_c$ is much smaller than that of $\rho_b$ the $c$ cross section is measured in more $y$--$Q^2$ intervals than the $b$ cross section. Therefore two sets of fits are performed, one with a fine binning of $29$ bins and the other with a coarse binning of $12$ bins. The two sets of fits are performed in an identical manner. The results of the two sets of fits are listed in tables~\ref{tab:rhotabledeffine} and~\ref{tab:rhotabledefcoarse}. Also included in the tables are the $\chi^2/$n.d.f. evaluated using statistical errors only. Acceptable values are obtained for all bins. The fit is also performed to the complete data sample and shown in figure~\ref{fig:s1s2nnnegsub}. The stability of the method is checked by repeating the fit to the complete data in a variety of ways: fitting the $e^+p$ and $e^-p$ data separately; using $S_3$ instead of the NN; using $S_L$ instead of the NN; using the NN alone without $S_1$ and $S_2$; using $S_1$ and $S_2$ without the NN; using the NN alone without negative subtraction. All give consistent results within statistical and systematic errors. The results of the fit in each $y$--$Q^2$ interval are converted to a measurement of the `reduced $c$ cross section' defined from the differential cross section as \begin{equation} \tilde{\sigma}^{c\bar{c}} (x, Q^2) = \frac{{\rm d}^2\sigma^{c\bar{c}} }{{\rm d} x\,{\rm d} Q^2} \frac {x Q^4 } {2 \pi \alpha^2 (1+ (1-y)^2)}, \end{equation} \noindent using: \begin{equation} \tilde{\sigma}^{c\bar{c}} (x, Q^2) = \tilde{\sigma} (x, Q^2) \frac{\rho_c N^{\rm MC gen}_c}{\rho_l N^{\rm MC gen}_l+\rho_c N^{\rm MC gen}_c+\rho_b N^{\rm MC gen}_b} \delta_{\rm BCC}, \label{eqn:cross} \end{equation} where $\alpha$ is the fine structure constant evaluated with scale $Q^2$, and $N^{\rm MC gen}_l$, $N^{\rm MC gen}_c$, and $N^{\rm MC gen}_b$ are the number of light, $c$, and $b$ quark events generated from the Monte Carlo simulation in each bin. The inclusive reduced cross section $\tilde{\sigma} (x, Q^2)$ is taken from H1 measurements: Tables~17 and 19 from \cite{h1newf2paper1}, Tables~10 and 11 from \cite{h1newf2paper2} and Table~11 from \cite{H19900NCCC}. A bin centre correction $\delta_{\rm BCC}$ is applied using a NLO QCD expectation for $\tilde{\sigma}^{c\bar{c}}$ and $\tilde{\sigma}$ to convert the bin averaged measurement into a measurement at a given $x$--$Q^2$ point. This NLO QCD expectation is calculated from the results of a fit similar to that performed in~\cite{Adloff:1999ah} but using the massive scheme to generate heavy flavours. The reduced cross section is corrected using the Monte Carlo simulation for pure QED radiative effects. The photoproduction background is not subtracted, which, due to the method used to calculate the reduced cross sections, means that the fraction of $c$ and $b$ events in the photoproduction background is assumed to be the same as in the DIS data. In most of the $y$ range the background from photoproduction events is negligible and does not exceed $4\%$ in any $y$--$Q^2$ interval used in this analysis. Events that contain $c$ hadrons via the decay of $b$ hadrons are not included in the definition of the reduced $c$ cross section. The reduced $b$ cross section is evaluated in the same manner. \section{Systematic Uncertainties} The following uncertainties are taken into account in order to evaluate the systematic error. \label{sec:systematics} \begin{itemize} \item An uncertainty in the $\delta$ resolution of the tracks is estimated by varying the resolution by an amount that encompasses any difference between the data and simulation evaluated from figure~\ref{fig:dca}. This was achieved by applying an additional Gaussian smearing in the Monte Carlo simulation of $200$~$\mu{\rm m}$ to $5\%$ of randomly selected tracks and $12$~$\mu{\rm m}$ to the rest. \item A track efficiency uncertainty is assigned of $1\%$ due to the CTD and of $2\%$ due to the CST. \item The uncertainties on the various $D$ and $B$ meson lifetimes, decay branching fractions and mean charge multiplicities are estimated by varying the input values of the Monte Carlo simulation by the errors on the world average measurements. For the branching fractions of $b$ quarks to hadrons and the lifetimes of the $D$ and $B$ mesons the central values and errors on the world averages are taken from\cite{pdg2006}. For the branching fractions of $c$ quarks to hadrons the values and uncertainties are taken from\cite{Gladilin:1999pj}, which are consistent with measurements made in DIS at HERA\cite{Aktas:2004ka}. For the mean charged track multiplicities the values and uncertainties for $c$ and $b$ quarks are taken from MarkIII\cite{Coffman:1991ud} and LEP/SLD\cite{lepjetmulti} measurements, respectively. \item The uncertainty on the fragmentation function of the heavy quarks is estimated by re\-weighting the events according to the longitudinal string momentum fraction $z$ carried by the heavy hadron in the Lund model using weights of $(1 \mp 0.7)\cdot(1-z) + z \cdot (1 \pm 0.7)$ for charm quarks and of $(1 \mp 0.5)\cdot(1-z) + z \cdot (1 \pm 0.5)$ for beauty quarks. The variations for the charm fragmentation are motivated by comparison of the Monte Carlo simulation with H1 data~\cite{h1frag}. \item An uncertainty on the QCD model of heavy quark production is estimated by reweighting the jet transverse momentum and pseudorapidity by $(P^{\rm jet}_T/(10 ~{\rm GeV}))^{\pm 0.2}$ and $(1 \pm \eta^{\rm jet})^{\pm 0.15}$ for charm jets and $(P^{\rm jet}_T/(10 ~{\rm GeV}))^{\pm 0.3}$ and $(1 \pm \eta^{\rm jet})^{\pm 0.3}$ for bottom jets. These values are obtained by comparing these variations with the measured reduced cross section for $b$ and $c$ jets\cite{rahmatthesis}. The effects of each of these uncertainties on the subtracted reconstructed distributions of $P_T^q$ and $\eta_q$ are shown in figure~\ref{fig:controlneg}, where the data is seen to be consistent with the Monte Carlo simulation within the uncertainties. \item The uncertainty on the asymmetry of the light quark $\delta$ distribution is estimated by repeating the fits with the subtracted light quark distributions (figure~\ref{fig:s1s2nnnegsub}) changed by $\pm30\%$. The light quark asymmetry was checked to be within this uncertainty by comparing the asymmetry of Monte Carlo simulation events to that of the data for $K^0$ candidates, in the region $0.1<|\delta|<0.5~{\rm cm}$, where the light quark asymmetry is enhanced. \item An error on $\phi^q$ is estimated by shifting $\phi^q$ by $2^\circ$($5^\circ$) for events with (without) a reconstructed jet. These shifts were estimated by comparing the difference between $\phi^q$ and the track azimuthal angle in data and Monte Carlo simulation. \item The uncertainty on the hadronic energy scale is estimated by changing the hadronic energy by $\pm 2\%$. \item The uncertainty in the photoproduction background is taken as $100\%$ of the fraction of photoproduction events in each bin, for events with $N_{\rm track} \ge 1$. \item Uncertainties on the acceptance and bin centre correction due to the input heavy quark structure functions used are estimated by reweighting the input $\tilde{\sigma}^{c\bar{c}}$ distribution by $x^{\pm0.1}$ and $1 \pm 0.2 \ln [Q^2/(10~{\rm GeV}^2)]$ and $\tilde{\sigma}^{b\bar{b}}$ by $x^{\pm 0.3}$ and $1 \pm 0.4 \ln [Q^2/(10~{\rm GeV}^2)]$. The range of variation of the input structure functions was estimated by comparing to the measured values obtained in this analysis. The effects of each of these uncertainties on the subtracted distributions of $\log Q^2$ and $\log x$ are shown in figure~\ref{fig:controlneg}, where the data is seen to be consistent with the Monte Carlo simulation within the uncertainties. \end{itemize} The above systematic uncertainties are evaluated by making the changes described above to the Monte Carlo simulation and repeating the procedure to evaluate the reduced $c$ and $b$ cross sections, including the fits. The uncertainties are evaluated separately for each $x$--$Q^2$ measurement. Additional contributions to the systematic error due to the inclusive cross section are taken from the corresponding $x$--$Q^2$ bin in \cite{h1newf2paper1, h1newf2paper2, H19900NCCC}, since the measurements are normalised to the inclusive cross section measurements (see equation~\ref{eqn:cross}). The error due to the inclusive DIS selection includes a $1.1$--$1.5\%$ uncertainty on the luminosity measurement; an uncertainty on the scattered electron polar angle of $0.2$--$3.0~{\rm mrad}$ and energy of $0.2$--$2.0\%$ depending on the energy and angle; typically $<1\%$ combined error due to trigger and scattered electron identification efficiency; and a $0.5$--$1.0\%$ uncertainty on the reduced cross section evaluation due to QED radiative corrections. A detailed list of the systematic effect on each reduced cross section measurement is given in tables~\ref{tab:sig} and~\ref{tab:sysmod}. The errors of $\delta$ resolution and track efficiency are considered uncorrelated between the HERA~I and HERA~II datasets. All other errors are assumed $100\%$ correlated. \section{Results} \label{results} \subsection{Comparison and Combination of Data} The measurements of $\tilde{\sigma}^{c\bar{c}}$ and $\tilde{\sigma}^{b\bar{b}}$ are shown as a function of $x$ for fixed values of $Q^2$ in figures~\ref{fig:f2cc} and \ref{fig:f2bb} respectively. Also shown in these figures are the HERA~I data extracted using measurements based on the displacement of tracks\cite{Aktas:2004az,Aktas:2005iw}. The HERA~I measurements for $Q^2 \le 60 \ {\rm GeV^2}$ are normalised to the H1 inclusive cross sections measurements from~\cite{h1newf2paper1} and~\cite{h1newf2paper2}. The $\tilde{\sigma}^{c\bar{c}}$ and $\tilde{\sigma}^{b\bar{b}}$ data from HERA~I and HERA~II show good agreement for all measured $x$ and $Q^2$ values. In figure~\ref{fig:f2cc} the $\tilde{\sigma}^{c\bar{c}}$ data are also compared with those extracted from $D^*$ meson measurements by H1~\cite{H1Dstar}, which were obtained using a NLO program~\cite{Harris:1997zq} based on DGLAP evolution to extrapolate the measurements outside the visible $D^*$ range. The $D^*$ data agree well with the measurements from the present analysis. The HERA~I and HERA~II datasets are combined for each $x$-$Q^2$ point where there are two measurements. The combination is performed using a weighted mean using those errors considered uncorrelated between the two data sets (the statistical errors and the systematic errors of $\delta$ resolution and track efficiency, see section~\ref{sec:systematics}): \begin{eqnarray} \sigma_{\rm comb}= w_I\sigma_{\rm I}+ w_{II}\sigma_{\rm II}, \end{eqnarray} \begin{eqnarray} {\rm with} \,\,\,\, w_I=\frac{\delta^2_{\rm II}}{\delta^2_{\rm I}+\delta^2_{\rm II}} \,\,\,\, {\rm and} \,\,\,\, w_{II}=\frac{\delta^2_{\rm I}}{\delta^2_{\rm I}+\delta^2_{\rm II}}, \end{eqnarray} where $\sigma_{\rm comb}$ is the combined measurement and $\sigma_{\rm I}$ and $\sigma_{\rm II}$ are the HERA~I and HERA~II measurements respectively, with their respective errors of $\delta_{\rm I}=\sum_i \delta^i_{\rm I}$ and $\delta_{\rm II}=\sum_i \delta^i_{\rm II}$, where $i$ denotes the statistical error and each source of uncorrelated error between the two data sets. These contributions to the combined error are evaluated as: \begin{eqnarray} (\delta^i_{\rm comb})^2= (w_I \delta^i_{\rm I})^2+ (w_{II} \delta^i_{\rm II})^2. \end{eqnarray} For those systematics $j$ considered correlated between the two datasets (all apart from the errors of $\delta$ resolution and track efficiency) the systematic error as evaluated from HERA~II is taken, $\delta^j_{\rm II}$. The total error on the combined measurement is therefore evaluated from: \begin{eqnarray} (\delta^{\rm tot}_{\rm comb})^2=\sum_i (\delta^i_{\rm comb})^2+\sum_j (\delta^j_{\rm II})^2. \end{eqnarray} All data that is subsequently shown is from the combined dataset. Tables~\ref{tab:sig} and~\ref{tab:sysmod} list the combined results for $c$ and $b$. The results listed in these tables supersede the HERA~I measurements published in~\cite{Aktas:2004az} and~\cite{Aktas:2005iw}. If fits are performed with the data listed in these tables it should be noted that, although each $x$--$Q^2$ measurement is statistically independent, there is a correlation between the $c$ and $b$ measurements. Since the $c$ and $b$ binning is very different this correlation may be neglected for most bins. However, for the first and last bin the binning is identical so the correlation coefficient $C_{bc}$ listed in table~\ref{tab:rhotabledeffine} should be used. \subsection{Comparison with QCD} The leading contribution to heavy flavour production in the region $Q^2 \lapprox M^2$ is given by the massive boson-gluon fusion matrix element~\cite{massive,massivenlo} convoluted with the gluon density of the proton. In the region where $Q^2$ is much larger than $M^2$ the massive approach may be a poor approximation due to the large logarithms $\log Q^2/M^2$ which are not resummed~\cite{VFNS1}. Here, the heavy quarks can be treated as massless partons with the leading order contribution coming from the quark parton model and the heavy quark parton densities. In QCD fits to global hard-scattering data the parton density functions are usually extracted using the general mass variable flavour number scheme (GM VFNS)~\cite{collins,VFNS1,VFNS2,VFNS3} for heavy quarks which interpolates from the massive approach at low scales to the massless approach at high scales. The data are compared with recent QCD predictions based on the GM VFNS from MSTW\cite{mstw} (at NLO and NNLO), CTEQ\cite{cteqvfns} (at NLO) and from H1\cite{h1newf2paper2} (at NLO). The MSTW predictions use the MSTW08 PDFs which have $m_c=1.4~{\rm GeV}$ and $m_b=4.75~{\rm GeV}$, and the renormalization and factorization scales are set to $\mu_r = \mu_f = Q$. The CTEQ predictions use the CTEQ6.6 PDFs where $m_c=1.3~{\rm GeV}$, $m_b=4.5~{\rm GeV}$ and $\mu_r = \mu_f = \sqrt{Q^2+M^2}$. The predictions from H1 use the H1PDF 2009 PDFs and the same heavy flavour treatment, including the quark masses and perturbative scales, as for the MSTW08 NLO predictions~\cite{mstw}. The data are also compared with predictions based on CCFM\cite{ccfm2} parton evolution and %CASCADE is an implementation of the CCFM\cite{ccfm2} %evolution equation and uses off shell matrix elements convoluted with %$k_t$ unintegrated proton parton distributions. massive heavy flavour production. The CCFM predictions use the A0 PDF set~\cite{a0} with $m_c=1.4~{\rm GeV}$, $m_b=4.75~{\rm GeV}$ and $\mu_r = \mu_f = \sqrt{\hat{s}+Q_{T}^2}$, where $\hat{s}$ is the square of the partonic centre of mass energy and $Q_{T}$ is the transverse momentum of the heavy quark pair. The $\tilde{\sigma}^{c\bar{c}}$ data as a function of $x$ for fixed values of $Q^2$ are compared in figure~\ref{fig:f2ccav} with the QCD predictions from CCFM, CTEQ and MSTW at NNLO, and in figure~\ref{fig:f2ccav2} with the predictions from H1 and MSTW at NLO. In figure~\ref{fig:f2ccav} the GM VFNS predictions from CTEQ and MSTW at NNLO are observed to be similar with the size of the largest differences between the two being at the level of the total experimental errors on the data. The CTEQ and MSTW predictions provide a reasonable description of the rise of the data with decreasing $x$ across the whole of the measured kinematic range thus supporting the validity of PDFs extracted using the GM VFNS. The predictions based on CCFM evolution tend to undershoot the data at the lowest values of $Q^2$ and $x$ but also provide a reasonable description for the rest of the measured phase space. In figure~\ref{fig:f2ccav2} the GM VFNS predictions from H1 and MSTW at NLO, which implement the same heavy flavour treatment~\cite{mstw}, are similar and also provide a reasonable description of the data. The H1 predictions are shown with uncertainty bands representing the experimental and theoretical uncertainties~\cite{h1newf2paper2}. The inner error band describes the experimental fit uncertainty, the middle error band represents the experimental and model uncertainties added in quadrature and the outer error band represents the fit parameterisation uncertainty added in quadrature with all the other uncertainties. The largest contribution to the uncertainty comes from the model which is dominated by the variation of the charm quark mass, which is varied from $1.38$ to $1.47\ {\rm GeV}$. The total uncertainties on the H1PDF 2009 reduced charm cross section predictions are generally smaller than those on the data. The $\tilde{\sigma}^{b\bar{b}}$ data as a function of $x$ for fixed values of $Q^2$ are compared with the QCD predictions in figures~\ref{fig:f2bbav} and~\ref{fig:f2bbav2}. In figure~\ref{fig:f2bbav} the CTEQ and CCFM predictions are seen to be very similar across the whole range of the measurements. The MSTW NNLO predictions are around $35\%$ higher than CTEQ and CCFM at low values of $Q^2$, with the difference decreasing with increasing values of $Q^2$. The differences between the theory predictions for the $b$ cross section at low $Q^2$ are much reduced compared with the theoretical status at the time of the HERA~I publication where there was a factor $2$ difference at $Q^2 = 12 \ {\rm GeV^2}$~\cite{Aktas:2005iw}. In figure~\ref{fig:f2bbav2} the MSTW and H1 NLO QCD predictions for $\tilde{\sigma}^{b\bar{b}}$ are observed, as for the case of $\tilde{\sigma}^{c\bar{c}}$, to be very similar. The uncertainty on the H1 predictions is again dominated by the model uncertainty due to the variation of the $b$ quark mass, which is varied from $4.3$ to $5.0\ {\rm GeV}$. At lower values of $Q^2$ the uncertainties on the H1 PDF predictions are larger than those on the data. The reduced $b$ cross section data, including the points in the newly measured regions, are well described by all the present QCD predictions. The structure function $F_2^{c\bar{c}}$ is evaluated from the reduced cross section \begin{equation} \tilde{\sigma}^{c\bar{c}} = F_2^{c\bar{c}} - \frac{y^2 }{1+ (1-y)^2} F_L^{c\bar{c}}, \label{eq:sigcc} \end{equation} where the longitudinal structure function $F_L^{c\bar{c}}$ is estimated from the same NLO QCD expectation as used for the bin centre correction. The correction due to $F_L^{c\bar{c}}$ is negligible for most bins but contributes up to $6.7$\% of the reduced cross section at the highest value of $y$. The structure function $F_2^{b\bar{b}}$ is evaluated in the same manner. The measurements of $F_2^{c\bar{c}}$ and $F_2^{b\bar{b}}$ are presented in table~\ref{tab:sig} and shown as a function of $Q^2$ in figure~\ref{fig:f2ccq2} and figure~\ref{fig:f2bbq2} respectively. The data are compared with the GM VFNS QCD predictions from CTEQ\cite{cteqvfns} at NLO and from MSTW at NLO and NNLO\cite{mstw}. The description of the charm data by the MSTW QCD calculations is reasonable, with the NNLO being somewhat better than NLO. The CTEQ NLO prediction also gives a reasonable description of the data. The measurements are presented in figure~\ref{fig:frac} in the form of the fractional contribution to the total $ep$ cross section \begin{equation} f^{c\bar{c}} = \frac{{\rm d}^2 \sigma^{c\bar{c}}} {{\rm d} x\, {\rm d} Q^2} / \frac{ {\rm d}^2 \sigma}{ {\rm d} x\, {\rm d} Q^2 }. \end{equation} The $b$ fraction $f^{b\bar{b}}$ is defined in the same manner. In the present kinematic range the value of $f^{c\bar{c}}$ is around $17\%$ on average and increases slightly with increasing $Q^2$ and decreasing $x$. The value of $f^{b\bar{b}}$ increases rapidly with $Q^2$ from about $0.2\%$ at $Q^2 = 5 \ {\rm GeV^2}$ to around $1\%$ for $Q^2 \gapprox 60 \ {\rm GeV^2}$. The NNLO QCD predictions of MSTW shown in figure~\ref{fig:frac} are found to describe the data well. \section{Conclusion} The reduced charm and beauty cross sections in deep inelastic scattering are measured for a wide range of $Q^2$ and Bjorken $x$ using the HERA~II data. The analysis was performed using several variables including the significance (the impact parameter divided by its error) and the position of the secondary vertex as reconstructed from the vertex detector. For selected track multiplicities of $1$ or $2$ the highest and second highest significance distributions are used to evaluate the charm and beauty content of the data. For selected track multiplicities $\ge 3$ several variables are combined using an artificial neural network. The reduced cross sections agree with previous measurements using a similar technique, but have reduced errors and cover an extended $Q^2$ range. HERA~I and HERA~II data are combined resulting in more precise reduced cross section and structure function measurements. The charm and beauty fractional contributions to the total $ep$ cross section are also measured. In this kinematic range the charm cross section contributes on average $17\%$ and the beauty fraction increases from about $0.2\%$ at $Q^2 = 5 \ {\rm GeV^2}$ to $1.0\%$ for $Q^2 \gapprox 60 \ {\rm GeV^2}$. The measurements are described by predictions using perturbative QCD in the general mass variable flavour number scheme at NLO and NNLO. \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. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{thebibliography}{99} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{collins} J.~C.~Collins, %``Hard-scattering factorization with heavy quarks: A general treatment,'' Phys.\ Rev.\ D {\bf 58} (1998) 094002 \ifpdf \href{http://arXiv.org/abs/hep-ph/9806259}{[hep-ph/9806259].} \else [hep-ph/9806259]. \fi %%CITATION = PHRVA,D58,094002;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % ACOT \bibitem{VFNS1} % SACOT %\bibitem{Aivazis:1993kh} M.~A.~G.~Aivazis, F.~I.~Olness and W.~K.~Tung, %``Leptoproduction of heavy quarks. 1. General formalism and kinematics of %charged current and neutral current production processes,'' Phys.\ Rev.\ D {\bf 50} (1994) 3085 \ifpdf \href{http://arXiv.org/abs/hep-ph/9312318}{[hep-ph/9312318];} \else [hep-ph/9312318]; \fi %%CITATION = HEP-PH 9312318;%% % SACOT %\bibitem{Aivazis:1993pi} M.~A.~G.~Aivazis, J.~C.~Collins, F.~I.~Olness and W.~K.~Tung, %``Leptoproduction of heavy quarks. 2. A Unified QCD formulation of charged and %neutral current processes from fixed target to collider energies,'' Phys.\ Rev.\ D {\bf 50} (1994) 3102 \ifpdf \href{http://arXiv.org/abs/hep-ph/9312319}{[hep-ph/9312319];} \else [hep-ph/9312319]; \fi %%CITATION = HEP-PH 9312319;%% % SACOT %\bibitem{Kramer:2000hn} M.~Kr\"{a}mer, F.~I.~Olness and D.~E.~Soper, %``Treatment of heavy quarks in deeply inelastic scattering,'' Phys.\ Rev.\ D {\bf 62} (2000) 096007 \ifpdf \href{http://arXiv.org/abs/hep-ph/0003035}{[hep-ph/0003035];} \else [hep-ph/0003035]; \fi %%CITATION = HEP-PH 0003035;%% %\bibitem{Tung:2001mv} W.~K.~Tung, S.~Kretzer and C.~Schmidt, %``Open heavy flavor production in QCD: Conceptual framework and %implementation issues,'' J.\ Phys.\ G {\bf 28} (2002) 983 [hep-ph/0110247]; %%CITATION = JPHGB,G28,983;%% %\bibitem{Tung:2006tb} W.~K.~Tung, H.~L.~Lai, A.~Belyaev, J.~Pumplin, D.~Stump and C.~P.~Yuan, %``Heavy quark mass effects in deep inelastic scattering and global QCD %analysis,'' JHEP {\bf 0702} (2007) 053 [hep-ph/0611254]. %%CITATION = JHEPA,0702,053;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{VFNS2} %\bibitem{Thorne:1997ga} R.~S.~Thorne and R.~G.~Roberts, %``An ordered analysis of heavy flavour production in deep inelastic %scattering,'' Phys.\ Rev.\ D {\bf 57} (1998) 6871 \ifpdf \href{http://arXiv.org/abs/hep-ph/9709442}{[hep-ph/9709442];} \else [hep-ph/9709442]; \fi %%CITATION = HEP-PH 9709442;%% %\bibitem{Thorne:1997uu} R.~S.~Thorne and R.~G.~Roberts, %``A practical procedure for evolving heavy flavour structure functions,'' Phys.\ Lett.\ B {\bf 421} (1998) 303 \ifpdf \href{http://arXiv.org/abs/hep-ph/9711223}{[hep-ph/9711223];} \else [hep-ph/9711223]; \fi %%CITATION = HEP-PH 9711223;%% %\bibitem{NNLO} %\cite{Thorne:2006qt} R.~S.~Thorne, %``A variable-flavour number scheme for NNLO,'' Phys.\ Rev.\ D {\bf 73} (2006) 054019 \ifpdf \href{http://arXiv.org/abs/hep-ph/0601245}{[hep-ph/0601245].} \else [hep-ph/0601245]. \fi %%CITATION = PHRVA,D73,054019;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{VFNS3} %%Chuvakin:2000zj below references 2 O(alpha_s^2) VFNS treatments %\bibitem{Buza:1997nv} M.~Buza, Y.~Matiounine, J.~Smith and W.~L.~van Neerven, %``Comparison between the various descriptions for charm electroproduction and %the HERA-data,'' Phys.\ Lett.\ B {\bf 411} (1997) 211 \ifpdf \href{http://arXiv.org/abs/hep-ph/9707263}{[hep-ph/9707263];} \else [hep-ph/9707263]; \fi %%CITATION = HEP-PH 9707263;%% %\bibitem{Chuvakin:1999nx} A.~Chuvakin, J.~Smith and W.~L.~van Neerven, %``Comparison between variable flavor number schemes for charm quark %electroproduction,'' Phys.\ Rev.\ D {\bf 61} (2000) 096004 \ifpdf \href{http://arXiv.org/abs/hep-ph/9910250}{[hep-ph/9910250];} \else [hep-ph/9910250]; \fi %%CITATION = HEP-PH 9910250;%% %\bibitem{Chuvakin:2000jm} A.~Chuvakin, J.~Smith and W.~L.~van Neerven, %``Bottom quark electroproduction in variable flavor number schemes,'' Phys.\ Rev.\ D {\bf 62} (2000) 036004 \ifpdf \href{http://arXiv.org/abs/hep-ph/0002011}{[hep-ph/0002011].} \else [hep-ph/0002011]. \fi %%CITATION = HEP-PH 0002011;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{VFNS4} %\cite{Kretzer:2003it} S.~Kretzer, H.~L.~Lai, F.~I.~Olness and W.~K.~Tung, %``CTEQ6 parton distributions with heavy quark mass effects,'' Phys.\ Rev.\ D {\bf 69} (2004) 114005 [hep-ph/0307022]; %%CITATION = PHRVA,D69,114005;%% %\bibitem{Martin:2007bv} A.~D.~Martin, W.~J.~Stirling, R.~S.~Thorne and G.~Watt, %``Update of parton distributions at NNLO,'' Phys.\ Lett.\ B {\bf 652} (2007) 292 [arXiv:0706.0459]; %%CITATION = PHLTA,B652,292;%% S.~Alekhin and S.~Moch, Phys.\ Lett.\ B {\bf 672} (2009) 166 %``Higher order QCD corrections to charged-lepton deep-inelastic scattering %and global fits of parton distributions,'' \ifpdf \href{http://arXiv.org/abs/arXiv:0811.1412}{[arXiv:0811.1412];} \else [arXiv:0811.1412]; \fi %%CITATION = ARXIV:0811.1412;%% %\bibitem{Bierenbaum:2009zt} I.~Bierenbaum, J.~Blumlein and S.~Klein, %``The Gluonic Operator Matrix Elements at O(\alpha_s^2) for DIS Heavy Flavor %Production,'' Phys.\ Lett.\ B {\bf 672} (2009) 401 [arXiv:0901.0669]; %%CITATION = PHLTA,B672,401;%% %\bibitem{Bierenbaum:2009mv} I.~Bierenbaum, J.~Blumlein and S.~Klein, %``Mellin Moments of the {$O(\alpha_s^3$)} Heavy Flavor Contributions to %unpolarized Deep-Inelastic Scattering at $Q^2 \gg m^2$ and Anomalous %Dimensions,'' Nucl.\ Phys.\ B {\bf 820} (2009) 417 [arXiv:0904.3563]; %%CITATION = NUPHA,B820,417;%% %\bibitem{Nadolsky:2009ge} P.~M.~Nadolsky and W.~K.~Tung, %``Improved Formulation of Global QCD Analysis with Zero-mass Matrix %Elements,'' Phys.\ Rev.\ D {\bf 79} (2009) 113014 [arXiv:0903.2667]; %%CITATION = PHRVA,D79,113014;%% %\bibitem{ThorneTung} R.~S.~Thorne and W.~K.~Tung, %``PQCD Formulations with Heavy Quark Masses and Global Analysis'', proceedings of the HERA and the LHC workshop series on the implications of HERA for LHC physics (2009) 332 \ifpdf \href{http://arXiv.org/abs/arXiv:0903.3861}{[arXiv:0903.3861].} \else [arXiv:0903.3861]. \fi %%CITATION = ARXIV:0903.3861;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{cteqvfns} %\cite{Nadolsky:2008zw} %\bibitem{Nadolsky:2008zw} P.~M.~Nadolsky {\it et al.}, %``Implications of CTEQ global analysis for collider observables,'' Phys.\ Rev.\ D {\bf 78} (2008) 013004 \ifpdf \href{http://arXiv.org/abs/arXiv:0802.0007}{[arXiv:0802.0007].} \else [arXiv:0802.0007]. \fi %%CITATION = ARXIV:0802.0007;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{mstw} A.~D.~Martin, W.~J.~Stirling, R.~S.~Thorne and G.~Watt, %``Parton distributions for the LHC,'' \ifpdf \href{http://arXiv.org/abs/arXiv:0901.0002}{arXiv:0901.0002.} \else arXiv:0901.0002. \fi %%CITATION = ARXIV:0901.0002;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{bbh} D.~Dicus, T.~Stelzer, Z.~Sullivan, and S.~Willenbrock, Phys. Rev. D~{\bf 59} (1999) 094016; C.~S.~Huang and S.~H.~Zhu, Phys. Rev. D~{\bf 60} (1999) 075012; C.~Balazs, H.~J.~He, and C.~P.~Yuan, Phys. Rev. D~{\bf 60} (1999) 114001; J.~Campbell, R.~K.~Ellis, F.~Maltoni, and S.~Willenbrock, Phys. Rev. D~{\bf 67} (2003) 095002; F.~Maltoni, Z.~Sullivan, and S.~Willenbrock, Phys. Rev. D~{\bf 67} (2003) 093005. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{Aktas:2004az} A.~Aktas {\it et al.} [H1 Collaboration], %``Measurement of F2(c anti-c) and F2(b anti-b) at high Q**2 using the H1 %vertex detector at HERA,'' Eur.\ Phys.\ J.\ C {\bf 40} (2005) 349 \ifpdf \href{http://arXiv.org/abs/hep-ex/0411046}{[hep-ex/0411046].} \else [hep-ex/0411046]. \fi %%CITATION = HEP-EX 0411046;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{Aktas:2005iw} A.~Aktas {\it et al.} [H1 Collaboration], %``Measurement of F2(c anti-c) and F2(b anti-b) at low Q**2 and x using the %H1 vertex detector at HERA,'' Eur.\ Phys.\ J.\ C {\bf 45} (2006) 23 \ifpdf \href{http://arXiv.org/abs/hep-ex/0507081}{[hep-ex/0507081].} \else [hep-ex/0507081]. \fi %%CITATION = EPHJA,C45,23;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%% OTHER MEASUREMENTS %\cite{Adloff:1996xq} \bibitem{H1ZEUSDstar} C.~Adloff {\it et al.} [H1 Collaboration], %``Inclusive D0 and D*+- production in neutral current deep inelastic e p %scattering at HERA,'' Z.\ Phys.\ C {\bf 72} (1996) 593 \ifpdf \href{http://arXiv.org/abs/hep-ex/9607012}{[hep-ex/9607012];} \else [hep-ex/9607012]; \fi %%CITATION = HEP-EX 9607012;%% %%\cite{Breitweg:1997mj} J.~Breitweg {\it et al.} [ZEUS Collaboration], %``D* production in deep inelastic scattering at HERA,'' Phys.\ Lett.\ B {\bf 407} (1997) 402 \ifpdf \href{http://arXiv.org/abs/hep-ex/9706009}{[hep-ex/9706009];} \else [hep-ex/9706009]; \fi %%CITATION = HEP-EX 9706009;%% %%\cite{Adloff:1998vb} C.~Adloff {\it et al.} [H1 Collaboration], %``Measurement of D* meson cross sections at HERA and determination of the %gluon density in the proton using NLO QCD,'' Nucl.\ Phys.\ B {\bf 545} (1999) 21 \ifpdf \href{http://arXiv.org/abs/hep-ex/9812023}{[hep-ex/9812023];} \else [hep-ex/9812023]; \fi %%CITATION = HEP-EX 9812023;%% %%\cite{Breitweg:1999ad} J.~Breitweg {\it et al.} [ZEUS Collaboration], %``Measurement of D*+- production and the charm contribution to F2 in deep %inelastic scattering at HERA,'' Eur.\ Phys.\ J.\ C {\bf 12} (2000) 35 \ifpdf \href{http://arXiv.org/abs/hep-ex/9908012}{[hep-ex/9908012];} \else [hep-ex/9908012]; \fi %%CITATION = HEP-EX 9908012;%% %\bibitem{ZEUSDstar} S.~Chekanov {\it et al.} [ZEUS Collaboration], %``Measurement of D*+- production in deep inelastic e+- p scattering at HERA,'' Phys.\ Rev.\ D {\bf 69} (2004) 012004 \ifpdf \href{http://arXiv.org/abs/hep-ex/0308068}{[hep-ex/0308068];} \else [hep-ex/0308068]; \fi %%CITATION = HEP-EX 0308068;%% %\bibitem{Chekanov:2008yd} S.~Chekanov {\it et al.} [ZEUS Collaboration], %``Measurement of D+- and D0 production in deep inelastic scattering using a %lifetime tag at HERA,'' DESY-08-201, \ifpdf \href{http://arXiv.org/abs/arXiv:0812.3775}{arXiv:0812.3775.} \else arXiv:0812.3775. \fi %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{H1Dstar} C.~Adloff {\it et al.} [H1 Collaboration], %``Measurement of D*+- meson production and F2(c) in deep inelastic scattering %at HERA,'' Phys.\ Lett.\ B {\bf 528} (2002) 199 \ifpdf \href{http://arXiv.org/abs/hep-ex/0108039}{[hep-ex/0108039].} \else [hep-ex/0108039]. \fi %%CITATION = HEP-EX 0108039;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{H1ZEUSmuons} %\bibitem{Chekanov:2004tk} S.~Chekanov {\it et al.} [ZEUS Collaboration], %``Measurement of beauty production in deep inelastic scattering at HERA,'' Phys.\ Lett.\ B {\bf 599} (2004) 173 \ifpdf \href{http://arXiv.org/abs/hep-ex/0405069}{[hep-ex/0405069];} \else [hep-ex/0405069]; \fi %%CITATION = PHLTA,B599,173;%% %\bibitem{Aktas:2005zc} A.~Aktas {\it et al.} [H1 Collaboration], %``Measurement of beauty production at HERA using events with muons and %jets,'' Eur.\ Phys.\ J.\ C {\bf 41} (2005) 453 \ifpdf \href{http://arXiv.org/abs/hep-ex/0502010}{[hep-ex/0502010];} \else [hep-ex/0502010]; \fi %%CITATION = EPHJA,C41,453;%% %\bibitem{Chekanov:2006sg} S.~Chekanov {\it et al.} [ZEUS Collaboration], %``Measurement of open beauty production at HERA in the $D^{*}$ muon final %state,'' Eur.\ Phys.\ J.\ C {\bf 50} (2007) 299 \ifpdf \href{http://arXiv.org/abs/hep-ex/0609050}{[hep-ex/0609050];} \else [hep-ex/0609050]; \fi %%CITATION = EPHJA,C50,299;%% %\bibitem{Chekanov:2009kj} S.~Chekanov {\it et al.} [ZEUS Collaboration], %``Measurement of charm and beauty production in deep inelastic ep scattering %from decays into muons at HERA,'' DESY-09-056, \ifpdf \href{http://arXiv.org/abs/arXiv:0904.3487}{arXiv:0904.3487.} \else arXiv:0904.3487. \fi %%CITATION = ARXIV:0904.3487;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%% RAPGAP MC \bibitem{Jung:1993gf} H.~Jung, %``Hard diffractive scattering in high-energy e p collisions and the Monte Carlo %generation RAPGAP,'' RAPGAP 3.1, Comput.\ Phys.\ Commun.\ {\bf 86} (1995) 147\\ %%CITATION = CPHCB,86,147;%% \ifpdf (see also \href{http://www.desy.de/~jung/rapgap/} {http://www.desy.de/$\sim$jung/rapgap/}). \else (see also http://www.desy.de/$\sim$jung/rapgap/). \fi %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % FFNS \bibitem{massive} % leading order %\bibitem{Ingelman:1988qn} G.~Ingelman and G.~A.~Schuler, %``CHARACTERISTICS OF HEAVY FLAVOR PRODUCTION IN e p COLLISIONS,'' Z.\ Phys.\ C {\bf 40} (1988) 299; %%CITATION = ZEPYA,C40,299;%% %\bibitem{Schuler:1987wj} G.~A.~Schuler, %``HEAVY FLAVOR PRODUCTION AT HERA,'' Nucl.\ Phys.\ B {\bf 299} (1988) 21. %%CITATION = NUPHA,B299,21;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%% fixed flavour \bibitem{Martin:2006qz} A.~D.~Martin, W.~J.~Stirling and R.~S.~Thorne, %``MRST partons generated in a fixed-flavour scheme,'' Phys.\ Lett.\ B {\bf 636} (2006) 259 \ifpdf \href{http://arXiv.org/abs/hep-ph/0603143}{[hep-ph/0603143].} \else [hep-ph/0603143]. \fi %%CITATION = PHLTA,B636,259;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Lund String model \bibitem{Andersson:1983ia} B.~Andersson, G.~Gustafson, G.~Ingelman and T.~Sj\"{o}strand, %``Parton Fragmentation And String Dynamics,'' Phys.\ Rept.\ {\bf 97} (1983) 31. %%CITATION = PRPLC,97,31;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%% pythia 6.2 %\cite{Sjostrand:2001yu} \bibitem{Sjostrand:2001yu} T.~Sj\"{o}strand {\it et al.} Comput.\ Phys.\ Commun.\ {\bf 135} (2001) 238; T.~Sj\"{o}strand, L.~L\"{o}nnblad and S.~Mrenna, ``PYTHIA 6.2: Physics and manual,'' \ifpdf \href{http://arXiv.org/abs/hep-ph/0108264}{hep-ph/0108264.} \else hep-ph/0108264. \fi %%CITATION = HEP-PH/0108264;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{bowler} M.~G.~Bowler, Z.\ Phys.C {\bf 11} (1981) 169. %%CITATION=ZEPYA,C11,169;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%% Heracles \bibitem{Kwiatkowski:1990es} A.~Kwiatkowski, H.~Spiesberger and H.~J.~M\"{o}hring, %``Heracles: An Event Generator For E P Interactions At Hera Energies Including %Radiative Processes: Version 1.0,'' Comput.\ Phys.\ Commun.\ {\bf 69} (1992) 155. %%CITATION = CPHCB,69,155;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %\cite{Engel:1995yd} \bibitem{Engel:1995yd} R.~Engel and J.~Ranft, %``Hadronic photon-photon interactions at high energies,'' Phys.\ Rev.\ D {\bf 54} (1996) 4244 \ifpdf \href{http://arXiv.org/abs/hep-ph/9509373}{[hep-ph/9509373].} \else [hep-ph/9509373]. \fi %%CITATION = HEP-PH 9509373;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%% GEANT \bibitem{Brun:1978fy} %R.~Brun, R.~Hagelberg, M.~Hansroul and J.~C.~Lassalle, R.~Brun {\it et al.}, GEANT 3, %``Geant: Simulation Program For Particle Physics Experiments. User Guide And %Reference Manual,'' CERN-DD-78-2-REV. %\href{http://www.slac.stanford.edu/spires/find/hep/www?r=cern-dd-78-2-rev}{SPIRES entry} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{Abt:1997xv} I.~Abt {\it et al.} [H1 Collaboration], %``The Tracking, calorimeter and muon detectors of the H1 experiment at HERA ,'' Nucl.\ Instrum.\ Meth.\ A {\bf 386} (1997) 310; %CITATION = NUIMA,A386,348;%% I.~Abt {\it et al.} [H1 Collaboration], Nucl.\ Instrum.\ Meth.\ A {\bf 386} (1997) 348. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{cst} % CST paper D.~Pitzl {\it et al.}, %``The H1 silicon vertex detector,'' Nucl.\ Instrum.\ Meth.\ A {\bf 454} (2000) 334 \ifpdf \href{http://arXiv.org/abs/hep-ex/0002044}{[hep-ex/0002044];} \else [hep-ex/0002044]; \fi %%CITATION = HEP-EX 0002044;%% B.~List, %``The H1 central silicon tracker,'' Nucl.\ Instrum.\ Meth.\ A {\bf 501} (2001) 49. %%CITATION = NUIMA,A501,49;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{Nicholls:1996di} T.~Nicholls {\it et al.} [H1 SPACAL Group], %``Performance of an electromagnetic lead / scintillating %fiber calorimeter for the H1 detector,'' Nucl.\ Instrum.\ Meth.\ A {\bf 374} (1996) 149. %%CITATION = NUIMA,A374,149;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{peezportheault} M.~Peez, ``Recherche de d\'eviations au Mod\`ele Standard dans les processus de grande \'energie transverse sur le collisionneur \'electron-proton HERA'', Ph.D. thesis. Universit\'e de Lyon (2003), DESY-THESIS-2003-023 \ifpdf (available at \href{http://www-h1.desy.de/publications/theses_list.html} {http://www-h1.desy.de/publications/theses\_list.html}); \else (available at http://www-h1.desy.de/publications/theses\_list.html); \fi B.~Portheault, ``Premi\`ere mesure des sections efficaces de courant charg\'e et neutre avec le faisceau de positrons polaris\'e \'a HERA~II et analyses QCD-\'electrofaibles'', Ph.D. thesis. Universit\'e Paris XI (2005), LAL-05-05 \ifpdf (available at \href{http://www-h1.desy.de/publications/theses_list.html} {http://www-h1.desy.de/publications/theses\_list.html}. \else (available at http://www-h1.desy.de/publications/theses\_list.html). \fi %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %\cite{Bassler:1994uq} \bibitem{Bassler:1994uq} U.~Bassler and G.~Bernardi, %``On the kinematic reconstruction of deep inelastic scattering at HERA: The %Sigma method,'' Nucl.\ Instrum.\ Meth.\ A {\bf 361} (1995) 197 \ifpdf \href{http://arXiv.org/abs/hep-ex/9412004}{[hep-ex/9412004].} \else [hep-ex/9412004]. \fi %%CITATION = HEP-EX 9412004;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{eflow} T.~Ahmed {\it et al.} [H1 Collaboration], Z.\ Phys.\ C {\bf 63} (1994) 377. %\cite{Ellis:tq} \bibitem{KTJET} S.~D.~Ellis and D.~E.~Soper, %``Successive Combination Jet Algorithm For Hadron Collisions,'' Phys.\ Rev.\ D {\bf 48} (1993) 3160 \ifpdf \href{http://arXiv.org/abs/hep-ph/9305266}{[hep-ph/9305266];} \else [hep-ph/9305266]; \fi %%CITATION = HEP-PH 9305266;%% S.~Catani, Y.~L.~Dokshitzer, M.~H.~Seymour and B.~R.~Webber, %``Longitudinally invariant K(t) clustering algorithms for hadron hadron %collisions,'' Nucl.\ Phys.\ B {\bf 406} (1993) 187. %%CITATION = NUPHA,B406,187;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{svref1} R.~Fr\"{u}hwirth {\it et al.}, Nucl.\ Instrum.\ Meth.\ A {\bf 502} (2003) 699. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{svref2} R.~Fr\"{u}hwirth and A.~Strandlie, Comput.\ Phys.\ Commun.\ {\bf 120} (1999) 197. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{jetnet} C.~Peterson, T.~R\"ognvaldsson and L.~L\"onnblad, %``JETNET 3.0: A Versatile artificial neural network package,'' Comput.\ Phys.\ Commun.\ {\bf 81} (1994) 185. %%CITATION = CPHCB,81,185;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% 99/00 low Q2 for F2 points %\cite{h1newf2paper} \bibitem{h1newf2paper1} F.D.~Aaron {\it et al.} [H1 Collaboration], submitted to Eur.\ Phys.\ J.C, %DESY-08-171 \ifpdf \href{http://arXiv.org/abs/arXiv:0904.0929}{arXiv:0904.0929.} \else arXiv:0904.0929. \fi %%CITATION = ARXIV:0904.0929;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{h1newf2paper2} F.D.~Aaron {\it et al.} [H1 Collaboration], submitted to Eur.\ Phys.\ J.C, %DESY-09-005 \ifpdf \href{http://arXiv.org/abs/arXiv:0904.3513}{arXiv:0904.3513.} \else arXiv:0904.3513. \fi %%CITATION = ARXIV:0904.3513;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%% high Q2 results for F2 \bibitem{H19900NCCC} 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. C {\bf 30} (2003) 1 \ifpdf \href{http://arXiv.org/abs/hep-ex/0304003}{[hep-ex/0304003].} \else [hep-ex/0304003]. \fi %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%% fit of which we use a FFNS version \bibitem{Adloff:1999ah} C.~Adloff {\it et al.} [H1 Collaboration], %``Measurement of neutral and charged current cross-sections in positron %proton collisions at large momentum transfer,'' Eur.\ Phys.\ J.\ C {\bf 13} (2000) 609 \ifpdf \href{http://arXiv.org/abs/hep-ex/9908059}{[hep-ex/9908059].} \else [hep-ex/9908059]. \fi %%CITATION = HEP-EX 9908059;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{pdg2006} W.~M.~Yao {\it et al.} [Particle Data Group], %``Review of particle physics,'' J.\ Phys.\ G {\bf 33} (2006) 1. %%CITATION = JPHGB,G33,1;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{Gladilin:1999pj} L.~Gladilin, ``Charm hadron production fractions'', \ifpdf \href{http://arXiv.org/abs/hep-ex/9912064}{hep-ex/9912064.} \else hep-ex/9912064. \fi %%CITATION = HEP-EX 9912064;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{Aktas:2004ka} A.~Aktas {\it et al.} [H1 Collaboration], %``Inclusive production of D+, D0, D/s+ and D*+ mesons in deep inelastic %scattering at HERA,'' Eur.\ Phys.\ J.\ C {\bf 38} (2005) 447 \ifpdf \href{http://arXiv.org/abs/hep-ex/0408149}{[hep-ex/0408149].} \else [hep-ex/0408149]. \fi %%CITATION = HEP-EX 0408149;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{Coffman:1991ud} D.~Coffman {\it et al.} [MARK-III Collaboration], %``Measurement of the inclusive decay properties of charmed mesons,'' Phys.\ Lett.\ B {\bf 263} (1991) 135. %%CITATION = PHLTA,B263,135;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{lepjetmulti} D.~Abbaneo {\it et al.} [LEP/SLD Heavy Flavour Working Group], LEPHF 2001-01\\ \ifpdf (available from \href{http://lepewwg.web.cern.ch/LEPEWWG/heavy/} {http://lepewwg.web.cern.ch/LEPEWWG/heavy/}). \else (available from http://lepewwg.web.cern.ch/LEPEWWG/heavy/). \fi %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{h1frag} F.~D.~Aaron {\it et al.} [H1 Collaboration], %``Study of Charm Fragmentation into D^{*\pm} Mesons in Deep-Inelastic %Scattering at HERA,'' Eur.\ Phys.\ J.\ C {\bf 59} (2009) 1 \ifpdf \href{http://arXiv.org/abs/arXiv:0808.1003}{[arXiv:0808.1003].} \else [arXiv:0808.1003]. \fi %%CITATION = ARXIV:0808.1003;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{rahmatthesis} A.J.\ Rahmat, ``Measurement of the Charm and Beauty Jet Cross Sections in Neutral Current Deep Inelastic Scattering at HERA'', Ph.D. Thesis, University of Liverpool, UK, 2008 \ifpdf (available at \href{http://www-h1.desy.de/publications/theses_list.html} {http://www-h1.desy.de/publications/theses\_list.html}). \else (available at http://www-h1.desy.de/publications/theses\_list.html). \fi %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{Harris:1997zq} B.~W.~Harris and J.~Smith, %``Charm quark and D*+- cross sections in deeply inelastic scattering at %HERA,'' Phys.\ Rev.\ D {\bf 57} (1998) 2806 \ifpdf \href{http://arXiv.org/abs/hep-ph/9706334}{[hep-ph/9706334].} \else [hep-ph/9706334]. \fi %%CITATION = HEP-PH 9706334;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %\bibitem{Laenen:1992zk} \bibitem{massivenlo} E.~Laenen, S.~Riemersma, J.~Smith and W.~L.~van Neerven, %``Complete O (alpha-s) corrections to heavy flavor structure functions in %electroproduction,'' Nucl.\ Phys.\ B {\bf 392} (1993) 162; %CITATION = NUPHA,B392,162;%% %%\bibitem{Laenen:1992xs} E.~Laenen, S.~Riemersma, J.~Smith and W.~L.~van Neerven, %``O(alpha-s) corrections to heavy flavor inclusive distributions in %electroproduction,'' Nucl.\ Phys.\ B {\bf 392} (1993) 229. %%CITATION = NUPHA,B392,229;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{ccfm2} M.~Ciafaloni, %%``Coherence Effects In Initial Jets At Small Q**2 / S,'' Nucl.\ Phys.\ B {\bf 296} (1988) 49; %%%CITATION = NUPHA,B296,49;%% S.~Catani, F.~Fiorani and G.~Marchesini, %%%``QCD Coherence In Initial State Radiation,'' Phys.\ Lett.\ B {\bf 234} (1990) 339; %%%CITATION = PHLTA,B234,339;%% S.~Catani, F.~Fiorani and G.~Marchesini, %%%``Small X Behavior Of Initial State Radiation In Perturbative QCD,'' Nucl.\ Phys.\ B {\bf 336} (1990) 18; %%%CITATION = NUPHA,B336,18;%% G.~Marchesini, %%%``QCD coherence in the structure function and associated distributions at %%%small x,'' Nucl.\ Phys.\ B {\bf 445} (1995) 49 \ifpdf \href{http://arXiv.org/abs/hep-ph/9412327}{[hep-ph/9412327].} \else [hep-ph/9412327]. \fi %%%CITATION = HEP-PH 9412327;%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \bibitem{a0} H.\ Jung, %``Unintegrated PDFs in CCFM'', proceedings of the $12^{\rm th}$ International Workshop on Deep Inelastic Scattering (DIS 2004), Strbske Pleso, Slovakia (2004) 299 \ifpdf \href{http://arXiv.org/abs/hep-ph/0411287}{[hep-ph/0411287].} \else [hep-ph/0411287]. \fi %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \end{thebibliography} \newpage \begin{table}[htb] \tiny \begin{center} \renewcommand{\arraystretch}{1.15} \begin{tabular}{|r|r|c|c|c|c|c|c|c|c|c|} \hline \multicolumn{1}{|c|}{bin} & \multicolumn{1}{c|}{$Q^2$} & $x$ & $y$ & $\rho_l$ & $\rho_c$ & $\rho_b$ & $\chi^2/{\rm n.d.f.}$ & $C_{lc}$ & $C_{lb}$ & $C_{bc}$ \bigstrut[t] \\ \hline $ 1 $ & $ 5.0 $ & $ 0.00020 $ & $ 0.246 $ & $ 1.27 \pm 0.02 $ & $ 1.39 \pm 0.13 $ & $ 1.72 \pm 0.58 $ & $ 19.9 / 24 $ & $ -0.99 $ & $ 0.56 $ & $ -0.61 $ \\ $ 2 $ & $ 8.5 $ & $ 0.00050 $ & $ 0.167 $ & $ 1.15 \pm 0.01 $ & $ 1.27 \pm 0.08 $ & $ 0.58 \pm 0.56 $ & $ 26.4 / 24 $ & $ -0.99 $ & $ 0.57 $ & $ -0.64 $ \\ $ 3 $ & $ 8.5 $ & $ 0.00032 $ & $ 0.262 $ & $ 1.15 \pm 0.01 $ & $ 1.07 \pm 0.07 $ & $ 0.86 \pm 0.22 $ & $ 35.9 / 39 $ & $ -0.99 $ & $ 0.52 $ & $ -0.57 $ \\ $ 4 $ & $ 12.0 $ & $ 0.00130 $ & $ 0.091 $ & $ 1.11 \pm 0.01 $ & $ 1.13 \pm 0.08 $ & $ 0.94 \pm 0.92 $ & $ 25.6 / 24 $ & $ -0.98 $ & $ 0.62 $ & $ -0.73 $ \\ $ 5 $ & $ 12.0 $ & $ 0.00080 $ & $ 0.148 $ & $ 1.14 \pm 0.01 $ & $ 0.99 \pm 0.06 $ & $ 1.29 \pm 0.41 $ & $ 67.3 / 40 $ & $ -0.98 $ & $ 0.59 $ & $ -0.67 $ \\ $ 6 $ & $ 12.0 $ & $ 0.00050 $ & $ 0.236 $ & $ 1.09 \pm 0.01 $ & $ 1.17 \pm 0.06 $ & $ 0.48 \pm 0.24 $ & $ 34.9 / 38 $ & $ -0.99 $ & $ 0.52 $ & $ -0.59 $ \\ $ 7 $ & $ 12.0 $ & $ 0.00032 $ & $ 0.369 $ & $ 1.15 \pm 0.02 $ & $ 1.11 \pm 0.07 $ & $ 0.88 \pm 0.21 $ & $ 64.3 / 39 $ & $ -0.99 $ & $ 0.55 $ & $ -0.60 $ \\ $ 8 $ & $ 20.0 $ & $ 0.00200 $ & $ 0.098 $ & $ 1.12 \pm 0.01 $ & $ 1.10 \pm 0.06 $ & $ 1.77 \pm 0.72 $ & $ 23.3 / 26 $ & $ -0.97 $ & $ 0.62 $ & $ -0.77 $ \\ $ 9 $ & $ 20.0 $ & $ 0.00130 $ & $ 0.151 $ & $ 1.14 \pm 0.01 $ & $ 1.10 \pm 0.05 $ & $ 0.97 \pm 0.30 $ & $ 52.3 / 39 $ & $ -0.98 $ & $ 0.57 $ & $ -0.67 $ \\ $ 10 $ & $ 20.0 $ & $ 0.00080 $ & $ 0.246 $ & $ 1.15 \pm 0.01 $ & $ 1.14 \pm 0.05 $ & $ 0.60 \pm 0.21 $ & $ 36.5 / 38 $ & $ -0.98 $ & $ 0.52 $ & $ -0.61 $ \\ $ 11 $ & $ 20.0 $ & $ 0.00050 $ & $ 0.394 $ & $ 1.25 \pm 0.02 $ & $ 0.98 \pm 0.06 $ & $ 1.00 \pm 0.15 $ & $ 33.2 / 40 $ & $ -0.99 $ & $ 0.56 $ & $ -0.61 $ \\ $ 12 $ & $ 35.0 $ & $ 0.00320 $ & $ 0.108 $ & $ 1.12 \pm 0.02 $ & $ 1.09 \pm 0.07 $ & $ 0.93 \pm 0.72 $ & $ 27.2 / 26 $ & $ -0.96 $ & $ 0.62 $ & $ -0.78 $ \\ $ 13 $ & $ 35.0 $ & $ 0.00200 $ & $ 0.172 $ & $ 1.18 \pm 0.02 $ & $ 0.97 \pm 0.06 $ & $ 1.52 \pm 0.35 $ & $ 43.5 / 39 $ & $ -0.97 $ & $ 0.58 $ & $ -0.70 $ \\ $ 14 $ & $ 35.0 $ & $ 0.00130 $ & $ 0.265 $ & $ 1.19 \pm 0.02 $ & $ 1.08 \pm 0.06 $ & $ 0.67 \pm 0.22 $ & $ 58.2 / 39 $ & $ -0.98 $ & $ 0.54 $ & $ -0.63 $ \\ $ 15 $ & $ 35.0 $ & $ 0.00080 $ & $ 0.431 $ & $ 1.25 \pm 0.02 $ & $ 1.04 \pm 0.06 $ & $ 1.02 \pm 0.17 $ & $ 46.5 / 39 $ & $ -0.99 $ & $ 0.57 $ & $ -0.64 $ \\ $ 16 $ & $ 60.0 $ & $ 0.00500 $ & $ 0.118 $ & $ 1.15 \pm 0.02 $ & $ 0.99 \pm 0.06 $ & $ 2.17 \pm 0.49 $ & $ 50.6 / 40 $ & $ -0.95 $ & $ 0.59 $ & $ -0.76 $ \\ $ 17 $ & $ 60.0 $ & $ 0.00320 $ & $ 0.185 $ & $ 1.16 \pm 0.02 $ & $ 1.10 \pm 0.06 $ & $ 0.66 \pm 0.25 $ & $ 61.6 / 37 $ & $ -0.96 $ & $ 0.52 $ & $ -0.66 $ \\ $ 18 $ & $ 60.0 $ & $ 0.00200 $ & $ 0.295 $ & $ 1.15 \pm 0.02 $ & $ 1.00 \pm 0.06 $ & $ 0.90 \pm 0.20 $ & $ 38.6 / 39 $ & $ -0.97 $ & $ 0.52 $ & $ -0.64 $ \\ $ 19 $ & $ 60.0 $ & $ 0.00130 $ & $ 0.454 $ & $ 1.32 \pm 0.02 $ & $ 0.80 \pm 0.07 $ & $ 1.06 \pm 0.18 $ & $ 76.8 / 38 $ & $ -0.98 $ & $ 0.57 $ & $ -0.66 $ \\ $ 20 $ & $ 120.0 $ & $ 0.01300 $ & $ 0.091 $ & $ 1.13 \pm 0.02 $ & $ 0.91 \pm 0.09 $ & $ 1.37 \pm 0.46 $ & $ 32.4 / 23 $ & $ -0.95 $ & $ 0.61 $ & $ -0.77 $ \\ $ 21 $ & $ 120.0 $ & $ 0.00500 $ & $ 0.236 $ & $ 1.20 \pm 0.02 $ & $ 0.88 \pm 0.06 $ & $ 1.27 \pm 0.20 $ & $ 37.0 / 38 $ & $ -0.96 $ & $ 0.55 $ & $ -0.69 $ \\ $ 22 $ & $ 120.0 $ & $ 0.00200 $ & $ 0.591 $ & $ 1.25 \pm 0.03 $ & $ 0.92 \pm 0.08 $ & $ 0.94 \pm 0.17 $ & $ 42.5 / 37 $ & $ -0.97 $ & $ 0.54 $ & $ -0.65 $ \\ $ 23 $ & $ 200.0 $ & $ 0.01300 $ & $ 0.151 $ & $ 1.14 \pm 0.02 $ & $ 1.03 \pm 0.08 $ & $ 0.65 \pm 0.21 $ & $ 40.7 / 36 $ & $ -0.95 $ & $ 0.53 $ & $ -0.67 $ \\ $ 24 $ & $ 200.0 $ & $ 0.00500 $ & $ 0.394 $ & $ 1.19 \pm 0.03 $ & $ 0.83 \pm 0.08 $ & $ 0.75 \pm 0.16 $ & $ 41.1 / 36 $ & $ -0.97 $ & $ 0.53 $ & $ -0.65 $ \\ $ 25 $ & $ 300.0 $ & $ 0.02000 $ & $ 0.148 $ & $ 1.08 \pm 0.03 $ & $ 0.83 \pm 0.12 $ & $ 1.16 \pm 0.27 $ & $ 32.7 / 33 $ & $ -0.95 $ & $ 0.52 $ & $ -0.67 $ \\ $ 26 $ & $ 300.0 $ & $ 0.00800 $ & $ 0.369 $ & $ 1.12 \pm 0.03 $ & $ 0.99 \pm 0.09 $ & $ 0.54 \pm 0.17 $ & $ 35.8 / 34 $ & $ -0.96 $ & $ 0.51 $ & $ -0.64 $ \\ $ 27 $ & $ 650.0 $ & $ 0.03200 $ & $ 0.200 $ & $ 1.08 \pm 0.03 $ & $ 0.72 \pm 0.15 $ & $ 0.75 \pm 0.29 $ & $ 28.2 / 32 $ & $ -0.96 $ & $ 0.52 $ & $ -0.65 $ \\ $ 28 $ & $ 650.0 $ & $ 0.01300 $ & $ 0.492 $ & $ 1.09 \pm 0.03 $ & $ 0.85 \pm 0.10 $ & $ 0.82 \pm 0.16 $ & $ 37.0 / 35 $ & $ -0.96 $ & $ 0.54 $ & $ -0.66 $ \\ $ 29 $ & $ 2000.0 $ & $ 0.05000 $ & $ 0.394 $ & $ 1.07 \pm 0.04 $ & $ 0.73 \pm 0.20 $ & $ 0.67 \pm 0.37 $ & $ 25.5 / 25 $ & $ -0.96 $ & $ 0.55 $ & $ -0.67 $ \\ \hline \end{tabular} \end{center} \normalsize \caption{The fit parameters $\rho_l$, $\rho_c$ and $\rho_b$, along with their errors, the $\chi^2$ per degree of freedom and the correlation coefficients of the fit parameters for each bin in $Q^2$ and $x$. The parameters are shown for the fine binning scheme used to evaluate the reduced $c$ cross section.} \label{tab:rhotabledeffine} \end{table} \begin{table}[htb] \tiny \begin{center} \renewcommand{\arraystretch}{1.15} \begin{tabular}{|r|r|c|c|c|c|c|c|c|c|c|} \hline \multicolumn{1}{|c|}{bin} & \multicolumn{1}{c|}{$Q^2$} & $x$ & $y$ & $\rho_l$ & $\rho_c$ & $\rho_b$ & $\chi^2/{\rm n.d.f.}$ & $C_{lc}$ & $C_{lb}$ & $C_{bc}$ \bigstrut[t] \\ \hline $ 1 $ & $ 5.0 $ & $ 0.00020 $ & $ 0.246 $ & $ 1.27 \pm 0.02 $ & $ 1.39 \pm 0.13 $ & $ 1.72 \pm 0.58 $ & $ 19.9 / 24 $ & $ -0.99 $ & $ 0.56 $ & $ -0.61 $ \\ $ 2 $ & $ 12.0 $ & $ 0.00032 $ & $ 0.369 $ & $ 1.17 \pm 0.01 $ & $ 1.08 \pm 0.06 $ & $ 1.06 \pm 0.17 $ & $ 49.4 / 40 $ & $ -0.99 $ & $ 0.55 $ & $ -0.60 $ \\ $ 3 $ & $ 12.0 $ & $ 0.00080 $ & $ 0.148 $ & $ 1.12 \pm 0.01 $ & $ 1.12 \pm 0.03 $ & $ 0.75 \pm 0.15 $ & $ 61.0 / 46 $ & $ -0.99 $ & $ 0.54 $ & $ -0.61 $ \\ $ 4 $ & $ 25.0 $ & $ 0.00050 $ & $ 0.492 $ & $ 1.29 \pm 0.02 $ & $ 0.95 \pm 0.05 $ & $ 1.15 \pm 0.13 $ & $ 56.3 / 42 $ & $ -0.99 $ & $ 0.57 $ & $ -0.63 $ \\ $ 5 $ & $ 25.0 $ & $ 0.00130 $ & $ 0.189 $ & $ 1.15 \pm 0.01 $ & $ 1.08 \pm 0.02 $ & $ 0.93 \pm 0.13 $ & $ 67.3 / 45 $ & $ -0.98 $ & $ 0.55 $ & $ -0.65 $ \\ $ 6 $ & $ 60.0 $ & $ 0.00130 $ & $ 0.454 $ & $ 1.20 \pm 0.01 $ & $ 0.91 \pm 0.03 $ & $ 1.00 \pm 0.10 $ & $ 61.1 / 43 $ & $ -0.97 $ & $ 0.54 $ & $ -0.65 $ \\ $ 7 $ & $ 60.0 $ & $ 0.00500 $ & $ 0.118 $ & $ 1.13 \pm 0.01 $ & $ 1.01 \pm 0.04 $ & $ 1.16 \pm 0.20 $ & $ 51.3 / 44 $ & $ -0.96 $ & $ 0.55 $ & $ -0.71 $ \\ $ 8 $ & $ 200.0 $ & $ 0.00500 $ & $ 0.394 $ & $ 1.14 \pm 0.02 $ & $ 0.83 \pm 0.04 $ & $ 0.88 \pm 0.09 $ & $ 35.6 / 40 $ & $ -0.97 $ & $ 0.55 $ & $ -0.67 $ \\ $ 9 $ & $ 200.0 $ & $ 0.01300 $ & $ 0.151 $ & $ 1.19 \pm 0.02 $ & $ 1.06 \pm 0.07 $ & $ 1.01 \pm 0.24 $ & $ 35.3 / 39 $ & $ -0.95 $ & $ 0.53 $ & $ -0.68 $ \\ $ 10 $ & $ 650.0 $ & $ 0.01300 $ & $ 0.492 $ & $ 1.07 \pm 0.03 $ & $ 0.91 \pm 0.09 $ & $ 0.50 \pm 0.13 $ & $ 42.9 / 36 $ & $ -0.97 $ & $ 0.55 $ & $ -0.66 $ \\ $ 11 $ & $ 650.0 $ & $ 0.03200 $ & $ 0.200 $ & $ 1.08 \pm 0.02 $ & $ 0.74 \pm 0.09 $ & $ 1.12 \pm 0.19 $ & $ 48.8 / 37 $ & $ -0.96 $ & $ 0.52 $ & $ -0.65 $ \\ $ 12 $ & $ 2000.0 $ & $ 0.05000 $ & $ 0.394 $ & $ 1.07 \pm 0.04 $ & $ 0.73 \pm 0.20 $ & $ 0.67 \pm 0.37 $ & $ 25.5 / 25 $ & $ -0.96 $ & $ 0.55 $ & $ -0.67 $ \\ \hline \end{tabular} \end{center} \normalsize \caption{The fit parameters $\rho_l$, $\rho_c$ and $\rho_b$, along with their errors, the $\chi^2$ per degree of freedom and the correlation coefficients of the fit parameters for each bin in $Q^2$ and $x$. The parameters are shown for the coarse binning scheme used to evaluate the reduced $b$ cross section.} \label{tab:rhotabledefcoarse} \end{table} \begin{sidewaystable}[h!] \tiny \begin{center} \renewcommand{\arraystretch}{1.15} \begin{tabular}{|c|d|c|c|c|c|d|d|d|d|d|d|d|d|d|d|d|d|d|d|d|d|d|} \hline bin & Q^2 & $x$ & $y$ & $\tilde{\sigma}^{q\bar{q}}$ & $F^{q\bar{q}}_2$ & \multicolumn{1}{r|}{$\delta_{\rm stat}$} & \multicolumn{1}{r|}{$\delta_{\rm sys}$} & \multicolumn{1}{r|}{$\delta_{\rm tot}$} & \multicolumn{1}{r|}{$\delta_{\rm unc}$} & \multicolumn{1}{r|}{$\delta_{\rm res}$} & \multicolumn{1}{r|}{$\delta_{\rm eff CJC}$} & \multicolumn{1}{r|}{$\delta_{\rm eff CST}$} & \multicolumn{1}{r|}{$\delta_{\rm frag C}$} & \multicolumn{1}{r|}{$\delta_{\rm frag B}$} & \multicolumn{1}{r|}{$\delta_{uds}$} & \multicolumn{1}{r|}{$\delta_{\phi^q}$} & \multicolumn{1}{r|}{$\delta_{\rm had E}$} & \multicolumn{1}{r|}{$\delta_{\rm gp}$} & \multicolumn{1}{r|}{$\delta_{\rm F_2}$} \bigstrut[t] \\ & \multicolumn{1}{r|}{(GeV$^2$)} & & & & & \multicolumn{1}{r|}{($\%$)} & \multicolumn{1}{r|}{($\%$)} & \multicolumn{1}{r|}{($\%$)} & \multicolumn{1}{r|}{($\%$)} & \multicolumn{1}{r|}{($\%$)} & \multicolumn{1}{r|}{($\%$)} & \multicolumn{1}{r|}{($\%$)} & \multicolumn{1}{r|}{($\%$)} & \multicolumn{1}{r|}{($\%$)} & \multicolumn{1}{r|}{($\%$)} & \multicolumn{1}{r|}{($\%$)} & \multicolumn{1}{r|}{($\%$)} & \multicolumn{1}{r|}{($\%$)} & \multicolumn{1}{r|}{($\%$)} \\ \hline $c$ 1 & $ 5.0 $ & $ 0.00020 $ & $ 0.246 $ & $ 0.148 $ & $ 0.149 $ & 9.8 & 14.6 & 17.6 & 3.3 & 1.9 & -1.9 & -2.8 & -3.5 & 0.3 & -7.7 & 4.6 & -1.2 & 3.1 & 1.0 \\ $c$ 2 & $ 8.5 $ & $ 0.00050 $ & $ 0.167 $ & $ 0.176 $ & $ 0.176 $ & 6.5 & 13.3 & 14.8 & 2.0 & 1.9 & -1.3 & -2.0 & -2.8 & 0.0 & -10.0 & 1.8 & -1.2 & 0.4 & 1.0 \\ $c$ 3 & $ 8.5 $ & $ 0.00032 $ & $ 0.262 $ & $ 0.186 $ & $ 0.187 $ & 6.4 & 14.1 & 15.5 & 3.8 & 1.1 & -1.0 & -1.5 & -5.8 & 0.0 & -7.2 & 2.1 & -0.2 & 3.6 & 1.0 \\ $c$ 4 & $ 12.0 $ & $ 0.00130 $ & $ 0.091 $ & $ 0.150 $ & $ 0.150 $ & 7.3 & 17.2 & 18.7 & 1.2 & 1.5 & -0.6 & -0.9 & -2.5 & 0.1 & -15.7 & 1.7 & 1.7 & 0.0 & 1.0 \\ $c$ 5 & $ 12.0 $ & $ 0.00080 $ & $ 0.148 $ & $ 0.177 $ & $ 0.177 $ & 5.2 & 15.0 & 15.9 & 1.3 & 1.0 & -0.8 & -1.3 & -2.5 & 0.0 & -13.5 & 1.2 & -0.7 & 0.2 & 1.1 \\ $c$ 6 & $ 12.0 $ & $ 0.00050 $ & $ 0.236 $ & $ 0.240 $ & $ 0.242 $ & 4.9 & 10.1 & 11.2 & 1.3 & 2.0 & -1.0 & -1.4 & -3.0 & 0.0 & -5.6 & 1.7 & -1.7 & 0.6 & 1.0 \\ $c$ 7 & $ 12.0 $ & $ 0.00032 $ & $ 0.369 $ & $ 0.273 $ & $ 0.277 $ & 5.6 & 12.6 & 13.8 & 1.4 & 2.1 & -0.8 & -1.3 & -5.6 & 0.1 & -5.8 & 1.6 & -0.3 & 3.1 & 1.1 \\ $c$ 8 & $ 20.0 $ & $ 0.00200 $ & $ 0.098 $ & $ 0.187 $ & $ 0.188 $ & 4.0 & 12.0 & 12.7 & 1.6 & 1.4 & -0.4 & -0.8 & -2.2 & 0.2 & -9.9 & 1.2 & 2.2 & 0.0 & 1.1 \\ $c$ 9 & $ 20.0 $ & $ 0.00130 $ & $ 0.151 $ & $ 0.219 $ & $ 0.219 $ & 4.6 & 11.0 & 11.9 & 1.2 & 1.2 & -0.9 & -1.4 & -2.2 & 0.0 & -8.7 & 1.6 & 0.7 & 0.2 & 1.1 \\ $c$ 10 & $ 20.0 $ & $ 0.00080 $ & $ 0.246 $ & $ 0.274 $ & $ 0.276 $ & 4.5 & 9.2 & 10.2 & 1.2 & 1.5 & -1.0 & -1.4 & -2.6 & 0.0 & -5.8 & 1.4 & -0.8 & 0.5 & 1.0 \\ $c$ 11 & $ 20.0 $ & $ 0.00050 $ & $ 0.394 $ & $ 0.281 $ & $ 0.287 $ & 4.9 & 12.9 & 13.8 & 1.5 & 1.6 & -0.5 & -0.7 & -5.9 & 0.2 & -6.7 & 1.7 & -1.3 & 3.2 & 1.1 \\ $c$ 12 & $ 35.0 $ & $ 0.00320 $ & $ 0.108 $ & $ 0.200 $ & $ 0.200 $ & 6.9 & 10.7 & 12.7 & 1.9 & 2.4 & -0.6 & -0.8 & -1.8 & 0.3 & -8.8 & 1.2 & 3.0 & 0.1 & 1.1 \\ $c$ 13 & $ 35.0 $ & $ 0.00200 $ & $ 0.172 $ & $ 0.220 $ & $ 0.220 $ & 6.1 & 10.1 & 11.8 & 1.5 & 1.9 & -0.7 & -1.0 & -1.9 & 0.0 & -7.9 & 2.5 & -0.3 & 0.2 & 1.0 \\ $c$ 14 & $ 35.0 $ & $ 0.00130 $ & $ 0.265 $ & $ 0.295 $ & $ 0.297 $ & 5.4 & 8.1 & 9.7 & 1.2 & 2.2 & -0.7 & -1.0 & -2.0 & 0.0 & -5.4 & 1.4 & -0.4 & 0.3 & 1.0 \\ $c$ 15 & $ 35.0 $ & $ 0.00080 $ & $ 0.431 $ & $ 0.349 $ & $ 0.360 $ & 6.1 & 11.1 & 12.7 & 1.3 & 0.6 & -0.5 & -0.8 & -4.6 & 0.1 & -6.4 & 2.0 & -1.0 & 1.8 & 1.1 \\ $c$ 16 & $ 60.0 $ & $ 0.00500 $ & $ 0.118 $ & $ 0.198 $ & $ 0.199 $ & 5.1 & 9.5 & 10.8 & 1.8 & 1.9 & -0.5 & -0.8 & -1.7 & 0.7 & -6.8 & 2.3 & 1.2 & 0.1 & 1.1 \\ $c$ 17 & $ 60.0 $ & $ 0.00320 $ & $ 0.185 $ & $ 0.263 $ & $ 0.264 $ & 5.5 & 6.3 & 8.4 & 1.1 & 1.5 & -0.7 & -1.1 & -1.5 & 0.0 & -3.3 & 2.6 & 1.3 & 0.0 & 1.0 \\ $c$ 18 & $ 60.0 $ & $ 0.00200 $ & $ 0.295 $ & $ 0.335 $ & $ 0.339 $ & 4.3 & 7.7 & 8.8 & 1.1 & 1.8 & -0.5 & -0.8 & -1.5 & 0.1 & -5.8 & 1.7 & -0.3 & 0.1 & 1.0 \\ $c$ 19 & $ 60.0 $ & $ 0.00130 $ & $ 0.454 $ & $ 0.296 $ & $ 0.307 $ & 8.3 & 12.6 & 15.1 & 1.2 & 1.6 & -0.1 & -0.2 & -3.7 & 0.2 & -9.6 & 1.5 & -1.3 & 1.0 & 1.0 \\ $c$ 20 & $ 120.0 $ & $ 0.01300 $ & $ 0.091 $ & $ 0.133 $ & $ 0.133 $ & 9.7 & 10.2 & 14.1 & 2.7 & 2.6 & -0.7 & -1.0 & -2.2 & 0.9 & -7.7 & 3.2 & 1.2 & 0.0 & 1.2 \\ $c$ 21 & $ 120.0 $ & $ 0.00500 $ & $ 0.236 $ & $ 0.218 $ & $ 0.220 $ & 6.8 & 8.8 & 11.1 & 1.6 & 2.5 & -0.3 & -0.5 & -1.2 & 0.1 & -6.6 & 3.1 & 0.7 & 0.0 & 1.1 \\ $c$ 22 & $ 120.0 $ & $ 0.00200 $ & $ 0.591 $ & $ 0.351 $ & $ 0.375 $ & 8.6 & 9.5 & 12.8 & 3.2 & 3.3 & -0.3 & -0.4 & -2.7 & 0.2 & -5.1 & 1.9 & 0.0 & 0.9 & 2.9 \\ $c$ 23 & $ 200.0 $ & $ 0.01300 $ & $ 0.151 $ & $ 0.160 $ & $ 0.160 $ & 7.1 & 9.5 & 11.9 & 2.8 & 2.1 & -1.0 & -1.5 & -1.8 & 0.2 & -6.5 & 3.1 & 1.7 & 0.0 & 2.7 \\ $c$ 24 & $ 200.0 $ & $ 0.00500 $ & $ 0.394 $ & $ 0.237 $ & $ 0.243 $ & 8.8 & 10.2 & 13.5 & 3.4 & 3.3 & -0.1 & -0.2 & -2.0 & 0.0 & -6.2 & 3.7 & 0.3 & 0.5 & 2.9 \\ $c$ 25 & $ 300.0 $ & $ 0.02000 $ & $ 0.148 $ & $ 0.117 $ & $ 0.117 $ & 14.1 & 12.0 & 18.5 & 3.3 & 2.6 & -0.4 & -0.6 & -1.3 & 0.4 & -9.3 & 3.8 & 1.9 & 0.0 & 2.9 \\ $c$ 26 & $ 300.0 $ & $ 0.00800 $ & $ 0.369 $ & $ 0.273 $ & $ 0.278 $ & 9.5 & 8.4 & 12.7 & 3.3 & 3.6 & -0.2 & -0.3 & -1.1 & 0.0 & -3.4 & 4.8 & -0.1 & 0.5 & 2.9 \\ $c$ 27 & $ 650.0 $ & $ 0.03200 $ & $ 0.200 $ & $ 0.084 $ & $ 0.085 $ & 16.7 & 26.0 & 30.9 & 3.8 & 5.2 & -2.2 & -3.3 & -3.4 & 0.2 & -20.8 & 10.7 & 4.5 & 0.0 & 3.4 \\ $c$ 28 & $ 650.0 $ & $ 0.01300 $ & $ 0.492 $ & $ 0.195 $ & $ 0.203 $ & 10.8 & 12.0 & 16.2 & 3.7 & 3.0 & -0.1 & -0.2 & -1.0 & 0.2 & -6.9 & 6.9 & -0.2 & 0.4 & 3.3 \\ $c$ 29 & $ 2000.0 $ & $ 0.05000 $ & $ 0.394 $ & $ 0.059 $ & $ 0.060 $ & 28.3 & 22.9 & 36.4 & 6.7 & 4.0 & -1.7 & -2.6 & -1.2 & 0.1 & -14.1 & 15.6 & 0.0 & 0.2 & 4.3 \\ \hline $b$ 1 & $ 5.0 $ & $ 0.00020 $ & $ 0.246 $ & $ 0.00244 $ & $ 0.00244 $ & 33.8 & 31.3 & 46.1 & 8.1 & -2.4 & -2.8 & -5.6 & -3.7 & -6.6 & 6.0 & -3.7 & 3.1 & 3.1 & 1.0 \\ $b$ 2 & $ 12.0 $ & $ 0.00032 $ & $ 0.369 $ & $ 0.00487 $ & $ 0.00490 $ & 16.0 & 27.4 & 31.8 & 8.4 & -3.4 & -4.1 & -8.2 & -1.1 & -4.5 & 5.1 & 2.2 & -0.8 & 4.1 & 1.1 \\ $b$ 3 & $ 12.0 $ & $ 0.00080 $ & $ 0.148 $ & $ 0.00247 $ & $ 0.00248 $ & 16.7 & 40.2 & 43.5 & 9.6 & -8.1 & -7.4 & -14.8 & -1.6 & -5.3 & 17.8 & 2.2 & 6.6 & 0.6 & 1.1 \\ $b$ 4 & $ 25.0 $ & $ 0.00050 $ & $ 0.492 $ & $ 0.01189 $ & $ 0.01206 $ & 10.1 & 23.0 & 25.1 & 5.0 & -4.0 & -3.8 & -7.5 & 0.2 & -4.5 & 3.5 & 0.8 & 0.3 & 3.1 & 1.3 \\ $b$ 5 & $ 25.0 $ & $ 0.00130 $ & $ 0.189 $ & $ 0.00586 $ & $ 0.00587 $ & 12.1 & 31.9 & 34.1 & 7.8 & -8.7 & -5.6 & -11.2 & -0.8 & -4.9 & 10.5 & 0.9 & 3.5 & 0.2 & 1.0 \\ $b$ 6 & $ 60.0 $ & $ 0.00130 $ & $ 0.454 $ & $ 0.01928 $ & $ 0.01969 $ & 9.5 & 23.1 & 25.0 & 10.9 & -4.6 & -4.0 & -8.1 & -0.1 & -2.9 & 4.1 & 1.4 & 0.7 & 0.8 & 1.0 \\ $b$ 7 & $ 60.0 $ & $ 0.00500 $ & $ 0.118 $ & $ 0.00964 $ & $ 0.00965 $ & 14.6 & 29.1 & 32.6 & 4.3 & -7.0 & -5.2 & -10.4 & -0.5 & -7.9 & 2.3 & -2.8 & 8.6 & 0.0 & 1.1 \\ $b$ 8 & $ 200.0 $ & $ 0.00500 $ & $ 0.394 $ & $ 0.02365 $ & $ 0.02422 $ & 9.4 & 21.2 & 23.2 & 7.9 & -4.6 & -4.2 & -8.5 & 0.2 & -2.7 & 1.1 & 1.0 & 0.5 & 1.9 & 2.9 \\ $b$ 9 & $ 200.0 $ & $ 0.01300 $ & $ 0.151 $ & $ 0.01139 $ & $ 0.01142 $ & 19.5 & 28.4 & 34.4 & 4.6 & -4.4 & -4.5 & -8.9 & 0.2 & -10.3 & 5.7 & 0.0 & 10.4 & 0.0 & 2.7 \\ $b$ 10 & $ 650.0 $ & $ 0.01300 $ & $ 0.492 $ & $ 0.01331 $ & $ 0.01394 $ & 23.6 & 25.4 & 34.7 & 13.4 & -4.6 & -5.0 & -10.1 & -1.1 & -2.2 & 5.8 & 4.5 & 1.5 & 0.8 & 3.3 \\ $b$ 11 & $ 650.0 $ & $ 0.03200 $ & $ 0.200 $ & $ 0.01018 $ & $ 0.01024 $ & 16.0 & 25.5 & 30.1 & 10.8 & -3.3 & -3.3 & -6.5 & -0.2 & -9.0 & 3.6 & 2.1 & 3.1 & 0.0 & 3.4 \\ $b$ 12 & $ 2000.0 $ & $ 0.05000 $ & $ 0.394 $ & $ 0.00499 $ & $ 0.00511 $ & 55.7 & 25.1 & 61.1 & 15.0 & -2.0 & -4.4 & -8.7 & 1.7 & -2.7 & 9.4 & -3.0 & 1.7 & 0.2 & 4.3 \\ \hline \end{tabular} \end{center} \normalsize \caption{ The measured values and relative errors for the reduced DIS cross section ($\tilde{\sigma}^{q\bar{q}}$) and structure function $F_2^{q\bar{q}}$ for charm ($c$) and beauty ($b$) quarks for the combined datasets. The values for $F_2^{q\bar{q}}$ are obtained from the measured reduced cross sections using the NLO QCD fit to correct for the contributions from $F_L^{q\bar{q}}$. The table shows the statistical error ($\delta_{\rm stat}$), the systematic error ($\delta_{\rm sys}$), the total error ($\delta_{\rm tot}$) and the uncorrelated systematic error ($\delta_{\rm unc}$) on the reduced cross section. The next ten columns show the effect of a $+ 1 \sigma$ shift for the correlated systematic error contributions to the reduced cross section from: track impact parameter resolution, CJC track efficiency, CST track efficiency, $c$ fragmentation, $b$ fragmentation, light quark contribution, struck quark angle $\phi^q$, hadronic energy scale, photoproduction background and the DIS event selection. The $-1 \sigma$ errors are taken as the negative of the $+1 \sigma$ errors. The correlated systematic errors are continued in table~\ref{tab:sysmod}. } \label{tab:sig} \end{sidewaystable} \begin{sidewaystable} \tiny \begin{center} \renewcommand{\arraystretch}{1.15} \begin{tabular}{|c|d|d|d|d|d|d|d|d|d|d|d|d|d|d|} \hline bin & \multicolumn{1}{r|}{$\delta_{{\rm mod} Q^2 c}$} & \multicolumn{1}{r|}{$\delta_{{\rm mod} Q^2 b}$} & \multicolumn{1}{r|}{$\delta_{{\rm mod} x c}$} & \multicolumn{1}{r|}{$\delta_{{\rm mod} x b}$} & \multicolumn{1}{r|}{$\delta_{{\rm mod} P_T c}$} & \multicolumn{1}{r|}{$\delta_{{\rm mod} P_T b}$} & \multicolumn{1}{r|}{$\delta_{{\rm mod} \eta c}$} & \multicolumn{1}{r|}{$\delta_{{\rm mod} \eta b}$} & \multicolumn{1}{r|}{$\delta_{{\rm BF} D^+}$} & \multicolumn{1}{r|}{$\delta_{{\rm BF} D^0}$} & \multicolumn{1}{r|}{$\delta_{{\rm Mult} D^+}$} & \multicolumn{1}{r|}{$\delta_{{\rm Mult} D^0}$} & \multicolumn{1}{r|}{$\delta_{{\rm Mult} D_s}$} & \multicolumn{1}{r|}{$\delta_{{\rm Mult} B}$} \\ & \multicolumn{1}{r|}{$(\%)$} & \multicolumn{1}{r|}{$(\%)$} & \multicolumn{1}{r|}{$(\%)$} & \multicolumn{1}{r|}{$(\%)$} & \multicolumn{1}{r|}{$(\%)$} & \multicolumn{1}{r|}{$(\%)$} & \multicolumn{1}{r|}{$(\%)$} & \multicolumn{1}{r|}{$(\%)$} & \multicolumn{1}{r|}{$(\%)$} & \multicolumn{1}{r|}{$(\%)$} & \multicolumn{1}{r|}{$(\%)$} & \multicolumn{1}{r|}{$(\%)$} & \multicolumn{1}{r|}{$(\%)$} & \multicolumn{1}{r|}{$(\%)$} \\ \hline $ c 1 $ & $ -0.1 $ & $ 0.0 $ & $ 2.3 $ & $ -0.4 $ & $ -7.8 $ & $ 0.9 $ & $ 2.8 $ & $ 0.4 $ & $ -0.6 $ & $ 0.2 $ & $ -2.0 $ & $ -1.4 $ & $ -1.7 $ & $ -0.7 $ \\ $ c 2 $ & $ -0.7 $ & $ 0.0 $ & $ 1.4 $ & $ -0.1 $ & $ -5.1 $ & $ 0.3 $ & $ 3.6 $ & $ -0.1 $ & $ -0.6 $ & $ 0.1 $ & $ -1.6 $ & $ -1.6 $ & $ -1.4 $ & $ -0.5 $ \\ $ c 3 $ & $ -0.5 $ & $ 0.0 $ & $ -0.1 $ & $ 0.0 $ & $ -8.3 $ & $ 0.8 $ & $ 0.6 $ & $ 0.4 $ & $ -0.6 $ & $ 0.3 $ & $ -1.8 $ & $ -1.4 $ & $ -1.4 $ & $ -0.9 $ \\ $ c 4 $ & $ -0.3 $ & $ 0.0 $ & $ 1.3 $ & $ 0.0 $ & $ -3.5 $ & $ -0.1 $ & $ 3.5 $ & $ 0.0 $ & $ -0.6 $ & $ 0.0 $ & $ -1.7 $ & $ -1.3 $ & $ -1.7 $ & $ 0.4 $ \\ $ c 5 $ & $ -0.2 $ & $ 0.0 $ & $ 0.8 $ & $ 0.0 $ & $ -4.5 $ & $ 0.4 $ & $ 0.0 $ & $ 0.2 $ & $ -0.6 $ & $ 0.2 $ & $ -1.7 $ & $ -1.3 $ & $ -1.6 $ & $ -1.0 $ \\ $ c 6 $ & $ -0.2 $ & $ 0.0 $ & $ 0.4 $ & $ 0.0 $ & $ -5.9 $ & $ 0.3 $ & $ 1.5 $ & $ 0.1 $ & $ -0.5 $ & $ 0.2 $ & $ -1.7 $ & $ -1.4 $ & $ -1.6 $ & $ -0.3 $ \\ $ c 7 $ & $ -0.2 $ & $ 0.0 $ & $ -0.3 $ & $ 0.0 $ & $ -7.8 $ & $ 1.0 $ & $ 0.3 $ & $ 0.4 $ & $ -0.5 $ & $ 0.3 $ & $ -1.7 $ & $ -1.4 $ & $ -1.7 $ & $ -0.7 $ \\ $ c 8 $ & $ -0.1 $ & $ 0.0 $ & $ 1.5 $ & $ -0.1 $ & $ -2.6 $ & $ 0.3 $ & $ 3.6 $ & $ 0.0 $ & $ -0.3 $ & $ -0.2 $ & $ -1.2 $ & $ -1.6 $ & $ -1.7 $ & $ -0.6 $ \\ $ c 9 $ & $ -0.1 $ & $ 0.0 $ & $ 0.9 $ & $ 0.0 $ & $ -3.7 $ & $ 0.3 $ & $ 3.1 $ & $ 0.1 $ & $ -0.5 $ & $ 0.1 $ & $ -1.6 $ & $ -1.4 $ & $ -1.7 $ & $ -0.5 $ \\ $ c 10 $ & $ -0.1 $ & $ 0.0 $ & $ 0.3 $ & $ 0.0 $ & $ -4.8 $ & $ 0.3 $ & $ 1.6 $ & $ 0.1 $ & $ -0.5 $ & $ 0.2 $ & $ -1.7 $ & $ -1.4 $ & $ -1.8 $ & $ -0.5 $ \\ $ c 11 $ & $ -0.2 $ & $ 0.0 $ & $ -0.8 $ & $ -0.1 $ & $ -7.2 $ & $ 1.3 $ & $ 0.6 $ & $ 0.6 $ & $ -0.4 $ & $ 0.2 $ & $ -1.5 $ & $ -1.5 $ & $ -1.5 $ & $ -1.5 $ \\ $ c 12 $ & $ 0.1 $ & $ 0.0 $ & $ 1.6 $ & $ 0.0 $ & $ -1.6 $ & $ 0.4 $ & $ 0.0 $ & $ -0.1 $ & $ -0.1 $ & $ -0.2 $ & $ -1.0 $ & $ -1.4 $ & $ -1.8 $ & $ -1.1 $ \\ $ c 13 $ & $ 0.0 $ & $ 0.0 $ & $ 0.9 $ & $ 0.0 $ & $ -2.3 $ & $ 0.6 $ & $ 2.9 $ & $ 0.2 $ & $ -0.5 $ & $ 0.2 $ & $ -1.2 $ & $ -1.4 $ & $ -1.4 $ & $ -1.0 $ \\ $ c 14 $ & $ -0.1 $ & $ 0.0 $ & $ 0.3 $ & $ 0.0 $ & $ -3.4 $ & $ 0.5 $ & $ 1.8 $ & $ 0.2 $ & $ -0.5 $ & $ 0.2 $ & $ -1.4 $ & $ -1.5 $ & $ -1.4 $ & $ -0.6 $ \\ $ c 15 $ & $ -0.2 $ & $ 0.0 $ & $ -1.0 $ & $ 0.0 $ & $ -6.3 $ & $ 1.0 $ & $ 0.5 $ & $ 0.5 $ & $ -0.5 $ & $ 0.3 $ & $ -1.8 $ & $ -1.3 $ & $ -1.7 $ & $ -1.1 $ \\ $ c 16 $ & $ 0.0 $ & $ -0.1 $ & $ 1.4 $ & $ -0.3 $ & $ -1.1 $ & $ 0.7 $ & $ 3.5 $ & $ -0.4 $ & $ -0.3 $ & $ 0.0 $ & $ -1.1 $ & $ -1.1 $ & $ -1.8 $ & $ -1.7 $ \\ $ c 17 $ & $ -0.1 $ & $ 0.0 $ & $ 0.6 $ & $ 0.0 $ & $ -1.6 $ & $ 0.3 $ & $ 2.4 $ & $ 0.0 $ & $ -0.4 $ & $ 0.1 $ & $ -1.2 $ & $ -1.4 $ & $ -1.2 $ & $ -0.6 $ \\ $ c 18 $ & $ -0.2 $ & $ 0.0 $ & $ -0.4 $ & $ -0.1 $ & $ -2.4 $ & $ 0.6 $ & $ 1.0 $ & $ 0.1 $ & $ -0.4 $ & $ 0.1 $ & $ -1.4 $ & $ -1.2 $ & $ -1.5 $ & $ -1.1 $ \\ $ c 19 $ & $ -0.6 $ & $ 0.0 $ & $ -1.8 $ & $ 0.0 $ & $ -5.3 $ & $ 1.3 $ & $ 0.0 $ & $ 0.5 $ & $ -0.3 $ & $ 0.2 $ & $ -1.5 $ & $ -1.3 $ & $ -1.6 $ & $ -1.8 $ \\ $ c 20 $ & $ 0.3 $ & $ 0.1 $ & $ 1.5 $ & $ -0.2 $ & $ -0.1 $ & $ 0.5 $ & $ 0.0 $ & $ -0.6 $ & $ -0.2 $ & $ -0.2 $ & $ -1.1 $ & $ -0.9 $ & $ -1.3 $ & $ -2.2 $ \\ $ c 21 $ & $ 0.3 $ & $ 0.0 $ & $ -0.2 $ & $ -0.3 $ & $ -0.9 $ & $ 0.7 $ & $ 0.9 $ & $ -0.2 $ & $ -0.3 $ & $ 0.1 $ & $ -1.2 $ & $ -1.0 $ & $ -1.2 $ & $ -2.4 $ \\ $ c 22 $ & $ -0.9 $ & $ 0.1 $ & $ -1.9 $ & $ 0.2 $ & $ -4.3 $ & $ 0.8 $ & $ 0.2 $ & $ 0.4 $ & $ -0.4 $ & $ 0.3 $ & $ -1.4 $ & $ -1.2 $ & $ -1.1 $ & $ -1.6 $ \\ $ c 23 $ & $ -0.4 $ & $ 0.0 $ & $ 1.6 $ & $ -0.2 $ & $ -0.7 $ & $ 0.2 $ & $ 0.0 $ & $ -0.1 $ & $ -0.4 $ & $ 0.0 $ & $ -1.3 $ & $ -1.1 $ & $ -1.3 $ & $ -1.1 $ \\ $ c 24 $ & $ -0.1 $ & $ 0.1 $ & $ -0.8 $ & $ 0.1 $ & $ -2.3 $ & $ 0.7 $ & $ 0.3 $ & $ 0.3 $ & $ -0.3 $ & $ 0.1 $ & $ -1.2 $ & $ -1.2 $ & $ -1.5 $ & $ -2.2 $ \\ $ c 25 $ & $ -0.1 $ & $ 0.1 $ & $ 0.9 $ & $ -0.2 $ & $ -0.4 $ & $ 0.6 $ & $ 0.0 $ & $ -0.4 $ & $ -0.5 $ & $ 0.1 $ & $ -1.5 $ & $ -1.0 $ & $ -2.0 $ & $ -3.2 $ \\ $ c 26 $ & $ 0.0 $ & $ 0.0 $ & $ -0.5 $ & $ 0.1 $ & $ -1.4 $ & $ 0.3 $ & $ 0.8 $ & $ 0.1 $ & $ -0.3 $ & $ 0.3 $ & $ -1.1 $ & $ -0.8 $ & $ -1.6 $ & $ -1.6 $ \\ $ c 27 $ & $ 0.1 $ & $ 0.1 $ & $ 1.7 $ & $ -0.2 $ & $ 0.1 $ & $ 0.2 $ & $ 3.5 $ & $ -0.1 $ & $ -0.6 $ & $ 0.0 $ & $ -1.3 $ & $ -1.4 $ & $ -2.3 $ & $ -2.5 $ \\ $ c 28 $ & $ 0.1 $ & $ 0.3 $ & $ -0.7 $ & $ 0.3 $ & $ -2.0 $ & $ 0.4 $ & $ -0.1 $ & $ 0.1 $ & $ -0.3 $ & $ 0.2 $ & $ -0.7 $ & $ -0.6 $ & $ -1.1 $ & $ -2.9 $ \\ $ c 29 $ & $ 0.6 $ & $ 0.6 $ & $ 0.3 $ & $ 0.1 $ & $ -0.9 $ & $ 0.5 $ & $ 0.0 $ & $ 0.2 $ & $ -0.5 $ & $ 0.3 $ & $ -1.6 $ & $ -1.0 $ & $ 0.9 $ & $ -2.1 $ \\ \hline $ b 1 $ & $ 0.0 $ & $ -0.4 $ & $ 0.6 $ & $ 14.0 $ & $ -7.2 $ & $ -17.8 $ & $ -2.2 $ & $ 9.0 $ & $ -2.2 $ & $ 0.6 $ & $ -3.0 $ & $ -0.8 $ & $ -1.8 $ & $ 9.2 $ \\ $ b 2 $ & $ -0.2 $ & $ -1.1 $ & $ -0.4 $ & $ -1.3 $ & $ -7.3 $ & $ -17.3 $ & $ -3.8 $ & $ -1.0 $ & $ -1.2 $ & $ 0.9 $ & $ -3.4 $ & $ -0.3 $ & $ -1.1 $ & $ 11.4 $ \\ $ b 3 $ & $ -0.3 $ & $ -2.6 $ & $ 0.8 $ & $ 4.7 $ & $ -13.7 $ & $ -14.5 $ & $ 0.0 $ & $ 8.2 $ & $ -3.9 $ & $ 2.6 $ & $ -10.5 $ & $ -1.3 $ & $ -3.0 $ & $ 12.9 $ \\ $ b 4 $ & $ -0.2 $ & $ -0.5 $ & $ -0.1 $ & $ -1.0 $ & $ -3.8 $ & $ -15.2 $ & $ -1.8 $ & $ -1.4 $ & $ -0.9 $ & $ 0.4 $ & $ -2.3 $ & $ -0.5 $ & $ -0.5 $ & $ 10.7 $ \\ $ b 5 $ & $ -0.6 $ & $ -0.3 $ & $ 0.3 $ & $ 7.1 $ & $ -8.1 $ & $ -11.9 $ & $ 0.0 $ & $ 7.7 $ & $ -3.4 $ & $ 2.0 $ & $ -8.6 $ & $ -1.2 $ & $ -2.8 $ & $ 13.3 $ \\ $ b 6 $ & $ -0.3 $ & $ -0.8 $ & $ -0.4 $ & $ -2.9 $ & $ -2.8 $ & $ -10.2 $ & $ -1.4 $ & $ -0.6 $ & $ -1.4 $ & $ 0.8 $ & $ -3.3 $ & $ -1.1 $ & $ -1.5 $ & $ 12.1 $ \\ $ b 7 $ & $ -0.4 $ & $ -0.5 $ & $ 0.3 $ & $ 7.2 $ & $ -3.0 $ & $ -7.4 $ & $ 0.0 $ & $ 10.1 $ & $ -2.9 $ & $ 1.3 $ & $ -6.8 $ & $ -1.4 $ & $ -2.8 $ & $ 14.8 $ \\ $ b 8 $ & $ -0.4 $ & $ -1.8 $ & $ -0.6 $ & $ -2.1 $ & $ -1.3 $ & $ -8.7 $ & $ 0.0 $ & $ 0.8 $ & $ -0.9 $ & $ 0.3 $ & $ -2.5 $ & $ -1.0 $ & $ -1.2 $ & $ 12.6 $ \\ $ b 9 $ & $ -0.3 $ & $ -2.4 $ & $ 0.4 $ & $ 4.7 $ & $ -1.4 $ & $ -5.1 $ & $ 0.0 $ & $ 10.3 $ & $ -2.0 $ & $ 0.5 $ & $ -4.5 $ & $ -1.3 $ & $ -1.0 $ & $ 15.1 $ \\ $ b 10 $ & $ -0.8 $ & $ -1.5 $ & $ -1.0 $ & $ -2.7 $ & $ -1.0 $ & $ -5.4 $ & $ -1.9 $ & $ 0.3 $ & $ -1.1 $ & $ 0.3 $ & $ -4.0 $ & $ -2.2 $ & $ -2.8 $ & $ 13.0 $ \\ $ b 11 $ & $ -0.2 $ & $ -0.6 $ & $ 0.1 $ & $ 7.9 $ & $ -0.3 $ & $ -1.8 $ & $ 0.0 $ & $ 10.6 $ & $ -0.6 $ & $ -0.5 $ & $ -1.3 $ & $ -0.9 $ & $ -0.1 $ & $ 13.0 $ \\ $ b 12 $ & $ -0.9 $ & $ -0.2 $ & $ 0.2 $ & $ 0.8 $ & $ 0.5 $ & $ -2.2 $ & $ 2.4 $ & $ 2.3 $ & $ -0.5 $ & $ -0.2 $ & $ -1.8 $ & $ -1.2 $ & $ -6.0 $ & $ 11.8 $ \\ \hline \end{tabular} \end{center} \normalsize \caption{ The correlated errors continued from table~\ref{tab:sig}. The first eight errors represent a $+ 1 \sigma$ shift for the correlated systematic error contributions from: reweighting the $Q^2$ distribution, reweighting the $x$ distribution, reweighting the jet transverse momentum distribution $P^{\rm jet}_T$, and reweighting the jet pseudorapidity $\eta^{\rm jet}$ distribution for $c$ and $b$ events. The remaining six columns show the contributions from: $c$ hadron branching fractions and multiplicities, and the $b$ quark multiplicity. Only those uncertainties where there is an effect of $>1\%$ for any $x$--$Q^2$ point are listed; the remaining uncertainties are included in the uncorrelated error.} \label{tab:sysmod} \end{sidewaystable} \clearpage %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[h!] \begin{center} \includegraphics[width=0.75\textwidth]{d09-096f1a} \vspace{0.6cm} \includegraphics[width=0.72\textwidth]{d09-096f1b} \caption{Diagrams of a track in the $X$--$Y$ plane. If the angle $\alpha$ is less than $90^\circ$, $\delta$ is defined as positive otherwise $\delta$ is defined as negative. } \label{fig:alpha} \end{center} \end{figure} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[h!] \begin{center} \includegraphics[width=0.9\textwidth]{d09-096f2} \caption{The distribution of the signed impact parameter $\delta$ of a track to the primary vertex in the $X$--$Y$ plane. Included in the figure is the expectation from the Monte Carlo simulation for light, $c$ and $b$ quarks. The contributions from the various quark flavours in the Monte Carlo are shown after applying the scale factors $\rho_l$, $\rho_c$ and $\rho_b$ obtained from the fit to the complete data sample (see section~\ref{quarkflavourseparation}). The arrows indicate the range over which tracks are selected for analysis. } \label{fig:dca} \end{center} \end{figure} \newpage %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[htb] \includegraphics[width=0.5\textwidth]{d09-096f3a} \includegraphics[width=0.5\textwidth]{d09-096f3b} % \put(40,10){ \includegraphics[width=0.5\textwidth]{H1prelim-07-171.fig4}} \caption{The significance $\delta /\sigma(\delta)$ distribution (a) of the highest absolute significance track ($S_1$) and (b) of the track with the second highest absolute significance ($S_2$). Included in the figure is the expectation from the Monte Carlo simulation for light, $c$ and $b$ quarks. The contributions from the various quark flavours in the Monte Carlo are shown after applying the scale factors $\rho_l$, $\rho_c$ and $\rho_b$ obtained from the fit to the complete data sample.} \label{fig:s1s2} \end{figure} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[htb] \begin{center} \includegraphics[width=0.494\textwidth]{d09-096f4a} \includegraphics[width=0.494\textwidth]{d09-096f4b} \includegraphics[width=0.494\textwidth]{d09-096f4c} \includegraphics[width=0.494\textwidth]{d09-096f4d} \end{center} \caption{Inputs to the NN for events with at least 3 CST tracks: $|S_1|$, $|S_2|$, $|S_3|$ and $S_L$. Included in the figure is the expectation from the Monte Carlo simulation for light, $c$ and $b$ quarks. The contributions from the various quark flavours in the Monte Carlo are shown after applying the scale factors $\rho_l$, $\rho_c$ and $\rho_b$ obtained from the fit to the complete data sample.} \label{fig:nninputsa} \end{figure} \begin{figure}[htb] \begin{center} \includegraphics[width=0.494\textwidth]{d09-096f5a} \includegraphics[width=0.494\textwidth]{d09-096f5b} \includegraphics[width=0.494\textwidth]{d09-096f5c} \includegraphics[width=0.494\textwidth]{d09-096f5d} \end{center} \caption{ Inputs to the NN for events with at least 3 CST tracks: $P_T^{\rm track}$ of the highest and second highest transverse momentum track, $N_{\rm track}$ and $N^{\rm SV}_{\rm track}$. Included in the figure is the expectation from the Monte Carlo simulation for light, $c$ and $b$ quarks. The contributions from the various quark flavours in the Monte Carlo are shown after applying the scale factors $\rho_l$, $\rho_c$ and $\rho_b$ obtained from the fit to the complete data sample.} \label{fig:nninputsb} \end{figure} \begin{figure}[htb] \begin{center} \includegraphics[width=0.494\textwidth]{d09-096f6a} \includegraphics[width=0.494\textwidth]{d09-096f6b} \includegraphics[width=0.494\textwidth]{d09-096f6c} \includegraphics[width=0.494\textwidth]{d09-096f6d} \end{center} \caption{The subtracted distributions for $|S_1|$, $|S_2|$, $|S_3|$ and $S_L$. Each plot shows the difference between the distributions for those events with $S_2>3$ and the corresponding distribution with $S_2<-3$. Included in the figure is the expectation from the Monte Carlo simulation for light, $c$ and $b$ quarks. The contributions from the various quark flavours in the Monte Carlo are shown after applying the scale factors $\rho_l$, $\rho_c$ and $\rho_b$ obtained from the fit to the complete data sample.} \label{fig:nninputsnega} \end{figure} \begin{figure}[htb] \begin{center} \includegraphics[width=0.494\textwidth]{d09-096f7a} \includegraphics[width=0.494\textwidth]{d09-096f7b} \includegraphics[width=0.494\textwidth]{d09-096f7c} \includegraphics[width=0.494\textwidth]{d09-096f7d} \end{center} \caption{The subtracted distributions for $P_T^{\rm track}$ of the highest and second highest transverse momentum track, $N_{\rm track}$ and $N^{\rm SV}_{\rm track}$. Each plot shows the difference between the distributions for those events with $S_2>3$ and the corresponding distribution with $S_2<-3$. Included in the figure is the expectation from the Monte Carlo simulation for light, $c$ and $b$ quarks. The contributions from the various quark flavours in the Monte Carlo are shown after applying the scale factors $\rho_l$, $\rho_c$ and $\rho_b$ obtained from the fit to the complete data sample.} \label{fig:nninputsnegb} \end{figure} \begin{figure}[htb] \begin{center} \includegraphics[width=0.494\textwidth]{d09-096f8a} \includegraphics[width=0.494\textwidth]{d09-096f8b} \includegraphics[width=0.494\textwidth]{d09-096f8c} \includegraphics[width=0.494\textwidth]{d09-096f8d} \end{center} \caption{The subtracted distributions for $P_T^q$, $\eta^q$, $\log Q^2$, $\log x$. Each plot shows the difference between the distributions for those events with $S_2>3$ and the corresponding distribution with $S_2<-3$. Included in the figure is the expectation from the Monte Carlo simulation for light, $c$ and $b$ quarks. The contributions from the various quark flavours in the Monte Carlo are shown after applying the scale factors $\rho_l$, $\rho_c$ and $\rho_b$ obtained from the fit to the complete data sample. Also shown are the variations in $P_T^q$ and $\eta^q$ resulting from uncertainties in the QCD model of heavy quark production, as well as the variations in $\log Q^2$ and $\log x$, resulting from uncertainties on the input heavy quark structure functions.} \label{fig:controlneg} \end{figure} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[htb] \includegraphics[width=0.9\textwidth]{d09-096f9} \caption{The output of the NN. Included in the figure is the expectation from the Monte Carlo simulation for light, $c$ and $b$ quarks. The contributions from the various quark flavours in the Monte Carlo are shown after applying the scale factors $\rho_l$, $\rho_c$ and $\rho_b$ obtained from the fit to the complete data sample.} \label{fig:nnoutput} \end{figure} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[htb] { \includegraphics[width=0.494\textwidth]{d09-096f10a}} { \includegraphics[width=0.494\textwidth]{d09-096f10b}} { \includegraphics[width=0.494\textwidth]{d09-096f10c}} \caption{The subtracted distributions of (a) $S_1$ and (b) $S_2$ and (c) NN output. Included in the figure is the result from the fit to the complete data sample of the Monte Carlo distributions of the various quark flavours to obtain the scale factors $\rho_l$, $\rho_c$ and $\rho_b$.} \label{fig:s1s2nnnegsub} \end{figure} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[htb] \begin{center} \includegraphics[width=0.85\textwidth]{d09-096f11} \caption{The measured reduced cross section $\tilde{\sigma}^{c\bar{c}}$ shown as a function of $x$ for different $Q^2$ values. The inner error bars show the statistical error, the outer error bars represent the statistical and systematic errors added in quadrature. The \mbox{HERA II} measurements are compared with those from \mbox{HERA I}. The measurements obtained from $D^*$ mesons by H1~\cite{H1Dstar} using \mbox{HERA I} data are also shown. The $x$ values of the \mbox{HERA I} data are shifted for visual clarity. } \label{fig:f2cc} \end{center} \end{figure} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[htb] \begin{center} \includegraphics[width=0.75\textwidth]{d09-096f12} \caption{The measured reduced cross section $\tilde{\sigma}^{b\bar{b}}$ shown as a function of $x$ for different $Q^2$ values. The inner error bars show the statistical error, the outer error bars represent the statistical and systematic errors added in quadrature. The \mbox{HERA II} measurements are compared with those from \mbox{HERA I}. The combined H1 data are also shown. The $x$ values of the \mbox{HERA I} and \mbox{HERA II} data are shifted for visual clarity.} \label{fig:f2bb} \end{center} \end{figure} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[htb] \begin{center} \includegraphics[width=0.89\textwidth]{d09-096f13} \caption{ The combined reduced cross section $\tilde{\sigma}^{c\bar{c}}$ shown as a function of $x$ for different $Q^2$ values. The inner error bars show the statistical error, the outer error bars represent the statistical and systematic errors added in quadrature. The predictions of QCD calculations are also shown.} \label{fig:f2ccav} \end{center} \end{figure} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[htb] \begin{center} \includegraphics[width=0.89\textwidth]{d09-096f14} \caption{ The combined reduced cross section $\tilde{\sigma}^{c\bar{c}}$ (as in figure~\ref{fig:f2ccav}). The predictions of the H1PDF 2009 and MSTW08 NLO QCD fits are also shown. For the H1 QCD fit the inner error bands show the experimental uncertainty, the middle error bands include the theoretical model uncertainties of the fit assumptions, and the outer error band represents the total uncertainty including the parameterisation uncertainty.} \label{fig:f2ccav2} \end{center} \end{figure} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[htb] \begin{center} \includegraphics[width=0.79\textwidth]{d09-096f15} \caption{The combined reduced cross section $\tilde{\sigma}^{b\bar{b}}$ shown as a function of $x$ for different $Q^2$ values. The inner error bars show the statistical error, the outer error bars represent the statistical and systematic errors added in quadrature. The predictions of QCD calculations are also shown.} \label{fig:f2bbav} \end{center} \end{figure} \clearpage %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[htb] \begin{center} \includegraphics[width=0.79\textwidth]{d09-096f16} \caption{ The combined reduced cross section $\tilde{\sigma}^{b\bar{b}}$ (as in figure~\ref{fig:f2bbav}). The predictions of the H1PDF 2009 and MSTW08 NLO QCD fits are also shown. For the H1 QCD fit the inner error bands show the experimental uncertainty, the middle error bands include the theoretical model uncertainties of the fit assumptions, and the outer error band represents the total uncertainty including the parameterisation uncertainty. } \label{fig:f2bbav2} \end{center} \end{figure} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[htb] \begin{center} \includegraphics[width=0.85\textwidth]{d09-096f17} \caption{The combined $F_2^{c\bar{c}}$ shown as a function of $Q^2$ for various $x$ values. The inner error bars show the statistical error, the outer error bars represent the statistical and systematic errors added in quadrature. The predictions of QCD calculations are also shown.} \label{fig:f2ccq2} \end{center} \end{figure} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[htb] \begin{center} \includegraphics[width=0.9\textwidth]{d09-096f18} \caption{The combined $F_2^{b\bar{b}}$ shown as a function of $Q^2$ for various $x$ values. The inner error bars show the statistical error, the outer error bars represent the statistical and systematic errors added in quadrature. The predictions of QCD calculations are also shown. Note that some points have been interpolated in $x$ for visual clarity.} \label{fig:f2bbq2} \end{center} \end{figure} %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \begin{figure}[htb] \begin{center} \includegraphics[width=0.88\textwidth]{d09-096f19} \caption{The charm and beauty contributions to the total cross section, $f^{c\bar{c}}$ and $f^{b\bar{b}}$, shown as a function of $Q^2$ for different $x$ values. The inner error bars show the statistical error, the outer error bars represent the statistical and systematic errors added in quadrature. A prediction of NNLO QCD is also shown. The charm data point at $x=0.005$ and $Q^2=300 \ {\rm GeV^2}$ has been interpolated from $x=0.008$.} \label{fig:frac} \end{center} \end{figure} \end{document}