\documentclass[12pt]{article}
\usepackage{graphicx}
\usepackage{epstopdf}

\usepackage[T1]{fontenc}
\usepackage{amsmath}
\usepackage{amssymb}
\usepackage{times}
\usepackage{txfonts}
\usepackage[mathlines,displaymath]{lineno}
\usepackage{pdflscape}
\usepackage{afterpage}
\usepackage{rotating}
\DeclareMathAlphabet{\mathbold}{OML}{txr}{b}{it}


%\documentclass[12pt]{article}
\usepackage{epstopdf}
\usepackage{epsfig}
%\usepackage{amsmath}
\usepackage{hhline}
%\usepackage{amssymb}
%\usepackage{times}
\usepackage{cite}
\usepackage{subfig} 
\usepackage{multirow}
%%%%%%%%%%%% Comment the next two lines to remove the line numbering
%\usepackage[]{lineno}
%\linenumbers
%%%%%%%%%%%%
\RequirePackage{filecontents} 
\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{\pom}{{I\!\!P}}
\newcommand{\reg}{{I\!\!R}}
\newcommand{\qsq}{\ensuremath{Q^2} }
\newcommand{\mxq}{\ensuremath{M_{X}^{2}} }
\newcommand{\gevsq}{\ensuremath{\mathrm{GeV}^2} }
\newcommand{\xpom}{x_{\pom}}
\newcommand{\zpom}{z_{\pom}}
\newcommand{\muf}{\ensuremath{\mu_{f}}}
\newcommand{\mur}{\ensuremath{\mu_{r}}}
\newcommand{\as}{\ensuremath{\alpha_s}}
\newcommand{\asmz}{\ensuremath{\alpha_s(M_Z)}}
\newcommand{\asmur}{\ensuremath{\alpha_s(\mur)}}
\newcommand{\gp}{\gamma\ -\ \rm p}

\renewcommand{\arraystretch}{1.2}

\begin{document}

\begin{titlepage}

\noindent
\begin{flushleft}
{\tt DESY 14-200    \hfill    ISSN 0418-9833} \\
{\tt November 2014}                  \\
\vspace{\baselineskip}
%{Date: November 14, 2014} \\
%{Version: 4.3, after final reading} \\
%{Editors: Boris Pokorny, Daniel Britzger} \\
%{Referees: Daniel Britzger, Sergey Levonian} \\
\end{flushleft}

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

{\bf Measurement of Dijet Production in Diffractive Deep-Inelastic $\mathbold{ep}$ Scattering at HERA}

\vspace{2cm}

H1 Collaboration

\end{Large}
\end{center}

\vspace{2cm}

\begin{abstract}
A measurement is presented of single- and double-differential dijet cross sections
in diffractive deep-inelastic $ep$ scattering at HERA 
using data collected by the H1 experiment corresponding to an integrated luminosity
of $290\ \rm{pb^{-1}}$.
The investigated phase space is spanned by the photon virtuality
in the range of $4<Q^{2}<100\ \rm GeV^{2}$ and by the fractional proton
longitudinal momentum loss $\xpom<0.03$.
The resulting cross sections are compared with next-to-leading order QCD predictions
based on diffractive parton distribution functions and 
the value of the strong coupling constant is extracted.
%the measured cross sections as $\alpha_{s}(M_{Z})=0.119\ (4)_{\rm exp}\ (12)_{\rm theo}$.
\end{abstract}

\vspace{1.5cm}

\begin{center}
%To be submitted to \it{Eur. Phys. J.} \bf{C} 
Submitted to \it{JHEP} 
\end{center}

\end{titlepage}

%\include{h1auts}
 
%-- H1AUTS Author list by names 
%-- Status: Wed Sep 17 11:48:12 CEST 2014  Number of authors = 143 

V.~Andreev$^{21}$,             %LPI -PD        8/88            Andreev             
A.~Baghdasaryan$^{33}$,        %YERE-PD        09/03           Baghdasaryana       
K.~Begzsuren$^{30}$,           %ULBA-PD        04/06           Begzsuren           
A.~Belousov$^{21}$,            %LPI -HON       12/12           Belousov            
V.~Boudry$^{24}$,              %ECPL-PD        1/93            Boudry              
G.~Brandt$^{45}$,              %OXFU-PD        03/07           Brandt              
V.~Brisson$^{23}$,             %ORSA-HON       01/12           Brisson             
D.~Britzger$^{10}$,            %DESY-PD        06/13           Britzger            
A.~Buniatyan$^{2}$,            %BIRM-HON       06/14           Buniatyan           
A.~Bylinkin$^{20,42}$,         %ITEP-ST        09/12           Bylinkin            
L.~Bystritskaya$^{20}$,        %ITEP-PD        05/99           Bystritskaya        
A.J.~Campbell$^{10}$,          %DESY-PD        10/84           Campbella           
K.B.~Cantun~Avila$^{19}$,      %MEX1-PD        05/13           Cantunavila         
F.~Ceccopieri$^{3}$,           %BRUX-LEFT      10/13           Ceccopieri          
K.~Cerny$^{27}$,               %PRG2-PD        11/08           Cernyk              
V.~Chekelian$^{22}$,           %MPIM-PD        01/90           Chekelian           
J.G.~Contreras$^{19}$,         %MEX1-PROF      03/98           Contreras           
J.~Cvach$^{26}$,               %PRAG-HON       01/10           Cvach
J.B.~Dainton$^{16}$,           %LIVE-PROF      01/94           Dainton             
K.~Daum$^{32,37}$,             %WUPP-PROF      06/90           Daum                
C.~Diaconu$^{18}$,             %MARS-PD        09/96           Diaconu             
M.~Dobre$^{4}$,                %BUCH-PD        12/12           Dobre               
V.~Dodonov$^{10}$,             %DESY-PD        04/98           Dodonov             
G.~Eckerlin$^{10}$,            %DESY-PD        8/88            Eckerlin            
S.~Egli$^{31}$,                %PSI -PD        01/10           Egli                
E.~Elsen$^{10}$,               %DESY-PROF      01/06           Elsen               
L.~Favart$^{3}$,               %BRUX-PROF      10/99           Favart              
A.~Fedotov$^{20}$,             %ITEP-PD        8/88            Fedotov             
J.~Feltesse$^{9}$,             %SACL-HON       11/06           Feltesse            
J.~Ferencei$^{14}$,            %KOSI-PD        08/88           Ferencei            
M.~Fleischer$^{10}$,           %DESY-PD        07/95           Fleischer           
A.~Fomenko$^{21}$,             %LPI -PD        02/89           Fomenko             
E.~Gabathuler$^{16}$,          %LIVE-HON       01/12           Gabathulere         
J.~Gayler$^{10}$,              %DESY-HON       02/05           Gayler              
S.~Ghazaryan$^{10}$,           %DFLC-PD        09/09           Ghazaryan           
A.~Glazov$^{10}$,              %DESY-PD        01/04           Glazov              
L.~Goerlich$^{6}$,             %CRAC-PD        8/88            Goerlich            
N.~Gogitidze$^{21}$,           %LPI -PD        05/91           Gogitidze           
M.~Gouzevitch$^{10,38}$,       %DESY-PD        10/11           Gouzevitch          
C.~Grab$^{35}$,                %ZUTH-PROF      08/07           Grab                
A.~Grebenyuk$^{3}$,            %BRUX-PD        04/12           Grebenyuk           
T.~Greenshaw$^{16}$,           %LIVE-PD        8/88            Greenshaw           
G.~Grindhammer$^{22}$,         %MPIM-HON       09/11           Grindhammer         
D.~Haidt$^{10}$,               %DESY-HON       03/04           Haidt               
R.C.W.~Henderson$^{15}$,       %LANC-PROF      10/04           Henderson           
M.~Herbst$^{13}$,              %HDB2-PD        11/11           Herbst              
J.~Hladk\`y$^{26}$,            %PRAG-HON       12/02           Hladky              
D.~Hoffmann$^{18}$,            %MARS-PD        07/00           Hoffmann            
R.~Horisberger$^{31}$,         %PSI -PD        01/10           Horisberger         
T.~Hreus$^{3}$,                %BRUX-PD        10/08           Hreus               
F.~Huber$^{12}$,               %HDB1-PD        01/13           Huberf              
M.~Jacquet$^{23}$,             %ORSA-PD        09/96           Jacquet             
X.~Janssen$^{3}$,              %ANTW-PD        02/03           Janssenx            
H.~Jung$^{10,3}$,              %DESY-PROF      12/99           Jungh               
M.~Kapichine$^{8}$,            %JINR-PD        3/97            Kapichine           
C.~Kiesling$^{22}$,            %MPIM-HON       09/11           Kiesling            
M.~Klein$^{16}$,               %LIVE-PROF      12/06           Klein               
C.~Kleinwort$^{10}$,           %DESY-PD        8/88            Kleinwort           
R.~Kogler$^{11}$,              %HAM2-PD        12/10           Kogler              
P.~Kostka$^{16}$,              %ZEUT-PD        10/87           Kostka              
J.~Kretzschmar$^{16}$,         %LIVE-PD        03/08           Kretzschmar         
K.~Kr\"uger$^{10}$,            %DESY-PD        01/04           Kruegerk            
M.P.J.~Landon$^{17}$,          %QMWC-PD        10/83           Landon              
W.~Lange$^{34}$,               %ZEUT-PD        8/88            Lange               
P.~Laycock$^{16}$,             %LIVE-PD        11/03           Laycock             
A.~Lebedev$^{21}$,             %LPI -HON       12/12           Lebedev             
S.~Levonian$^{10}$,            %DESY-PD        08/88           Levonian            
K.~Lipka$^{10,41}$,            %DESY-PD        06/01           Lipka               
B.~List$^{10}$,                %DESY-PD        11/99           Listb               
J.~List$^{10}$,                %DFLC-PD        10/00           Listj               
B.~Lobodzinski$^{22}$,         %MPIM-PD        01/06           Lobodzinski         
E.~Malinovski$^{21}$,          %LPI -PD        01/89           Malinovskie         
H.-U.~Martyn$^{1}$,            %AAC1-HON       01/12           Martyn              
S.J.~Maxfield$^{16}$,          %LIVE-PD        8/88            Maxfield            
A.~Mehta$^{16}$,               %LIVE-PROF      06/00           Mehta               
A.B.~Meyer$^{10}$,             %DESY-PD        10/97           Meyeran             
H.~Meyer$^{32}$,               %WUPP-HON       12/00           Meyerhi             
J.~Meyer$^{10}$,               %DESY-PROF      10/84           Meyerj              
S.~Mikocki$^{6}$,              %CRAC-PROF      10/10           Mikocki             
A.~Morozov$^{8}$,              %JINR-PD        06/99           Morozova            
K.~M\"uller$^{36}$,            %ZUER-PD        11/94           Muellerk            
Th.~Naumann$^{34}$,            %ZEUT-PROF      01/05           Naumannt            
P.R.~Newman$^{2}$,             %BIRM-PROF      10/09           Newman              
C.~Niebuhr$^{10}$,             %DESY-PD        3/93            Niebuhr             
G.~Nowak$^{6}$,                %CRAC-PROF      02/12           Nowakg              
J.E.~Olsson$^{10}$,            %DESY-HON       09/10           Olsson              
D.~Ozerov$^{10}$,              %DESY-PD        09/08           Ozerov              
P.~Pahl$^{10}$,                %DESY-LEFT      05/14           Pahl                
C.~Pascaud$^{23}$,             %ORSA-HON       09/05           Pascaud             
G.D.~Patel$^{16}$,             %LIVE-PD        8/88            Patel               
E.~Perez$^{9,39}$,             %SACL-PD        10/07           Perez               
A.~Petrukhin$^{10}$,           %DESY-PD        09/09           Petrukhin           
I.~Picuric$^{25}$,             %PODG-PROF      12/07           Picuric             
H.~Pirumov$^{10}$,             %DESY-PD        04/13           Pirumov             
D.~Pitzl$^{10}$,               %DESY-PD        8/88            Pitzl               
R.~Pla\v{c}akyt\.{e}$^{10,41}$, %DESY-PD        10/06           Placakyte           
B.~Pokorny$^{27}$,             %PRG2-ST        10/09           Pokorny             
R.~Polifka$^{27,43}$,          %PRG2-PD        05/11           Polifka             
V.~Radescu$^{10,41}$,          %DESY-PD        10/06           Radescu             
N.~Raicevic$^{25}$,            %PODG-PD        06/02           Raicevic            
T.~Ravdandorj$^{30}$,          %ULBA-PROF      09/00           Ravdandorj          
P.~Reimer$^{26}$,              %PRAG-PD        07/90           Reimer              
E.~Rizvi$^{17}$,               %QMWC-PROF      03/04           Rizvi               
P.~Robmann$^{36}$,             %ZUER-PD        8/88            Robmann             
R.~Roosen$^{3}$,               %BRUX-HON       03/12           Roosen              
A.~Rostovtsev$^{20}$,          %ITEP-PROF      02/11           Rostovtsev          
M.~Rotaru$^{4}$,               %BUCH-PD        06/11           Rotaru              
S.~Rusakov$^{21}$,             %LPI -HON       01/12           Rusakov             
D.~\v S\'alek$^{27}$,          %PRG2-PD        10/10           Salek               
D.P.C.~Sankey$^{5}$,           %RAL -PD        01/90           Sankey              
M.~Sauter$^{12}$,              %HDB1-PD        10/09           Sauter              
E.~Sauvan$^{18,44}$,           %MARS-PD        11/1            Sauvan              
S.~Schmitt$^{10}$,             %DESY-PD        09/99           Schmittst           
L.~Schoeffel$^{9}$,            %SACL-PROF      10/10           Schoeffel           
A.~Sch\"oning$^{12}$,          %HDB1-PROF      01/09           Schoening           
H.-C.~Schultz-Coulon$^{13}$,   %HDB2-PROF      01/04           Schultzcoulon       
F.~Sefkow$^{10}$,              %DFLC-PD        09/99           Sefkow              
S.~Shushkevich$^{10}$,         %DESY-PD        08/11           Shushkevich         
Y.~Soloviev$^{10,21}$,         %LPI -PD        08/89           Soloviev            
P.~Sopicki$^{6}$,              %CRAC-PD        03/14           Sopicki             
D.~South$^{10}$,               %DESY-PD        05/03           South               
V.~Spaskov$^{8}$,              %JINR-PD        12/97           Spaskov             
A.~Specka$^{24}$,              %ECPL-PROF      09/05           Specka              
M.~Steder$^{10}$,              %DESY-PD        09/08           Steder              
B.~Stella$^{28}$,              %ROME-HON       11/10           Stella              
U.~Straumann$^{36}$,           %ZUER-PROF      09/95           Straumann           
T.~Sykora$^{3,27}$,            %ANTW-PD        01/06           Sykora              
P.D.~Thompson$^{2}$,           %BIRM-PD        08/99           Thompsonp           
D.~Traynor$^{17}$,             %QMWC-PD        10/01           Traynor             
P.~Tru\"ol$^{36}$,             %ZUER-HON       10/06           Truoel              
I.~Tsakov$^{29}$,              %SOFI-PROF      03/12           Tsakov              
B.~Tseepeldorj$^{30,40}$,      %ULBA-PROF      06/06           Tseepeldorj         
J.~Turnau$^{6}$,               %CRAC-PROF      09/95           Turnau              
A.~Valk\'arov\'a$^{27}$,       %PRG2-PD        8/88            Valkarova           
C.~Vall\'ee$^{18}$,            %MARS-PD        06/86           Vallee              
P.~Van~Mechelen$^{3}$,         %ANTW-PROF      01/06           Vanmechelen         
Y.~Vazdik$^{21}$,              %LPI -PD        8/88            Vazdik              
D.~Wegener$^{7}$,              %DORT-PROF      01/06           Wegener             
E.~W\"unsch$^{10}$,            %DESY-PD        8/88            Wuensch             
J.~\v{Z}\'a\v{c}ek$^{27}$,     %PRG2-PROF      10/05           Zacek               
Z.~Zhang$^{23}$,               %ORSA-PD        10/92           Zhang               
R.~\v{Z}leb\v{c}\'{i}k$^{27}$, %PRG2-ST        03/12           Zlebcik             
H.~Zohrabyan$^{33}$,           %YERE-PD        11/02           Zohrabyan           
and
F.~Zomer$^{23}$                %ORSA-PROF      09/06           Zomer          

%\newpage

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

\bigskip
 $ ^{37}$ Also at Rechenzentrum, Universit\"at Wuppertal,
          Wuppertal, Germany \\
 $ ^{38}$ Also at IPNL, Universit\'e Claude Bernard Lyon 1, CNRS/IN2P3,
          Villeurbanne, France \\
 $ ^{39}$ Also at CERN, Geneva, Switzerland \\
 $ ^{40}$ Also at Ulaanbaatar University, Ulaanbaatar, Mongolia \\
 $ ^{41}$ Supported by the Initiative and Networking Fund of the
          Helmholtz Association (HGF) under the contract VH-NG-401 and S0-072 \\
 $ ^{42}$ Also at Moscow Institute of Physics and Technology, Moscow, Russia \\
 $ ^{43}$ Also at  Department of Physics, University of Toronto,
          Toronto, Ontario, Canada M5S 1A7 \\
 $ ^{44}$ Also at LAPP, Universit\'e de Savoie, CNRS/IN2P3,
          Annecy-le-Vieux, France \\
 $ ^{45}$ Department of Physics, Oxford University,
          Oxford, UK$^{ b}$ \\


\bigskip
 $ ^a$ Supported by the Bundesministerium f\"ur Bildung und Forschung, FRG,
      under contract numbers 05H09GUF, 05H09VHC, 05H09VHF,  05H16PEA \\
 $ ^b$ Supported by the UK Science and Technology Facilities Council,
      and formerly by the UK Particle Physics and
      Astronomy Research Council \\
 $ ^c$ Supported by FNRS-FWO-Vlaanderen, IISN-IIKW and IWT
      and  by Interuniversity
Attraction Poles Programme,
      Belgian Science Policy \\
 $ ^d$ Partially Supported by Polish Ministry of Science and Higher
      Education, grant  DPN/N168/DESY/2009 \\
 $ ^e$ Supported by VEGA SR grant no. 2/7062/ 27 \\
 $ ^f$ Supported by the Ministry of Education of the Czech Republic
      under the projects  LC527, INGO-LA09042 and
      MSM0021620859 \\
 $ ^g$ Supported by the Swiss National Science Foundation \\
 $ ^h$ Supported by  CONACYT,
      M\'exico, grant 48778-F \\
 $ ^i$ Russian Foundation for Basic Research (RFBR), grant no 1329.2008.2
      and Rosatom \\
 $ ^j$ Supported by the Romanian National Authority for Scientific Research
      under the contract PN 09370101 \\
 $ ^k$ Partially Supported by Ministry of Science of Montenegro,
      no. 05-1/3-3352 \\
}
%\tableofcontents


%---------------------
%---------------------
\section{Introduction}

In deep-inelastic scattering (DIS),
diffractive reactions of the type $ep\to eXY$, where $X$ is a high-mass hadronic final state and
$Y$ is either the elastically scattered proton or its low-mass excitation,
represent about $10\%$ of the events at HERA and provide rich experimental input for testing
quantum chromodynamics (QCD) in the diffractive regime. These processes can be understood 
as probing by a virtual photon
emitted from the beam lepton a net colour singlet carrying vacuum quantum numbers (a pomeron)~\cite{pomeranchuk,goulianos}.
Due to the colourless exchange the systems $X$ and $Y$ are separated by a rapidity
interval free of hadronic activities.
In these processes at least one hard scale is involved such that perturbative QCD (pQCD) can be applied.

According to the QCD collinear factorisation theorem~\cite{trentadue},  
calculations of diffractive cross sections factorise into process dependent hard scattering coefficient functions
and a set of process independent diffractive parton distribution functions (DPDFs).  
While the hard scattering coefficient functions are calculable in pQCD,
the DPDFs have to be determined from QCD fits to the measured inclusive diffractive cross sections.
In such QCD fits~\cite{fitb}, DGLAP evolution~\cite{gribov,dokschitzer,altarelli} of the DPDFs is assumed.
The QCD factorisation theorem is proven to hold for inclusive and dijet diffractive processes~\cite{collins},
assuming high enough photon virtuality such that higher twist effects can be neglected.
The DPDFs are experimentally determined by assuming an additional
factorisation of the DPDFs dependence on the scattered proton momentum from the dependence on the other variables, ascribed to the structure of the colourless exchange.
This assumption is known as proton vertex factorisation.
A pomeron flux in the proton is introduced and universal parton densities
are attributed to the diffractively exchanged object.
Many measurements of diffraction in DIS suggest the validity 
of the proton vertex factorisation assumption in DIS~\cite{rfact_1,rfact_2,rfact_3,fitb}.

In leading order the inclusive diffractive cross section in $ep$ scattering
is proportional to the charge-squared weighted sum of the quark distribution functions in the pomeron,
while %the gluon component 
its gluon content can be determined only indirectly via scaling violations.
%As heavy quarks or jets are dominantly produced in gluon-induced processes~\cite{bgf_1,fitb},
As events with two jets (dijets) are readily produced in gluon-induced processes,
measurements of diffractive dijet cross sections are 
sensitive to the value of the strong coupling $\alpha_s$  and
to the gluon content of the pomeron.
The production of dijets in diffractive DIS has previously been studied
at HERA using either the large rapidity gap (LRG) method~\cite{jetdiff_1,jetdiff_3,jetdiff_zeus} or
by direct detection of the outgoing proton~\cite{jetdiff_2}.
 

%A diffractive dijet production cross section measurement, based on $ep$ collision data 
%collected in the years 2005-2007 by the H1 detector, is presented in this paper.
In this paper cross section measurements of dijet production in diffractive $ep$ scattering
are presented, based on data collected in the years 2005-2007 with the H1 detector at HERA.
Diffractive events are selected
by means of the LRG method, requiring a clear separation in rapidity of the final state systems
$X$ and $Y$. % to be recognised by the H1 detector. 
The measured cross sections are compared to next-to-leading order (NLO) QCD predictions evaluated with input DPDFs
determined in previous inclusive diffractive measurements by the H1 collaboration~\cite{fitb}.

The present analysis is based on the full HERA-II data sample resulting in significantly increased
statistics with respect to previous analyses. Furthermore, the cross sections are determined
using a regularised unfolding procedure which fully accounts for efficiencies, migrations and  correlations among the measurements. 
The measured dijet cross sections are used to extract the strong coupling constant $\alpha_{s}$ 
in diffractive DIS processes for the first time.
%The presented analysis profits from the use of a data sample of unprecedentedly high luminosity
%and incorporates a regularised unfolding procedure for treating the effects of acceptance, migrations and of statistical correlations among different
%observables. 

%Jet production is a key process for a precise determination of the strong coupling constant $\alpha_{s}$.
%In $ep$ collisions, the value of $\alpha_{s}$ was extracted using data covering 
%a wide range in the photon virtuality~\cite{Aaron:10:1,HighQ2Paper} as well as photoproduction~\cite{Abramowicz:2012jz}.   
%The achieved experimental precision of the measurement presented here allows
%for the first determination of $\alpha_{s}$ in diffractive DIS processes.

%---------------------
%---------------------
\section{Kinematics}

\begin{figure}[ht]
\centering
\includegraphics[width=.5\textwidth]{d14-200f1.eps}
\caption{Leading order diagram for the production of dijets in diffractive DIS.}
\label{fig:feyn}
\end{figure}

A leading order (LO) diagram of boson-gluon fusion, which is the dominant process for the production of two jets in diffractive DIS,
is depicted in Figure~\ref{fig:feyn}. 
The incoming electron\footnote{In this paper the term "electron" is used generically to refer to both electrons and positrons.} 
of four-momentum $k$ interacts with %undergoes scattering on 
the incoming proton of four-momentum $p$ via the exchange
of a virtual photon of four-momentum $q=k-k'$. The outgoing proton or its low-mass dissociation state carries
four-momentum $p'$. 
The DIS kinematics is described
by the following set of variables:
\begin{equation}
\label{_kine_dis}
\qsq=-q^{2}=(k-k')^{2}, \qquad x=\frac{\qsq}{2p\cdot q}, \qquad y=\frac{p\cdot q}{p\cdot k}, 
\end{equation}
where $\qsq$, $x$ and $y$ denote the photon virtuality,
the Bjorken-$x$ variable and the inelasticity of the process, respectively.
%with $\qsq$ denoting the photon virtuality, $x$ the longitudinal momentum fraction of the incoming proton
%carried by the quark coupling to the virtual photon
%and $y$ the inelasticity of the process.
Conservation laws stipulate the relation $\qsq=xys$, where $s$ stands for the $ep$ centre-of-mass energy squared.

The kinematics of the diffractive exchange is described in terms of the additional quantities
\begin{equation}
\label{_kine_diff}
%\xpom=\frac{q\cdot(p-p')}{q\cdot p}, \qquad \beta=\frac{\qsq}{2q\cdot(p-p')}=\frac{x}{\xpom}, \qquad t=(p-p')^{2}
\xpom=\frac{q\cdot(p-p')}{q\cdot p}, \qquad t=(p-p')^{2}
\end{equation}
with $\xpom$ and $t$ being the longitudinal momentum fraction of the incoming proton carried by the pomeron and 
the squared four-momentum transfer at the proton vertex, respectively. 
%The variable $\beta$ is interpreted as the longitudinal momentum fraction of the pomeron carried by the struck parton.
The fractional longitudinal momentum of the pomeron transferred to the dijet system is given by
\begin{equation}
\label{_kine_jet}
\zpom=\frac{q\cdot v}{q\cdot(p-p')}=\frac{x}{\xpom},
\end{equation}
where $v$ is the four-momentum of the parton entering the hard interaction. %in the LO picture.

%---------------------
%---------------------
\section{Monte Carlo Models and Fixed Order QCD Calculations}
\label{sec:model}
%elastic
The RAPGAP event generator~\cite{rapgap} allows for the simulation of processes $ep\to eXY$ 
including both leading (pomeron) and sub-leading (reggeon) exchanges.
Assuming the proton vertex factorisation, the parton densities obtained
in the previous QCD analysis of inclusive diffractive data (H12006 Fit-B)~\cite{fitb}
are convoluted with leading order QCD matrix elements. 
%The DPDF evolution is based on the DGLAP equations using QCD scales  $\mu_{r,f}=\sqrt{\hat{p}^{\ast2}_{\rm T}+Q^{2}}$,
%with $\hat{p}^{\ast}_{\rm T}$ being the transverse momentum of the outgoing partons with respect to the photon axis
%in the $\gamma^{\ast}$-$p$ centre-of-mass frame. 
Higher order QCD radiation effects are modelled via
initial and final state parton showers in the leading-$\log$ approximation~\cite{llog}.
Hadronisation is accounted for by making use of the Lund string model~\cite{lund} as implemented in PYTHIA~\cite{pythia}.

%dissociative
Within the diffractive selection based on the LRG method, the system $Y$ may also be a low mass dissociative system.
Proton dissociation events are simulated in the the range of $M_Y<20\ \mathrm{GeV}$ using the RAPGAP event generator, where  
$M_Y$ is the mass of the system $Y$.
Resonant contributions together with the continuum part of the $M_Y$ distribution are modelled similarly 
to the DIFFVM event generator~\cite{diffvm}. 
A small admixture of resolved $\gamma^{*}p$ scattering is included in fixed LO mode of jet production
in the low $\qsq$  region~\cite{resgamma}. The resolved photon contribution is simulated
with the RAPGAP event generator using the SaS-G PDF set~\cite{sasg2dfit} as the input PDF of the photon. 
% qed radiation
QED radiation effects are simulated with the HERACLES~\cite{heracles} program interfaced to RAPGAP.
Besides the Born level contribution, the simulated cross sections
include contributions from initial and final state emission of real photons from the electron, from vertex corrections
as well as from self energy diagrams.
%reweighting
As the H12006 Fit-B DPDF set has previously been observed to underestimate the data in the low $\qsq$ region,
a weighting is applied for $\qsq<7\ \mathrm{GeV^{2}}$,
parametrised as the ratio of the data in~\cite{fitb} to the Monte Carlo expectation based on the H12006 Fit-B DPDF set.

% non-diffr. background
Background arising from non-diffractive DIS processes
is also simulated with the RAPGAP event generator
using its inclusive mode together with the CTEQ6L PDF set~\cite{cteq}.

% simulation
The MC simulation is used to correct the data for detector effects. 
The generated events undergo the full GEANT~\cite{geant} simulation of the H1 detector
and are analysed in the same way as the real data.
In order to describe the measured distributions, the diffractive MC
is reweighted in several variables as discussed in~\ref{subsec:data}. 

QCD predictions of the dijet cross sections at the parton level are evaluated at NLO
using the NLOJET++ program~\cite{nlojetprog,nlojet}. 
The NLO pQCD predictions are calculated in the $\overline{\rm MS}$-scheme with five active flavors.
The two-loop approximation of the renormalisation group equation is used for the running of the strong coupling constant
with a coupling strength of $\alpha_s(M_Z)=0.118$.
The cross sections are evaluated in intervals of $\xpom$, effectively
replacing the beam proton by a pomeron (slicing method).
The H12006 Fit-B DPDF set is used in the calculation. 
The renormalisation and factorisation scales $\mu_{r}$ and $\mu_{f}$ are provided by the photon virtuality and  
the average transverse momentum of the leading and sub-leading jet, $\langle p^{\ast}_{\rm T}\rangle$, 
in the $\gamma^{\ast}$-$p$ centre-of-mass frame and are defined as $\mu_{r}=\mu_{f}=\sqrt{\langle p^{\ast}_{\rm T}\rangle^{2}+\qsq}$.
The uncertainty on the prediction due to missing higher orders is estimated by simultaneous variation of the renormalisation
and factorisation scales by factors of $0.5$ or $2$. 
An uncertainty on the NLO prediction from the experimental uncertainties on
the DPDF set is obtained using the eigenvector decomposition
of the uncertainties of the H12006 Fit-B DPDF set.
This uncertainty is propagated to the NLO prediction
using the sign-improved formulae for error propagation~\cite{dpdferr}.
A significant contribution to the uncertainty of the H12006 Fit-B set
originates from the restriction of the input data to $\zpom<0.8$ and the extrapolation of the DPDF to $\zpom>0.8$.

Whereas the measured cross sections are compared to the predictions obtained by the slicing method,
an alternative method of adapting the NLO calculations
for diffractive DIS is used in the $\alpha_s$ extraction.
In order to provide theory predictions  
with different values of $\asmz$, the fastNLO method~\cite{H1fastNLO,fastNLO,fnlo} is used.
%
Cross section predictions are obtained by folding tabulated matrix
elements obtained from NLOJET++~\cite{nlojetprog,nlojet} with the DPDF
parametrisation.
%
The matrix elements are determined as a function of the observable of
interest, the factorisation scale $\mu_F$ and the convolution variable $x$.
%
The relation $x=\xpom\zpom$ is used when folding with the DPDF.
%
This way predictions can be obtained for different choices of DPDFs, of $\alpha_s$ and of
%This method allows thus for a change of the DPDFs, the running coupling, 
the renormalisation and the factorisation scales
without having to calculate the matrix elements all over again.
%
Settings identical to the slicing method are used for parameters such
as renormalisation and factorisation scales or DPDF set and 
very good numerical agreement with the slicing method is found.
%
The uncertainty on the prediction due to missing higher orders
is estimated by varying the scales by a factor $f$, where $0.5<f<2$.

Since the measured cross sections are given at the level of stable hadrons, the QCD predicted cross sections
have to be corrected for effects of initial and final state parton showers, hadronisation and fragmentation.
These corrections are determined for each of the measured cross sections as the ratio of hadron to parton level cross sections,
predicted with the RAPGAP event generator. 
Two distinct models of parton showers, 
the leading-$\log$ approximation and the colour dipole model as implemented in the ARIADNE program~\cite{ariadne},
are used in this calculation. 
In each measurement interval the resulting correction is
taken as the average of the values predicted by the two
models and the uncertainties on the correction factors
are taken as half the difference of the two predictions.
The hadron level cross sections are on average about $5\%$  higher than the parton level cross sections.
The total uncertainty on the NLO QCD predictions is obtained as the quadratic sum of the uncertainties
from scale variation, DPDF fit and hadronisation uncertainties.

%---------------------
%---------------------
\section{Experimental Technique}
\subsection{H1 detector}

A detailed description of the detector can be found elsewhere~\cite{det_1}. 
Here only those detector components 
relevant for the present analysis are briefly described.  
A right-handed coordinate system with the origin at the nominal interaction point
and with the $z$-axis pointing in the proton beam direction is conventionally chosen
as the laboratory frame. The polar angle $\theta$ is measured with respect to the $z$-axis,
while the direction in the $x$-$y$ plane is defined by the azimuthal angle $\phi$.
The pseudorapidity is defined as $\eta=-\ln\tan(\theta/2)$. 

The liquid argon (LAr) sampling calorimeter~\cite{lar} is located inside a $1.15\ \rm T$ solenoidal field 
and covers the polar angular range $4^\circ\!<\!\theta\!<\!154^\circ$.
The energy resolutions for electromagnetic and hadronic showers as determined in test beam measurements~\cite{testbeam_1,testbeam_2}
are $\sigma(E)/E \propto 11\%/\sqrt{E/\rm GeV}\oplus 1\%$ and
$\sigma(E)/E \propto 50\%/\sqrt{E/\rm GeV}\oplus 2\%$, respectively.
The energy and scattering angle of the scattered electron is measured in a scintillating fibre calorimeter SpaCal~\cite{spac1,spac2}
with a resolution of $\sigma(E)/E \propto 7\%/\sqrt{E/\rm GeV}\oplus 1\%$.
The precision of the energy scale is $1\%$ covering the polar angular range $154^\circ\!<\theta_{e'}<\!174^\circ$. %for electrons and $7\%$ hadrons.
The measurement of the polar angle of the scattered electron $\theta_{e'}$ is improved by means of a backward proportional chamber (BPC).
%The BPC reconstruction efficiency is $80\%$.
The precision of the polar angle measurement is $1\ \rm mrad$.

Trajectories of charged particles are measured with the central tracking detector (CTD) 
%The central tracking detector (CTD) is 
located inside the LAr calorimeter 
%and covers the range in polar angle $20^\circ\!<\!\theta\!<\!160^\circ$. 
%in the polar angular range of $15^\circ\!<\!\theta\!<\!165^\circ$
%The measurement of charged particle transverse momenta is possible.
%The magnetic field enables the measurement 
%The trajectories of charged particles are measured 
with a transverse momentum resolution of 
$\sigma_{p_T}/p_T\simeq 0.2\ \%\cdot\ p_{T}/\rm GeV\oplus1.5\%$
in the polar angular range of $15^\circ\!<\!\theta\!<\!165^\circ$.

The information from CTD and LAr is used for the reconstruction
of the system $X$.
The interaction vertex position is determined event-by-event 
using the particle trajectories measured in CTD. 


The following H1 forward detectors are used in the LRG selection of diffractive events.
The forward muon detector (FMD) consists of six proportional chambers which are grouped into
two three-layer sections separated by a toroidal magnet. Although the nominal
coverage of FMD is $1.9<\eta<3.7$, particles with pseudorapidity up to $\eta\! \sim\! 6.5$
can be detected indirectly through their interactions with the beam transport system and detector support structures. 
The lead-scintillator Plug calorimeter is located at $z=4.9\ \rm m$ and covers the range $3.5<\eta<5.5$.
The very forward region is covered by the forward tagging system (FTS) comprising
scintillators surrounding the beam pipe. Only one station of FTS, situated at $z=28\ \rm m$
and covering the range $6.0<\eta<7.5$, is included in the present analysis. 

The instantaneous luminosity is monitored based on the rate of the Bethe-Heitler process $ep\to ep\gamma$.
The final state photon is detected by a photon detector located close to the beam pipe at $z=-103\ \rm m$. 
The precision of the integrated luminosity measurement is improved in a dedicated analysis of the QED Compton
process~\cite{qedclumi}.

\subsection{Reconstruction of observables}

The DIS observables $\qsq$, $x$ and $y$ are reconstructed using the electron-$\Sigma$ method ~\cite{esigma}.
Within this method, the photon virtuality $\qsq$ is reconstructed based on the measured four-momentum of the scattered electron,
while the inelasticity $y$ and Bjorken-$x$ are determined making use
of combined information from the  
hadronic final state (HFS) and the scattered electron.  

The four-momenta of the particles attributed to HFS 
are reconstructed using an algorithm which combines information
provided by the tracking system and the LAr calorimeter
by avoiding double counting of hadronic energy~\cite{peez:03,hellwig:04}. 
The calibration of the HFS energy scale derived in~\cite{kogler} is applied.
The performance of the calibration was studied by comparing 
the transverse momentum balance in data and MC in the kinematic domain of this analysis.
  
Jets are reconstructed in the $\gamma^{\ast}$-$p$ centre-of-mass frame
using the inclusive 
$k_{\rm T}$ jet algorithm ~\cite{kt} with the $p_{\rm T}$ recombination scheme as implemented in the FastJet program ~\cite{fastjet}. 
The jet distance parameter is set to $R=1.0$.
The transverse momenta and pseudorapidities of the leading
and sub-leading jets are denoted as $p^{\ast}_{\rm T,1}$, $\eta^{\ast}_{1}$ and $p^{\ast}_{\rm T,2}$, $\eta^{\ast}_{2}$, respectively\footnote{Observables in the $\gamma^{\ast}$-$p$ centre-of-mass frame are labelled with an asterisk.}.  
The invariant mass of the final state system $X$ is reconstructed as:
\begin{equation}
M_{X}=c(\eta_{max})\sqrt{P_{X}^{2}}, 
\end{equation}
where $P_{X}$ is the four-momentum of the system $X$ obtained as a vector sum of all particles contained 
in the HFS.
The MC simulation is used in order to derive the average correction for detector losses $c(\eta_{max})$, 
where $\eta_{max}$ is the pseudorapidity of
the most forward energy deposition above $800\ \rm MeV$ in the LAr calorimeter.
The momentum fractions $\xpom$ and $\zpom$ are reconstructed as: 
\begin{equation}
\label{_reco_xpom}
\xpom=\frac{\qsq + \mxq}{ys}
\end{equation}
and
\begin{equation}
\label{_reco_zpom}
\zpom=\frac{\qsq + M_{12}^{2}}{\qsq + \mxq},
\end{equation}
where $M_{12}$ is the invariant mass of the dijet system.  

%The variable $x_{\gamma}$, which is interpreted as the fractional momentum of the photon in resolved photon interactions~\cite{karel},
%is reconstructed as
%\begin{equation}
%\label{_reco_xgamma}
%x_{\gamma}=\frac{\sum_{1,2}(E^{\ast jet}-p^{\ast jet}_{z})_{i}}{\sum_{HFS}(E^{\ast}-p^{\ast}_{z})_{i}},
%\end{equation}
%where the sum in the numerator counts the two jets and the sum in the denominator runs over all HFS objects.

Cross sections for dijet production in diffractive DIS are measured differentially with respect to the variables $\qsq$, $y$, $\xpom$, $\zpom$,
$p^{\ast}_{\rm T,1}$, $p^{\ast}_{\rm T,2}$, $\langle p^{\ast}_{\rm T}\rangle=(p^{*}_{\rm T,1}+p^{\ast}_{\rm T,2})/2$ and 
$\Delta\eta^{\ast}=|\eta^{\ast}_{1}-\eta^{\ast}_{2}|$. 

\subsection{Event selection}

The measurement is based on the H1 data collected
in the years $2005$ to $2007$ with a total integrated luminosity of $290\ \rm pb^{-1}$. 
The nominal beam energies of the protons and electrons are $E_{p}=920\ \rm GeV$ and $E_{e}=27.6\ \rm GeV$, respectively.

The longitudinal position of the reconstructed event vertex 
is restricted to the range $-35<z_{vtx}<35\ \rm cm$. 
DIS events are selected by the identification of the scattered electron in the backward calorimeter SpaCal.
The isolated energy deposit of electromagnetic structure with the highest transverse momentum is identified as scattered electron and 
has to have a measured energy of at least $9.5\ \rm GeV$. 

Only events accepted by a trigger combining signals induced by the scattered electron in the SpaCal 
with minimum track information of  the CTD  are used in the analysis.
The trigger efficiency related to the CTD condition is found to be
$98\%$-$99\%$, depending on the detector configuration
and is reproduced by the MC simulation within $2\%$.
The trigger efficiency related to the SPACAL condition is better than $99\%$.

Residual non-DIS background is dominated by photoproduction processes, 
where a hadron is misidentified as the scattered electron, whereas the true
scattered electron escapes detection due to its small scattering angle.
This background is reduced to a negligible level by demanding
$35<\sum_{i}(E-p_{z})_{i}<75\ \rm GeV$, where the sum runs over all HFS particles and the scattered electron candidate.
Elastic QED Compton scattering $ep\to e\gamma p$ introduces another background contribution
which is suppressed by rejecting configurations with two back-to-back clusters in SpaCal.
  
Diffractive events are identified with the LRG method 
which requires an empty interval in rapidity between the systems $X$ and $Y$. 
The low-mass system $Y$ is produced at very large pseudorapidities
and escapes detection. The diffractive signature is thus defined by
the systems $X$ (in the main detector) and $Y$ (undetected).
The energy of any cluster in the forward region
of the LAr calorimeter is required to be below the noise level of $800\ \rm MeV$,
which is ensured by demanding $\eta_{max}<3.2$.
The variable $\eta_{max}$ corresponds to the LAr cluster above
the noise threshold which has the largest pseudorapidity.
Information provided by the forward detectors FMD, FTS and the Plug calorimeter
is used in order to extend the gap to rapidities beyond the LAr acceptance
and in order to suppress the proton dissociation contribution.
These detectors are required to show no signal above noise level~\cite{mythesis}.
At high momentum fractions $\xpom$,
the system $X$ tends to extend into the direction of the outgoing system $Y$
and the experimental separation of the systems $X$ and $Y$ is not possible.
The LRG selection method is thus applicable only in the region of $\xpom\lesssim0.03$.
The sample of DIS events satisfying the LRG criteria is dominated by the diffractive exchange, 
as the system $X$ is isolated
in the main part of the H1 detector, while the system $Y$ escapes undetected
down the beam pipe. The signal is dominated by proton-elastic processes, $ep\to eXp$, however,
a small fraction of proton dissociation events is also accepted by the LRG selection. 
The LRG requirements impose restrictions on the mass and scattering angle
of the hadronic system $Y$. These correspond approximately to the requirements
$M_{Y}<1.6\ \rm GeV$ and $|t|<1\ \rm GeV^{2}$. 
Migrations in these variables are modelled using MC simulations.

\begin{table}[h!]
\begin{center}
\begin{tabular}{|c|c|c|}
\hline
& \bf{Extended Analysis Phase Space} & \bf{Measurement Cross Section Phase Space} \\
\hline
\multirow{2}{*}{DIS} & $3 < Q^{2} < 100 \ \mathrm{GeV^{2}}$ & $4 < Q^{2} < 100 \ \mathrm{GeV^{2}}$ \\
& $y < 0.7$ & $0.1 < y < 0.7$ \\
\hline
\multirow{3}{*}{Diffraction} & $x_{\pom} < 0.04$ & $x_{\pom} < 0.03$ \\
& LRG requirements & $|t|<1 \ \mathrm{GeV^{2}}$ \\
& & $M_{Y} < 1.6 \ \mathrm{GeV} $ \\
\hline
\multirow{3}{*}{Dijets} & $p_{\rm T,1}^{*} > 3.0 \ \mathrm{GeV}$ & $p_{\rm T,1}^{*} > 5.5 \ \mathrm{GeV} $\\
& $p_{\rm T,2}^{*} > 3.0 \ \mathrm{GeV} $ & $p_{\rm T,2}^{*} > 4.0 \ \mathrm{GeV} $ \\
& $-2 < \eta^{\rm lab}_{1,2} < 2$ & $-1 < \eta^{\rm lab}_{1,2} < 2$ \\
\hline
\end{tabular}
\end{center}
\caption{Summary of the extended analysis phase space and the phase space for the dijet cross sections measurements.}
\label{tab:ps}
\end{table}

Events are selected in a phase space which is extended compared to the measurement phase space
in order to improve the precision of the measurement 
by accounting for migrations at the phase space boundaries.
Events within the DIS phase space of $y<0.7$ and $3<\qsq<100\ \rm GeV^2$ are selected.
The events are required to have at least two jets in the pseudorapidity range $-2<\eta^{\rm lab}_{1,2}<2$
and transverse momenta greater than $3\ \rm GeV$ in the $\gamma^{\ast}$-$p$ centre-of-mass frame.

The measurement phase is defined by the DIS requirements of $0.1<y<0.7$ and $4<\qsq<100\ \rm GeV^2$.
The pseudorapidity of jets is restricted in the laboratory frame to $-1<\eta^{\rm lab}_{1,2}<2$ to
ensure the jets to be contained well within the central detector. 
The transverse momenta
of the leading and sub-leading jets are required to be larger than $5.5\ \rm GeV$ and $4.0\ \rm GeV$, respectively.
The extended phase space and the measurement phase space definitions are summarised in table~\ref{tab:ps}.
The total number of events accepted by the LRG selection criteria together with the DIS and jet requirements is $\sim\! 50000$
and $\sim\! 15000$ for the extended and measurement phase space, respectively.

\subsection{Corrections to the data}
\label{subsec:data}

Cross sections at the level of stable hadrons are
obtained from the measured event rates in data
by applying corrections determined using the MC simulation.
In figure~\ref{fig:figctr} kinematic distributions
of the observables $\qsq$, $p^{\ast}_{\rm T,1}$, $\xpom$ and $\zpom$ as observed in the detector are shown
in comparison to the expectations from the reweighted MC simulation. The overall good description
of the data is achieved after applying 
a dedicated weighting of the MC simulation
in the variables $\zpom$, $\xpom$ and
%
 $x_{dijet}=
\sum_{1,2}(E^{\ast jet}-p^{\ast
  jet}_{z})_{i}/\sum_{HFS}(E^{\ast}-p^{\ast}_{z})_{i}$.
%
Weights are obtained from the reconstructed kinematic distributions 
and are applied at the hadron level.
This procedure is iterated until a good description of the shapes of the observables is achieved.

The data are corrected for detector inefficiencies, acceptance and 
finite resolution using the regularised unfolding procedure as implemented in TUnfold~\cite{unf_3}. 
A detector response matrix $A$, with elements $a_{ij}$ expressing the probability for an observable
originating in the generated MC sample from an interval $i$ to be measured in an interval $j$, 
is determined using the MC simulation. Migrations from outside the measurement 
phase space are included by additional rows of the detector response matrix.
The domains of jets with $3.0<p^{*}_{\rm T,1}<5.5\ \rm GeV$ and of events
with $0.03<\xpom<0.04$ are found to be the dominating sources of these migrations.
The MC simulation is reweighted in order to describe
the data also in these regions beyond the nominal phase space.

Two sources of background are considered in this analysis
and are subtracted from the data using Monte Carlo simulations prior to unfolding:
diffractive dijet events with $M_{Y}>1.6\ \rm GeV$ and $|t|<1\ \rm GeV^2$
and background from non-diffractive DIS.

For a background subtracted measurement $y_{j}$, the corresponding
number of events in the truth bin $i$, $x_{i}$,
is found by solving a minimisation problem %~\cite{unf_1,unf_2} 
for a $\chi^{2}$ function
\begin{equation}
\chi^{2}=(y-Ax)^{T}V_{yy}^{-1}(y-Ax)+\tau^{2}x^{2}, 
\label{eq:unf}
\end{equation} 
where $x$ and $y$ are vectors defined by $y_{j}$ and $x_{i}$, respectively,
$V_{yy}$ is the covariance matrix accounting for the statistical uncertainties of $y_{j}$ and 
$\tau$ is a regularisation parameter introduced in order to damp statistical fluctuations of the solution.
The regularisation parameter $\tau$ is determined using the L-Curve scan~\cite{unf_3}.

The cross section %is evaluated 
in each measurement interval $i$ is given by
%as:
\begin{equation}
\sigma_{i}(ep\to ep'X)=\frac{x_{i}}{\cal L}(1+\delta_{i, rad}),
\end{equation}
where $\cal L$ is the integrated luminosity of the data sample and $(1+\delta_{i, rad})$
is the correction for QED radiation effects in the interval $i$.   
These corrections are calculated as a ratio of RAPGAP predictions
with and without QED radiation simulated.
The differential cross section is determined by dividing $\sigma_{i}$ by the area of the corresponding interval.

\subsection{Systematic uncertainties}

The systematic uncertainties induced by experimental effects
and by model adequateness are propagated to each measurement interval 
in the unfolding procedure (eq.~\ref{eq:unf}). A dedicated detector response matrix is constructed
for each variation related to particular sources of uncertainties:  
\begin{itemize}
\item{The energy of the scattered electron is varied by $\pm 1\%$ with a resulting uncertainty on the integrated 
dijet cross section of $1\%.$}
\item{The polar angle of the scattered electron is varied by $\pm 1\ \rm mrad$ with a resulting uncertainty 
on the integrated dijet cross section of $1\%.$}
\item{The energy of each particle contained in HFS is varied by $\pm 1\%$~\cite{kogler} which translates into an uncertainty 
on the integrated dijet cross section of $4\%$.}
\item{
Uncertainties related to the model dependent corrections of the data
are accounted for by varying the shape of the kinematic distributions in
$\qsq$, $\xpom$, $\beta$, $p^{\ast}_{\rm T,1}$, $\zpom$, $x_{dijet}$ and $\Delta\eta^{\ast}$
in the MC such that the data are still described  within the statistical uncertainties.
For this purpose, the multiplicative weights
$(\log\qsq)^{\pm 0.2}$, $\xpom^{\pm 0.05}$, $\beta^{\pm 0.01}(1-\beta^{\pm 0.01})$, $p^{\ast \pm 0.04}_{\rm T,1}$, $\zpom^{\pm 0.15}$,
$x^{\pm 0.15}_{dijet}$ and $(1.5+\Delta\eta^{\ast})^{\pm 0.5}$ are applied, respectively.
The largest resulting uncertainty of $3\%$ arises from the variation of the shape in $p^{\ast}_{\rm T,1}$.
The shape of the distribution in $t$ is varied within the experimental uncertainty on the $t$-slope~\cite{tslope}
by applying a weight of $e^{\pm t}$ in MC, which translates into an uncertainty on the integrated dijet cross section of $1\,\%$.
The integrated cross section uncertainty due to the model dependence of the measurement is of the order of $5\,\%$.
}
\end{itemize}
 
The following uncertainties on the global normalisation are considered:
\begin{itemize}
\item{The luminosity of the data is measured with a precision $\pm 2.7\ \%$~\cite{qedclumi}.}
\item{The trigger efficiency related to the tracking and  SpaCal condition 
	induces an uncertainty of $2\%$ and $1\%$, respectively.}
\item{The uncertainty accounting for the LRG selection efficiency is $7\%$~\cite{incldiff_lrg}.}
\item The normalisation of the non-diffractive DIS background modelled by RAPGAP is varied by $\pm 50\,\%$
and the normalisation of the diffractive background is varied by $\pm 100\,\%$,
yielding a resulting uncertainty on the integrated dijet cross section below $1\%$ in both cases.

\end{itemize}

The total systematic uncertainty is obtained by adding the individual contributions in quadrature.
%The uncertainties originating from the energy and the polar angle of the scattered electron variation are treated as correlated,
%whereas the uncertainties originating from the HFS energy variation are treated as $50\%$ correlated and $50\%$ uncorrelated.
%The model uncertainties originating from the shape variations are treated as fully correlated.


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

%\subsection{Dijet cross sections}

The integrated cross section in the measurement phase space specified in table~\ref{tab:ps} is found to be
\begin{equation}
\sigma^{dijet}_{meas}(ep\to eXY) = 73\pm{2}\ (\rm stat.)\ \pm{7}\ (\rm syst.)\ \rm pb~.
\end{equation}
The NLO QCD prediction of the total diffractive dijet cross section is
\begin{equation}
\sigma^{dijet}_{theo}(ep\to eXY) = 77\ ^{+25}_{-20}\ (\rm scale)\ ^{+4}_{-14}\ (\rm DPDF)\ \pm{3}\ (\rm had)\ \rm pb~,
\end{equation}
in very good agreement with the measurement.
The uncertainty on the NLO prediction is found to be significantly larger than the experimental uncertainty.

Single differential cross sections are given in tables~\ref{tab:tab1} and~\ref{tab:tab2} and are shown in 
figures~\ref{fig:fig1}-\ref{fig:fig4}. The statistical correlations between measurements
in different bins are given in tables~\ref{tab:tabrho1} and~\ref{tab:tabrho2}.  
The differential cross sections as a function of the DIS variables $\qsq$ and $y$
are shown in figure~\ref{fig:fig1}, as a function of the momentum fractions $\xpom$ and $\zpom$
are shown in figure~\ref{fig:fig2} and as a function of the jet variables $p^{\ast}_{T,1}$, 
$p^{\ast}_{T,2}$, $\langle p^{\ast}_{T}\rangle$ and $\Delta\eta^{\ast}$ 
are shown in Figure~\ref{fig:fig3} and~\ref{fig:fig4}.
For the majority of the measurements, the data precision is limited by systematic effects.
The statistical correlations are small for the inclusive kinematic variables $\qsq$ and $y$
and moderate ($|\rho|< 0.6$) for the other variables.
The figures also include the NLO QCD predictions
which describe within their large uncertainties the data well.

The dynamics of dijet production is further studied in terms of double differential cross sections
in bins of $\zpom$ and of the QCD scale defining observables $\qsq$ and $p^{\ast}_{\rm T,1}$.
The double differential cross sections are listed in tables~\ref{tab:tab3}-\ref{tab:tab5}
and are shown in figures~\ref{fig:fig5}-\ref{fig:fig10}.
The corresponding statistical correlations between measurements
in different bins are given in tables~\ref{tab:tabrho3}-\ref{tab:tabrho5}.
Figure~\ref{fig:fig5} shows the double differential cross section measured in bins of $\zpom$ and $\qsq$.
The ratio of the data to the theory prediction is shown in figure~\ref{fig:fig6}. 
The data are well described by the NLO prediction in most of the phase space. 
The double differential cross section measured in bins of $p^{\ast}_{\rm T,1}$ and $\qsq$ is shown in figure~\ref{fig:fig9}
and the corresponding ratios of the measurements to the NLO predictions are shown in figure~\ref{fig:fig10}.

The present measurement is based on a %$\sim6$ 
six times increased
luminosity as compared to the previous H1 measurement of dijet
production with LRG~\cite{jetdiff_3} and is using a more sophisticated data
correction method. 
A direct comparison of the present data to other measurements of dijet production in diffractive DIS is not possible because of different phase space definitions.
Measurements based on the direct detection of a
forward proton~\cite{jetdiff_2} are limited in statistical precision due
to the restricted geometrical acceptance of the proton taggers.

The experimental uncertainties on both single- and double-differential cross sections
are in general smaller than the theory uncertainties.
The data thus have the power to constrain QCD in diffractive DIS.
%\subsection{Extraction of the strong coupling constant $\alpha_s(M_Z)$}
Here, the double-differential dijet cross sections as a function of $\qsq$ and $p_{\rm T,1}^*$ are used 
to determine the value of the strong coupling constant $\asmz$ at the scale 
of the mass of the $Z$-boson, $M_Z$. % in the framework of pQCD.
The value of $\asmz$ is determined by an iterative $\chi^2$-minimisation 
procedure using NLO calculations, corrected for hadronisation effects
following the method~\cite{HighQ2Paper}.
In the fit, the uncertainties on the HFS energy scale are treated as $50\%$ correlated and $50\%$ uncorrelated. 
All other experimental uncertainties are treated as correlated.
Scale uncertainties, hadronisation uncertainties and  DPDF uncertainties of the NLO calculation
are propagated to the fit result as described in~\cite{HighQ2Paper}.

The fit yields a value of $\chi^2/n_{\rm dof} = 16.7/14$, with $n_{\rm dof}$ being the number of degrees of freedom,
thus indicating good agreement 
of theory to data. The nuisance parameters of the correlated systematic uncertainties 
are equally distributed around zero with absolute values below one.
The value of $\asmz$ determined in the fit to the dijet cross sections is
\begin{flalign}
 \asmz &= 0.119\pm0.004\,({\rm exp})\pm0.002\,({\rm had})\pm0.005\,({\rm DPDF})\pm0.010\,({\mu_r})\pm0.004\,({\mu_f}) \\
       & =0.119\pm0.004\,({\rm exp})\pm0.012\,({\rm DPDF,theo}) \nonumber
\end{flalign}
The largest uncertainties arise from the estimate of the contributions from  orders beyond NLO
and from the poor knowledge of the DPDF.
The largest contribution to the experimental uncertainty of $0.003$
arises from the global normalisation uncertainty.

The result for $\asmz$ is consistent within the uncertainties with the world average% value of $\asmz=0.1185\pm0.0006$
~\cite{Beringer:1900zz,Bethke:2012jm} and with values from other jet data 
in DIS and photoproduction~\cite{Aaron:10:1,HighQ2Paper,Abramowicz:2012jz} as well as values 
of $\asmz$ determined from jet data at 
the Tevatron~\cite{Abazov:2009nc,Abazov:2012lua} and at the LHC~\cite{Chatrchyan:2013txa,alps_cms}.
Although the uncertainty of this $\asmz$ extraction is not competitive
with measurements in other processes the agreement with the other
measurements supports the underlying concept of treating dijet
production in diffractive DIS with perturbative QCD calculations. 
%---------------------
%---------------------
\section{Conclusions}

Integrated, single- and double-differential cross sections of diffractive DIS dijet production
are measured with the H1 experiment in $ep$ collisions at HERA
and compared with NLO QCD predictions.%theory predictions calculated in pQCD up to NLO.

The integrated diffractive dijet cross section is found to be well
described by the NLO QCD prediction using the H12006 Fit-B DPDF set.
Both shapes and normalisation of the single-differential cross sections are 
reproduced by the theory within the experimental and theory uncertainties,
confirming at improved precision the conclusions from previous H1 measurements. %~\cite{jetdiff_1,jetdiff_2},
Good agreement of the theory with the measurement is also
found for the shapes and normalisation of the double differential cross sections.
The cross section measurements presented here show 
experimental uncertainties significantly smaller than the uncertainties of the
theory predictions. 
From a fit of the NLO prediction to the double differential cross sections in $\qsq$ and $p^{*}_{T,1}$,  
the strong coupling constant has been determined to be $\alpha_{s}(M_{Z})=0.119\ (4)_{\rm exp}\ (12)_{\rm theo}$.
%Although the precision of $\alpha_s$ is not competitive with measurements from other processes, the agreement with other measurements supports the underlying concept of treating diffractive jet production in QCD.

%---------------------
%---------------------
\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
the hospitality which they extend to
the non-DESY members of the collaboration.
We would like to give credit to all partners contributing to the EGI computing infrastructure for their support for the H1 Collaboration.

\begin{flushleft}
%\bibliographystyle{is-unsrt} %The style you want to use for references.
%\bibliography{cite} %The files containing all the articles and books
                    %you ever referenced.
                    %you ever referenced.
\begin{thebibliography}{10}
\ifx \showCODEN  \undefined \def \showCODEN #1{CODEN #1}  \fi
\ifx \showISBN   \undefined \def \showISBN  #1{ISBN #1}   \fi
\ifx \showISSN   \undefined \def \showISSN  #1{ISSN #1}   \fi
\ifx \showLCCN   \undefined \def \showLCCN  #1{LCCN #1}   \fi
\ifx \showPRICE  \undefined \def \showPRICE #1{#1}        \fi
\ifx \showURL    \undefined \def \showURL {URL }          \fi
\ifx \path       \undefined \input path.sty               \fi
\ifx \ifshowURL \undefined
     \newif \ifshowURL
     \showURLtrue
\fi

\bibitem{pomeranchuk}
E.~Feinberg and I.~Pomeranchuk,
\newblock Suppl. Nuovo Cimento {\bf3 } (1956)  652.

\bibitem{goulianos}
K.~Goulianos,
\newblock Phys. Rep. {\bf 101} (1983)  171.

\bibitem{trentadue}
L.~Trentadue and G.~Veneziano,
\newblock Phys. Lett. B {\bf 323} (1994)  201.

\bibitem{fitb}
A.~Aktas et~al. [H1~Collaboration],
\newblock Eur. Phys. J. C {\bf 48} (2006)  715, [hep--ex/0606004].

\bibitem{gribov}
V.~N. Gribov and L.~N. Lipatov,
\newblock Sov. J. Nucl. Phys. {\bf 15} (1972)  438.

\bibitem{dokschitzer}
Y.~L. Dokshitzer,
\newblock Sov. Phys. JETP {\bf 46} (1977)  641.

\bibitem{altarelli}
G.~Altarelli and G.~Parisi,
\newblock Nucl. Phys. B {\bf 126} (1977)  298.

\bibitem{collins}
J.~C. Collins,
\newblock Phys. Rev. D {\bf 57} (1998)  3051 [Erratum--ibid. D 61 (2000)
  019902] [hep--ph/9709499].

\bibitem{rfact_1}
G.~Ingelman and P.~Schlein,
\newblock Phys. Lett. B {\bf 152} (1985)  256.

\bibitem{rfact_2}
A.~Donnachie and P.~Landshoff,
\newblock Phys. Lett. B {\bf 191} (1987) 309.

\bibitem{rfact_3}
S.~Chekanov et~al. [ZEUS~Collaboration],
\newblock Eur. Phys. J. C {\bf 38} (2004)  43, [hep--ex/0408009].

\bibitem{jetdiff_1}
C.~Adloff et~al. [H1~Collaboration],
\newblock Eur. Phys. J. C {\bf 20} (2001)  29, [hep--ex/0012051].

\bibitem{jetdiff_3}
A.~Aktas et~al. [H1~Collaboration],
\newblock JHEP 0710 (2007)  42, [arxiv:0708.3217].

\bibitem{jetdiff_zeus}
S.~Chekanov et~al. [ZEUS~Collaboration],
\newblock Eur. Phys. J. C {\bf 52} (2007)  813.

\bibitem{jetdiff_2}
F.~D.~Aaron et~al. [H1~Collaboration],
\newblock Eur. Phys. J. C {\bf 72} (2011)  1970, [arXiv:1111.0584].

\bibitem{rapgap}
H.~Jung,
\newblock Comput. Phys. Commun. {\bf 86} (1995)  147.

\bibitem{llog}
M.~Bengtsson and T.~Sj\"{o}strand,
\newblock Z. Phys. C {\bf 37} (1988)  465.

\bibitem{lund}
B.~Andersson et~al.,
\newblock Phys. Rep. {\bf 97} (1983) ~31.

\bibitem{pythia}
T.~Sj\"{o}strand,
\newblock Comput. Phys. Commun. {\bf 82} (1994) ~72.

\bibitem{diffvm}
B.~List and A.~Mastroberardino,
\newblock In {\em Proceedings of the Monte Carlo Generators for HERA Physics
  Workshop 1998/99}. Eds. A.T. Doyle, G. Grindhammer, G. Ingelman, H. Jung,
  Hamburg, Germany (1999)  396, [DESY--PROC--1999--02].

\bibitem{resgamma}
C.~Adloff et.~al [H1~Collaboration],
\newblock Eur. Phys. J. C {\bf 13} (2000)  397, [hep--ex/9812024].

\bibitem{sasg2dfit}
G.~A. Schuller and T.~Sj\"{o}strand,
\newblock Phys. lett. B {\bf 376} (1996)  193.

\bibitem{heracles}
A.~Kwiatkowski et~al.,
\newblock Comput. Phys. Commun. {\bf 69} (1992)  155.

\bibitem{cteq}
J.~Pumplin et~al.,
\newblock JHEP {\bf 207} (2002)  12, [hep--ph/0201195].

\bibitem{geant}
R.~Brun et~al.,
\newblock (1987), [CERN-DD/EE-84-1].

\bibitem{nlojetprog}
Z.~Nagy and Z.~Trocsanyi,
\newblock Phys. Rev. D {\bf 59} (1999)  014020, [arXiv:hep--ph/9806317],
  Erratum--ibid. D 62 (2000) 099902.

\bibitem{nlojet}
Z.~Nagy and Z.~Trocsanyi,
\newblock Phys. Rev. Lett. {\bf 87} (2001)  082001, [hep--ph/0104315].

\bibitem{dpdferr}
W.~Giele et~al.,
\newblock In {\em Proceedings of the Workshop on Physics at TeV Colliders
  (2001)}. Les Houches, France (2002)  275, [hep--ph/0204316].

\bibitem{H1fastNLO}
C.~Adloff et~al. [H1~Collaboration],
\newblock Eur. Phys. J. C {\bf 19} (2001)  289, [arXiv:hep--ex/0010054].

\bibitem{fastNLO}
T.~Kluge et~al.,
\newblock In {\em {Proceedings of the 14th International Workshop on Deep
  Inelastic Scattering (DIS 2006)}}. Eds. M. Kuze, K. Nagano and K. Tokushuku,
  Tsukuba, Japan (2007)  483, [hep--ph/0609285].

\bibitem{fnlo}
D.~Britzger et~al.,
\newblock In {\em Proceedings of the 20th International Workshop on
  Deep-Inelastic Scattering and Related Subjects (DIS 2012)}. Ed. I. Brock,
  Bonn, Germany, (2012)  217, [arXiv:1208.3641].

\bibitem{ariadne}
L.~L\"{o}nnblad,
\newblock Comput. Phys. Commun. {\bf 71} (1992) ~15.

\bibitem{det_1}
I.~Abt et~al. [H1~Collaboration],
\newblock Nucl. Instrum. Meth. A {\bf 386} (1997)  310.

\bibitem{lar}
B.~Andrieu et~al. [H1 Calorimeter~Group],
\newblock Nucl. Instrum. Meth. A {\bf 336} (1993)  460.

\bibitem{testbeam_1}
B.~Andrieu et~al. [H1 Calorimeter~Group],
\newblock Nucl. Instrum. Meth. A {\bf 350} (1994) ~57.

\bibitem{testbeam_2}
B.~Andrieu et~al. [H1 Calorimeter~Group],
\newblock Nucl. Instrum. Meth. A {\bf 336} (1993)  499.

\bibitem{spac1}
R.~Appuhn et~al.[H1 SPACAL~Group],
\newblock Nucl. Instrum. Meth. A {\bf 386} (1997)  397.

\bibitem{spac2}
R.~Appuhn et~al.[H1 SPACAL~Group],
\newblock Nucl. Instrum. Meth. A {\bf 382} (1996)  395.

\bibitem{qedclumi}
F.~D.~Aaron et~al. [H1~Collaboration],
\newblock Eur. Phys. J. C {\bf 72} (2012)  2163, erratum ibid. C {\bf 74}
  (2014) 2733, [arXiv:1205.2448].

\bibitem{esigma}
C.~Adloff et~al. [H1~Collaboration],
\newblock Z. Phys. C {\bf 76} (1997)  613, [hep--ex/9708016].

\bibitem{peez:03}
M.~Peez,
\newblock {\em {Search for deviations from the standard model in high
  transverse energy processes at the electron proton collider HERA (in
  French)}},
\newblock Dissertation, Univ. Claude Bernard, Lyon, DESY-THESIS-2003-023,
  CPPM-T-2003-04 (2003), Available at
  {http://www-h1.desy.de/psfiles/theses/h1th-317.pdf}.

\bibitem{hellwig:04}
S.~Hellwig,
\newblock {\em {Investigation of the $D^{*} - \pi_{slow}$ double tagging method
  for the analysis of charm (in German)}},
\newblock Diploma thesis, Hamburg University (2004), Available at
  {http://www-h1.desy.de/psfiles/theses/h1th-341.pdf}.

\bibitem{kogler}
R.~Kogler,
\newblock {\em {Measurement of jet production in deep-inelastic ep scattering
  at HERA}},
\newblock Dissertation, Universit\"at Hamburg, DESY-THESIS-2011-003,
  MPP-2010-175 (2010), Available at
  {http://www-h1.desy.de/psfiles/theses/h1th-590.pdf}.

\bibitem{kt}
Y.~Dokshitzer S.~Catani and B.~Weber,
\newblock Phys. Lett. B {\bf 285} (1992)  291.

\bibitem{fastjet}
M.~Cacciari and G.~P. Salam,
\newblock Phys. Lett. B {\bf 641} (2006)  57, [hep--ph/0512210].

\bibitem{mythesis}
B.~Pokorny,
\newblock {\em {Measurement of Diffractive Dijet Production in Deep Inelastic
  Scattering at HERA Collider}},
\newblock Dissertation, Charles University in Prague (2014) Available at
  {http://www-h1.desy.de/psfiles/theses/h1th-850.pdf}.

\bibitem{unf_3}
S.~Schmitt,
\newblock JINST {\bf 7} (2012), [arXiv:1205.6201] T10003.

\bibitem{tslope}
F.~D.~Aaron et~al. [H1~Collaboration],
\newblock Eur. Phys. J. C {\bf 71} (2011)  1578, [arXiv:1010.1476].

\bibitem{incldiff_lrg}
C.~Adloff et~al. [H1~Collaboration],
\newblock Eur. Phys. J. C {\bf 72} (2007)  2074, [arXiv:1203.4495].

\bibitem{HighQ2Paper}
V.~Andreev et~al. [H1~Collaboration],
\newblock Eur. Phys. J. C (2014)  [arXiv:1406.4709].

\bibitem{Beringer:1900zz}
J.~Beringer et~al.,
\newblock Phys. Rev. D {\bf 86} (2012)  010001,
\newblock and 2013 partial update for the 2014 edition.

\bibitem{Bethke:2012jm}
S.~Bethke,
\newblock Nucl. Phys. Proc. Suppl. {\bf 234} (2013)  229, [arXiv:1210.0325].

\bibitem{Aaron:10:1}
F.~D.~Aaron et~al. [H1~Collaboration],
\newblock Eur. Phys. J. C {\bf67} (2010)  1, [arXiv:0911.5678].

\bibitem{Abramowicz:2012jz}
H.~Abramowicz et~al. [ZEUS~Collaboration],
\newblock Nucl. Phys. B {\bf 864} (2012)  1, [arXiv:1205.6153].

\bibitem{Abazov:2009nc}
V.~M. Abazov et~al.,
\newblock Phys. Rev. D {\bf 80} (2009)  111107, [arXiv:0911.2710].

\bibitem{Abazov:2012lua}
V.~M. Abazov et~al.,
\newblock Phys. Lett. B {\bf{718}} (2012)  56, [arXiv:1207.4957].

\bibitem{Chatrchyan:2013txa}
S.~Chatrchyan et~al. [CMS~Collaboration],
\newblock Eur. Phys. J. C {\bf 73} (2013)  2604, [arXiv:1304.7498].

\bibitem{alps_cms}
V.~Khachatryan et~al. [CMS~Collaboration],
\newblock Eur. Phys. J. C (2014), [arXiv:1410.6765].

\end{thebibliography}
\end{flushleft}

\newpage
\clearpage
\begin{landscape}
\begin{table}
\resizebox{.8\textwidth}{!}{\begin{minipage}{\textwidth}
\begin{tabular}{|c|c|c|c|c||c|c|c|c|c|c|c|c|c|c|c|c||c||c|}
\hline
$\qsq$ & $d\sigma/d\qsq$ & $\delta_{tot}$ & $\delta_{stat}$ & $\delta_{sys}$ & $\delta_{\theta}$ & $\delta_{E}$ & $\delta_{HFS}$
        & $\delta_{\qsq}$ & $\delta_{\xpom}$ & $\delta_{\beta}$ & $\delta_{p*_{T,1}}$ 
        & $\delta_{\zpom}$ & $\delta_{x_{dijet}}$ & $\delta_{\Delta\eta*}$ & $\delta_{t}$
        & $\delta_{bgr}$ & $1+\delta_{had}$ & $1+\delta_{rad}$\\
$\rm [GeV^2]$ & $\rm [pb/GeV^2]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &
        $[\%]$ &  $[\%]$ &  $[\%]$ & $[\%]$ & & \\
\hline
$4 \div 6$ &$8.20$ &$13.2$ & $5.7$ & $11.9$ &$1.0$ & $4.5$ & $-3.9$ &$2.1$ & $1.1$ &$2.8$ & $-5.0$ &$1.9$ & $0.7$ &$1.2$ & $-0.9$ &$0.1$ &$1.05 \pm 0.05$ &$1.05$ \\
$6 \div 10$ &$4.23$ &$11.8$ & $4.0$ & $11.0$ &$2.6$ & $1.7$ & $-5.0$ &$-0.5$ & $0.5$ &$-0.5$ & $-3.1$ &$-3.2$ & $-1.6$ &$1.6$ & $-1.5$ &$0.3$ &$1.05 \pm 0.04$ &$1.03$ \\
$10 \div 18$ &$1.92$ &$11.4$ & $4.0$ & $10.7$ &$1.0$ & $1.9$ & $-4.6$ &$-0.9$ & $0.6$ &$-0.8$ & $-3.1$ &$-3.1$ & $-1.5$ &$1.7$ & $-0.9$ &$0.4$ &$1.05 \pm 0.04$ &$1.03$ \\
$18 \div 34$ &$0.797$ &$11.6$ & $4.8$ & $10.5$ &$1.1$ & $2.1$ & $-5.1$ &$0.1$ & $0.6$ &$-0.1$ & $-2.9$ &$-2.5$ & $-1.3$ &$1.4$ & $-0.6$ &$0.2$ &$1.06 \pm 0.04$ &$1.03$ \\
$34 \div 100$ &$0.164$ &$12.3$ & $6.2$ & $10.6$ &$0.9$ & $2.3$ & $-5.0$ &$-0.2$ & $0.5$ &$-0.6$ & $-2.7$ &$-2.9$ & $-1.5$ &$1.6$ & $-0.8$ &$0.1$ &$1.06 \pm 0.04$ &$1.03$ \\
\hline
$y$ & $d\sigma/dy$ & $\delta_{tot}$ & $\delta_{stat}$ & $\delta_{sys}$ & $\delta_{\theta}$ & $\delta_{E}$ & $\delta_{HFS}$
        & $\delta_{\qsq}$ & $\delta_{\xpom}$ & $\delta_{\beta}$ & $\delta_{p*_{T,1}}$
        & $\delta_{\zpom}$ & $\delta_{x_{dijet}}$ & $\delta_{\Delta\eta*}$ & $\delta_{t}$
        & $\delta_{bgr}$ & $1+\delta_{had}$ & $1+\delta_{rad}$\\
& $\rm [pb]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &
        $[\%]$ &  $[\%]$ &  $[\%]$ & $[\%]$ & & \\
\hline
$0.10 \div 0.22$ &$113$ &$18.4$ & $6.5$ & $17.2$ &$2.1$ & $0.2$ & $-8.7$ &$-3.6$ & $-0.2$ &$-4.2$ & $-3.5$ &$-8.9$ & $-3.6$ &$3.9$ & $-1.5$ &$0.6$ &$1.01 \pm 0.06$ &$1.07$ \\
$0.22 \div 0.34$ &$163$ &$12.7$ & $4.5$ & $11.9$ &$2.0$ & $1.1$ & $-5.9$ &$-2.0$ & $0.5$ &$-1.5$ & $-3.2$ &$-4.1$ & $-2.1$ &$1.4$ & $-0.9$ &$0.6$ &$1.02 \pm 0.04$ &$1.05$ \\
$0.34 \div 0.46$ &$144$ &$11.2$ & $4.3$ & $10.4$ &$1.6$ & $2.8$ & $-4.2$ &$-0.4$ & $0.8$ &$-0.1$ & $-3.1$ &$-2.3$ & $-1.3$ &$1.0$ & $-1.1$ &$0.3$ &$1.06 \pm 0.04$ &$1.04$ \\
$0.46 \div 0.58$ &$106$ &$11.2$ & $5.0$ & $10.0$ &$1.2$ & $3.2$ & $-3.2$ &$0.7$ & $0.8$ &$0.9$ & $-3.1$ &$-1.0$ & $-0.6$ &$1.9$ & $-0.3$ &$0.4$ &$1.13 \pm 0.03$ &$1.02$ \\
$0.58 \div 0.70$ &$76.5$ &$12.4$ & $7.0$ & $10.2$ &$0.7$ & $4.3$ & $-2.3$ &$1.0$ & $0.6$ &$1.6$ & $-3.3$ &$0.3$ & $0.4$ &$1.2$ & $-1.5$ &$0.2$ &$1.17 \pm 0.02$ &$0.97$ \\
\hline
$\xpom$ & $d\sigma/d\xpom$ & $\delta_{tot}$ & $\delta_{stat}$ & $\delta_{sys}$ & $\delta_{\theta}$ & $\delta_{E}$ & $\delta_{HFS}$
        & $\delta_{\qsq}$ & $\delta_{\xpom}$ & $\delta_{\beta}$ & $\delta_{p*_{T,1}}$
        & $\delta_{\zpom}$ & $\delta_{x_{dijet}}$ & $\delta_{\Delta\eta*}$ & $\delta_{t}$
        & $\delta_{bgr}$ & $1+\delta_{had}$ & $1+\delta_{rad}$\\
& $\rm [pb]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &
        $[\%]$ &  $[\%]$ &  $[\%]$ & $[\%]$ & & \\
\hline
$-2.30 \div -2.10$ &$14.2$ &$42.0$ & $36.2$ & $21.1$ &$1.8$ & $3.9$ & $-9.3$ &$-2.7$ & $4.0$ &$-3.7$ & $-5.8$ &$-11.4$ & $-4.7$ &$7.6$ & $-1.2$ &$0.7$ &$1.17 \pm 0.13$ &$1.06$ \\
$-2.10 \div -1.90$ &$53.5$ &$14.7$ & $8.9$ & $11.7$ &$1.6$ & $2.4$ & $-5.6$ &$-0.6$ & $1.2$ &$-0.8$ & $-3.2$ &$-3.7$ & $-1.8$ &$2.3$ & $-1.4$ &$0.0$ &$1.10 \pm 0.08$ &$1.04$ \\
$-1.90 \div -1.70$ &$111$ &$11.6$ & $5.5$ & $10.2$ &$1.5$ & $1.3$ & $-4.5$ &$-1.1$ & $0.1$ &$-0.2$ & $-3.5$ &$-1.5$ & $-1.0$ &$1.4$ & $-1.1$ &$0.0$ &$1.06 \pm 0.04$ &$1.04$ \\
$-1.70 \div -1.52$ &$196$ &$10.9$ & $4.9$ & $9.8$ &$1.3$ & $2.5$ & $-3.6$ &$-1.0$ & $-0.5$ &$0.5$ & $-3.4$ &$-0.0$ & $-0.3$ &$-0.5$ & $-0.5$ &$0.8$ &$1.03 \pm 0.03$ &$1.03$ \\
\hline
$\zpom$ & $d\sigma/d\zpom$ & $\delta_{tot}$ & $\delta_{stat}$ & $\delta_{sys}$ & $\delta_{\theta}$ & $\delta_{E}$ & $\delta_{HFS}$
        & $\delta_{\qsq}$ & $\delta_{\xpom}$ & $\delta_{\beta}$ & $\delta_{p*_{T,1}}$
        & $\delta_{\zpom}$ & $\delta_{x_{dijet}}$ & $\delta_{\Delta\eta*}$ & $\delta_{t}$
        & $\delta_{bgr}$ & $1+\delta_{had}$ & $1+\delta_{rad}$\\
& $\rm [pb]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &
        $[\%]$ &  $[\%]$ &  $[\%]$ & $[\%]$ & & \\
\hline
$0.00 \div 0.22$ &$70.4$ &$20.3$ & $9.3$ & $18.0$ &$1.4$ & $3.4$ & $-4.0$ &$1.2$ & $-0.5$ &$4.6$ & $-2.4$ &$12.0$ & $4.1$ &$6.5$ & $-0.5$ &$0.8$ &$1.10 \pm 0.03$ &$1.06$ \\
$0.22 \div 0.40$ &$132$ &$11.9$ & $6.3$ & $10.1$ &$1.5$ & $3.0$ & $-1.2$ &$-0.9$ & $0.3$ &$-0.1$ & $-3.9$ &$-2.3$ & $-2.2$ &$1.0$ & $-0.9$ &$0.4$ &$1.07 \pm 0.02$ &$1.04$ \\
$0.40 \div 0.60$ &$89.7$ &$14.9$ & $6.8$ & $13.3$ &$1.2$ & $1.6$ & $-9.1$ &$-1.3$ & $0.8$ &$-1.2$ & $-2.8$ &$-3.9$ & $-1.4$ &$0.5$ & $-0.6$ &$0.3$ &$1.10 \pm 0.03$ &$1.02$ \\
$0.60 \div 0.80$ &$54.8$ &$14.9$ & $7.5$ & $12.9$ &$2.5$ & $1.9$ & $7.6$ &$-1.4$ & $0.9$ &$-1.2$ & $-3.2$ &$-4.2$ & $-1.4$ &$0.2$ & $-2.0$ &$0.1$ &$1.10 \pm 0.10$ &$1.02$ \\
$0.80 \div 1.00$ &$19.9$ &$45.0$ & $11.4$ & $43.5$ &$0.8$ & $0.6$ & $-42.1$ &$-1.9$ & $1.3$ &$-2.4$ & $-2.5$ &$-5.1$ & $-2.0$ &$3.0$ & $-1.5$ &$0.6$ &$0.57 \pm 0.10$ &$1.00$ \\
\hline
\end{tabular}
\end{minipage}}
\caption{Diffractive DIS dijet cross section measured differentially as a function of $\qsq$, $y$, $\log\xpom$ and $\zpom$.
        The statistical $\delta_{stat}$ and systematic $\delta_{sys}$ uncertainties are given together with the total uncertainty $\delta_{tot}$. The next 12 columns represent $+1\sigma$ shifts for the systematic error contributions from:
%        Uncertainties corresponding to individual sources are also listed: 
        electron polar angle measurement $\delta_{\theta}$,
        electron energy scale $\delta_{E}$, HFS energy scale $\delta_{HFS}$,
        model uncertainties $\delta_{\qsq}$,  $\delta_{\xpom}$, $\delta_{\beta}$, $\delta_{p*_{T,1}}$,
	$\delta_{\zpom}$, $\delta_{x_{dijet}}$, $\delta_{\Delta\eta*}$ and $\delta_{t}$
        and the background normalisation uncertainty $\delta_{bgr}$.
%	A sign of each uncertainty corresponds to the positive variation of the source of the uncertainty.
	The global normalisation uncertainty of $7.8\%$ is not listed explicitly but is included in the total 
%        The 
	systematic uncertainty $\delta_{sys}$.
%includes also the global normalisation uncertainties of $7.8\%$.
        The %next-to-last and 
        last two column show the correction factors for hadronisation and QED radiation, respectively.
        }
\label{tab:tab1}
\end{table}
\end{landscape}


\newpage
\clearpage
\begin{landscape}
\begin{table}
\resizebox{.8\textwidth}{!}{\begin{minipage}{\textwidth}
\begin{tabular}{|c|c|c|c|c||c|c|c|c|c|c|c|c|c|c|c|c||c||c|}
\hline
$p^{\ast}_{T,1}$ & $d\sigma/dp^{\ast}_{T,1}$ & $\delta_{tot}$ & $\delta_{stat}$ & $\delta_{sys}$ & $\delta_{\theta}$ & $\delta_{E}$ & $\delta_{HFS}$
        & $\delta_{\qsq}$ & $\delta_{\xpom}$ & $\delta_{\beta}$ & $\delta_{p*_{T,1}}$
        & $\delta_{\zpom}$ & $\delta_{x_{dijet}}$ & $\delta_{\Delta\eta*}$ & $\delta_{t}$
        & $\delta_{bgr}$ & $1+\delta_{had}$ & $1+\delta_{rad}$\\
$[GeV]$ & $\rm [pb/GeV]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &
        $[\%]$ &  $[\%]$ &  $[\%]$ & $[\%]$ & & \\
\hline
$5.50 \div 7.00$ &$30.8$ &$9.6$ & $3.2$ & $9.0$ &$1.4$ & $1.4$ & $-3.1$ &$-0.9$ & $0.5$ &$-0.6$ & $-1.3$ &$-1.7$ & $-0.9$ &$0.6$ & $-0.8$ &$0.1$ &$1.05 \pm 0.05$ &$1.03$ \\
$7.00 \div 9.00$ &$10.5$ &$11.8$ & $6.1$ & $10.0$ &$1.3$ & $3.0$ & $-4.6$ &$-0.6$ & $0.8$ &$-0.4$ & $-1.4$ &$-1.5$ & $-0.9$ &$1.2$ & $-1.0$ &$0.7$ &$1.06 \pm 0.04$ &$1.04$ \\
$9.00 \div 15.00$ &$1.07$ &$19.6$ & $12.7$ & $14.9$ &$1.3$ & $2.3$ & $-9.8$ &$-0.1$ & $1.0$ &$-0.1$ & $-4.2$ &$-2.7$ & $-1.1$ &$5.3$ & $-1.5$ &$0.7$ &$1.04 \pm 0.03$ &$1.06$ \\
\hline
$p^{\ast}_{T,2}$ & $d\sigma/dp^{\ast}_{T,2}$ & $\delta_{tot}$ & $\delta_{stat}$ & $\delta_{sys}$ & $\delta_{\theta}$ & $\delta_{E}$ & $\delta_{HFS}$
        & $\delta_{\qsq}$ & $\delta_{\xpom}$ & $\delta_{\beta}$ & $\delta_{p*_{T,1}}$
        & $\delta_{\zpom}$ & $\delta_{x_{dijet}}$ & $\delta_{\Delta\eta*}$ & $\delta_{t}$
        & $\delta_{bgr}$ & $1+\delta_{had}$ & $1+\delta_{rad}$\\
$[GeV]$ & $\rm [pb/GeV]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &
        $[\%]$ &  $[\%]$ &  $[\%]$ & $[\%]$ & & \\
\hline
$4.00 \div 6.50$ &$22.3$ &$10.4$ & $3.7$ & $9.7$ &$1.5$ & $2.3$ & $-3.8$ &$-1.0$ & $0.7$ &$-0.8$ & $-1.2$ &$-2.0$ & $-1.2$ &$1.2$ & $-0.8$ &$0.1$ &$1.10 \pm 0.06$ &$1.03$ \\
$6.50 \div 9.00$ &$5.67$ &$12.2$ & $6.9$ & $10.1$ &$1.2$ & $2.0$ & $-4.9$ &$-0.6$ & $0.6$ &$-0.2$ & $-2.6$ &$-1.3$ & $-0.5$ &$1.2$ & $-1.0$ &$0.6$ &$0.97 \pm 0.02$ &$1.04$ \\
$9.00 \div 15.00$ &$0.539$ &$18.2$ & $12.3$ & $13.4$ &$1.1$ & $1.5$ & $-7.8$ &$0.5$ & $0.6$ &$0.8$ & $-6.2$ &$-2.4$ & $-0.9$ &$2.8$ & $-1.3$ &$0.1$ &$0.97 \pm 0.02$ &$1.06$ \\
\hline
$\langle p^{\ast}_{T}\rangle$ & $d\sigma/d\langle p^{\ast}_{T}\rangle$ & $\delta_{tot}$ & $\delta_{stat}$ & $\delta_{sys}$ & $\delta_{\theta}$ & $\delta_{E}$ & $\delta_{HFS}$
        & $\delta_{\qsq}$ & $\delta_{\xpom}$ & $\delta_{\beta}$ & $\delta_{p*_{T,1}}$
        & $\delta_{\zpom}$ & $\delta_{x_{dijet}}$ & $\delta_{\Delta\eta*}$ & $\delta_{t}$
        & $\delta_{bgr}$ & $1+\delta_{had}$ & $1+\delta_{rad}$\\
$[GeV]$ & $\rm [pb/GeV]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &
        $[\%]$ &  $[\%]$ &  $[\%]$ & $[\%]$ & & \\
\hline
$4.75 \div 6.50$ &$27.6$ &$9.9$ & $3.5$ & $9.3$ &$1.5$ & $2.0$ & $-3.3$ &$-1.1$ & $0.5$ &$-0.8$ & $-1.0$ &$-1.9$ & $-1.0$ &$0.8$ & $-0.8$ &$0.1$ &$1.09 \pm 0.06$ &$1.03$ \\
$6.50 \div 9.00$ &$8.52$ &$11.3$ & $5.2$ & $10.0$ &$1.4$ & $2.4$ & $-5.0$ &$-0.4$ & $0.8$ &$-0.1$ & $-1.7$ &$-0.8$ & $-0.4$ &$1.5$ & $-1.1$ &$0.5$ &$1.01 \pm 0.03$ &$1.04$ \\ 
$9.00 \div 15.00$ &$0.701$ &$19.7$ & $13.4$ & $14.4$ &$0.7$ & $1.2$ & $-9.2$ &$-0.3$ & $0.7$ &$-0.2$ & $-5.4$ &$-3.5$ & $-1.3$ &$3.8$ & $-0.9$ &$0.5$ &$1.01 \pm 0.03$ &$1.06$ \\
\hline
$\Delta\eta^{\ast}$ & $d\sigma/d\Delta\eta^{\ast}$ & $\delta_{tot}$ & $\delta_{stat}$ & $\delta_{sys}$ & $\delta_{\theta}$ & $\delta_{E}$ & $\delta_{HFS}$
        & $\delta_{\qsq}$ & $\delta_{\xpom}$ & $\delta_{\beta}$ & $\delta_{p*_{T,1}}$
        & $\delta_{\zpom}$ & $\delta_{x_{dijet}}$ & $\delta_{\Delta\eta*}$ & $\delta_{t}$
        & $\delta_{bgr}$ & $1+\delta_{had}$ & $1+\delta_{rad}$\\
& $\rm [pb]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &  $[\%]$ &
        $[\%]$ &  $[\%]$ &  $[\%]$ & $[\%]$ & & \\
\hline
$0.00 \div 0.15$ &$51.6$ &$17.9$ & $9.5$ & $15.1$ &$1.6$ & $2.8$ & $-4.4$ &$-1.0$ & $1.0$ &$-0.8$ & $-3.8$ &$-2.3$ & $-1.4$ &$10.6$ & $-1.2$ &$0.1$ &$1.04 \pm 0.03$ &$1.03$ \\
$0.15 \div 0.40$ &$57.8$ &$14.1$ & $7.3$ & $12.1$ &$1.2$ & $1.0$ & $-5.1$ &$-0.9$ & $0.9$ &$-0.7$ & $-3.0$ &$-2.1$ & $-1.5$ &$6.2$ & $-1.1$ &$0.2$ &$1.05 \pm 0.03$ &$1.04$ \\
$0.40 \div 0.80$ &$45.1$ &$12.5$ & $5.7$ & $11.1$ &$1.9$ & $2.5$ & $-4.3$ &$-0.9$ & $0.8$ &$-0.5$ & $-3.8$ &$-2.2$ & $-1.2$ &$3.2$ & $-1.3$ &$0.5$ &$1.06 \pm 0.04$ &$1.04$ \\
$0.80 \div 1.30$ &$33.9$ &$12.3$ & $5.5$ & $10.9$ &$1.7$ & $2.4$ & $-4.7$ &$-1.0$ & $0.5$ &$-0.3$ & $-3.7$ &$-2.4$ & $-1.0$ &$-2.5$ & $-0.6$ &$0.3$ &$1.07 \pm 0.05$ &$1.03$ \\
$1.30 \div 3.00$ &$9.29$ &$15.0$ & $6.7$ & $13.4$ &$1.2$ & $3.4$ & $-5.3$ &$-1.0$ & $0.2$ &$-0.0$ & $-2.8$ &$-3.4$ & $-1.2$ &$-7.4$ & $-1.1$ &$0.3$ &$1.04 \pm 0.06$ &$1.03$ \\
\hline
\end{tabular}
\end{minipage}}
\caption{Diffractive DIS dijet cross section measured differentially as a function of $p^{\ast}_{T,1}$, $p^{\ast}_{T,2}$, $\langle p^{\ast}_{T}\rangle$ 
	and $\Delta\eta^{\ast}$.
        The statistical $\delta_{stat}$ and systematic $\delta_{sys}$ uncertainties are given together with the total uncertainty $\delta_{tot}$.
	Further details are given in Table~\ref{tab:tab1}.
        }
\label{tab:tab2}
\end{table}
\end{landscape}



\newpage
\clearpage
\begin{landscape}
\begin{table}[!ht]
\resizebox{.8\textwidth}{!}{\begin{minipage}{\textwidth}
\begin{tabular}{|c|c|c|c|c|c||c|c|c|c|c|c|c|c|c|c|c|c||c||c|}
\hline
$\zpom$ & $\qsq$ & $\frac{d^{2}\sigma}{d\zpom d\qsq}$ & $\delta_{tot}$ & $\delta_{stat}$ & $\delta_{sys}$ & $\delta_{\theta}$ & $\delta_{E}$ & $\delta_{HFS}$
& $\delta_{\qsq}$ & $\delta_{\xpom}$ & $\delta_{\beta}$ & $\delta_{p*_{T,1}}$
& $\delta_{\zpom}$ & $\delta_{x_{dijet}}$ & $\delta_{\Delta\eta*}$ & $\delta_{t}$
& $\delta_{bgr}$ & $1+\delta_{had}$ & $1+\delta_{rad}$\\
& $\rm [GeV^{2}]$ & $\rm [pb/GeV^{2}]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$
& $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & & \\
\hline
$0.0 \div 0.3$ &$4 \div 10$ &$7.67$ &$14.8$ & $7.7$ & $12.7$ &$2.0$ & $3.5$ & $-1.0$ &$1.6$ & $-0.3$ &$3.4$ & $-3.3$ &$6.5$ & $2.2$ &$3.1$ & $-1.3$ &$0.5$ &$1.08 \pm 0.03$ &$1.05$ \\
 & $10 \div 20$ &$2.40$ &$15.6$ & $10.0$ & $12.1$ &$0.7$ & $2.9$ & $-2.1$ &$-2.0$ & $-0.2$ &$1.1$ & $-3.2$ &$5.7$ & $1.5$ &$4.6$ & $-0.4$ &$0.7$ &$1.08 \pm 0.02$ &$1.05$ \\
 & $20 \div 40$ &$0.544$ &$27.6$ & $20.8$ & $18.2$ &$2.2$ & $4.7$ & $-3.3$ &$0.9$ & $-0.2$ &$4.7$ & $-3.5$ &$11.6$ & $2.8$ &$7.3$ & $-0.4$ &$0.2$ &$1.09 \pm 0.02$ &$1.05$ \\
 & $40 \div 100$ &$0.0994$ &$41.6$ & $35.7$ & $21.3$ &$1.2$ & $7.3$ & $-3.3$ &$-0.4$ & $-0.1$ &$4.1$ & $-4.1$ &$12.9$ & $3.1$ &$10.0$ & $-3.7$ &$2.0$ &$1.09 \pm 0.02$ &$1.06$ \\
$0.3 \div 0.5$ &$4 \div 10$ &$8.80$ &$18.4$ & $9.3$ & $15.9$ &$2.0$ & $3.1$ & $-8.6$ &$-0.7$ & $1.2$ &$-1.0$ & $-5.3$ &$-7.4$ & $-3.9$ &$-0.6$ & $-1.5$ &$0.1$ &$1.08 \pm 0.02$ &$1.03$ \\
 & $10 \div 20$ &$2.31$ &$19.9$ & $13.6$ & $14.5$ &$2.0$ & $2.4$ & $-6.9$ &$-0.4$ & $0.7$ &$-1.5$ & $-3.3$ &$-7.7$ & $-4.1$ &$-0.0$ & $-0.4$ &$1.0$ &$1.08 \pm 0.02$ &$1.03$ \\
 & $20 \div 40$ &$1.12$ &$17.0$ & $12.6$ & $11.4$ &$-0.4$ & $3.2$ & $-3.7$ &$0.2$ & $0.5$ &$-0.7$ & $-3.0$ &$-5.3$ & $-2.8$ &$-0.3$ & $-0.0$ &$0.2$ &$1.08 \pm 0.03$ &$1.02$ \\
 & $40 \div 100$ &$0.264$ &$20.1$ & $17.1$ & $10.6$ &$0.7$ & $2.6$ & $-4.4$ &$-0.2$ & $0.4$ &$-0.9$ & $-2.1$ &$-3.4$ & $-2.1$ &$1.4$ & $1.3$ &$0.2$ &$1.07 \pm 0.03$ &$1.03$ \\
$0.5 \div 0.7$ &$4 \div 10$ &$4.50$ &$17.8$ & $13.3$ & $11.8$ &$3.3$ & $1.9$ & $6.5$ &$-1.4$ & $1.1$ &$-0.7$ & $-3.1$ &$-2.8$ & $-0.4$ &$0.9$ & $-0.7$ &$0.2$ &$1.14 \pm 0.06$ &$1.03$ \\
 & $10 \div 20$ &$1.86$ &$15.2$ & $11.8$ & $9.6$ &$0.8$ & $1.0$ & $3.1$ &$-0.5$ & $0.5$ &$-0.4$ & $-3.1$ &$-2.6$ & $-0.6$ &$0.2$ & $-1.8$ &$0.1$ &$1.12 \pm 0.06$ &$1.02$ \\
 & $20 \div 40$ &$0.703$ &$16.2$ & $13.5$ & $8.9$ &$2.0$ & $0.8$ & $-0.2$ &$-0.7$ & $0.5$ &$-0.6$ & $-2.2$ &$-2.3$ & $-0.7$ &$0.2$ & $-1.6$ &$0.2$ &$1.12 \pm 0.06$ &$1.02$ \\
 & $40 \div 100$ &$0.109$ &$31.9$ & $29.7$ & $11.4$ &$2.2$ & $-0.8$ & $3.2$ &$-0.8$ & $0.1$ &$-1.3$ & $-1.4$ &$-4.0$ & $-1.0$ &$-1.6$ & $-5.6$ &$0.1$ &$1.12 \pm 0.06$ &$1.01$ \\
$0.7 \div 1.0$ &$4 \div 10$ &$1.99$ &$27.8$ & $11.7$ & $25.2$ &$2.2$ & $2.9$ & $-21.9$ &$-1.6$ & $1.8$ &$-1.9$ & $-3.8$ &$-6.5$ & $-2.6$ &$1.7$ & $-2.9$ &$0.3$ &$0.79 \pm 0.11$ &$1.02$ \\
 & $10 \div 20$ &$0.639$ &$26.9$ & $11.2$ & $24.5$ &$1.4$ & $0.4$ & $-22.1$ &$-0.4$ & $1.1$ &$-1.6$ & $-1.7$ &$-5.2$ & $-2.0$ &$2.4$ & $-1.8$ &$0.3$ &$0.81 \pm 0.11$ &$1.01$ \\
 & $20 \div 40$ &$0.248$ &$22.4$ & $13.0$ & $18.2$ &$0.9$ & $1.6$ & $-15.3$ &$-0.3$ & $0.9$ &$-1.1$ & $-2.6$ &$-4.3$ & $-1.6$ &$1.1$ & $-0.3$ &$0.0$ &$0.85 \pm 0.11$ &$1.00$ \\
 & $40 \div 100$ &$0.0968$ &$18.5$ & $13.3$ & $13.0$ &$0.3$ & $2.1$ & $-9.0$ &$-0.5$ & $0.4$ &$-1.1$ & $-2.2$ &$-3.4$ & $-1.5$ &$1.3$ & $0.4$ &$0.5$ &$0.89 \pm 0.10$ &$1.01$ \\
\hline
\end{tabular}
\end{minipage}}
\caption{Diffractive DIS dijet cross section measured differentially
        as a function of $\zpom$ and $\qsq$.
        The statistical $\delta_{stat}$ and systematic $\delta_{sys}$ uncertainties are given together with the total uncertainty $\delta_{tot}$.
        Further details are given in Table~\ref{tab:tab1}.
	}
\label{tab:tab3}
\end{table}
\end{landscape}


\newpage
\clearpage
\begin{landscape}
\begin{table}[!ht]
\resizebox{.8\textwidth}{!}{\begin{minipage}{\textwidth}
\begin{tabular}{|c|c|c|c|c|c||c|c|c|c|c|c|c|c|c|c|c|c||c||c|}
\hline
$p^{\ast}_{\rm T,1}$ & $\qsq$ & $\frac{d^{2}\sigma}{dp^{\ast}_{\rm T,1}d\qsq}$  & $\delta_{tot}$ & $\delta_{stat}$ & $\delta_{sys}$
& $\delta_{\theta}$ & $\delta_{E}$ & $\delta_{HFS}$
& $\delta_{\qsq}$ & $\delta_{\xpom}$ & $\delta_{\beta}$ & $\delta_{p*_{T,1}}$
& $\delta_{\zpom}$ & $\delta_{x_{dijet}}$ & $\delta_{\Delta\eta*}$ & $\delta_{t}$
& $\delta_{bgr}$ & $1+\delta_{had}$ & $1+\delta_{rad}$\\
$\rm [GeV^{2}]$ & $\rm [GeV]$ & $\rm [pb/GeV^{3}]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$
& $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & $[\%]$ & & \\
\hline
$5.5 \div 7.0$ &$4 \div 6$ &$3.35$ &$15.6$ & $9.1$ & $12.7$ &$0.4$ & $5.9$ & $-1.2$ &$2.7$ & $1.1$ &$3.4$ & $-1.3$ &$5.8$ & $2.9$ &$0.6$ & $-1.0$ &$0.1$ &$1.05 \pm 0.05$ &$1.04$ \\
 & $6 \div 10$ &$1.84$ &$12.7$ & $7.1$ & $10.5$ &$3.0$ & $0.8$ & $-4.7$ &$-1.5$ & $0.3$ &$-1.0$ & $-1.8$ &$-2.9$ & $-1.3$ &$0.7$ & $-1.6$ &$0.1$ &$1.05 \pm 0.05$ &$1.02$ \\
 & $10 \div 18$ &$0.834$ &$12.3$ & $7.2$ & $9.9$ &$1.2$ & $0.9$ & $-3.2$ &$-1.1$ & $0.5$ &$-1.4$ & $-1.6$ &$-3.9$ & $-1.7$ &$0.6$ & $-1.1$ &$0.2$ &$1.05 \pm 0.05$ &$1.02$ \\
 & $18 \div 34$ &$0.344$ &$13.3$ & $8.6$ & $10.1$ &$1.5$ & $0.4$ & $-5.1$ &$0.3$ & $0.5$ &$-0.4$ & $-1.4$ &$-2.6$ & $-1.7$ &$1.1$ & $-0.0$ &$0.0$ &$1.06 \pm 0.05$ &$1.03$ \\
 & $34 \div 100$ &$0.0613$ &$15.8$ & $11.7$ & $10.6$ &$1.5$ & $0.7$ & $-5.0$ &$-0.4$ & $0.5$ &$-1.1$ & $-1.7$ &$-3.4$ & $-2.0$ &$0.9$ & $-1.6$ &$0.1$ &$1.07 \pm 0.04$ &$1.02$ \\
$7.0 \div 9.0$ &$4 \div 6$ &$1.23$ &$18.6$ & $15.1$ & $10.9$ &$-0.1$ & $3.0$ & $-6.1$ &$1.3$ & $1.2$ &$1.5$ & $-2.5$ &$0.1$ & $-0.7$ &$0.0$ & $-0.3$ &$0.5$ &$1.06 \pm 0.04$ &$1.05$ \\
 & $6 \div 10$ &$0.578$ &$16.4$ & $12.9$ & $10.1$ &$1.9$ & $3.3$ & $-3.9$ &$1.1$ & $0.8$ &$0.0$ & $-0.6$ &$-2.1$ & $-1.3$ &$1.5$ & $-0.8$ &$0.4$ &$1.06 \pm 0.04$ &$1.05$ \\
 & $10 \div 18$ &$0.287$ &$16.6$ & $12.6$ & $10.7$ &$0.4$ & $3.1$ & $-6.0$ &$-0.3$ & $0.7$ &$-0.1$ & $-1.3$ &$-0.3$ & $-0.4$ &$2.3$ & $0.7$ &$0.6$ &$1.06 \pm 0.05$ &$1.04$ \\
 & $18 \div 34$ &$0.100$ &$20.4$ & $17.6$ & $10.3$ &$0.3$ & $5.2$ & $-3.5$ &$-0.2$ & $0.8$ &$-0.2$ & $-0.5$ &$-1.0$ & $-0.4$ &$0.0$ & $-1.5$ &$0.7$ &$1.07 \pm 0.04$ &$1.04$ \\
 & $34 \div 100$ &$0.0276$ &$19.9$ & $17.4$ & $9.6$ &$-0.6$ & $3.7$ & $-3.1$ &$-0.6$ & $0.5$ &$-1.0$ & $-0.6$ &$-1.6$ & $-0.4$ &$1.0$ & $1.5$ &$0.6$ &$1.06 \pm 0.06$ &$1.04$ \\
$9.0 \div 15.0$ &$4 \div 6$ &$0.122$ &$30.1$ & $26.6$ & $14.2$ &$7.8$ & $0.5$ & $-5.5$ &$-0.1$ & $0.7$ &$0.3$ & $-5.5$ &$-1.2$ & $-0.5$ &$3.3$ & $-2.2$ &$0.8$ &$1.04 \pm 0.03$ &$1.06$ \\
 & $6 \div 10$ &$0.0511$ &$30.4$ & $24.7$ & $17.8$ &$1.9$ & $1.1$ & $-12.4$ &$-0.6$ & $1.3$ &$-0.0$ & $-6.3$ &$-2.8$ & $-0.9$ &$6.1$ & $-3.0$ &$0.4$ &$1.03 \pm 0.03$ &$1.05$ \\
 & $10 \div 18$ &$0.0207$ &$35.5$ & $30.0$ & $19.0$ &$1.4$ & $1.7$ & $-11.6$ &$-1.1$ & $1.0$ &$-1.4$ & $-6.5$ &$-6.4$ & $-2.9$ &$5.6$ & $-6.0$ &$0.5$ &$1.03 \pm 0.02$ &$1.05$ \\
 & $18 \div 34$ &$0.0160$ &$24.5$ & $20.1$ & $14.0$ &$1.6$ & $1.8$ & $-8.2$ &$0.0$ & $0.6$ &$-0.1$ & $-4.6$ &$-2.6$ & $-0.6$ &$5.2$ & $-2.5$ &$0.4$ &$1.04 \pm 0.04$ &$1.06$ \\
 & $34 \div 100$ &$0.0034$ &$31.9$ & $27.8$ & $15.7$ &$3.6$ & $0.3$ & $-9.0$ &$1.8$ & $0.7$ &$1.0$ & $-3.5$ &$-3.1$ & $-1.5$ &$6.0$ & $-4.9$ &$1.0$ &$1.05 \pm 0.03$ &$1.07$ \\
\hline
\end{tabular}
\end{minipage}}
\caption{Diffractive DIS dijet cross section measured differentially
        as a function of $p^{\ast}_{\rm T,1}$ and $\qsq$.
        The statistical $\delta_{stat}$ and systematic $\delta_{sys}$ uncertainties are given together with the total uncertainty $\delta_{tot.}$.
        Further details are given in Table~\ref{tab:tab1}.
	}
\label{tab:tab5}
\end{table}
\end{landscape}

\clearpage
\begin{table}[!ht]
\resizebox{.8\textwidth}{!}{\begin{minipage}{\textwidth}
\centering
\begin{tabular}{|c|c|ccccc|}
\hline
$\qsq\ \rm [GeV]$ & $\#Bin$ & $1$ & $2$ & $3$ & $4$ & $5$\\
\hline
$4 \div 6$ & $1$ &$100$ & $-5$ & $5$ &  & \\
$6 \div 10$ & $2$ & & $100$ & $1$ & $1$ & \\
$10 \div 18$ & $3$ & &  & $100$ & $-2$ & $1$\\
$18 \div 34$ & $4$ & &  &  & $100$ & $8$\\
$34 \div 100$ & $5$ & &  &  &  & $100$\\
\hline
$y$ & $\#Bin$ & $1$ & $2$ & $3$ & $4$ & $5$\\
\hline
$0.1 \div 0.2$ & $1$ &$100$ & $-7$ & $8$ & $5$ & $4$\\
$0.2 \div 0.3$ & $2$ & & $100$ & $-6$ & $8$ & $4$\\
$0.3 \div 0.5$ & $3$ & &  & $100$ & $-4$ & $7$\\
$0.5 \div 0.6$ & $4$ & &  &  & $100$ & $-10$\\
$0.6 \div 0.7$ & $5$ & &  &  &  & $100$\\
\hline
$\xpom$ & $\#Bin$ & $1$ & $2$ & $3$ & $4$ &\\
\hline
$-2.30 \div -2.10$ & $1$ &$100$ & $-55$ & $17$ & $-2$ &\\
$-2.10 \div -1.90$ & $2$ & & $100$ & $-41$ & $11$ &\\
$-1.90 \div -1.70$ & $3$ & &  & $100$ & $-31$ &\\
$-1.70 \div -1.52$ & $4$ & &  &  & $100$ &\\
\hline
$\zpom$ & $\#Bin$ & $1$ & $2$ & $3$ & $4$ & $5$\\
\hline
$0.0 \div 0.2$ & $1$ &$100$ & $-24$ & $8$ & $1$ & \\
$0.2 \div 0.4$ & $2$ & & $100$ & $-31$ & $10$ & $-2$\\
$0.4 \div 0.6$ & $3$ & &  & $100$ & $-45$ & $17$\\
$0.6 \div 0.8$ & $4$ & &  &  & $100$ & $-52$\\
$0.8 \div 1.0$ & $5$ & &  &  &  & $100$\\
\hline
\end{tabular}
\end{minipage}}
\caption{Correlation coefficients between data points for the single-differential measurements in 
	$\qsq$, $y$, $\xpom$ and $\zpom$.
	The values are given in per cent.}
\label{tab:tabrho1}
\end{table}

\clearpage
\begin{table}[!ht]
\resizebox{.8\textwidth}{!}{\begin{minipage}{\textwidth}
\centering
\begin{tabular}{|c|c|ccccc|}
\hline
$p^{\ast}_{\rm T,1}\ \rm [GeV]$ & $\#Bin$
 & $1$ & $2$ & $3$ & &\\
\hline
$5.5 \div 7.0$ & $1$ &$100$ & $-26$ & $1$ &&\\
$7.0 \div 9.0$ & $2$ & & $100$ & $-54$ &&\\
$9.0 \div 15.0$ & $3$ & &  & $100$ &&\\
\hline
$p^{\ast}_{\rm T,2}\ \rm [GeV]$ & $\#Bin$ 
& $1$ & $2$ & $3$ & & \\
\hline
$4.0 \div 6.5$ & $1$ &$100$ & $-36$ & $13$ &&\\
$6.5 \div 9.0$ & $2$ & & $100$ & $-46$ &&\\
$9.0 \div 15.0$ & $3$ & &  & $100$ &&\\
\hline
$\langle p^{\ast}_{\rm T}\rangle\ \rm [GeV]$ & $\#Bin$
        & $1$ & $2$ & $3$ & &\\
\hline
$4.75 \div 6.50$ & $1$ &$100$ & $-33$ & $12$ &&\\
$6.50 \div 9.00$ & $2$ & & $100$ & $-49$ &&\\
$9.00 \div 15.00$ & $3$ & &  & $100$ &&\\
\hline
$\Delta\eta^{\ast}$ & $\#Bin$ 
& $1$ & $2$ & $3$ & $4$ & $5$ \\
\hline
$0.00 \div 0.15$ & $1$ &$100$ & $-49$ & $13$ & $1$ & $2$\\
$0.15 \div 0.40$ & $2$ & & $100$ & $-29$ & $9$ & $1$\\
$0.40 \div 0.80$ & $3$ & &  & $100$ & $-19$ & $7$\\
$0.80 \div 1.30$ & $4$ & &  &  & $100$ & $-20$\\
$1.30 \div 3.00$ & $5$ & &  &  &  & $100$\\
\hline
\end{tabular}
\end{minipage}}
\caption{Correlation coefficients between data points for the single-differential measurements in $p^{\ast}_{\rm T,1}$, $p^{\ast}_{\rm T,2}$, 
		  $\langle p^{\ast}_{\rm T}\rangle$ and $\Delta\eta^{\ast}$.
        The values are given in per cent.}
\label{tab:tabrho2}
\end{table}

\clearpage
\begin{table}[!ht]
%\resizebox{.75\textwidth}{!}{\begin{minipage}{\textwidth}
\resizebox{\columnwidth}{!}{
\centering
\begin{tabular}{|c|c|c|cccccccccccccccc|}
\hline
$\zpom$ & $\qsq\ \rm [GeV^{2}]$ & $\#Bin$
 & $1$ & $2$ & $3$ & $4$ & $5$ & $6$ & $7$ & $8$ & $9$ & $10$ & $11$ & $12$ & $13$ & $14$ & $15$ & $16$\\
\hline
$0.0 \div 0.3$ & $4 \div 10$ & $1$ &$100$ & $1$ &  &  & $-32$ & $3$ &  &  & $11$ &  &  &  & $-3$ &  &  & \\
 & $10 \div 20$ & $2$ & & $100$ & $-2$ & $1$ & $3$ & $-41$ & $2$ &  &  & $19$ &  &  &  & $-5$ &  & \\
 & $20 \div 40$ & $3$ & &  & $100$ & $4$ &  & $2$ & $-37$ & $4$ &  &  & $20$ & $1$ &  &  & $-4$ & $1$\\
 & $40 \div 100$ & $4$ & &  &  & $100$ &  &  & $3$ & $-37$ &  &  & $1$ & $22$ &  &  & $1$ & $-5$\\
$0.3 \div 0.5$ & $4 \div 10$ & $5$ & &  &  &  & $100$ & $-3$ &  &  & $-46$ & $2$ &  &  & $15$ &  &  & \\
 & $10 \div 20$ & $6$ & &  &  &  &  & $100$ & $-3$ &  & $3$ & $-53$ & $2$ &  &  & $19$ &  & \\
 & $20 \div 40$ & $7$ & &  &  &  &  &  & $100$ & $-3$ &  & $2$ & $-51$ & $2$ &  & $-1$ & $17$ & \\
 & $40 \div 100$ & $8$ & &  &  &  &  &  &  & $100$ &  &  & $2$ & $-51$ &  &  &  & $21$\\
$0.5 \div 0.7$ & $4 \div 10$ & $9$ & &  &  &  &  &  &  &  & $100$ & $-5$ &  &  & $-47$ & $2$ &  & \\
 & $10 \div 20$ & $10$ & &  &  &  &  &  &  &  &  & $100$ & $-3$ &  & $1$ & $-46$ & $1$ & \\
 & $20 \div 40$ & $11$ & &  &  &  &  &  &  &  &  &  & $100$ & $-2$ &  & $2$ & $-44$ & $1$\\
 & $40 \div 100$ & $12$ & &  &  &  &  &  &  &  &  &  &  & $100$ &  &  & $1$ & $-51$\\
$0.7 \div 1.0$ & $4 \div 10$ & $13$ & &  &  &  &  &  &  &  &  &  &  &  & $100$ & $-4$ &  & \\
 & $10 \div 20$ & $14$ & &  &  &  &  &  &  &  &  &  &  &  &  & $100$ & $-5$ & \\
 & $20 \div 40$ & $15$ & &  &  &  &  &  &  &  &  &  &  &  &  &  & $100$ & $-2$\\
 & $40 \div 100$ & $16$ & &  &  &  &  &  &  &  &  &  &  &  &  &  &  & $100$\\
\hline
\end{tabular}
%\end{minipage}
}
\caption{Correlation coefficients between data points for the double-differential measurement in $\zpom$ and $\qsq$.
        The values are given in per cent.}
\label{tab:tabrho3}
\end{table}

%\clearpage
\begin{table}[!ht]
%\resizebox{.75\textwidth}{!}{\begin{minipage}{\textwidth}
\resizebox{\columnwidth}{!}{
\begin{tabular}{|c|c|c|ccccccccccccccc|}
\hline
$p^{\ast}_{\rm T,1}\ \rm [GeV]$ & $\qsq\ \rm [GeV^{2}]$ & $\#Bin$
 & $1$ & $2$ & $3$ & $4$ & $5$ & $6$ & $7$ & $8$ & $9$ & $10$ & $11$ & $12$ & $13$ & $14$ & $15$\\
\hline
$5.5 \div 7.0$ & $4 \div 6$ & $1$ &$100$ & $-7$ & $1$ &  &  & $-44$ & $2$ &  &  &  & $13$ & $1$ & $1$ & $1$ & \\
 & $6 \div 10$ & $2$ & & $100$ & $-3$ &  &  & $3$ & $-57$ & $3$ & $-1$ &  &  & $17$ & $1$ & $1$ & $1$\\
 & $10 \div 18$ & $3$ & &  & $100$ & $-2$ & $1$ & $1$ & $3$ & $-59$ & $1$ & $-1$ & $2$ & $1$ & $22$ & $1$ & $1$\\
 & $18 \div 34$ & $4$ & &  &  & $100$ & $3$ & $-1$ &  &  & $-58$ & $1$ & $2$ & $1$ & $2$ & $25$ & $2$\\
 & $34 \div 100$ & $5$ & &  &  &  & $100$ &  &  &  & $1$ & $-56$ &  &  & $1$ & $2$ & $27$\\
$7.0 \div 9.0$ & $4 \div 6$ & $6$ & &  &  &  &  & $100$ & $-7$ & $3$ & $3$ & $1$ & $-60$ & $2$ & $-5$ & $-6$ & $-3$\\
 & $6 \div 10$ & $7$ & &  &  &  &  &  & $100$ & $-4$ & $2$ & $1$ & $3$ & $-57$ &  & $-3$ & $-2$\\
 & $10 \div 18$ & $8$ & &  &  &  &  &  &  & $100$ &  & $2$ & $-6$ &  & $-60$ & $-4$ & $-4$\\
 & $18 \div 34$ & $9$ & &  &  &  &  &  &  &  & $100$ & $1$ & $-7$ & $-3$ & $-6$ & $-62$ & $-3$\\
 & $34 \div 100$ & $10$ & &  &  &  &  &  &  &  &  & $100$ & $-4$ & $-2$ & $-5$ & $-5$ & $-64$\\
$9.0 \div 15.0$ & $4 \div 6$ & $11$ & &  &  &  &  &  &  &  &  &  & $100$ & $-5$ & $13$ & $14$ & $7$\\
 & $6 \div 10$ & $12$ & &  &  &  &  &  &  &  &  &  &  & $100$ &  & $6$ & $3$\\
 & $10 \div 18$ & $13$ & &  &  &  &  &  &  &  &  &  &  &  & $100$ & $10$ & $9$\\
 & $18 \div 34$ & $14$ & &  &  &  &  &  &  &  &  &  &  &  &  & $100$ & $8$\\
 & $34 \div 100$ & $15$ & &  &  &  &  &  &  &  &  &  &  &  &  &  & $100$\\
\hline
\end{tabular}
%\end{minipage}
}
\caption{Correlation coefficients between data points for the double-differential measurement in $p^{\ast}_{\rm T,1}$ and $\qsq$.
        The values are given in per cent.}
\label{tab:tabrho5}
%\vspace{20cm}
\end{table}

\newpage
\begin{figure}[ht]
\centering
\subfloat{
\includegraphics[width=.48\textwidth]{d14-200f2a.eps}
}\hfill
\subfloat{
\includegraphics[width=.48\textwidth]{d14-200f2b.eps}
}\hfill
\subfloat{
\includegraphics[width=.48\textwidth]{d14-200f2c.eps}
}\hfill
\subfloat{
\includegraphics[width=.48\textwidth]{d14-200f2d.eps}
}\hfill
\caption{Distributions of the kinematic quantities $\qsq$, $p^{\ast}_{\rm T,1}$, $\xpom$
		and $\zpom$. The data are shown as black points compared to the sum of MC simulation estimates.
		The filled area shows the contribution of non-diffractive DIS,
		the dotted line shows the diffractive contribution with the elastically scattered proton added to the non-diffractive DIS
		and the dashed line displays the proton dissociation contribution added to the diffractive contribution with the elastically 
		scattered proton and the non-diffractive DIS contribution.
		The sum of all contributions including the resolved photon processes is given by the full line.
		The MC is reweighted to the data. The ratio of data to the MC prediction is shown in the lower part of %each plot.
		of the individual figures.
		}
\label{fig:figctr}
\end{figure}

\newpage
\begin{figure}[b]
\subfloat{
\includegraphics[bb=80bp  100bp 380bp 400bp, width=.4\textwidth]{d14-200f3a.eps}
}\hfill
\subfloat{
\includegraphics[bb=30bp  100bp 330bp 400bp, width=.4\textwidth]{d14-200f3b.eps}
}\hfill
\vspace{2cm}
\caption{Diffractive dijet differential cross section as a function of $\qsq$ and $y$.
        The inner error bars on the data points represent the statistical uncertainties, while the outer error bars include
        the systematic uncertainties added in quadrature. The NLO QCD prediction based on the H12006 Fit-B DPDF set is displayed as a white line. 
        The light shaded band indicates the uncertainty arising from hadronisation and the DPDF fit added
        in quadrature. The outer dark band shows the full theory uncertainty including the QCD scale uncertainty added in quadrature. 
	The ratio of the single-differential cross section to the NLO prediction is shown in the lower part of %the plot.
	the individual figures.
	}
\label{fig:fig1}
\end{figure}

\newpage
\begin{figure}[b]
\subfloat{
\includegraphics[bb=80bp  100bp 380bp 400bp, width=.4\textwidth]{d14-200f4a.eps}
}\hfill
\subfloat{
\includegraphics[bb=30bp  100bp 330bp 400bp, width=.4\textwidth]{d14-200f4b.eps}
}\hfill
\vspace{2cm}
\caption{Diffractive dijet differential cross section as a function of $\log\xpom$ and $\zpom$.
        The inner error bars on the data points represent the statistical uncertainties, while the outer error bars include
        the systematic uncertainties added in quadrature.
	Further details are given in figure~\ref{fig:fig1}.} 
\label{fig:fig2}
\end{figure}

\newpage
\clearpage
\begin{figure}[b]
\subfloat{
\includegraphics[bb=80bp  100bp 380bp 400bp, width=.4\textwidth]{d14-200f5a.eps}
}\hfill
\subfloat{
\includegraphics[bb=30bp  100bp 330bp 400bp, width=.4\textwidth]{d14-200f5b.eps}
}\hfill
\vspace{2cm}
\caption{Diffractive dijet differential cross section as a function of $p^{\ast}_{\rm T,1}$ and $p^{\ast}_{\rm T,2}$.
        The inner error bars on the data points represent the statistical uncertainties, while the outer error bars include
        the systematic uncertainties added in quadrature. 
        Further details are given in figure~\ref{fig:fig1}.}
\label{fig:fig3}
\end{figure}

\newpage
\clearpage
\begin{figure}[b]
\subfloat{
\includegraphics[bb=80bp  100bp 380bp 400bp, width=.4\textwidth]{d14-200f6a.eps}
}\hfill
\subfloat{
\includegraphics[bb=30bp  100bp 330bp 400bp, width=.4\textwidth]{d14-200f6b.eps}
}\hfill
\vspace{2cm}
\caption{Diffractive dijet differential cross section as a function of $\langle p^{\ast}_{\rm T}\rangle$ and $\Delta\eta^{\ast}$.
        The inner error bars on the data points represent the statistical uncertainties, while the outer error bars include
        the systematic uncertainties added in quadrature. 
        Further details are given in figure~\ref{fig:fig1}.}
\label{fig:fig4}
\end{figure}

\newpage
\clearpage
\begin{figure}[b!]
\includegraphics[bb=100bp  140bp 440bp 440bp, width=.64\textwidth]{d14-200f7.eps}
\vspace{4cm}
\caption{Double-differential cross section as a function of $\zpom$ and $\qsq$.
        The inner error bars on the data points represent the statistical uncertainties, while the outer error bars include
        the systematic uncertainties added in quadrature. 
        Further details are given in figure~\ref{fig:fig1}.}
\label{fig:fig5}
\end{figure}

\newpage
\begin{figure}[!b]
\includegraphics[bb=100bp  140bp 440bp 440bp, width=.64\textwidth]{d14-200f8.eps}
\vspace{4cm}
\caption{Ratio of the double-differential cross section to the NLO prediction as a function of $\zpom$  and $\qsq$.
%        The ratio of the double-differential cross section to the NLO prediction is displayed.
        The inner error bars on the data points represent the statistical uncertainties, while the outer error bars include
        the systematic uncertainties added in quadrature. 
        Further details are given in figure~\ref{fig:fig1}.}
\label{fig:fig6}
\end{figure}

\newpage
\clearpage
\begin{figure}[!b]
\includegraphics[bb=40bp  180bp 380bp 480bp, width=.64\textwidth]{d14-200f9.eps}
\vspace{6cm}
\caption{Double-differential cross section as a function of $p^{\ast}_{\rm T,1}$ and $\qsq$.
        The inner error bars on the data points represent the statistical uncertainties, while the outer error bars include
        the systematic uncertainties added in quadrature. 
        Further details are given in figure~\ref{fig:fig1}.}
\label{fig:fig9}
\end{figure}

\newpage
\clearpage
\begin{figure}[!b]
\includegraphics[bb=40bp  180bp 380bp 480bp, width=.64\textwidth]{d14-200f10.eps}
\vspace{6cm}
\caption{Ratio of the double-differential cross section to the NLO prediction as a function of $p^{\ast}_{\rm T,1}$ and $\qsq$.
%        The ratio of the double-differential cross section to the NLO prediction is displayed.
        The inner error bars on the data points represent the statistical uncertainties, while the outer error bars include
        the systematic uncertainties added in quadrature. 
        Further details are given in figure~\ref{fig:fig1}.}
\label{fig:fig10}
\end{figure}

%\subfloat{
%\includegraphics[width=0.3\textwidth]{fig1a.eps}
%}\hfill

\end{document}
