%================================================================ % LaTeX file with preferred layout for H1 paper drafts % use: dvips -D600 file-name %================================================================ %%%%%%%%%%%%% LATEX HEADER %%%%%%%%%%%%% DO NOT DELETE NOT CHANGE ANYTHING IN THE HEADER %%%%%%%%%%%%% UNLESS CLEARLY RECOMMENDED \documentclass[12pt]{article} \usepackage{color} \usepackage{epsfig} \usepackage{amsmath} \usepackage{hhline} \usepackage{amssymb} \usepackage{times} \usepackage{cite} %%%%%%%%%%%% Comment the next two lines to remove the line numbering %\usepackage[]{lineno} %\linenumbers %%%%%%%%%%%% \renewcommand{\topfraction}{1.0} \renewcommand{\bottomfraction}{1.0} \renewcommand{\textfraction}{0.0} \newlength{\dinwidth} \newlength{\dinmargin} \setlength{\dinwidth}{21.0cm} \textheight23.5cm \textwidth16.0cm \setlength{\dinmargin}{\dinwidth} \setlength{\unitlength}{1mm} \addtolength{\dinmargin}{-\textwidth} \setlength{\dinmargin}{0.5\dinmargin} \oddsidemargin -1.0in \addtolength{\oddsidemargin}{\dinmargin} \setlength{\evensidemargin}{\oddsidemargin} \setlength{\marginparwidth}{0.9\dinmargin} \marginparsep 8pt \marginparpush 5pt \topmargin -42pt \headheight 12pt \headsep 30pt \footskip 24pt \parskip 3mm plus 2mm minus 2mm %%%%%%%%%%%%%%%% END OF LATEX HEADER %===============================title page============================= \begin{document} %%%%%%%%%%%%%%%% Pre-defined commands, you can use for the most obvious notations % numbers which show up several times in the text \newcommand{\bhluminosity}{\ensuremath{338.9\pm10.2\,\text{pb}^{-1}}} \newcommand{\totalsystheor}{\ensuremath{1.1\%}} % 1.0677 \newcommand{\totalsysbgr}{\ensuremath{1.2\%}} % 1.1705 \newcommand{\totalsysexp}{\ensuremath{1.4\%}} %1.4265 \newcommand{\totalerrstat}{\ensuremath{0.8\%}} % 0.8470 \newcommand{\totalerrsys}{\ensuremath{2.1\%}} % 2.132 \newcommand{\totalerr}{\ensuremath{2.3\%}} % 2.2941 \newcommand{\triggererror}{\ensuremath{0.2\%}} \newcommand{\isrerror}{\ensuremath{1.0\%}} \newcommand{\tpeerror}{\ensuremath{0.1\%}} \newcommand{\gegmerror}{\ensuremath{0.2\%}} \newcommand{\mcstaterror}{\ensuremath{0.3\%}} \newcommand{\escaleerror}{\ensuremath{0.6\%}} \newcommand{\ereserror}{\ensuremath{1.1\%}} \newcommand{\posreserror}{\ensuremath{0.3\%}} \newcommand{\cipefferror}{\ensuremath{0.2\%}} \newcommand{\cipconverror}{\ensuremath{0.3\%}} \newcommand{\alignmenterror}{\ensuremath{0.4\%}} \newcommand{\zvtxerror}{\ensuremath{0.1\%}} \newcommand{\ctdefferror}{\ensuremath{0.03\%}} \newcommand{\clustererror}{\ensuremath{0.04\%}} \newcommand{\larvetoerror}{\ensuremath{0.05\%}} \newlength{\figwidth} \newcommand{\textchanged}{\color{black}} \setlength{\figwidth}{0.9\textwidth} % Journal macro \def\Journal#1#2#3#4{{#1}~{\bf #2} (#3) #4} \def\NCA{\em Nuovo Cimento} \def\NIM{\em Nucl. Instrum. Methods} \def\NIMA{{\em Nucl. Instrum. Methods} {\bf A}} \def\NPB{{\em Nucl. Phys.} {\bf B}} \def\PLB{{\em Phys. Lett.} {\bf B}} \def\PRL{\em Phys. Rev. Lett.} \def\PRD{{\em Phys. Rev.} {\bf D}} \def\ZPC{{\em Z. Phys.} {\bf C}} \def\EJC{{\em Eur. Phys. J.} {\bf C}} \def\CPC{\em Comp. Phys. Commun.} \begin{titlepage} \noindent \begin{flushleft} {\tt DESY 12-062 \hfill ISSN 0418-9833} \\ {\tt April 2012} \\ \end{flushleft} \noindent %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%% For conference papers %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%% comment the header and fill the right conference %%%%% {\it {\large version of \today}} \\[.3em] \vspace{2cm} \begin{center} \begin{Large} {\bf Determination of the Integrated Luminosity at HERA using Elastic QED Compton Events} \vspace{2cm} H1 Collaboration \end{Large} \end{center} \vspace{2cm} \begin{abstract} A measurement of the integrated luminosity at the $ep$ collider HERA is presented, exploiting the elastic QED Compton process $ep\to e\gamma p$. % The electron and the photon are detected in the backward calorimeter of the H1 experiment. % The integrated luminosity of the data recorded in 2003 to 2007 is determined with a precision of \totalerr{}. % The measurement is found to be compatible with the corresponding result obtained using the Bethe-Heitler process. \end{abstract} \vspace{1.5cm} \begin{center} Accepted by \EJC{} \end{center} \end{titlepage} % THE PAPER DRAFTS HAVE NO AUTHOR-LIST % % FOR PAPER ISSUED FOR THE FINAL READING % COPY THE AUTHOR AND INSTITUTE LISTS % INTO YOUR AREA % % from /h1/iww/ipublications/h1auts.tex % AND UNCOMMENT THE NEXT THREE LINES % %\begin{flushleft} % \input{h1auts} %\end{flushleft} % % Please not that the author list may need re-formatting. %-- H1AUTS Author list by names %-- Status: Tue Apr 3 13:47:10 CEST 2012 Number of authors = 197 \noindent F.D.~Aaron$^{5,48}$, %BUCH-LEFT 01/12 Aaron C.~Alexa$^{5}$, %BUCH-PD 06/06 Alexa V.~Andreev$^{25}$, %LPI -PD 8/88 Andreev S.~Backovic$^{30}$, %PODG-LEFT 08/11 Backovic A.~Baghdasaryan$^{38}$, %YERE-PD 09/03 Baghdasaryana S.~Baghdasaryan$^{38}$, %YERE-ST 02/10 Baghdasaryans E.~Barrelet$^{29}$, %PARI-HON 01/10 Barrelet W.~Bartel$^{11}$, %DESY-HON 02/03 Bartel K.~Begzsuren$^{35}$, %ULBA-PD 04/06 Begzsuren A.~Belousov$^{25}$, %LPI -HON 01/12 Belousov P.~Belov$^{11}$, %DESY-ST 07/10 Belov J.C.~Bizot$^{27}$, %ORSA-LEFT 01/12 Bizot V.~Boudry$^{28}$, %ECPL-PD 1/93 Boudry I.~Bozovic-Jelisavcic$^{2}$, %BEOG-PROF 04/07 Bozovicjelisavcic J.~Bracinik$^{3}$, %BIRM-LEFT 01/12 Bracinik G.~Brandt$^{11}$, %DESY-LEFT 06/11 Brandt M.~Brinkmann$^{11}$, %DESY-PD 03/10 Brinkmann V.~Brisson$^{27}$, %ORSA-HON 01/12 Brisson D.~Britzger$^{11}$, %DESY-ST 10/09 Britzger D.~Bruncko$^{16}$, %KOSI-LEFT 01/12 Bruncko A.~Bunyatyan$^{13,38}$, %MPIH-PD 12/95 Bunyatyan A.~Bylinkin$^{24}$, %ITEP-DIP 09/10 Bylinkin L.~Bystritskaya$^{24}$, %ITEP-PD 05/99 Bystritskaya A.J.~Campbell$^{11}$, %DESY-PD 10/84 Campbella K.B.~Cantun~Avila$^{22}$, %MEX1-ST 04/06 Cantunavila F.~Ceccopieri$^{4}$, %BRUX-PD 10/09 Ceccopieri K.~Cerny$^{32}$, %PRG2-PD 11/08 Cernyk V.~Cerny$^{16,47}$, %KOSI-LEFT 01/12 Cernyv V.~Chekelian$^{26}$, %MPIM-PD 01/90 Chekelian J.G.~Contreras$^{22}$, %MEX1-PROF 03/98 Contreras J.A.~Coughlan$^{6}$, %RAL -LEFT 01/12 Coughlan J.~Cvach$^{31}$, %PRAG-HON 01/10 Cvach J.B.~Dainton$^{18}$, %LIVE-PROF 01/94 Dainton K.~Daum$^{37,43}$, %WUPP-PROF 06/90 Daum B.~Delcourt$^{27}$, %ORSA-LEFT 01/12 Delcourt J.~Delvax$^{4}$, %BRUX-LEFT 10/11 Delvax E.A.~De~Wolf$^{4}$, %ANTW-HON 10/10 Dewolf C.~Diaconu$^{21}$, %MARS-PD 09/96 Diaconu M.~Dobre$^{12,50,51}$, %HAM2-ST 07/09 Dobre V.~Dodonov$^{13}$, %MPIH-PD 04/98 Dodonov A.~Dossanov$^{12,26}$, %MPIM-ST 05/07 Dossanov A.~Dubak$^{30,46}$, %PODG-LEFT 07/11 Dubak G.~Eckerlin$^{11}$, %DESY-PD 8/88 Eckerlin S.~Egli$^{36}$, %PSI -PD 01/10 Egli A.~Eliseev$^{25}$, %LPI -LEFT 12/11 Eliseev E.~Elsen$^{11}$, %DESY-PROF 01/06 Elsen L.~Favart$^{4}$, %BRUX-PROF 10/99 Favart A.~Fedotov$^{24}$, %ITEP-PD 8/88 Fedotov R.~Felst$^{11}$, %DESY-HON 06/02 Felst J.~Feltesse$^{10}$, %SACL-HON 11/06 Feltesse J.~Ferencei$^{16}$, %KOSI-PD 08/88 Ferencei D.-J.~Fischer$^{11}$, %DESY-PD 07/11 Fischer M.~Fleischer$^{11}$, %DESY-PD 07/95 Fleischer A.~Fomenko$^{25}$, %LPI -PD 02/89 Fomenko E.~Gabathuler$^{18}$, %LIVE-HON 01/12 Gabathulere J.~Gayler$^{11}$, %DESY-HON 02/05 Gayler S.~Ghazaryan$^{11}$, %DFLC-PD 09/09 Ghazaryan A.~Glazov$^{11}$, %DESY-PD 01/04 Glazov L.~Goerlich$^{7}$, %CRAC-PD 8/88 Goerlich N.~Gogitidze$^{25}$, %LPI -PD 05/91 Gogitidze M.~Gouzevitch$^{11,44}$, %DESY-PD 10/11 Gouzevitch C.~Grab$^{40}$, %ZUTH-PROF 08/07 Grab A.~Grebenyuk$^{11}$, %DESY-PD 03/12 Grebenyuk T.~Greenshaw$^{18}$, %LIVE-PD 8/88 Greenshaw G.~Grindhammer$^{26}$, %MPIM-HON 09/11 Grindhammer S.~Habib$^{11}$, %DESY-PD 10/09 Habib D.~Haidt$^{11}$, %DESY-HON 03/04 Haidt R.C.W.~Henderson$^{17}$, %LANC-PROF 10/04 Henderson E.~Hennekemper$^{15}$, %HDB2-PD 07/11 Hennekemper H.~Henschel$^{39}$, %ZEUT-LEFT 02/12 Henschel M.~Herbst$^{15}$, %HDB2-PD 11/11 Herbst G.~Herrera$^{23}$, %MEX2-PD 07/98 Herrera M.~Hildebrandt$^{36}$, %PSI -PD 01/10 Hildebrandtm K.H.~Hiller$^{39}$, %ZEUT-PD 01/92 Hiller D.~Hoffmann$^{21}$, %MARS-PD 10/0 Hoffmann R.~Horisberger$^{36}$, %PSI -PD 01/10 Horisberger T.~Hreus$^{4}$, %BRUX-PD 10/08 Hreus F.~Huber$^{14}$, %HDB1-ST 09/09 Huberf M.~Jacquet$^{27}$, %ORSA-PD 09/96 Jacquet X.~Janssen$^{4}$, %ANTW-PD 02/03 Janssenx L.~J\"onsson$^{20}$, %LUND-PROF 01/01 Joensson H.~Jung$^{11,4}$, %DESY-PROF 12/99 Jungh M.~Kapichine$^{9}$, %JINR-PD 3/97 Kapichine I.R.~Kenyon$^{3}$, %BIRM-LEFT 01/12 Kenyon C.~Kiesling$^{26}$, %MPIM-HON 09/11 Kiesling M.~Klein$^{18}$, %LIVE-PROF 12/06 Klein C.~Kleinwort$^{11}$, %DESY-PD 8/88 Kleinwort T.~Kluge$^{18}$, %LIVE-PD 05/04 Kluge R.~Kogler$^{12}$, %HAM2-PD 12/10 Kogler P.~Kostka$^{39}$, %ZEUT-PD 10/87 Kostka M.~Kr\"{a}mer$^{11}$, %DESY-PD 10/09 Kraemer J.~Kretzschmar$^{18}$, %LIVE-PD 03/08 Kretzschmar K.~Kr\"uger$^{15}$, %HDB2-PD 01/04 Kruegerk M.P.J.~Landon$^{19}$, %QMWC-PD 10/83 Landon W.~Lange$^{39}$, %ZEUT-PD 8/88 Lange G.~La\v{s}tovi\v{c}ka-Medin$^{30}$, %PODG-LEFT 01/12 Lastovickamedin P.~Laycock$^{18}$, %LIVE-PD 11/03 Laycock A.~Lebedev$^{25}$, %LPI -HON 01/12 Lebedev V.~Lendermann$^{15}$, %HDB2-LEFT 01/12 Lendermann S.~Levonian$^{11}$, %DESY-PD 08/88 Levonian K.~Lipka$^{11,50}$, %DESY-PD 06/01 Lipka B.~List$^{11}$, %DESY-PD 11/99 Listb J.~List$^{11}$, %DFLC-PD 10/00 Listj B.~Lobodzinski$^{11}$, %DESY-PD 01/06 Lobodzinski R.~Lopez-Fernandez$^{23}$, %MEX2-PD 03/2 Lopezfernandez V.~Lubimov$^{24}$, %ITEP-HON 01/12 Lubimov E.~Malinovski$^{25}$, %LPI -PD 01/89 Malinovskie H.-U.~Martyn$^{1}$, %AAC1-HON 01/12 Martyn S.J.~Maxfield$^{18}$, %LIVE-PD 8/88 Maxfield A.~Mehta$^{18}$, %LIVE-PROF 06/00 Mehta A.B.~Meyer$^{11}$, %DESY-PD 10/97 Meyeran H.~Meyer$^{37}$, %WUPP-HON 01/12 Meyerhi J.~Meyer$^{11}$, %DESY-PROF 10/84 Meyerj S.~Mikocki$^{7}$, %CRAC-PROF 10/10 Mikocki I.~Milcewicz-Mika$^{7}$, %CRAC-PD 12/11 Milcewicz F.~Moreau$^{28}$, %ECPL-PD 01/90 Moreau A.~Morozov$^{9}$, %JINR-PD 06/99 Morozova J.V.~Morris$^{6}$, %RAL -HON 01/10 Morris K.~M\"uller$^{41}$, %ZUER-PD 11/94 Muellerk Th.~Naumann$^{39}$, %ZEUT-PROF 01/05 Naumannt P.R.~Newman$^{3}$, %BIRM-PROF 10/09 Newman C.~Niebuhr$^{11}$, %DESY-PD 3/93 Niebuhr D.~Nikitin$^{9}$, %JINR-PD 06/08 Nikitin G.~Nowak$^{7}$, %CRAC-PROF 02/12 Nowakg K.~Nowak$^{12}$, %HAM2-PD 03/10 Nowakk J.E.~Olsson$^{11}$, %DESY-HON 09/10 Olsson D.~Ozerov$^{11}$, %DESY-PD 09/08 Ozerov P.~Pahl$^{11}$, %DESY-ST 10/08 Pahl V.~Palichik$^{9}$, %JINR-PD 01/04 Palichik I.~Panagoulias$^{l,}$$^{11,42}$, %DESY-LEFT 05/11 Panagoulias M.~Pandurovic$^{2}$, %BEOG-PD 03/11 Pandurovic Th.~Papadopoulou$^{l,}$$^{11,42}$, %DESY-LEFT 05/11 Papadopoulou C.~Pascaud$^{27}$, %ORSA-HON 09/05 Pascaud G.D.~Patel$^{18}$, %LIVE-PD 8/88 Patel E.~Perez$^{10,45}$, %SACL-PD 10/07 Perez A.~Petrukhin$^{11}$, %DESY-PD 10/09 Petrukhin I.~Picuric$^{30}$, %PODG-PROF 12/07 Picuric H.~Pirumov$^{14}$, %HDB1-ST 09/09 Pirumov D.~Pitzl$^{11}$, %DESY-PD 8/88 Pitzl R.~Pla\v{c}akyt\.{e}$^{11}$, %DESY-PD 10/06 Placakyte B.~Pokorny$^{32}$, %PRG2-ST 10/09 Pokorny R.~Polifka$^{32,52}$, %PRG2-PD 05/11 Polifka B.~Povh$^{13}$, %MPIH-PROF 09/95 Povh V.~Radescu$^{11}$, %DESY-PD 10/06 Radescu N.~Raicevic$^{30}$, %PODG-PD 06/02 Raicevic T.~Ravdandorj$^{35}$, %ULBA-PROF 09/00 Ravdandorj P.~Reimer$^{31}$, %PRAG-PD 07/90 Reimer E.~Rizvi$^{19}$, %QMWC-PROF 03/04 Rizvi P.~Robmann$^{41}$, %ZUER-PD 8/88 Robmann R.~Roosen$^{4}$, %BRUX-HON 03/12 Roosen A.~Rostovtsev$^{24}$, %ITEP-PROF 02/11 Rostovtsev M.~Rotaru$^{5}$, %BUCH-PD 06/11 Rotaru J.E.~Ruiz~Tabasco$^{22}$, %MEX1-PD 05/10 Ruiztabascojuliaelis S.~Rusakov$^{25}$, %LPI -HON 01/12 Rusakov D.~\v S\'alek$^{32}$, %PRG2-PD 10/10 Salek D.P.C.~Sankey$^{6}$, %RAL -PD 01/90 Sankey M.~Sauter$^{14}$, %HDB1-PD 10/09 Sauter E.~Sauvan$^{21,53}$, %MARS-PD 11/1 Sauvan S.~Schmitt$^{11}$, %DESY-PD 09/99 Schmittst L.~Schoeffel$^{10}$, %SACL-PROF 10/10 Schoeffel A.~Sch\"oning$^{14}$, %HDB1-PROF 01/09 Schoening H.-C.~Schultz-Coulon$^{15}$, %HDB2-PD 01/04 Schultzcoulon F.~Sefkow$^{11}$, %DFLC-PD 09/99 Sefkow L.N.~Shtarkov$^{25}$, %LPI -LEFT 01/12 Shtarkov S.~Shushkevich$^{11}$, %DESY-PD 08/11 Shushkevich T.~Sloan$^{17}$, %LANC-LEFT 01/12 Sloan Y.~Soloviev$^{11,25}$, %LPI -PD 08/89 Soloviev P.~Sopicki$^{7}$, %CRAC-ST 09/07 Sopicki D.~South$^{11}$, %DESY-PD 05/03 South V.~Spaskov$^{9}$, %JINR-PD 12/97 Spaskov A.~Specka$^{28}$, %ECPL-PROF 09/05 Specka Z.~Staykova$^{4}$, %ANTW-PD 10/10 Staykova M.~Steder$^{11}$, %DESY-PD 09/08 Steder B.~Stella$^{33}$, %ROME-HON 11/10 Stella G.~Stoicea$^{5}$, %BUCH-PD 02/08 Stoicea U.~Straumann$^{41}$, %ZUER-PROF 09/95 Straumann T.~Sykora$^{4,32}$, %ANTW-PD 01/06 Sykora P.D.~Thompson$^{3}$, %BIRM-PD 08/99 Thompsonp T.H.~Tran$^{27}$, %ORSA-LEFT 10/11 Tran D.~Traynor$^{19}$, %QMWC-PD 10/01 Traynor P.~Tru\"ol$^{41}$, %ZUER-HON 10/06 Truoel I.~Tsakov$^{34}$, %SOFI-PROF 03/12 Tsakov B.~Tseepeldorj$^{35,49}$, %ULBA-PD 06/06 Tseepeldorj J.~Turnau$^{7}$, %CRAC-PROF 09/95 Turnau A.~Valk\'arov\'a$^{32}$, %PRG2-PD 8/88 Valkarova C.~Vall\'ee$^{21}$, %MARS-PD 06/86 Vallee P.~Van~Mechelen$^{4}$, %ANTW-PROF 01/06 Vanmechelen Y.~Vazdik$^{25}$, %LPI -PD 8/88 Vazdik D.~Wegener$^{8}$, %DORT-PROF 01/06 Wegener E.~W\"unsch$^{11}$, %DESY-PD 8/88 Wuensch J.~\v{Z}\'a\v{c}ek$^{32}$, %PRG2-PD 8/88 Zacek J.~Z\'ale\v{s}\'ak$^{31}$, %PRAG-PD 01/05 Zalesak Z.~Zhang$^{27}$, %ORSA-PD 10/92 Zhang A.~Zhokin$^{24}$, %ITEP-LEFT 06/11 Zhokine R.~\v{Z}leb\v{c}\'{i}k$^{32}$, %PRG2-ST 03/12 Zlebcik H.~Zohrabyan$^{38}$, %YERE-PD 11/02 Zohrabyan and F.~Zomer$^{27}$ %ORSA-PROF 09/06 Zomer %\newpage %-- H1 Institutes \bigskip{\it \noindent $ ^{1}$ I. Physikalisches Institut der RWTH, Aachen, Germany \\ $ ^{2}$ Vinca Institute of Nuclear Sciences, University of Belgrade, 1100 Belgrade, Serbia \\ $ ^{3}$ School of Physics and Astronomy, University of Birmingham, Birmingham, UK$^{ b}$ \\ $ ^{4}$ Inter-University Institute for High Energies ULB-VUB, Brussels and Universiteit Antwerpen, Antwerpen, Belgium$^{ c}$ \\ $ ^{5}$ National Institute for Physics and Nuclear Engineering (NIPNE) , Bucharest, Romania$^{ m}$ \\ $ ^{6}$ STFC, Rutherford Appleton Laboratory, Didcot, Oxfordshire, UK$^{ b}$ \\ $ ^{7}$ Institute for Nuclear Physics, Cracow, Poland$^{ d}$ \\ $ ^{8}$ Institut f\"ur Physik, TU Dortmund, Dortmund, Germany$^{ a}$ \\ $ ^{9}$ Joint Institute for Nuclear Research, Dubna, Russia \\ $ ^{10}$ CEA, DSM/Irfu, CE-Saclay, Gif-sur-Yvette, France \\ $ ^{11}$ DESY, Hamburg, Germany \\ $ ^{12}$ Institut f\"ur Experimentalphysik, Universit\"at Hamburg, Hamburg, Germany$^{ a}$ \\ $ ^{13}$ Max-Planck-Institut f\"ur Kernphysik, Heidelberg, Germany \\ $ ^{14}$ Physikalisches Institut, Universit\"at Heidelberg, Heidelberg, Germany$^{ a}$ \\ $ ^{15}$ Kirchhoff-Institut f\"ur Physik, Universit\"at Heidelberg, Heidelberg, Germany$^{ a}$ \\ $ ^{16}$ Institute of Experimental Physics, Slovak Academy of Sciences, Ko\v{s}ice, Slovak Republic$^{ f}$ \\ $ ^{17}$ Department of Physics, University of Lancaster, Lancaster, UK$^{ b}$ \\ $ ^{18}$ Department of Physics, University of Liverpool, Liverpool, UK$^{ b}$ \\ $ ^{19}$ School of Physics and Astronomy, Queen Mary, University of London, London, UK$^{ b}$ \\ $ ^{20}$ Physics Department, University of Lund, Lund, Sweden$^{ g}$ \\ $ ^{21}$ CPPM, Aix-Marseille Univ, CNRS/IN2P3, 13288 Marseille, France \\ $ ^{22}$ Departamento de Fisica Aplicada, CINVESTAV, M\'erida, Yucat\'an, M\'exico$^{ j}$ \\ $ ^{23}$ Departamento de Fisica, CINVESTAV IPN, M\'exico City, M\'exico$^{ j}$ \\ $ ^{24}$ Institute for Theoretical and Experimental Physics, Moscow, Russia$^{ k}$ \\ $ ^{25}$ Lebedev Physical Institute, Moscow, Russia \\ $ ^{26}$ Max-Planck-Institut f\"ur Physik, M\"unchen, Germany \\ $ ^{27}$ LAL, Universit\'e Paris-Sud, CNRS/IN2P3, Orsay, France \\ $ ^{28}$ LLR, Ecole Polytechnique, CNRS/IN2P3, Palaiseau, France \\ $ ^{29}$ LPNHE, Universit\'e Pierre et Marie Curie Paris 6, Universit\'e Denis Diderot Paris 7, CNRS/IN2P3, Paris, France \\ $ ^{30}$ Faculty of Science, University of Montenegro, Podgorica, Montenegro$^{ n}$ \\ $ ^{31}$ Institute of Physics, Academy of Sciences of the Czech Republic, Praha, Czech Republic$^{ h}$ \\ $ ^{32}$ Faculty of Mathematics and Physics, Charles University, Praha, Czech Republic$^{ h}$ \\ $ ^{33}$ Dipartimento di Fisica Universit\`a di Roma Tre and INFN Roma~3, Roma, Italy \\ $ ^{34}$ Institute for Nuclear Research and Nuclear Energy, Sofia, Bulgaria$^{ e}$ \\ $ ^{35}$ Institute of Physics and Technology of the Mongolian Academy of Sciences, Ulaanbaatar, Mongolia \\ $ ^{36}$ Paul Scherrer Institut, Villigen, Switzerland \\ $ ^{37}$ Fachbereich C, Universit\"at Wuppertal, Wuppertal, Germany \\ $ ^{38}$ Yerevan Physics Institute, Yerevan, Armenia \\ $ ^{39}$ DESY, Zeuthen, Germany \\ $ ^{40}$ Institut f\"ur Teilchenphysik, ETH, Z\"urich, Switzerland$^{ i}$ \\ $ ^{41}$ Physik-Institut der Universit\"at Z\"urich, Z\"urich, Switzerland$^{ i}$ \\ \bigskip \noindent $ ^{42}$ Also at Physics Department, National Technical University, Zografou Campus, GR-15773 Athens, Greece \\ $ ^{43}$ Also at Rechenzentrum, Universit\"at Wuppertal, Wuppertal, Germany \\ $ ^{44}$ Also at IPNL, Universit\'e Claude Bernard Lyon 1, CNRS/IN2P3, Villeurbanne, France \\ $ ^{45}$ Also at CERN, Geneva, Switzerland \\ $ ^{46}$ Also at Max-Planck-Institut f\"ur Physik, M\"unchen, Germany \\ $ ^{47}$ Also at Comenius University, Bratislava, Slovak Republic \\ $ ^{48}$ Also at Faculty of Physics, University of Bucharest, Bucharest, Romania \\ $ ^{49}$ Also at Ulaanbaatar University, Ulaanbaatar, Mongolia \\ $ ^{50}$ Supported by the Initiative and Networking Fund of the Helmholtz Association (HGF) under the contract VH-NG-401. \\ $ ^{51}$ Absent on leave from NIPNE-HH, Bucharest, Romania \\ $ ^{52}$ Also at Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S 1A7 \\ $ ^{53}$ Also at LAPP, Universit\'e de Savoie, CNRS/IN2P3, Annecy-le-Vieux, France \\ \bigskip \noindent $ ^a$ Supported by the Bundesministerium f\"ur Bildung und Forschung, FRG, under contract numbers 05H09GUF, 05H09VHC, 05H09VHF, 05H16PEA \\ $ ^b$ Supported by the UK Science and Technology Facilities Council, and formerly by the UK Particle Physics and Astronomy Research Council \\ $ ^c$ Supported by FNRS-FWO-Vlaanderen, IISN-IIKW and IWT and by Interuniversity Attraction Poles Programme, Belgian Science Policy \\ $ ^d$ Partially Supported by Polish Ministry of Science and Higher Education, grant DPN/N168/DESY/2009 \\ $ ^e$ Supported by the Deutsche Forschungsgemeinschaft \\ $ ^f$ Supported by VEGA SR grant no. 2/7062/ 27 \\ $ ^g$ Supported by the Swedish Natural Science Research Council \\ $ ^h$ Supported by the Ministry of Education of the Czech Republic under the projects LC527, INGO-LA09042 and MSM0021620859 \\ $ ^i$ Supported by the Swiss National Science Foundation \\ $ ^j$ Supported by CONACYT, M\'exico, grant 48778-F \\ $ ^k$ Russian Foundation for Basic Research (RFBR), grant no 1329.2008.2 and Rosatom \\ $ ^l$ This project is co-funded by the European Social Fund (75\%) and National Resources (25\%) - (EPEAEK II) - PYTHAGORAS II \\ $ ^m$ Supported by the Romanian National Authority for Scientific Research under the contract PN 09370101 \\ $ ^n$ Partially Supported by Ministry of Science of Montenegro, no. 05-1/3-3352 \\ } \newpage %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \section{Introduction} For particle collider experiments, the precise knowledge of the luminosity is essential for any type of cross section measurement. % %For a particle collider with beam particles $p_1$ and $p_2$, the The instantaneous luminosity is defined as \begin{equation} L=\frac{fnN_{1}N_{2}}{A}, \end{equation} where $f$ is the revolution frequency for the two colliding particles $p_1$ and $p_2$, %assumed to be identical for the two particle species, $n$ is the number of colliding bunches per revolution, and $N_1$ ($N_2)$ is the number of particles of type $p_1$ ($p_2$) per bunch. The effective cross section of the beams is $A$. % The time-integrated luminosity ${\cal L}$ relates the cross section $\sigma_{p_1p_2\to X}$ of the reaction $p_1p_2\to X$ to the number of events $N_{p_1p_2\to X}$ expected in the time interval $T$ by \begin{equation} \label{eqn:intlumi} {\cal L}=\int_T Ldt = \frac{N_{p_1p_2\to X}}{\sigma_{p_1p_2\to X}}. \end{equation} Since it is difficult to monitor all beam parameters with a per cent level precision, in particular those defining the effective beam cross section $A$, the integrated luminosity is often determined by counting the number of observed events for a specific reaction $p_1p_2\to X$ with a well-known cross section. At HERA, the colliding beams are protons and electrons\footnote{ In this paper the term ``electron'' is used generically to refer to both electrons and positrons.}. % For the data taking period studied in this analysis, the proton beam energy is $E^0_p=920\,\mathrm{GeV}$ and the electron beam energy is $E^0_e=27.6\,\text{GeV}$. % The reaction used to determine the integrated luminosity is the production of a radiative photon in elastic $ep$ scattering, $ep\to e\gamma p$. % Depending on the phase space considered, this process is referred to as Bethe-Heitler (BH) scattering or QED Compton (QEDC) scattering. % In the BH process\cite{Bethe:1934za}, both the electron and the photon are emitted almost collinearly to the incident electron. % The corresponding cross section is very large, ${\cal O}(100\,\text{mb})$. Dedicated small angle detectors are used to record BH events. % In contrast, for QEDC scattering\cite{qedctheory}, the particles have a sizable transverse momentum with respect to the incident electron and can be detected in the main detector.%, while The momentum transfer squared at the proton vertex, $t$, is generally small. % At very small momentum transfer $\vert t\vert \ll 1\,\text{GeV}^2$, elastic scattering dominates. % At $\vert t\vert \gtrsim 1\,\text{GeV}^2$, inelastic processes are relevant and the reaction is sensitive to the proton structure. % In addition, there are quasi-elastic contributions to the cross section, where the outgoing proton forms an excited state, like $\Delta^{+}$ or $N^{\star}$, which then decays to a low mass hadronic system. % Within the phase space considered in this analysis, the elastic QEDC cross section is ${\cal O}(50\,\text{pb})$. %The events are recorded using the %H1 main detector, also used for other physics analyses. At HERA, the integrated luminosity is usually measured in the BH process, using dedicated detectors located at small angles. % The advantage of this process is its very large cross section, thus negligible statistical uncertainties are achieved for small amounts of integrated luminosity. % However, there are various sources of possibly large systematic uncertainty. % For example, there may be inevitable acceptance limitations for the small angle detectors, caused by elements of the beam transport system which separates the BH photons and electrons from the circulating proton and electron beams. % The acceptances of the photon and electron detectors may exhibit complex spatial structures and can vary in time as well.%, not only in geometry but also in time variations. % Another complication originates from synchrotron radiation emitted by the electron beam as it passes the focusing magnets surrounding the interaction region. % Furthermore $ep$ collisions can happen also outside the nominal interaction region. These contribute to the BH measurement of the integrated luminosity, but must be corrected for when analysing cross sections with the H1 main detector, which has a more limited acceptance as a function of the collision point position. In this paper, a determination of the integrated luminosity is presented, based on the elastic QEDC process, which is measured in the H1 main detector. % This method is insensitive to details of the beam optics. However, the smallness of the cross section leads to limited statistical precision, thus time-dependencies can not be resolved with high resolution by counting elastic QEDC events alone. % Auxiliary measurements of other reactions with higher cross section may be used to monitor time-dependencies. %, where only the relative event yields at different points in time are of interest. Comparisons of BH and QEDC measurements at HERA have been performed previously \cite{Ahmed:1995cf}. % The inelastic QEDC process also has been measured at large $\vert t \vert$ \cite{Aktas:2004ek}. The data available for the elastic QEDC analysis described in the following were recorded with the H1 detector in the years 2003 to 2007. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \section{H1 Detector} In the following, only those components of the H1 detector are briefly introduced which are essential for the present analysis. % A detailed description of the whole detector in its original configuration can be found elsewhere \cite{H1detector}. % Components which were part of later upgrades are referred to here separately. % %%%A detailed description of the H1 detector can be found elsewhere %%%\cite{H1detector}. %%%% %%%Only the components essential for the present analysis are described %%%here. % The origin of the H1 coordinate system is the nominal $ep$ interaction point. % The direction of the proton beam defines the positive $z$--axis (forward direction). % Transverse momenta are measured in the $xy$ plane. % Polar~($\theta$) and~azimuthal~($\varphi$) angles are measured with respect to this reference system. % The pseudo-rapidity is defined as $\eta= -\ln{\tan (\theta/2)}$. % A schematic view of the H1 detector with signals from an elastic QEDC candidate event is shown in figure \ref{fig:event}. In the backward region $-4.0 < \eta < -1.4$, a lead-scintillating fibre calorimeter\cite{spacal} (SpaCal) is used for the identification and measurement of both the scattered electron and the scattered photon. % The energy resolution for electromagnetic showers is $\sigma(E)/E \simeq 7.1\%/\sqrt{E/{\rm GeV}} \oplus 1\%$ \cite{spacaltest}. % The electromagnetic section of the SpaCal is read out in cells of size $4\times 4\,\text{cm}$ in the $xy$ plane, where the Moli\`{e}re radius is $2.5\,\text{cm}$. % The $xy$ position of a shower is reconstructed as a weighted mean of the cell centres, the weights taken proportional to the logarithm of the cell energies \cite{Glazov:2010zza}. % After applying $xy$ dependent corrections, the position resolution is about $3.5\,\text{mm}$ for an electromagnetic shower in the energy range relevant to this analysis. The liquid argon (LAr) calorimeter covers the range $-1.5 < \eta < 3.4$. % Its electromagnetic (hadronic) section is equipped with absorbers made of lead (steel) plates. % An energy resolution of $\sigma(E)/E \simeq 11\%/\sqrt{E/{\rm GeV}}$ for electromagnetic showers and of $\sigma(E)/E \simeq 50\%/\sqrt{E/{\rm GeV}}$ for hadronic showers is obtained from test beam measurements.%%%~\cite{reflar}. The central region of the detector is equipped with a set of tracking detectors (CTD). % There are the two concentric central jet chambers (CJC), interleaved by a $z$ chamber, and the central silicon tracker (CST)~\cite{Pitzl:2000wz}. The CTD measures the momenta of charged particles in the angular range $20^{\circ} < \theta < 160^{\circ}$. % The central inner proportional chambers (CIP)~\cite{Becker:2007ms} are located between the inner CJC and the CST. % The five CIP chambers have a radial spacing of $9\,\text{mm}$, where the innermost layer is located at a radius of $15.7\,\text{cm}$. % In $\varphi$ there is a $16$-fold segmentation, whereas in $z$ the segments have variable size, ranging from $1.8\,\text{cm}$ in the innermost layer to $2.3\,\text{cm}$ in the outermost layer. %119+112+105+99+92=527 pads per phi sector % The CIP has an angular acceptance in the range $10^{\circ}<\theta<170^{\circ}$. % In the backward region, the tracking is complemented by the backward proportional chamber (BPC), located directly in front of the SpaCal. The calorimeters and tracking detectors are located inside a large superconducting solenoid, providing a uniform field of $1.16\,\text{T}$ strength. % The return yoke of the solenoid is instrumented and serves as a muon detector. % Upstream and downstream of the interaction point there are systems of scintillators (VETO), providing time-of-flight information. % Timing signals from the VETO and the SpaCal were used during data taking to reject particles originating from non-$ep$ interactions of the proton beam in the HERA tunnel. % The luminosity system for measuring the Bethe-Heitler process consists of an electron tagger located at $z=-5.4\,\text{m}$ and a photon calorimeter located at $z = -103 \ {\rm m}$. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \section{Signal and Background processes} Monte Carlo event generators (MC) are used to predict event yields of signal and background processes. % A GEANT3 \cite{geant} simulation of the H1 detector is performed for each generated event, where also the relevant time-dependencies such as changes to the detector setup and varying beam conditions are taken into account. % Electromagnetic showers are simulated using a shower library\cite{Glazov:2010zza}. % After detector simulation, the events are passed through the same reconstruction algorithms as were used for the data. The QEDC signal is simulated using the COMPTON22 event generator \cite{compton22}. % This generator produces elastic, quasi-elastic and inelastic events. % The elastic QEDC events are taken as signal, since their cross section only depends on QED theory and on the proton elastic form factors, thus having small uncertainties. {\textchanged Details are discussed in Section \ref{text:theoryerror}.} % The quasi-elastic events are treated as background and suppressed in the analysis, because their cross section depends on less precisely known parameters such as probabilities to produce excited nucleons. % Similarly, the inelastic events are treated as background, because their cross section depends, for example, on {\textchanged parameterisations} of the proton structure functions at very low momentum transfer. % In COMPTON22, the fragmentation of quasi-elastic events is modelled using the SOPHIA package \cite{Mucke:1999yb}, whereas for inelastic events string fragmentation as implemented in PYTHIA \cite{Sjostrand:2000wi} is used. % For the elastic QEDC signal, final state radiation from the electron has been included in the COMPTON22 event generator using the relevant PYTHIA routines. An important source of background is electron-positron pair production, $ep\to ep\,e^{-}e^{+}$, simulated using the GRAPE event generator \cite{Abe:2000cv}. % Other background events originate mainly from various diffractive processes, namely deeply virtual Compton scattering (DVCS), diffractive vector meson (VM) production and non-resonant diffraction. % DVCS is modelled using the MILOU event generator \cite{milou}. % Diffractive VM production is simulated using the DIFFVM event generator \cite{diffvm}, where the production of $\rho^0$, $\omega$, $\phi$, $J/\psi$, $\psi^{\prime}$ and $\Upsilon$ mesons is considered. % For $\rho^0$ production, DIFFVM is modified such that decays to $\pi^0\gamma$ and $\eta^0\gamma$ are included. % Non-resonant diffraction is simulated using the RAPGAP event generator \cite{rapgap}. % Background from non-diffractive deep-inelastic scattering (DIS) is simulated using DJANGO \cite{Schuler:1991yg} and is found to be negligible. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \section{Event selection} During data taking, events with electrons in the SpaCal were recorded if certain trigger criteria were fulfilled. % The trigger required the presence of at least one localised energy deposit in the SpaCal calorimeter with energy $E>6\,\text{GeV}$. In addition, there were veto conditions to reject beam related background not originating from $ep$ collisions. Data periods where the SpaCal trigger was not fully efficient % or where one of the components relevant for this analysis was not % operational, are removed from the analysis, corresponding to about $6\%$ of the total H1 data sample recorded in 2003 to 2007. The elastic QEDC events are selected offline by requiring two clusters in the electromagnetic section of the SpaCal. % A summary of the selection criteria is given in table \ref{tab:cuts}. % % The transverse sizes of the SpaCal clusters are restricted to $R_{\log{}}<6\,\text{cm}$, where $R_{\log{}}$ is calculated from the SpaCal cell centres using logarithmic energy weighting \cite{Glazov:2010zza}. The cluster energies are required to be larger than $2.2\,\text{GeV}$. % In the range $30\le R<72\,\text{cm}$ of the radial distance from the beam, $R$, exactly two such clusters are required, whereas for $20\le R< 30\,\text{cm}$ no cluster is allowed. % The restriction in the number of clusters removes background and events with hard radiative photons. % The restriction $R\ge 30\,\text{cm}$ on the two clusters ensures that the particles are within the CIP acceptance. % Electron trajectories are reconstructed using the SpaCal cluster position together with position information from the CIP chambers. % Hits in the CIP chambers are considered if they match the SpaCal cluster in $\varphi$. % Adjacent hits are merged to CIP clusters in $\varphi$ and $z$. % This merging is done separately for each layer. % A straight line fit of the CIP clusters and the SpaCal cluster in the $rz$ plane is performed, where outliers are rejected. % The coordinate $r$ is the radial distance from the $z$ axis, measured along the azimuthal direction given by the SpaCal cluster. % After outlier rejection, there are up to five accepted CIP clusters, corresponding to the five CIP layers. % At least two accepted CIP clusters are required. % Next, the centre-of-gravity of the CIP clusters in the $rz$ plane is calculated. % Finally, the CIP centre-of-gravity in $rz$ , together with the SpaCal energy and the SpaCal position are used to reconstruct a helix trajectory in three dimensions, pointing back to the origin of the interaction. % In the determination of this helix, the beam spot and beam tilt are also used\footnote{The beam spot is defined as the average $x$ and $y$ position of interactions which take place at $z=0$. The beam tilt is a slope correction for interactions at $z\ne 0$. These parameters were monitored regularly in short time intervals using the CTD.}. % The direction of bending in the magnetic field is chosen assuming that the particle charge is equal to the charge of the beam lepton. % The algorithm finally returns the origin of the interaction (CIP vertex) and the momentum vector at the CIP vertex. The electron and the photon are then identified, making use of the helix fit results. If there is no CIP vertex, the event is rejected. % If there is only one SpaCal cluster linked to a CIP vertex, that cluster is taken as electron while the other cluster is taken as photon. % The photon momentum vector is calculated from the photon cluster energy and a straight line trajectory pointing from the electron CIP vertex to the photon cluster position. % If both SpaCal clusters are linked to CIP vertices, it is assumed that the photon has converted into an electron-positron pair while passing the material in front of the CIP detector. % For that case, two hypotheses are checked in the reconstruction. % First, one of the clusters is taken as the electron, and the photon momentum is calculated using the energy and position of the other cluster as described above. % The difference in azimuth between the electron and photon candidate momenta, $\Delta\varphi_1$, is determined. % Next, the particle hypotheses are interchanged and the corresponding difference in azimuth, $\Delta\varphi_2$, is calculated. % The particle assignment is done according to the hypothesis yielding an azimuthal opening angle closer to $180^{\circ}$. % In the simulation, the mis-identification probability is $0.3\%$ ($16\%$) if one (both) SpaCal clusters are linked to CIP vertices. \begin{table}[t] \begin{center} \begin{tabular}{l} \hline %Electromagnetic SpaCal clusters with transverse sizes $R_{\log{}}<6\,\text{cm} \\ Electromagnetic SpaCal clusters with transverse sizes $R_{\log{}}<6\,\text{cm}$ and $E>2.2\,\text{GeV}$ \\ Exactly two clusters with radial distance to the beam axis $30\le R<72\,\text{cm}$ \\ No additional cluster with $20\le R<30\,\text{cm}$ with energy $E>2.2\,\text{GeV}$ \\ At least one of the two clusters with CIP vertex \\ Electron and photon identification \\ Energies $\min(E_e,E_\gamma)>7\,\mathrm{GeV}$ and $\max(E_e,E_\gamma)>10\,\mathrm{GeV}$ \\ Vertex position $\vert z_{\text{vtx}}\vert<35\,\text{cm}$ \\ Polar angles $155.9^{\circ}<\theta_e,\theta_\gamma<169.5^{\circ}$ \\ Difference in azimuth $170^{\circ}<\Delta\varphi<190^{\circ}$ \\ Transverse momentum balance $\vert \vec{P}_T^{\text{miss}}\vert <0.3\,\text{GeV}$ \\ Energy in LAr at $\theta<10^{\circ}$ is below $0.5\,\text{GeV}$ \\ No additional CTD track \\ \hline \end{tabular} \caption{\label{tab:cuts} Summary of selection criteria.} \end{center} \end{table} Once the electron and photon are identified, the $z$ position of the electron CIP vertex, $z_{\text{vtx}}$, is verified. % Only events with $\vert z_{\text{vtx}}\vert<35\,\text{cm}$ are selected in the analysis. % To further suppress background contributions the following cuts are applied: % energy of the most (least) energetic particle greater than $10$~($7$)~$\text{GeV}$, % polar angles $\theta_e,\theta_\gamma$ within $155.9^{\circ}$ and $169.5^{\circ}$, % difference in azimuth between $170^{\circ}$ and $190^{\circ}$, % modulus of the transverse component of the missing momentum smaller than $0.3\,\text{GeV}$. The missing momentum vector is calculated from the reconstructed electron and photon four-momenta. The inelastic background sources are further suppressed by using conditions on additional activity in the detector. % Events are rejected if the energy in the forward part of the LAr calorimeter, with polar angle $\theta<10^{\circ}$, exceeds $0.5\,\text{GeV}$, or if there are %good CTD tracks %\cite{leewest} which can not be attributed to either the electron or the photon. % The total number of elastic QEDC candidate events is $14277$, after correcting for trigger efficiency effects. The efficiency of the reconstruction is determined for the following phase space defined for the generated elastic QEDC signal: the polar angular range of the generated electron and photon is restricted to $155.9^{\circ}<\theta<169.5^{\circ}$, the maximum fraction of incident electron energy carried away by initial state radiation is $0.35$, and the momentum transfer squared at the proton vertex is limited to $\vert t\vert<0.09\,\text{GeV}^2$. For this generated phase space, the cross section is $\sigma_{\text{gen}}=54.8\,\text{pb}$. % A total fraction of $f_{\text{vtx}}=2.5\%$ of the generated events have a simulated vertex position along the $z$ coordinate, $z_{gen}$, outside $\pm 35\,\text{cm}$ around the nominal vertex position, in order to be able to describe the observed longitudinal vertex distribution in the H1 detector. % Within the region $\vert z_{gen}\vert<35\,\text{cm}$ the reconstruction efficiency is found to be $\epsilon_{\text{rec}}=64.7\%$. % The sources of {\textchanged inefficiency} are investigated in the following. % Losses of $12.4\%$ originate from the cut in $\vert\vec{P}_T^{\text{miss}}\vert$ due to the limited experimental resolution. % The requirement $300$ ($z_{\text{gen}}<0$) and scattered at polar angles near $155.9^{\circ}$ ($169.5^{\circ}$) are outside the allowed range in $R$. % The track and calorimeter veto conditions contribute to the inefficiency by $3.4\%$, dominated by the restriction in forward LAr energy and the veto on additional SpaCal clusters. % The LAr condition is not fully efficient due to electronic noise and overlap with non $ep$ background. % The inefficiency due to a third SpaCal cluster originates from events with hard final state radiation. % Losses due to the other selection criteria are small ($1.8\%$ total). % Finally, $10.5\%$ are rejected by more than one selection criterion, where combinations involving the $\vert\vec{P}_T^{\text{miss}}\vert$ condition dominate. \begin{table}[t] \begin{center} \begin{tabular}{l|c|c} \hline & no $\vert\vec{P}_T^{\text{miss}}\vert$ cut & $\vert\vec{P}_T^{\text{miss}}\vert<0.3\,\text{GeV}$ \\ \hline quasi-elastic QEDC & $6.84\%$ & $2.93\%$ \\ inelastic QEDC & $7.02\%$ & $1.51\%$ \\ elastic DVCS & $2.10\%$ & $1.26\%$ \\ quasi-elastic DVCS & $0.55\%$ & $0.16\%$ \\ $ep \to ep\, e^{+}e^{-}$ & $1.15\%$ & $1.31\%$ \\ diffractive DIS & $2.78\%$ & $0.53\%$ \\ non-diffractive DIS & $0.02\%$ & $0.01\%$ \\ diffractive $\rho^0$ & $2.05\%$ & $0.15\%$ \\ diffractive $\omega$ & $0.43\%$ & $0.03\%$ \\ diffractive $\phi$ & $0.29\%$ & $0.02\%$ \\ diffractive $J/\psi$ & $0.20\%$ & $0.05\%$ \\ diffractive $\psi^{\prime}$ & $0.17\%$ & $0.08\%$ \\ diffractive $\Upsilon$ & $0.02\%$ & $0.01\%$ \\ \hline \end{tabular} \caption{\label{tab:bgr} Background processes contributing to the elastic QEDC selection. The background fractions are given for both a selection with the $\vert\vec{P}_T^{\text{miss}}\vert$ cut released and for the final selection. } \end{center} \end{table} A detailed breakdown of the different background sources contributing to the elastic QEDC selection as defined in table \ref{tab:cuts} is given in table~\ref{tab:bgr} without and with applying the $\vert\vec{P}_T^{\text{miss}}\vert$ cut. %Although this cut lowers the selection efficiency for elastic QEDC events by $12.4\%$ More than half of the background originates from inelastic and quasi-elastic QEDC processes. The $\vert\vec{P}_T^{\text{miss}}\vert$ cut significantly reduces the background by about a factor of three. For measuring the integrated luminosity ${\cal L}^{\text{QEDC}}$ of collisions originating from the region $\vert z_{gen}\vert<35\,\text{cm}$, the following relation is used \begin{equation} {\cal L}^{\text{QEDC}} = (1-f_{\text{vtx}})\frac{N_{\text{event}}(1-f_{\text{bgr}})}{\sigma_{\text{vis}}}, \end{equation} where $N_{\text{event}}$ is the number of QEDC candidate events observed in the detector, $f_{\text{bgr}}$ is the background fraction predicted by the MC simulation and $\sigma_{\text{vis}}=36.4\,\text{pb}$ is the visible QEDC cross section. % for the phase space of the present analysis. % for reconstructing elastic QEDC events. The main contribution to $\sigma_{\text{vis}}$ originates from genuine QEDC production in the phase space region of this analysis and with an interaction vertex within the accepted region %is from $(1-f_{\text{vtx}})\sigma_{\text{gen}}\times\epsilon_{\text{rec}}=34.6\,\text{pb}$. %but there are Additional contributions from %elastic QEDC reactions events outside the defined phase space or with an interaction vertex beyond the defined limits amount to $1.8 \,\text{pb}$. %outside the vertex region defined above. Distributions of variables used in the selection procedure are shown in figure \ref{fig:selection}. Within uncertainties, the data are described by the prediction. % Note that the prediction is normalised to the integrated luminosity ${\cal L}^{\text{QEDC}}$ as determined in the present analysis. % The $\vert \vec{P}_T^{\text{miss}}\vert$ distribution, figure \ref{fig:selection}e, is of particular interest as it shows a clear separation of background and signal. % The analysis cut of $0.3\,\text{GeV}$ is a compromise between inevitable systematic uncertainties due to the limited detector resolution, dominating at small $\vert \vec{P}_T^{\text{miss}}\vert$, and background contributions, increasing at large $\vert \vec{P}_T^{\text{miss}}\vert$. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \section{Systematic uncertainties} Systematic uncertainties on the elastic QEDC measurement may be categorised as follows: experimental uncertainties, background uncertainties and QEDC theory uncertainties. % The experimental uncertainties originate from two sources, trigger and reconstruction efficiencies. % A summary of the systematic uncertainties is given in table \ref{tab:syst}. % The individual contributions are discussed below. % % \begin{table}[t] \begin{center} \begin{tabular}{lr} \hline\hline \multicolumn{2}{c}{\bf Experimental uncertainties} \\ \hline trigger inefficiency & \triggererror{} \\ SpaCal energy scale & \escaleerror{} \\ SpaCal energy resolution & \ereserror{} \\ SpaCal position resolution & \posreserror{} \\ CIP efficiency & \cipefferror{} \\ conversion probability & \cipconverror{} \\ alignment & \alignmenterror{} \\ $z$-vertex distribution & \zvtxerror{} \\ SpaCal cluster finder & \clustererror{} \\ CTD efficiency & \ctdefferror{} \\ LAr energy veto & \larvetoerror{} \\ \hline %& \begin{boldmath}\totalsysexp\end{boldmath} \\ & \totalsysexp \\ \hline\hline \multicolumn{2}{c}{\bf Background uncertainties} \\ \hline non-elastic QEDC & $1.1\%$ \\ % 120311 elastic DVCS & $0.3\%$ \\ quasi-elastic DVCS & $0.1\%$ \\ diffractive VM production & $0.1\%$ \\ non-resonant diffractive DIS & $0.2\%$ \\ $ep\to ep\,e^{+}e^{-}$ & $0.1\%$ \\ \hline % & \begin{boldmath}\totalsysbgr\end{boldmath}\\ & \totalsysbgr\\ \hline\hline \multicolumn{2}{c}{\bf QEDC theory uncertainties} \\ \hline higher order corrections & \isrerror{} \\ {\textchanged proton form factor (TPE parametrisation)} & {\textchanged \tpeerror{}} \\ {\textchanged proton form factor (experimental)} & {\textchanged \gegmerror{}} \\ {\textchanged size of generated signal sample} & {\textchanged \mcstaterror{}} \\ \hline % & \begin{boldmath}\totalsystheor\end{boldmath} \\ & \totalsystheor \\ \hline\hline \end{tabular} \caption{\label{tab:syst} Systematic uncertainties on the determination of the integrated luminosity using elastic QEDC events. The different error sources are grouped into three categories: experimental, background related and theory related uncertainties. The error sources are described in detail in the text.} \end{center} \end{table} % Additional time-dependent uncertainties are present in cases where the integrated luminosity determined in the present analysis is applied to subsets of the H1 data. \subsection{Trigger uncertainties} The main trigger condition is based on calorimetric information in the SpaCal. % It has an efficiency of more than $95\%$ for clusters with energies $E>6\,\text{GeV}$, rising to above $99.8\%$ for energies $E>10\,\text{GeV}$. These efficiencies are verified using independent trigger conditions for a selection of DIS events with the scattered electron in the SpaCal. % Both the electron and the photon from the elastic QEDC reaction % were efficient to fire the trigger create clear signals above the trigger condition thresholds, hence the trigger inefficiency on the QEDC selection is negligible. % However, for certain time periods there were small regions opposite in $\varphi$ with reduced trigger efficiency. % This leads to an uncertainty of $0.02\%$. % The other trigger conditions are related to timing signals from the VETO system and from the SpaCal calorimeter, designed to veto non $ep$ background. % These trigger conditions in conjunction with the varying HERA beam conditions cause inefficiencies of typically $1\%$ for data taken up to the year 2005 and of typically $0.2\%$ afterwards. These inefficiencies mainly originate from particles from beam related background recorded within %overlay events of genuine $ep$ collision events. % and beam related background. % The veto inefficiencies are corrected for by applying time-dependent weights to the data events. % The corresponding systematic uncertainty is \triggererror{}. \subsection{Reconstruction uncertainties} Reconstruction uncertainties originate mainly from the understanding of the SpaCal response to electrons and photons.%, when reconstructing their energy. % The primary SpaCal energy calibration is done using electrons in DIS events \cite{Collaboration:2010ry}. % It corrects for time-dependent or spatial non-uniformities of the calorimeter response. % However, the response of the SpaCal is slightly different for electrons and photons, mostly due to the presence of dead material in front of the calorimeter and due to final state radiation of the electrons. % Furthermore, it is found that the primary calibration can be improved by correcting the energy response as a function of the transverse cluster size, $R_{\log{}}$. % For the QEDC analysis, multiplicative calibration factors are applied to the SpaCal cluster energies for electrons, non-converted photons and converted photons, respectively. % These factors are taken to have a linear dependence on $R_{\log{}}$. % The corresponding parameters are determined by applying the double-angle calibration method to both the photon and the electron differentially in $R_{\log{}}$. % Distributions of $P_T/P_{T,DA}$, where $P_T$ is the measured transverse momentum and $P_{T,DA}$ is the predicted transverse momentum, are investigated for the selection of QEDC events with the cut on the momentum balance, $\vert \vec{P}_T^{\text{miss}}\vert <0.3\,\text{GeV}$, relaxed. % The predicted transverse momentum is given by \cite{Bentvelsen:1992fu,Hoeger:1991wj} \begin{equation} P_{T,DA}=2E^0_e\left(\frac{1-\cos\theta_e}{\sin\theta_e}+\frac{1-\cos\theta_\gamma}{\sin\theta_\gamma}\right)^{-1}, \end{equation} where $\theta_e$ and $\theta_\gamma$ are the polar angles of the electron and the photon, respectively. % For each $P_T/P_{T,DA}$ distribution, the position of the maximum is determined in a fit. The calibration parameters are finally determined from a linear fit as a function of $R_{\log{}}$. % In data (MC), the energy response of the calorimeter to non-converted photons is found to be on average $3.5\%$ ($2.2\%$) higher than the response to electrons. % For the event sample of converted photons, the energy response to photon candidates is $0.4\%$ lower than the response to electrons, both for data and MC. % In order to determine the systematic uncertainty, the energy scale of electrons and photons is varied separately by $0.5\%$ each. % The size of this variation covers possible systematic effects originating from the calibration procedure described above. % Furthermore, a simultaneous variation of the electron and photon energy scale by another $0.5\%$ is considered as systematic uncertainty, originating from the primary energy calibration\cite{Collaboration:2010ry}. % In total, all energy scale variations together contribute to the uncertainty on ${\cal L}^{\text{QEDC}}$ by \escaleerror{}. In addition to the calibration factors, the energy resolution is determined from fits to the $P_T/P_{T,DA}$ distributions, however without dividing the sample into bins of $R_{\log{}}$. % Figure \ref{fig:calibration} shows the distributions of $P_{T,e}/P_{T,DA}$ and $P_{T,\gamma}/P_{T,DA}$ for electrons and photons, respectively. % The distributions are peaked at $1$, as expected after calibration. % Near the peak, the distribution is more asymmetric towards smaller energies for electrons as compared to photons. % This is attributed to final state radiation and energy losses in the material located in front of the calorimeter. % The original MC simulation (not shown in figure \ref{fig:calibration}) has deficits to describe both widths and tails towards lower transverse momenta. The effect exists for both electrons and photons. % It is corrected for by applying an extra smearing of the reconstructed energies in the MC simulation. % An energy offset $\Delta E=(\delta-\tau)\,E^0_e/2$ is subtracted, where $\delta$ is a random number drawn from an exponential distribution, i.e.~the probability density to find $\delta>0$ is given by $f(\delta)=1/\tau\exp[-\delta/\tau]$. % By construction, $\Delta E\sim (\delta-\tau)$ has an expectation value of $\langle\Delta E\rangle=0$. % This has the desired effect that the peak position of $P_T/P_{T,DA}$ is affected only little by the smearing. % Two independent parameters $\tau_e$ and $\tau_\gamma$ are foreseen to describe the expectation values of the exponential probability distributions for electrons and photons, respectively. % It turns out that both $\tau_e$ and $\tau_\gamma$ take the same central value, $\tau_e=\tau_\gamma=0.010$. % Figure \ref{fig:calibration}a shows that the distribution of $P_{T,e}/P_{T,DA}$ is described by the smeared simulation within a variation of the smearing parameter $\tau_e=0.010\pm0.005$. % Similarly, for photons, $P_{T,\gamma}/P_{T,DA}$ (Figure \ref{fig:calibration}b) is described within the variation $\tau_\gamma=0.010\pm0.005$. % The $\tau_e$ and $\tau_\gamma$ variations together cause an uncertainty on ${\cal L}^{\text{QEDC}}$ of \ereserror{}. The SpaCal cluster position resolution in MC is worse than in data. % This effect has been identified using the difference in azimuth of the electron and the photon, $\Delta\varphi$, shown in figure \ref{fig:posres}. % The original simulation has a deficit at $\Delta\varphi$ near $180^{\circ}$, corresponding to a resolution worse than in data. % In order to improve the description of data by MC, the reconstructed cluster positions in MC, $\vec{x}_{\text{rec}}$, are modified such that they are closer to the extrapolated SpaCal positions of the corresponding generated particles $\vec{x}_{\text{gen}}$. % For the analysis, the positions $\vec{x}_{\text{MC}}=(1-f)\vec{x}_{\text{rec}}+f\vec{x}_{\text{gen}}$ are used, where the constant is found to be $f=0.14$. The uncertainty of $f$ is taken as $0.05$, resulting in an uncertainty on ${\cal L}^{\text{QEDC}}$ of \posreserror{}. % The data are described by the prediction within that systematic variation, as demonstrated in figure \ref{fig:posres}. The CIP efficiency for electrons is determined in data and in the simulation using DIS events. It is found to be near $99\%$ in data and near $99.5\%$ for MC, varying as a function of the radial distance of the SpaCal cluster from the beam, $R$. % A correction as a function of $R$ is made by dropping a fraction of CIP vertices in the simulation. % The CIP spatial resolution is adjusted using elastic $\rho^0$ production events in DIS. The CIP vertex, reconstructed from the scattered electron, is compared to the CTD vertex, reconstructed from the $\pi^{+}\pi^{-}$ pair. % The conversion rate of photons in front of the CIP is underestimated in the simulation. % In data, the conversion probability is around $32\%$, whereas the MC predicts $23\%$. % This is corrected by mimicking conversion effects for a fraction of MC events with non-converted photons. % For these events, extra CIP clusters near the expected position are added, and the energy response is scaled to match the expectation for converted photons. % For estimating systematic effects, the three CIP related corrections described above are switched off one by one, and the resulting differences on ${\cal L}^{\text{QEDC}}$ are taken as uncertainties. For the CIP efficiency correction the uncertainty is \cipefferror{}. The CIP resolution tuning has negligible effect and the conversion probability leads to an uncertainty of \cipconverror{}. The alignment of the SpaCal and CIP detectors is done using DIS events. % The interaction vertex is reconstructed using tracks in the CTD, originating from the hadronic final state. % Using hits in the BPC detector and the energy measured in the SpaCal, the electron trajectory is extrapolated to the CIP and SpaCal detectors. % The alignment uncertainties are dominated by the uncertainty on the SpaCal $z$ position. % Systematic effects are estimated by varying the SpaCal $z$ position by $\pm 5\,\text{mm}$, resulting in an uncertainty of \alignmenterror{} on ${\cal L}^{\text{QEDC}}$. The longitudinal vertex distribution is dominated by a Gaussian near $z=0$ with a width of approximately $10\,\text{cm}$, as can be seen in figure \ref{fig:selection}f. % The longitudinal proton beam profile also exhibits prominent satellite peaks of similar width, leading to collisions in the H1 detector near $\pm70\,\text{cm}$. % In addition, there is an excess of collisions near $40\,\text{cm}$, as compared to a simple model which includes only collisions from the main bunch and from the satellites. % For this analysis, the simulated vertex distribution is re-weighted such that the full interaction region is described. The difference of \zvtxerror{} in ${\cal L}^{\text{QEDC}}$, obtained when using the simple beam profile model, is taken as a systematic uncertainty. % The reconstructed $z$-vertex distribution after re-weighting is compared to the data in figure \ref{fig:selection}f. Good agreement is found. % The regions of sizable systematic uncertainty due to the vertex re-weighting are visible. The identification of clusters in the SpaCal is checked by relaxing the $R_{\log{}}<6\,\text{cm}$ condition. Removing the $R_{\log{}}$ condition results in a somewhat smaller number of selected QEDC events, because a third SpaCal cluster is accepted more often, leading to the rejection of the corresponding events. %such that these events are rejected. % This procedure leads to a change in %The resulting shift on ${\cal L}^{\text{QEDC}}$ of \clustererror{} which is considered as systematic uncertainty related to the cut in $R_{\log{}}$. % The uncertainty on the CTD track reconstruction efficiency of typically $2\%$ per track affects the analysis through the track veto, resulting in an uncertainty of \ctdefferror{} on ${\cal L}^{\text{QEDC}}$. % The veto on the energy in the forward part of the LAr calorimeter is checked by relaxing the veto condition to $E<1\,\text{GeV}$, resulting in an uncertainty of \larvetoerror{} on ${\cal L}^{\text{QEDC}}$. \subsection{Background uncertainties} The normalisation of quasi-elastic and inelastic QEDC events predicted by the COMPTON22 event generator depends mainly on parameters related to the nucleon excitation and on the proton structure function parametrisation at low momentum transfer, respectively. % These parameters are not known very well. % For this reason, the normalisation of the sum of these contributions, referred to as ``non-elastic QEDC'', is verified by investigating the vector sum of the electron and photon transverse momenta, $\vec{P}_T^{\text{sum}}=-\vec{P}_T^{\text{miss}}$. % The vector $\vec{P}_T^{\text{sum}}$ is decomposed into components parallel to ($P_T^{\parallel}$) and perpendicular to ($P_T^{\perp}$) the electron transverse momentum. % The distributions of $P_T^{\parallel}$ and $P_T^{\perp}$ are shown in figure \ref{fig:ptcomponents} inside the analysis phase space as well as for $\vert\vec{P}_T^{\text{miss}}\vert>0.3\,\text{GeV}$ with all other selection criteria applied. % Both the parallel and the perpendicular components are described well. % Outside the nominal analysis phase space the non-elastic QEDC contributions dominate at large $P_T^{\parallel}$ or at large $P_T^{\perp}$. % The $P_T^{\parallel}$ distribution is somewhat asymmetric for $\vert\vec{P}_T^{\text{miss}}\vert>0.3\,\text{GeV}$, because in contrast to photons the SpaCal response to electrons has tails towards low energies, as discussed in the previous section. % The normalisation of the non-elastic QEDC contribution is tested by performing fits to either $P_T^{\parallel}$ or $P_T^{\perp}$ for $\vert\vec{P}_T^{\text{miss}}\vert>0.3\,\text{GeV}$. % The normalisation factors observed in these fits are compatible with the COMPTON22 prediction within $25\%$, which is used as normalisation uncertainty for the non-elastic QEDC processes. The DVCS cross section predictions obtained with the MILOU program are in agreement with recent H1 measurements \cite{:2009vda}. % Uncertainties of $20\%$ for the elastic DVCS process and $50\%$ for proton dissociative DVCS are considered. The elastic VM production rates are normalised using dedicated selections as close as possible to the QEDC analysis. % However, instead of requiring a photon in the SpaCal, a vector meson is reconstructed. % The decays $\rho^0\to\pi^{+}\pi^{-}$, $\phi\to K^{+}K^{-}$, $J/\psi\to \ell^{+}\ell^{-}$ ($\ell=e,\mu$) and $\psi^{\prime}\to \ell^{+}\ell^{-}$ are reconstructed from two oppositely charged tracks, detected in the CTD. % The decay $\omega\to\pi^{+}\pi^{-}\pi^{0}$ is reconstructed from two charged tracks and one or two neutral calorimeter clusters. % The decay $\Upsilon\to e^{+}e^{-}$ is reconstructed using a sample of photoproduction events, where an $e^{+}e^{-}$ pair from the $\Upsilon$ decay is reconstructed in the SpaCal, one of the SpaCal clusters matched with a central track, and the scattered electron is outside the acceptance of the H1 detector. The following normalisation uncertainties are found: $20\%$ on $\rho^0$ and $\phi$ production, $50\%$ on $J/\psi$, and $100\%$ on $\omega$, $\psi^{\prime}$ and $\Upsilon$. Possible contributions from $\rho(1450)$ production are covered by the $\rho^0$ normalisation uncertainty. The rate of non-resonant diffractive events, simulated using RAPGAP, is normalised using a selection of low multiplicity final states, where the electron is reconstructed in the SpaCal and one up to three additional charged or neutral particles are found. The uncertainty is estimated to be $30\%$. For the QED processes modelled by GRAPE, an uncertainty of $10\%$ is assumed, taking into account possible higher order effects. {\textchanged The uncertainties on the various background samples originating from the finite sizes of the generated event samples are negligible as compared to the uncertainties discussed above.} \subsection{QEDC theory uncertainties} \label{text:theoryerror} Uncertainties to the elastic QEDC cross section arise mainly from two sources: higher order corrections and the knowledge of the proton form factors. In the original COMPTON22 event generator, higher orders are simulated in the peaking approximation \cite{peakingisr}. % Improved higher order corrections have been calculated \cite{Anlauf:1991vi} %in the leading logarithmic approximation using a photon radiator \cite{Montagna:1991ku}. % For the purpose of this analysis, the COMPTON22 events are assigned weights, determined such that the cross section predicted by the photon radiator method is reproduced. % The difference {\textchanged of \isrerror{}} to the original COMPTON22 prediction is taken as systematic uncertainty due to higher order effects. % The elastic QEDC cross section also depends on the proton electric and magnetic form factors. % {\textchanged In the original COMPTON22 generator, only a simple parametrisation of the form factors, using one parameter, is implemented. % For this analysis, recent form-factor fits of elastic $ep$ scattering data \cite{Arrington:2007ux} are taken into account, using an event weighting technique. % The form-factor parameterisations \cite{Arrington:2007ux} are corrected for the effects of two-photon exchange (TPE), but parameterisations not including TPE corrections are also provided. % In the COMPTON22 generator, TPE contributions are not included when calculating the cross-section. % For this reason the form factor parameterisations not including TPE corrections are used to calculate elastic QEDC cross sections \cite{arringtonprivate}. % When computing the COMPTON22 cross section with the TPE corrected form factors, the analysis result changes by \tpeerror{}. This difference is included as systematic uncertainty. % Experimental uncertainties on the form factor parametrisation are also considered. % Such uncertainties are available with the parametrisation \cite{Venkat:2010by} which includes TPE corrections. It has been verified that the difference between \cite{Venkat:2010by} and \cite{Arrington:2007ux} is completely negligible for the purpose of this analysis if the TPE corrections are included. % The experimental uncertainties on the form factors $G_E(\vert t\vert)$ and $G_M(\vert t\vert)$ approach zero for $\vert t\vert\to 0$, because the parameterisations enforce $G_E(0)=1$ and $G_M(0)=\mu_p$. % At $\vert t\vert=0.007$, which is the average momentum transfer predicted for the events selected in this analysis, the uncertainties on $G_E$ ($G_M$) are $0.1\%$ ($0.2\%$). % The elastic QEDC cross section at fixed $\vert t\vert$ is given by a linear combination of $G_E^2(\vert t\vert)$ and $G_M^2(\vert t\vert)$, where the $G_E^2$ contribution dominates at small $\vert t \vert$. % When propagating the $\vert t\vert$ dependent parametrisation uncertainties on $G_E$ and $G_M$ to the luminosity measurement, an uncertainty of \gegmerror{} is found. % In addition to the theory uncertainties related to higher order corrections and proton form factors, discussed above, there is a statistical uncertainty on the theory prediction originating from the finite size of the generated signal event sample. % It amounts to \mcstaterror{}.} In figure \ref{fig:epz} the distribution of the variable $(E-p_z)/(2E^0_e)$ is studied. % This variable is calculated from the sum of the four-momenta of the electron and the photon. % The distribution of this variable is expected to peak at $1$. % The tail to small $(E-p_z)/(2E^0_e)$ originates from initial state radiation, whereas values larger than $1$ show up due to resolution effects. % The data are described within the systematic and statistical uncertainties for both small and large values of $(E-p_z)/(2E^0_e)$. As expected, the peak region is dominated by experimental uncertainties, whereas the region of small $(E-p_z)/(2E^0_e)$ is dominated by uncertainties of the QEDC cross section. \subsection{Time-dependent uncertainties} In order to apply the integrated luminosity ${\cal L}^{\text{QEDC}}$ to other physics analyses, possibly using time restricted H1 data sets, a luminosity calculation differential in time is required. % This is achieved using DIS events measured in the SpaCal. The DIS selection follows the selection described in \cite{Collaboration:2010ry} but is restricted in electron polar angle to the range $167^{\circ}<\theta_e<172^{\circ}$ such that the expected event yield is most insensitive to changing beam conditions, in particular to the average position of the interaction vertex in $z$. % In addition, the electron energy range is restricted to $150.3\,\text{GeV}$. The data are shown as black dots with the statistical uncertainties indicated as vertical bars. The simulation including background, normalised to the integrated luminosity as determined in this analysis, is indicated as a solid line, with the systematic uncertainties attached as shaded area. Also shown are the contributions from non-elastic QEDC and from other background sources.} \end{center} \end{figure} \begin{figure}[htbp] \begin{center} \includegraphics[width=0.5\figwidth]{d12-062f6a.eps} \includegraphics[width=0.5\figwidth]{d12-062f6b.eps} \caption{\label{fig:epz} Distribution of the variable $(E-p_z)/{2E^{e}_0}$ calculated from the sum of the electron and photon four-momenta. In (a) the event counts are shown, whereas in (b) the ratio of data to expectation is drawn. The data are shown as black dots with the statistical uncertainties indicated as vertical bars. The simulation including background, normalised to the integrated luminosity as determined in this analysis, is indicated as a solid line, with various components of the systematic uncertainties attached as shaded areas.} \end{center} \end{figure} \begin{figure}[htbp] \begin{center} \includegraphics[width=\figwidth]{d12-062f7.eps} \caption{\label{fig:yield} Integrated luminosity measured from elastic QEDC events in bins of approximately $25\,\text{pb}^{-1}$, divided by the integrated luminosity derived from the QEDC analysis on the full sample with time-dependent corrections applied. The statistical uncertainties of the binned QEDC analysis as well as the uncertainties of the time-dependent corrections, here applied to four data taking periods (I)--(IV), are indicated.} \end{center} \end{figure} \end{document}