DESY 09032 ISSN 04189833
November 2009
Jet Production in Collisions at High
and Determination of
H1 Collaboration
The production of jets is studied in deepinelastic scattering at large negative four momentum transfer squared GeV using HERA data taken in 19992007, corresponding to an integrated luminosity of . Inclusive jet, 2jet and 3jet cross sections, normalised to the neutral current deepinelastic scattering cross sections, are measured as functions of , jet transverse momentum and proton momentum fraction. The measurements are well described by perturbative QCD calculations at nexttoleading order corrected for hadronisation effects. The strong coupling as determined from these measurements is (pdf).
Accepted by Eur. Phys. J. C
F.D. Aaron, C. Alexa, K. Alimujiang, V. Andreev, B. Antunovic, A. Asmone, S. Backovic, A. Baghdasaryan, E. Barrelet, W. Bartel, K. Begzsuren, A. Belousov, J.C. Bizot, V. Boudry, I. BozovicJelisavcic, J. Bracinik, G. Brandt, M. Brinkmann, V. Brisson, D. Bruncko, A. Bunyatyan, G. Buschhorn, L. Bystritskaya, A.J. Campbell, K.B. Cantun Avila, F. CassolBrunner, K. Cerny, V. Cerny, V. Chekelian, A. Cholewa, J.G. Contreras, J.A. Coughlan, G. Cozzika, J. Cvach, J.B. Dainton, K. Daum, M. Deák, Y. de Boer, B. Delcourt, M. Del Degan, J. Delvax, A. De Roeck, E.A. De Wolf, C. Diaconu, V. Dodonov, A. Dossanov, A. Dubak, G. Eckerlin, V. Efremenko, S. Egli, A. Eliseev, E. Elsen, A. Falkiewicz, P.J.W. Faulkner, L. Favart, A. Fedotov, R. Felst, J. Feltesse, J. Ferencei, D.J. Fischer, M. Fleischer, A. Fomenko, E. Gabathuler, J. Gayler, S. Ghazaryan, A. Glazov, I. Glushkov, L. Goerlich, N. Gogitidze, M. Gouzevitch, C. Grab, T. Greenshaw, B.R. Grell, G. Grindhammer, S. Habib, D. Haidt, C. Helebrant, R.C.W. Henderson, E. Hennekemper, H. Henschel, M. Herbst, G. Herrera, M. Hildebrandt, K.H. Hiller, D. Hoffmann, R. Horisberger, T. Hreus, M. Jacquet, M.E. Janssen, X. Janssen, V. Jemanov, L. Jönsson, A.W. Jung, H. Jung, M. Kapichine, J. Katzy, I.R. Kenyon, C. Kiesling, M. Klein, C. Kleinwort, T. Kluge, A. Knutsson, R. Kogler, V. Korbel, P. Kostka, M. Kraemer, K. Krastev, J. Kretzschmar, A. Kropivnitskaya, K. Krüger, K. Kutak, M.P.J. Landon, W. Lange, G. LaštovičkaMedin, P. Laycock, A. Lebedev, G. Leibenguth, V. Lendermann, S. Levonian, G. Li, K. Lipka, A. Liptaj, B. List, J. List, N. Loktionova, R. LopezFernandez, V. Lubimov, L. Lytkin, A. Makankine, E. Malinovski, P. Marage, Ll. Marti, H.U. Martyn, S.J. Maxfield, A. Mehta, A.B. Meyer, H. Meyer, H. Meyer, J. Meyer, V. Michels, S. Mikocki, I. MilcewiczMika, F. Moreau, A. Morozov, J.V. Morris, M.U. Mozer, M. Mudrinic, K. Müller, P. Murín, B. Naroska, Th. Naumann, P.R. Newman, C. Niebuhr, A. Nikiforov, G. Nowak, K. Nowak, M. Nozicka, B. Olivier, J.E. Olsson, S. Osman, D. Ozerov, V. Palichik, I. Panagoulias, M. Pandurovic, Th. Papadopoulou, C. Pascaud, G.D. Patel, O. Pejchal, E. Perez, A. Petrukhin, I. Picuric, S. Piec, D. Pitzl, R. Plačakytė, B. Pokorny, R. Polifka, B. Povh, T. Preda, V. Radescu, A.J. Rahmat, N. Raicevic, A. Raspiareza, T. Ravdandorj, P. Reimer, E. Rizvi, P. Robmann, B. Roland, R. Roosen, A. Rostovtsev, M. Rotaru, J.E. Ruiz Tabasco, Z. Rurikova, S. Rusakov, D. Šálek, D.P.C. Sankey, M. Sauter, E. Sauvan, S. Schmitt, C. Schmitz, L. Schoeffel, A. Schöning, H.C. SchultzCoulon, F. Sefkow, R.N. ShawWest, I. Sheviakov, L.N. Shtarkov, S. Shushkevich, T. Sloan, I. Smiljanic, Y. Soloviev, P. Sopicki, D. South, V. Spaskov, A. Specka, Z. Staykova, M. Steder, B. Stella, G. Stoicea, U. Straumann, D. Sunar, T. Sykora, V. Tchoulakov, G. Thompson, P.D. Thompson, T. Toll, F. Tomasz, T.H. Tran, D. Traynor, T.N. Trinh, P. Truöl, I. Tsakov, B. Tseepeldorj, J. Turnau, K. Urban, A. Valkárová, C. Vallée, P. Van Mechelen, A. Vargas Trevino, Y. Vazdik, S. Vinokurova, V. Volchinski, M. von den Driesch, D. Wegener, Ch. Wissing, E. Wünsch, J. Žáček, J. Zálešák, Z. Zhang, A. Zhokin, T. Zimmermann, H. Zohrabyan, F. Zomer, and R. Zus
I. Physikalisches Institut der RWTH, Aachen, Germany
Vinca Institute of Nuclear Sciences, Belgrade, Serbia
School of Physics and Astronomy, University of Birmingham, Birmingham, UK
InterUniversity Institute for High Energies ULBVUB, Brussels; Universiteit Antwerpen, Antwerpen; Belgium
National Institute for Physics and Nuclear Engineering (NIPNE) , Bucharest, Romania
Rutherford Appleton Laboratory, Chilton, Didcot, UK
Institute for Nuclear Physics, Cracow, Poland
Institut für Physik, TU Dortmund, Dortmund, Germany
Joint Institute for Nuclear Research, Dubna, Russia
CEA, DSM/Irfu, CESaclay, GifsurYvette, France
DESY, Hamburg, Germany
Institut für Experimentalphysik, Universität Hamburg, Hamburg, Germany
MaxPlanckInstitut für Kernphysik, Heidelberg, Germany
Physikalisches Institut, Universität Heidelberg, Heidelberg, Germany
KirchhoffInstitut für Physik, Universität Heidelberg, Heidelberg, Germany
Institute of Experimental Physics, Slovak Academy of Sciences, Košice, Slovak Republic
Department of Physics, University of Lancaster, Lancaster, UK
Department of Physics, University of Liverpool, Liverpool, UK
Queen Mary and Westfield College, London, UK
Physics Department, University of Lund, Lund, Sweden
CPPM, CNRS/IN2P3  Univ. Mediterranee, Marseille, France
Departamento de Fisica Aplicada, CINVESTAV, Mérida, Yucatán, México
Departamento de Fisica, CINVESTAV, México
Institute for Theoretical and Experimental Physics, Moscow, Russia
Lebedev Physical Institute, Moscow, Russia
MaxPlanckInstitut für Physik, München, Germany
LAL, Univ ParisSud, CNRS/IN2P3, Orsay, France
LLR, Ecole Polytechnique, IN2P3CNRS, Palaiseau, France
LPNHE, Universités Paris VI and VII, IN2P3CNRS, Paris, France
Faculty of Science, University of Montenegro, Podgorica, Montenegro
Institute of Physics, Academy of Sciences of the Czech Republic, Praha, Czech Republic
Faculty of Mathematics and Physics, Charles University, Praha, Czech Republic
Dipartimento di Fisica Università di Roma Tre and INFN Roma 3, Roma, Italy
Institute for Nuclear Research and Nuclear Energy, Sofia, Bulgaria
Institute of Physics and Technology of the Mongolian Academy of Sciences , Ulaanbaatar, Mongolia
Paul Scherrer Institut, Villigen, Switzerland
Fachbereich C, Universität Wuppertal, Wuppertal, Germany
Yerevan Physics Institute, Yerevan, Armenia
DESY, Zeuthen, Germany
Institut für Teilchenphysik, ETH, Zürich, Switzerland
PhysikInstitut der Universität Zürich, Zürich, Switzerland
Also at Physics Department, National Technical University, Zografou Campus, GR15773 Athens, Greece
Also at Rechenzentrum, Universität Wuppertal, Wuppertal, Germany
Also at University of P.J. Šafárik, Košice, Slovak Republic
Also at CERN, Geneva, Switzerland
Also at MaxPlanckInstitut für Physik, München, Germany
Also at Comenius University, Bratislava, Slovak Republic
Also at DESY and University Hamburg, Helmholtz Humboldt Research Award
Also at Faculty of Physics, University of Bucharest, Bucharest, Romania
Supported by a scholarship of the World Laboratory Björn Wiik Research Project
Also at Ulaanbaatar University, Ulaanbaatar, Mongolia
Deceased
Supported by the Bundesministerium für Bildung und Forschung, FRG, under contract numbers 05 H1 1GUA /1, 05 H1 1PAA /1, 05 H1 1PAB /9, 05 H1 1PEA /6, 05 H1 1VHA /7 and 05 H1 1VHB /5
Supported by the UK Science and Technology Facilities Council, and formerly by the UK Particle Physics and Astronomy Research Council
Supported by FNRSFWOVlaanderen, IISNIIKW and IWT and by Interuniversity Attraction Poles Programme, Belgian Science Policy
Partially Supported by Polish Ministry of Science and Higher Education, grant PBS/DESY/70/2006
Supported by the Deutsche Forschungsgemeinschaft
Supported by VEGA SR grant no. 2/7062/ 27
Supported by the Swedish Natural Science Research Council
Supported by the Ministry of Education of the Czech Republic under the projects LC527, INGO1P05LA259 and MSM0021620859
Supported by the Swiss National Science Foundation
Supported by CONACYT, México, grant 48778F
Russian Foundation for Basic Research (RFBR), grant no 1329.2008.2
This project is cofunded by the European Social Fund (75%) and National Resources (25%)  (EPEAEK II)  PYTHAGORAS II
1 Introduction
Jet production in neutral current (NC) deepinelastic scattering (DIS) at HERA provides an important testing ground for Quantum Chromodynamics (QCD). While inclusive DIS gives only indirect information on the strong coupling via scaling violations of the proton structure functions, the production of jets allows a direct measurement of . The Born level contribution to DIS (figure 1a) generates no transverse momentum in the Breit frame, where the virtual boson and the proton collide head on [1]. Significant transverse momentum in the Breit frame is produced at leading order (LO) in the strong coupling by the QCDCompton (figure 1b) and bosongluon fusion (figure 1c) processes.
In leading order the proton’s momentum fraction carried by the emerging parton is given by . The variable denotes the Bjorken scaling variable, the invariant mass of two jets of highest and the negative four momentum transfer squared. In the kinematical regions of low , low and low , bosongluon fusion dominates the jet production and provides direct sensitivity to the gluon component of proton density functions (PDFs) [2].
Analyses of inclusive jet production in DIS at high were previously performed by the H1 [3] and ZEUS [4] collaborations at HERA. These analyses are based on data taken during 1999 and 2000 (HERAI) and use jet observables to test the running of the strong coupling and extract its value at the boson mass. In this paper an integrated luminosity six times larger than available in the previous H1 analysis[3] is used. The ratios of jet cross sections to the corresponding NC DIS cross sections, henceforth referred to as normalised jet cross sections, are measured. These ratios benefit from a partial cancellation of experimental and theoretical uncertainties. The measurements are compared with perturbative QCD (pQCD) predictions at nexttoleading order (NLO) corrected for hadronisation effects, and is extracted from a fit of the predictions to the data. The measurements presented in this paper supersede the previously published normalised jet cross sections in [3].
2 Experimental Method
The data sample was collected with the H1 detector at HERA in the years 1999 to 2007 when HERA collided electrons or positrons^{1}^{1}1Unless otherwise stated, the term ”electron” is used in the following to refer to both electron and positron. of energy GeV with protons of energy , providing a centreofmass energy . The data sample used in this analysis corresponds to an integrated luminosity of , comprising recorded in collisions and in collisions.
2.1 H1 detector
A detailed description of the H1 detector can be found in [5, 6]. H1 uses a righthanded coordinate system with the origin at the nominal interaction point and the axis along the beam direction. The positive direction, also called the forward direction, is given by the outgoing proton beam. Polar angles and azimuthal angles are defined with respect to this axis. The pseudorapidity is related to the polar angle by . The detector components important for this analysis are described below.
The electromagnetic and hadronic energies are measured using the Liquid Argon (LAr) calorimeter in the polar angular range and with full azimuthal coverage [7]. The LAr calorimeter consists of an electromagnetic section ( to radiation lengths) with lead absorbers and a hadronic section with steel absorbers. The total depth of the LAr calorimeter varies between and hadronic interaction lengths. The energy resolution is for electrons and for hadrons, as obtained from test beam measurements [8]. In the backward region () energy is measured by a lead/scintillating fibre Spaghettitype Calorimeter (SpaCal) composed of an electromagnetic and a hadronic section [6]. The central tracking system () is located inside the LAr calorimeter and consists of drift and proportional chambers, complemented by a silicon vertex detector covering the range [9]. The chambers and calorimeters are surrounded by a superconducting solenoid providing a uniform field of inside the tracking volume. The luminosity is determined by measuring the event rate of the BetheHeitler process (), where the photon is detected in a calorimeter close to the beam pipe at .
2.2 Event and jet selection
The NC DIS events are triggered and selected by requiring a compact energy deposit in the electromagnetic part of the LAr calorimeter. The scattered electron is identified as the isolated cluster of highest transverse momentum [10]. Its reconstructed energy is requested to exceed 11 GeV. Only the regions of the calorimeter where the trigger efficiency is greater than are used for the detection of the scattered electron. These requirements ensure that the overall trigger efficiency reaches . In the central region, , the cluster has to be associated with a track measured in the inner tracking chambers and matched to the primary event vertex. The coordinate of the primary event vertex is required to be within of the nominal position of the interaction point.
The remaining clusters in the calorimeters and charged tracks are attributed to the hadronic final state, which is reconstructed using an energy flow algorithm that avoids double counting of energy [11, 12]. Electromagnetic and hadronic energy calibration and the alignment of the H1 detector are performed following the same procedure as in [10]. The total longitudinal energy balance, calculated as the difference of the total energy and the longitudinal component of the total momentum , calculated from all detected particles including the scattered electron, must satisfy . This requirement reduces contributions of DIS events with hard initial state photon radiation. For the latter events, the undetected photons propagating in the negative direction lead to values of this observable significantly lower than the expected value of . The requirement together with the scattered electron selection also reduces contributions from photoproduction, estimated using Monte Carlo simulations. Cosmic muon and beam induced background is reduced to a negligible level after combining these cuts with the primary event vertex selection. Elastic QED Compton and lepton pair production processes are suppressed by rejecting events containing additional isolated electromagnetic deposits and low hadronic calorimeter activity.
The kinematical range of this analysis is defined by
GeV and ,
where quantifies the inelasticity of the interaction. These two variables are reconstructed from the four momenta of the scattered electron and the hadronic final state particles using the electronsigma method [13]. The selection of events passing all the above cuts is the NC DIS sample, which forms the basis of the subsequent analysis.
The jet finding is performed in the Breit frame, where the boost from the laboratory system is determined by , and by the azimuthal angle of the scattered electron. Particles of the hadronic final state are clustered into jets using the inclusive algorithm [14] with the massless recombination scheme and with the distance parameter in the plane. The cut , where is the jet pseudorapidity in the laboratory frame, ensures that jets are contained within the acceptance of the LAr calorimeter and are well calibrated.
Jets are ordered by decreasing transverse momentum in the Breit frame, which is identical to the transverse energy for massless jets. The jet with highest is referred to as the ”leading jet”. Every jet with the transverse momentum in the Breit frame satisfying contributes to the inclusive jet cross section. The upper cutoff is necessary for the integration of the NLO calculation. The steeply falling transverse momentum spectrum leaves almost no jets above 50 GeV. Events with at least two (three) jets with transverse momentum are considered as 2jet (3jet) events. In order to avoid regions of phasespace where fixed order perturbation theory is not reliable [15], 2jet events are accepted only if the invariant mass of the two leading jets exceeds . The same requirement, , is applied to the 3jet events so that the 3jet sample is a subset of the 2jet sample.
After this selection, the inclusive jet sample contains a total of 143811 jets in 104014 events. The 2jet sample contains 47278 events and the 3jet sample 7054 events.
2.3 Definition of the observables
The measurements presented in this paper refer to the phasespace given in table 1. Normalised inclusive jet cross sections are measured as functions of and double differentially as function of and the transverse jet momentum in the Breit frame. Normalised 2jet and 3jet cross sections are presented as a function of . In addition the 2jet cross sections are measured double differentially as function of and the average transverse momentum of the two leading jets or as function of and of the proton momentum fraction . The 3jet cross section normalised to the 2jet cross section as function of is also presented.
The normalised jet cross sections are defined as the ratio of the differential inclusive jet, 2jet and 3jet cross sections to the differential NC DIS cross section in a given bin, multiplied by the respective bin width in case of a double differential measurement as indicated by the following equations:
(1)  
(2)  
(3) 
The normalised inclusive jet cross section can be viewed as the average jet multiplicity in a given region and the normalised multijet cross sections as multijet event rates.
2.4 Determination of normalised cross sections
In each analysis bin the normalised jet cross section is determined as
(4) 
Here denotes the number of inclusive jets or the number of 2jet or 3jet events, respectively, while represents the number of NC DIS events in that bin. The bin dependent correction factor takes into account the limited detector acceptance and resolution. The correction factors are determined from Monte Carlo simulations as the ratio of the normalised jet cross sections obtained from particles at the hadron level to the normalised jet cross sections calculated using reconstructed particles.
The following LO Monte Carlo event generators are used for the correction procedure: DJANGOH[16], which uses the Color Dipole Model with QCD matrix element corrections as implemented in ARIADNE[17], and RAPGAP[18], based on QCD matrix elements matched with parton showers in leading log approximation. In both Monte Carlo generators the hadronisation is modelled with Lund string fragmentation [19]. All generated events are passed through a GEANT3 [20] based simulation of the H1 apparatus and are reconstructed using the same program chain as for the data. Both RAPGAP and DJANGOH provide a good overall description of the inclusive DIS sample. To further improve the agreement between Monte Carlo and data for the jet samples, the Monte Carlo events are weighted as a function of and and as function of and of the leading jet in the Breit frame. In addition, they are weighted as a function of of the second and third jets when present [21]. After weighting, the simulations provide a good description of the shapes of all data distributions, some of which are shown in figure 2.
The binnings in , and used to measure the jet observables are given in table 2. The associated bin purities, defined as the fraction of the events reconstructed in a particular bin that originate from that bin on the generator level, are typically and always greater than . The correction factors deviate typically by less than from unity, but reach difference from unity in the bin GeV for the 2jet cross section. Arithmetic means of the correction factors determined from the reweighted RAPGAP and DJANGOH event samples are used and half of the difference is assigned as a model uncertainty.
The above correction factors include QED radiation and electroweak effects. The effects of QED radiation, which are typically , are corrected for by means of the HERACLES [22] program. The LEPTO event generator [23] is used to correct the and data for their different electroweak effects which largely cancel in normalised jet cross sections leaving them below . The resulting pure photon exchange cross sections obtained from and data samples are then averaged.
2.5 Experimental uncertainties
The systematic uncertainties of the jet observables are determined by propagating the corresponding estimated measurement errors through the full analysis:

The relative uncertainty of the electron energy calibration is typically between and for most of the events and increases up to for electrons in the forward direction. The absolute uncertainty of the electron polar angle is mrad. Uncertainties in the electron reconstruction affect the event kinematics and thus the boost to the Breit frame. This in turn leads to a relative error of to on the normalised cross sections for each of the two sources, electron polar angle and energy.

The relative uncertainty on the energy of the total reconstructed hadronic final state as well as of jets is estimated to be [21]. It is dominated by the uncertainty of the hadronic energy scale of the calorimeter. This error is estimated using a procedure similar to that used in [10] based on the transverse momentum conservation in the laboratory frame between the hadronic final state and the electron . This systematic uncertainty is reduced with respect to the previous measurement [3] due to the restricted pseudorapidity range in which jets are reconstructed and due to the improved statistics in the calibration procedure. The hadronic energy scale uncertainty affects mainly the jet cross section through the calibration of and, to a lesser extent, the NC DIS cross section through the reconstruction of . The resulting errors range between and and increase up to when exceeds 30 GeV. The relative uncertainty due to the hadronic energy scale is reduced on average by about for the normalised jet cross sections compared to the jet cross sections.

The model dependence of the detector correction factors is estimated as described in section 2.4. It reflects the sensitivity of the detector simulation to the details of the model, especially the parton showering, and their impact on the migration between adjacent bins in . The model dependence ranges typically from to for below and to above, independently of .

The uncertainties of the luminosity measurements, the trigger efficiency and the electron identification efficiency cancel in the normalised cross section. In addition, the model dependence of the QED radiative corrections, which is estimated to be [10], is expected to cancel in the normalised cross sections.
The statistical errors for the normalised inclusive jet cross section take into account the statistical correlations which arise because there can be more than one jet per event [21]. The statistical errors are considerably smaller compared to the previous HERAI publication [3]. They are typically between and for the normalised inclusive and 2jet cross sections and do not exceed 10% in the regions of high transverse momentum or high boson virtuality .
The dominant experimental errors on the jet cross sections arise from the uncertainty on the hadronic energy scale. The second most important source of systematic errors is the model dependence of the data correction, which becomes comparable to or exceeds the former in regions of highest jet . The overall experimental error, calculated as the quadratic sum of all the contributions inventoried above, ranges typically between and , but increases up to in the regions of highest or , dominated there by statistical uncertainties. The experimental errors for normalised cross sections are reduced by up to compared to those for unnormalised cross sections.
3 NLO QCD prediction of jet cross sections
Reliable quantitative predictions of jet cross sections in DIS require the perturbative calculations to be performed at least to nexttoleading order in the strong coupling. By using the inclusive jet algorithm with radius parameter , the observables used in the present analysis are infrared and collinear safe and the nonperturbative effects are expected to be small [2]. In addition, applying this algorithm in the Breit frame has the advantage that initial state singularities can be absorbed in the definition of the proton parton densities [24].
Jet cross sections are predicted at the parton level using the same jet definition as in the data analysis. The QCD predictions for the jet cross sections are calculated using the NLOJET++ program at NLO in the strong coupling [25]. The NC DIS cross section is calculated at with the DISENT package [26]. The FastNLO program [27] provides an efficient method to calculate these cross sections based on matrix elements from NLOJET++ and DISENT, convoluted with the PDFs of the proton and as a function of . The program includes a coherent treatment of the renormalisation and factorisation scale dependences of all ingredients to the cross section calculation, namely the matrix elements, the PDFs and .
When comparing data and theory predictions the strong coupling at the boson mass is taken to be and is evolved as a function of the renormalisation scale with two loop precision. The calculations are performed in the scheme for five massless quark flavours. The PDFs of the proton are taken from the CTEQ6.5M set [28]. The factorisation scale is taken to be and the renormalisation scale to be for the NLO predictions, with denoting for the inclusive jet, for 2jet and for the 3jet cross sections. This choice of the renormalisation scale is motivated by the presence of two hard scales, and in the jet production in DIS. For the calculation of inclusive DIS cross sections, the renormalisation scale is used. No QED radiation or exchange is included in the calculations, but the running of the electromagnetic coupling with is taken into account.
Hadronisation corrections are calculated for each bin using Monte Carlo event generators. These corrections are determined as the ratio of the cross section at the hadron level to the cross section at the parton level after parton showers. They typically differ by less than from unity and are obtained using the event generators DJANGOH and RAPGAP which agree to within to . The arithmetic means of the two Monte Carlo hadronisation correction factors are used, while the full difference is considered as systematic error.
DJANGOH and RAPGAP both use the Lund string model of hadronisation. The analytic calculations carried out in [29] provide an alternative method to estimate the effects of hadronisation and to crosscheck the hadronisation correction procedure described above. They are based on soft gluon power corrections and result in a shift of the perturbatively calculated spectrum of the inclusive jets:
(5) 
The size of the nonperturbative shift can be calculated up to one single nonperturbative parameter , which is the first moment of the effective nonperturbative coupling matched to the strong coupling at the scale . The value of , expected to be universal [30], was measured to be using event shapes observables in DIS by the H1 Collaboration [31]. The hadronisation correction factors so calculated for the inclusive jet cross section differ in most of the bins by less than from the average correction factor obtained from DJANGOH and RAPGAP and the maximum difference in all bins does not exceed which is within the estimated uncertainty of the hadronisation correction.
The dominant theoretical error is due to the uncertainty related to the neglected higher orders in the perturbative calculation. The accuracy of the NLO calculation is conventionally estimated by separately varying the chosen scales for and by factors in the arbitrary range 0.5 to 2. At high transverse momentum, above GeV, the pQCD calculations do not depend monotonically on in some bins. This happens in the two highest bins for the inclusive jet cross section and in six bins for the 2jet cross section, where the largest deviation from the central value is found for factors well inside the range 0.5 to 2. In such cases the difference between maximum and minimum cross sections found in the variation interval is taken, in order not to underestimate the scale dependence. Renormalisation and factorisation scale uncertainties are added in quadrature, the former outweighing the latter by a factor of two on average. The uncertainties originating from the PDFs are estimated using the CTEQ6.5M set of parton densities.
Normalised jet cross sections are calculated by dividing the predicted jet cross sections by the NC DIS cross sections. The renormalisation scale uncertainties are assumed to be uncorrelated between NC DIS and jet cross sections, as well as between 3jet and 2jet cross sections for their ratio, whereas the factorisation scale and the parameterisation uncertainty of the PDFs are assumed to be fully correlated.
4 Results
In the following, the normalised differential cross sections are presented for inclusive jet, 2jet and 3jet production at the hadron level. Tables 3 to 6 and figures 3 to 6 present the measured observables together with their experimental uncertainties and hadronisation correction factors applied to the NLO predictions. These measurements are subsequently used to extract the strong coupling as shown in the table 9 and figures 7 to 12.
4.1 Cross section measurements compared to NLO predictions
The normalised inclusive jet cross sections as a function of are shown in figure 3a and table 3 together with the NLO predictions and previous measurements by H1 based on HERAI data [3]. For comparison, the HERAI data points were corrected for the phase space difference due to the slightly smaller jet pseudorapidity range of the present analysis. The double differential results as a function of in six ranges of are given in figure 4 and table 4. Normalised 2jet (3jet) cross sections as a function of and their comparison to NLO are also shown in figure 3b (3c) and table 3, while the ratio 3jet to 2jet is shown in figure 3d. Figures 5, 6 and tables 5, 6 present the normalised 2jet cross section as a function of and in six ranges of .
The new measurement of the normalised inclusive jet cross section is compatible with the previous H1 data. The precision is improved by typically a factor of two, as can be seen for example in figure 3a. The QCD NLO predictions for all normalised jet cross sections provide a good description of the data over the whole phase space. In almost all bins the theory error, dominated by the scale uncertainty, is significantly larger than the total experimental uncertainty, which is dominated by the hadronic energy scale uncertainty.
The normalised inclusive jet cross section, which may be interpreted as the average jet multiplicity produced in NC DIS, increases with as the available phase space opens (figure 3a) as do the 2jet and 3jet rates (figure 3b and 3c). As increases, the jet spectra become harder as can be seen in figure 4 and 5. The 3jet rate is observed to be nearly seven times smaller than the 2jet rate as shown in figure 3d. The 2jet rates measured as a function of and the momentum fraction are well described by the NLO calculations (figure 6). Kinematic constraints from the considered range and the restricted invariant mass of the jets lead to a reduction of the 2jet rate at low and a rise at large with increasing .
4.2 Extraction of the strong coupling
The QCD predictions for jet production depend on and on the parton density functions of the proton. The strong coupling is determined from the measured normalised jet cross sections using the parton density functions from global analyses, which include inclusive deepinelastic scattering and other data. The determination is performed from individual observables and also from their combination.
QCD predictions of the jet cross sections are calculated as a function of with the FastNLO package using the CTEQ6.5M proton PDFs and applying the hadronisation corrections as described in section 3. Measurements and theory predictions are used to calculate a with the Hessian method [33], where parameters representing systematic shifts of detector related observables are left free in the fit. The shifts in the electron energy scale, electron polar angle and the hadronic final state energy scale found by the fit are consistent with the a priori estimated uncertainties. This method takes into account correlations of experimental uncertainties and has also been used in global data analyses [33, 34] and in previous H1 publications [3, 35]. The experimental uncertainty of is defined by the change in which gives an increase in of one unit with respect to the minimal value.
The correlations of the experimental uncertainties between data points were estimated using Monte Carlo simulations:

The measurement of the electron polar angle is assumed to be fully correlated [10].

The model dependence of the experimental correction factors is considered as fully uncorrelated after the averaging procedure described in section 2.5.
The sharing of correlated and uncorrelated contributions between the different sources of uncertainty has the following impact on the determination: when going from uncorrelated to fully correlated error for each source, the fitted value of typically varies by half the total experimental error and the estimated uncertainty by less than of .
The theory error is estimated by the so called offset method as the difference between the value of from the nominal fit to the value when the fit is repeated with independent variations of different sources of theoretical uncertainties as described in section 3. The resulting uncertainties due to the different sources are summed in quadrature. The up (or down) variations are applied simultaneously to all bins in the fit. The impact of hadronisation corrections on is between and , while that of the factorisation scale amounts to . The sensitivity of to the renormalisation scale variation of the inclusive NC DIS cross section alone is typically . The largest uncertainty, of typically to , corresponds to the accuracy of the NLO approximation to the jet cross sections estimated by varying the renormalisation scale as described in section 3. An alternative method to estimate the impact of missing orders, called the band method, developed by Jones et al. [36] was tested and, for the present measurement, it leads to a smaller uncertainty on of typically .
The uncertainty due to PDFs is estimated by propagating the CTEQ6.5M errors. The typical size of the resulting error is for determined from the normalised inclusive jet or 2jet cross sections and when measured with the normalised 3jet cross sections. This uncertainty is twice as large as that estimated with the uncertainties given for the MSTW2008nlo90cl set [37] which in turn exceeds the difference between values extracted with the central sets of CTEQ6.5M and MSTW2008nlo. The PDFs also depend on the value of . Potential biases on the extraction from that source have been studied in detail previously [3]. For this analysis, the resulting uncertainty is found to be negligible.
Individual fits of are made to each of the 24 measurements of the normalised double differential inclusive jet cross section, as shown in figure 7a. These individual determinations show the expected scale dependence. Equivalently, the values at each scale can be related to the value of the strong coupling at the mass as shown in figure 8.
Then is determined by a common fit to the normalised inclusive jet cross section in four bins for each region in . The resulting six values are evolved from the scale to the average in that region (figure 10a). Finally, a central value is extracted from a common fit to all 24 measurements and given in table 9. The result of evolving this value together with its associated uncertainty is also shown as the curve and surrounding band in figure 10a.
The same fit procedure of successive combination steps is applied to the 24 points of the normalised 2jet cross section with GeV (figure 7b, 9 and 10b). The bins with GeV are not used for the extraction of the strong coupling since the theory uncertainty is significantly larger than in the other bins (figure 5). The fit procedure is also applied to the 6 points of the normalised 3jet cross section (figure 10c). The normalised 3jet cross section (figure 3c), which is , is preferred to the ratio of the 3jet cross section to the 2jet cross section (figure 3d), which is , due to better sensitivity to the strong coupling. The three values of determined from the normalised inclusive jet (24 points), 2jet (24 points) and 3jet (6 points) cross sections are given in table 9 with experimental and theoretical uncertainties. All obtained values are compatible with each other within two standard deviations of the experimental uncertainty.
The impact of the choice of renormalisation scale on the central value of is studied in the case of the normalised inclusive jet cross section by repeating the fit procedure with and instead of . In the first case the central value of the is found to be approximatively smaller and in the latter approximatively bigger with respect to the nominal fit, a difference which is well inside the estimated theoretical uncertainties. Similar deviations are observed for the normalised 2jet and 3jet cross sections when is used instead of . To get information on the description of the data by the NLO calculations as a function of the renormalisation scale, the of the fit is studied in the case of the normalised inclusive jet cross section for different values of the parameter , defined by . The results are shown in figure 11, where the fit is repeated for different choices of and the corresponding values are shown. The lowest value is obtained for while choices above and below are disfavoured.
The sensitivity of the strong coupling determination procedure to the choice of the jet definition is tested for the normalised inclusive jet and 2jet cross sections by repeating all the extraction procedure using the anti metric [38] instead of , but keeping the recombination scheme and the distance parameter unchanged. The resulting central value of differs in both cases by less than from the central value extracted using the metric.
In each region the values of from different observables are combined taking into account statistical and systematic correlations. The resulting values, evolved from the scale to the average of the measurements in each region, are shown in figure 12. This visualises the running of for scales between 10 and 100 GeV and the corresponding experimental and theory uncertainties. All 54 data points are used in a common fit of the strong coupling taking the correlations into account with a fit quality (see table 9), which is also shown in figure 12.
The values of obtained in this way are also consistent with the world averages [39] and [40], and with the previous H1 and ZEUS determinations from inclusive jet production measurements [3, 4] and multijet production [41]. The experimental error on measured with each observable typically amounts to . The combination of different observables, even though partially correlated, gives rise to additional constraints on the strong coupling and leads to an improved experimental uncertainty of . The experimental error on is independent of the choice of renormalisation scale within the variation used to determine the theoretical uncertainty. The total error is strongly dominated by the theoretical uncertainty due to missing higher orders in the perturbative calculation which is about .
5 Conclusion
Measurements of the normalised inclusive, 2jet and 3jet cross sections in the Breit frame in deepinelastic electronproton scattering in the range and using the H1 data taken in years 1999 to 2007 are presented. Calculations at NLO QCD, corrected for hadronisation effects, provide a good description of the single and double differential cross sections as functions of the jet transverse momentum , the boson virtuality as well as of the proton momentum fraction . The strong coupling is determined from a fit of the NLO prediction to the measured normalised jet cross sections. The normalisation leads to cancellations of systematic effects, resulting in improved experimental and PDF uncertainties. The experimentally most precise determination of is derived from a common fit to the normalised jet cross sections:
The dominating source of the uncertainty is due to the renormalisation scale dependence, which is used to estimate the effect of missing higher orders beyond NLO in the pQCD prediction. This measurement improves the experimental precision on determinations from other recent jet measurements at HERA [3, 4]. The result is competitive with those from data [40, 42] and is in good agreement with the world average [39, 40].
Acknowledgements
We are grateful to the HERA machine group whose outstanding efforts have made this experiment possible. We thank the engineers and technicians for their work in constructing and maintaining the H1 detector, our funding agencies for financial support, the DESY technical staff for continual assistance and the DESY directorate for support and for the hospitality which they extend to the non DESY members of the collaboration. Furthermore we thank Gavin Salam, Matteo Cacciari, Mrinal Dasgupta and Zoltan Nagy for fruitful discussions.
References
 [1] R. P. Feynman, ”PhotonHadron Interactions”, Benjamin, New York (1972).
 [2] C. Adloff et al. [H1 Collaboration], Eur. Phys. J. C 19 (2001) 289 [hepex/0010054].
 [3] A. Aktas et al. [H1 Collaboration], Phys. Lett. B 653 (2007) 134 [arXiv:0706.3722].
 [4] S. Chekanov et al. [ZEUS Collaboration], Phys. Lett. B 649 (2007) 12 [hepex/0701039].
 [5] I. Abt et al. [H1 Collaboration], Nucl. Instrum. Meth. A 386 (1997) 310, ibid., 348.
 [6] R. D. Appuhn et al. [H1 SPACAL Group], Nucl. Instrum. Meth. A 386 (1997) 397.
 [7] B. Andrieu et al. [H1 Calorimeter Group], Nucl. Instrum. Meth. A 336 (1993) 460.
 [8] B. Andrieu et al. [H1 Calorimeter Group], Nucl. Instrum. Meth. A 336 (1993) 499; B. Andrieu et al. [H1 Calorimeter Group], Nucl. Instrum. Meth. A 350 (1994) 57.
 [9] D. Pitzl et al., Nucl. Instrum. Meth. A 454 (2000) 334 [hepex/0002044]; B. List et al., Nucl. Instrum. Meth. A 549 (2005) 33.
 [10] C. Adloff et al. [H1 Collaboration], Eur. Phys. J. C 30 (2003) 1 [hepex/0304003].
 [11] M. Peez, ”Search for deviations from the standard model in high transverse energy processes at the electron proton collider HERA” (in French), CPPMT200304 (available at http://wwwh1.desy.de/psfiles/theses/).
 [12] B. Portheault, ”First measurement of charged and neutral current cross sections with the polarized positron beam at HERA II and QCDelectroweak analyses” (in French), LAL0505 (available at http://wwwh1.desy.de/psfiles/theses/).
 [13] U. Bassler and G. Bernardi, Nucl. Instrum. Meth. A 361 (1995) 197 [hepex/9412004]; U. Bassler and G. Bernardi, Nucl. Instrum. Meth. A 426 (1999) 583 [hepex/9801017].
 [14] S. D. Ellis and D. E. Soper, Phys. Rev. D 48 (1993) 3160 [hepph/9305266]; S. Catani et al., Nucl. Phys. B 406 (1993) 187.
 [15] S. Frixione and G. Ridolfi, Nucl. Phys. B 507 (1997) 315 [hepph/9707345].
 [16] K. Charchula, G. A. Schuler and H. Spiesberger, DJANGOH 1.4, Comput. Phys. Commun. 81 (1994) 381.
 [17] L. Lönnblad, ARIADNE 4.08, Comput. Phys. Commun. 71 (1992) 15.
 [18] H. Jung, RAPGAP 2.08, Comput. Phys. Commun. 86 (1995) 147.
 [19] B. Andersson, G. Gustafson, G. Ingelman and T. Sjostrand, JETSET 7.41, Phys. Rept. 97 (1983) 31.
 [20] R. Brun et al., ”GEANT3”, CERNDD/EE/841.
 [21] M. Gouzevitch, ”Measurement of the strong coupling constant alpha(s) with hadronic jets in Deep Inelastic Scattering” (in French), DESYTHESIS2008047 (available at http://wwwh1.desy.de/psfiles/theses/).
 [22] A. Kwiatkowski, H. Spiesberger and H. J. Möhring, HERACLES 4.63, Comput. Phys. Commun. 69 (1992) 155.
 [23] G. Ingelman, A. Edin and J. Rathsman, LEPTO 6.5, Comput. Phys. Commun. 101 (1997) 108 [hepph/9605286].
 [24] B.R. Webber, J. Phys. G19 (1993) 1567.
 [25] Z. Nagy and Z. Trocsanyi, NLOJET++ 4.0.1, Phys. Rev. Lett. 87 (2001) 082001 [hepph/0104315].
 [26] S. Catani and M. H. Seymour, Nucl. Phys. B 485 (1997) 291 [Erratumibid. B 510 (1998) 503] [hepph/9605323].
 [27] T. Kluge, K. Rabbertz and M. Wobisch, FastNLO 1.0, hepph/0609285, published in ”Deep inelastic scattering DIS2006, Proceedings of the 14th workshop”, eds. M. Kuze, K. Nagano, K. Tokushuka, 483.
 [28] W. K. Tung et al. JHEP 0702 (2007) 053 [hepph/0611254].
 [29] M. Dasgupta and Y. Delenda, JHEP 0907 (2009) 004 [arXiv:0903.2187]; M. Dasgupta, L. Magnea and G. P. Salam, JHEP 0802 (2008) 055 [arXiv:0712.3014]; M. Dasgupta and B. R. Webber, Eur. Phys. J. C 1 (1998) 539 [hepph/9704297]; M. Dasgupta and B. R. Webber, JHEP 9810 (1998) 001 [hepph/9809247].
 [30] Y. L. Dokshitzer, G. Marchesini and B. R. Webber, Nucl. Phys. B 469 (1996) 93 [hepph/9512336].
 [31] A. Aktas et al. [H1 Collaboration], Eur. Phys. J. C 46 (2006) 343 [hepex/0512014].
 [32] C. Pascaud and F. Zomer, preprint LAL 9505 (1995).
 [33] M. Botje, Eur. Phys. J. C 14 (2000) 285 [hepph/9912439].
 [34] V. Barone, C. Pascaud and F. Zomer, Eur. Phys. J. C 12 (2000) 243 [hepph/9907512].
 [35] C. Adloff et al. [H1 Collaboration], Phys. Lett. B 479 (2000) 358 [hepex/0003002].
 [36] R. W. L. Jones et al., JHEP 0312 (2003) 007 [hepph/0312016].
 [37] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, arXiv:0901.0002.
 [38] M. Cacciari, G. P. Salam and G. Soyez, JHEP 0804 (2008) 063 [arXiv:0802.1189].
 [39] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667 (2008) 1.
 [40] S. Bethke, Prog. Part. Nucl. Phys. 58 (2007) 351 [hepex/0606035].
 [41] S. Chekanov et al. [ZEUS Collaboration], Eur. Phys. J. C 44 (2005) 183 [hepex/0502007].
 [42] G. Dissertori et al., JHEP 0802 (2008) 040 [arXiv:0712.0327].
NC DIS Selection  GeV  

Inclusive jet  
2jet  
3jet 
bin number  corresponding range 

1  
2  
3  
4  
5  
6 
Normalised inclusive jet cross section in bins of  
total  total  single contributions to correlated uncertainty  hadronisation  hadronisation  
bin  normalised  statistical  total  uncorrelated  correlated  electron  electron  hadronic  correction  correction 
cross  uncert.  uncert.  uncertainty  uncert.  energy scale  polar angle  energy scale  factor  uncertainty  
section  (%)  (%)  (%)  (%)  (%)  (%)  (%)  (%)  
1  
2  
3  
4  
5  
6  
Normalised 2jet cross section in bins of  
1  
2  
3  
4  
5  
6  
Normalised 3jet cross section in bins of  
1  
2  
3  
4  
5  
6  
3jet cross section normalised to 2jet cross section in bins of  
1  
2  
3  
4  
5  
6 
Normalised inclusive jet cross section in bins of and  
total  total  single contributions to correlated uncertainty  hadronisation  hadronisation  
bin  normalised  statistical  total  uncorrelated  correlated  electron  electron  hadronic  correction  correction 
cross  uncert.  uncert.  uncertainty  uncert.  energy scale  polar angle  energy scale  factor  uncertainty  
section  (%)  (%)  (%)  (%)  (%)  (%)  (%)  (%)  
1 a  
1 b  
1 c  
1 d  
2 a  
2 b  
2 c  
2 d  
3 a  
3 b  
3 c  
3 d  
4 a  
4 b  
4 c  
4 d  
5 a  
5 b  
5 c  
5 d  
6 a 