Measurement of the production cross section for Z / gamma * in association with jets in pp collisions at root s = 7 TeV with the ATLAS detector

Aad, G.; Abbott, B.; Abdallah, J.; Abdelalim, A. A.; Abdesselam, A.; Abdinov, O.; Abi, B.; Abolins, M.; Abramowicz, H.; Abreu, H.; Acerbi, E.; Acharya, B. S.; Adams, D. L.; Addy, T. N.; Adelman, J.; Aderholz, M.; Adomeit, S.; Adragna, P.; Adye, T.; Aefsky, S.; Aguilar-Saavedra, J. A.; Aharrouche, M.; Ahlen, S. P.; Ahles, F.; Ahmad, A.; Ahsan, M.; Aielli, G.; Akdogan, T.; Åkesson, Torsten; Akimoto, G.; Akimov, A. V.; Akiyama, A.; Alam, M. S.; Alam, M. A.; Albert, J.; Albrand, S.; Aleksa, M.; Aleksandrov, I. N.; Alessandria, F.; Alexa, C.; Alexander, G.; Alexandre, G.; Alexopoulos, T.; Alhroob, M.; Aliev, M.; Alimonti, G.; Alison, J.; Aliyev, M.; Allport, P. P.; Allwood-Spiers, S. E.


I. INTRODUCTION
The study of the production of jets of particles in association with a Z/γ * boson in proton-proton collisions provides a stringent test of perturbative quantum chromodynamics (pQCD). In addition, the proper understanding of these processes in the Standard Model (SM) is a fundamental element of the LHC physics program, since they constitute backgrounds in searches for new physics. These SM background contributions are estimated using next-to-leading order (NLO) pQCD calculations, and Monte Carlo (MC) predictions that include leading-order (LO) matrix elements supplemented by parton showers. The latter are affected by large scale uncertainties and need to be tuned and validated using data. Measurements of Z/γ * +jets production have been previously reported in proton-antiproton collisions at √ s = 1.96 TeV [1] and in proton-proton collisions at √ s = 7 TeV [2]. This article presents measurements of jet production in events with a Z/γ * boson in the final state, using 36±1 pb −1 of data collected by the ATLAS experiment in 2010 at √ s = 7 TeV. In this period, the accelerator operated with a moderate instantaneous luminosity of up to 2.1 × 10 32 cm −2 s −1 , and a long spacing of 150 ns between proton bunches, leading to relatively low collision rates and low rates of multiple proton-proton interactions per bunch crossing (pileup) and out-of-time pileup, which makes this data sample especially suitable for cross section measurements at low jet transverse momentum p T [3]. Events are selected with a Z/γ * decaying into a pair of electrons (e + e − ) or muons (µ + µ − ), and the measurements are corrected for detector effects. Inclusive jet differential cross sections are measured as functions of jet transverse momentum, p T , and rapidity, |y|, and total cross sections as functions of jet multiplicity, N jet , in well-defined kinematic regions for the leptons and jets in the final state. Differential cross sections are also measured as functions of p T and |y| of the leading jet (highest p T ) and second leading jet in Z/γ * events with at least one and two jets in the final state, respectively. For the latter, the cross section is measured as a function of the invariant mass and the angular separation of the two leading jets. The data are compared to NLO pQCD predictions [4,5], including non-perturbative contributions, and to predictions from several MC programs.
The paper is organized as follows. The detector is described in the next Section. Section III discusses the event selection, while Section IV provides details of the simulations used in the measurements and Sections V and VI describe the reconstruction of jets and leptons, respectively. The estimation of background contributions is described in Section VII. Selected uncorrected distributions are presented in Section VIII, and the procedure used to correct the measurements for detector effects is explained in Section IX. The study of systematic uncertainties is discussed in Section X. The NLO pQCD predictions are described in Section XI. The measured cross sections are presented separately for the electron and muon channels in Section XII, where the combination of the electron and muon results is also discussed. Finally, Section XIII provides a summary.

II. EXPERIMENTAL SETUP
The ATLAS detector [6] covers almost the whole solid angle around the collision point with layers of tracking detectors, calorimeters and muon chambers. The ATLAS inner detector (ID) has full coverage in φ and covers the pseudorapidity range |η| < 2.5. It consists of a silicon pixel detector, a silicon microstrip detector (SCT), and a straw tube tracker (TRT) which also measures transition radiation for particle identification, all immersed in a 2 tesla axial magnetic field produced by a solenoid.
High-granularity liquid-argon (LAr) electromagnetic sampling calorimeters, with very good energy and position resolution [7], cover the pseudorapidity range |η| < 3.2. The hadronic calorimetry in the range |η| < 1.7 is provided by a scintillator-tile calorimeter, consisting of a large barrel and two smaller extended barrel cylinders, one on either side of the central barrel. In the end caps (|η| > 1. 5), LAr hadronic calorimeters match the outer |η| limits of the end cap electromagnetic calorimeters. The LAr forward calorimeters provide both electromagnetic and hadronic energy measurements, and extend the coverage to |η| < 4.9.
The muon spectrometer measures the deflection of muon tracks in the large superconducting air-core toroid magnets in the pseudorapidity range |η| < 2.7, instrumented with separate trigger and high-precision tracking chambers. Over most of the η-range, a precision measurement of the track coordinates in the principal bending direction of the magnetic field is provided by monitored drift tubes. At large pseudorapidities, cathode strip chambers with higher granularity are used in the innermost plane over 2.0 < |η| < 2.7. The muon trigger system, which covers the pseudorapidity range |η| < 2.4, consists of resistive plate chambers in the barrel (|η| < 1.05) and thin gap chambers in the end cap regions (1.05 < |η| < 2.4), with a small overlap in the |η| =1.05 region.
The data samples considered in this paper were collected with tracking detectors, calorimeters, muon chambers, and magnets fully operational, and correspond to a total integrated luminosity of 36 pb −1 .
In the case of the Z/γ * → e + e − analysis, events are selected online using a trigger that requires the presence of at least one identified electron candidate in the calorimeter with transverse energy above 15 GeV in the region |η| < 2.5.
The events are then selected to have two oppositely charged reconstructed electrons (medium quality electrons, as described in Ref. [8]) with transverse energy E e T > 20 GeV, pseudorapidity in the range |η e | < 2.47 (where the transition region between calorimeter sections 1.37 < |η e | < 1.52 is excluded), and a dilepton invariant mass in the range 66 GeV< m e + e − < 116 GeV, which optimizes the signal sensitivity.
The Z/γ * → µ + µ − sample is collected online using a trigger that requires the presence of at least one muon candidate reconstructed in the muon spectrometer, consistent with having originated from the interaction region with p T > 10 GeV or p T > 13 GeV, depending on the data period, and with the majority of the data taken with the higher threshold, and |η| < 2.4. The muon candidates are associated with track segments reconstructed in the inner detectors which, combined with the muon spectrometer information, define the final muon track. Combined muon tracks with p µ T > 20 GeV and |η µ | < 2.4 are selected. A number of quality requirements are applied to the muon candidates [9]: the associated inner detector track segment is required to have a minimum number of hits in the pixel, SCT and TRT detectors; and the muon transverse and longitudinal impact parameters, d 0 and z 0 , with respect to the reconstructed primary vertex are required to be d 0 /σ(d 0 ) < 3 and z 0 < 10 mm in the r − φ and r − z planes, respectively, where σ(d 0 ) denotes the d 0 resolution. The muons are required to be isolated: the scalar sum of the transverse momenta of the tracks in an η − φ cone of radius 0.2 around the muon candidate is required to be less than 10% of the muon p T . Events are selected with two oppositely charged muons and an invariant mass 66 GeV< m µ + µ − < 116 GeV.
In both analyses, events are required to have a reconstructed primary vertex of the interaction with at least 3 tracks associated to it, which suppresses beam-related background contributions and cosmic rays. The selected dilepton samples contain a total of 9705 and 12582 events for the electron and muon channels, respectively.

IV. MONTE CARLO SIMULATION
Monte Carlo event samples are used to compute detector acceptance and reconstruction efficiencies, determine background contributions, correct the measurements for detector effects, and estimate systematic uncertainties on the final results.
Samples of simulated Z/γ * (→ e + e − )+jets and Z/γ * (→ µ + µ − )+jets events with a dilepton invariant mass above 40 GeV are generated using ALPGEN v2.13 [10] (including LO matrix elements for up to 2 → 5 parton scatters) interfaced to HERWIG v6.510 [11] for parton shower and fragmentation into particles, and to JIMMY v4.31 [12] to model underlying event (UE) contributions. Similar samples are generated using Sherpa 1.2.3 [13] with an UE modeling according to Ref. [14]. For the ALPGEN samples CTEQ6L1 [15] parton density functions (PDFs) are employed, while for Sherpa CTEQ6.6 [16] is used. The ALPGEN and Sherpa samples are normalized to the next-tonext-to-leading order (NNLO) pQCD inclusive Drell-Yan prediction of 1.07±0.05 nb, as determined by the FEWZ [17] program using the MSTW2008 PDFs. In addition, Z/γ * +jets samples (qq → Z/γ * g and qg → Z/γ * q processes witĥ p T > 10 GeV, wherep T is the transverse momentum defined in the rest frame of the hard interaction) are produced using PYTHIA v6.423 [18] and HERWIG plus JIMMY with MRST2007LO * [19] PDFs. For the ALPGEN and HERWIG plus JIMMY MC samples the AUET1 [20] tuned set of parameters is used to model the UE activity in the final state. In the case of the PYTHIA samples, the AMBT1 [21] tune is employed.
The MC samples are generated with minimum bias interactions from PYTHIA overlaid on top of the hard-scattering event in order to account for the presence of the pileup experienced in the data. The number of minimum bias (MB) interactions follows a Poisson distribution with a mean of two, which is appropriate for the 2010 data. The MC generated samples are then passed through a full simulation [24] of the ATLAS detector and trigger system, based on GEANT4 [25]. The simulated events are reconstructed and analyzed with the same analysis chain as for the data, using the same trigger and event selection criteria, and re-weighted such that the distribution of the number of primary vertices matches that of the data.
The multi-jets background contributions in the electron and muon channels are determined using data, as discussed in Section VII.

V. JET RECONSTRUCTION
Jets are defined using the anti-k t jet algorithm [26] with the distance parameter set to R = 0.4. Energy depositions reconstructed as calorimeter clusters are the inputs to the jet algorithm in data and MC simulated events. The same jet algorithm is applied to final state particles in the MC generated events to define jets at particle level [27]. The jet kinematics in data and MC simulated events are corrected to account for the following effects: the presence of additional proton-proton interactions per bunch crossing, leading to an additional energy offset of (500 ± 160) MeV within the jet cone for each extra interaction [28]; the position of the primary vertex of the interaction; and the measurement biases induced by calorimeter non-compensation, additional dead material, and out-of-cone effects. The measured jet p T is corrected for detector effects back to the true jet energy [29] using an average correction, computed as a function of the jet transverse momentum and pseudorapidity, and extracted from inclusive jet MC samples. The measured jet p T is reconstructed with a resolution of about 10% at low p T which improves to 6% for p T about 200 GeV. The measured jet angular variables y and φ are reconstructed with no significant shift and a resolution better than 0.05, which improves as the jet transverse momentum increases.
In this analysis, jets are selected with corrected p T > 30 GeV and |y| < 4.4 to ensure full containment in the instrumented region. Events are required to have at least one jet well separated from the final state leptons from the Z/γ * decay. Jets within a cone of radius 0.5 around any selected lepton are not considered. Additional quality criteria are applied to ensure that jets are not produced by noisy calorimeter cells, and to avoid problematic detector regions.
The final sample for Z/γ * (→ e + e − )+jets contains 1514, 333, 62, and 15 events with at least one, two, three, and four jets in the final state, respectively. Similarly, the Z/γ * (→ µ + µ − )+jets sample contains 1885, 422, 93, and 20 events with at least one, two, three, and four jets in the final state, respectively.

VI. LEPTON RECONSTRUCTION
Samples of Z/γ * → e + e − and Z/γ * → µ + µ − events in data and MC simulation, together with the world average values for the Z boson mass and width, are used to determine the absolute scale and resolution of the energy/momentum of the leptons, to validate calibration-and alignment-related constants in data, and to check the MC description [30]. In addition, the trigger and offline lepton reconstruction efficiencies are studied using control samples in data, and the results are compared to the simulation. The differences observed between data and MC predictions define scale factors which are applied in the analysis to the simulated samples before they are used to correct the measurements for detector effects.
For the electron channel, the trigger and offline electron reconstruction and identification efficiencies for single electrons are estimated using W → eν and Z/γ * → e + e − events in data and compared to MC predictions. In the kinematic range for the electrons considered in the analysis (see Section III), the trigger and offline efficiencies per electron are above 99% and 93%, respectively. The study indicates a good agreement between data and simulated trigger efficiencies with a MC-to-data scale factor of 0.995 ± 0.005. The simulation tends to overestimate the offline efficiencies. Scale factors in the range between 0.901 ± 0.045 and 0.999 ± 0.016, depending on η e and E e T , for E e T > 20 GeV, are applied per lepton to the MC samples to account for this effect.
In the muon analysis, the trigger and offline muon reconstruction efficiencies are also estimated using the data and are compared to simulation. The measured average single muon trigger efficiency is about 85%, independent of p µ T , and varies from 80% for |η µ | < 0.63 and 73% for 0.63 < |η µ | < 1.05 to 94% for 1.05 < |η µ | < 2.4, limited mainly by the trigger chamber geometric acceptance. The measured average offline muon reconstruction efficiency is about 92% and approximately independent of p µ T . The MC simulation predicts efficiencies very similar to those in the data, but tends to overestimate the average offline reconstruction efficiency by about 1%. This originates from the transition region between the barrel part and the endcap wheels at |η| ∼ 1, where the simulation overestimates the offline reconstruction efficiency by about 6%. The latter is attributed to the limited accuracy of the magnetic field map used in this region which leads to a small mismeasurement of the stand-alone muon momentum and an overestimation in the simulated efficiency. Scale factors are applied in the analysis that take this effect into account.

VII. BACKGROUND ESTIMATION
The background contribution to the electron and muon analyses from SM processes is estimated using MC simulated samples, as discussed in Section IV, with the exception of the multi-jets background that is estimated using data.
The multi-jets background contribution in the Z/γ * (→ e + e − )+jets analysis is estimated using a control data sample with two electron candidates which pass a loose selection but fail to pass the medium identification requirements. This sample is dominated by jets faking electrons in the final state and is employed to determine the shape of the multi-jets background under each of the measured distributions. The normalization of the multi-jets background events in the signal region is extracted from a fit to the measured inclusive dilepton invariant mass spectrum with nominal lepton requirements, using as input the observed shape of the multi-jets contribution in data and the MC predictions for the shape of the signal and the rest of the SM background processes. The multi-jets background contribution to the measured inclusive jet multiplicity varies between 3.2 ± 0.5(stat.) +0.3 −0.2 (syst.)% for N jet ≥ 1 and 4.5 ± 1.9(stat.) +0. 4 −0.2 (syst.)% for N jet ≥ 4. The quoted total systematic uncertainty includes: uncertainties related to the details of the parameterization and the mass range used to fit the measured dilepton invariant mass spectrum; uncertainties on the shape of the dilepton invariant mass distribution, as determined in the control sample; and uncertainties on the shape of the simulated dilepton invariant mass distribution for the other SM processes.
In the Z/γ * (→ µ + µ − )+jets case, the multi-jets background mainly originates from heavy-flavour jet production processes, with muons from bottom and charm quark decays, as well as from the decay-in-flight of pions and kaons, which are highly suppressed by the isolation requirement applied to the muon candidates. The isolation criterion of the muon pair, defined as the isolation of the least-isolated muon candidate, is used together with the dimuon invariant mass to estimate the remaining multi-jets background contribution. The MC simulation indicates that, for multi-jet processes, the muon isolation is not correlated with the dimuon invariant mass, and so the ratio of isolated to non-isolated muon pairs (as defined with an inverted isolation criterion) does not depend on the dimuon mass. The multi-jets background with isolated muons with 66 GeV< m µ + µ − < 116 GeV is therefore extracted from data as the ratio between the number of isolated and non-isolated dimuon candidates in the region 40 GeV< m µ + µ − < 60 GeV multiplied by the number of non-isolated dimuon candidates in the range 66 GeV< m µ + µ − < 116 GeV. A small contribution from top pair production processes is subtracted from the data according to MC predictions. The multi-jets background contribution to the Z/γ * (→ µ + µ − )+jets analysis is of the order of one per mille and therefore negligible.
In the electron channel, the total background increases from 5% to 17% as the jet multiplicity increases and is dominated by multi-jet processes, followed by contributions from tt and diboson production at large jet multiplicities. In the muon channel, the SM background contribution increases from 2% to 10% as the jet multiplicity increases, dominated by tt and diboson processes. Table I shows, for the electron and muon analyses separately, the observed number of events for the different jet multiplicities in the final state compared to predictions for signal and background processes.

VIII. UNCORRECTED DISTRIBUTIONS
The uncorrected Z/γ * (→ e + e − )+jets and Z/γ * (→ µ + µ − )+jets data are compared to the predictions for signal and background contributions. For the signal, both ALPGEN and Sherpa predictions are considered. As an example, Fig. 1 shows, separately for the electron and muon channels, the measured dilepton invariant mass in events with at least one jet in the final state, as well as the measured uncorrected inclusive jet multiplicity. Other observables considered include: the uncorrected inclusive jet p T , y, and φ distributions; the corresponding p T , y, and φ distributions of the leading, second-leading and third-leading jet in events with at least one, two and three jets in the final state, respectively; the invariant mass of the two leading jets, m jj , and their rapidity difference, ∆y jj , their azimuthal  separation, ∆φ jj , and the angular separation in y − φ space, ∆R jj = (∆y jj ) 2 + (∆φ jj ) 2 , in events with at least two jets in the final state. In all cases, the data yields are described, within statistical uncertainties, by the MC predictions for the signal plus the estimated SM background contributions.

IX. CORRECTION FOR DETECTOR EFFECTS
The jet measurements are corrected for detector effects back to the particle level using a bin-by-bin correction procedure, based on MC simulated samples, that corrects for jet selection efficiency and resolution effects and also accounts for the efficiency of the Z/γ * selection.
The corrected measurements refer to particle level jets identified using the anti-k t algorithm with R = 0.4, for jets with p T > 30 GeV and |y| < 4.4. At particle level, the lepton kinematics in the MC generated samples include the contributions from the photons radiated within a cone of radius 0.1 around the lepton direction. The measured cross sections are defined in a limited kinematic range for the Z/γ * decay products.
• Similarly, in the muon case the measurements are presented in the region: 66 GeV< m µ + µ − < 116 GeV, p µ T > 20 GeV, |η µ | < 2.4, and ∆R(jet − muon) > 0.5. The ALPGEN samples for Z/γ * +jets processes provide a satisfactory description of both lepton and jet distributions in data and are employed to compute the correction factors. For each observable ξ the bin-by-bin correction factors U (ξ) are defined as the ratio between the simulated distribution, after all selection criteria are applied, and the corresponding distribution at the particle level defined in a limited fiducial kinematic region for the generated leptons and jets, as detailed above.
Correction factors are considered for the following measurements: the inclusive jet multiplicity, p T and |y| distributions; the p T and |y| distributions for the leading-and second-leading jets in events with at least one and two jets, respectively; and the invariant mass and angular separation distributions in the inclusive dijet sample. Typical correction factors are about 1.40 for the electron channel and about 1.15 for the muon channel (see below), where the difference is mainly attributed to the identification of the Z boson candidate in the final state.
The measured differential cross sections are defined as functions of a given ξ: where, for each bin in ξ, N data and N backg denote the number of entries (events or jets) observed in data and the background prediction, respectively, ∆ξ is the bin width, U (ξ) is the correction factor, and L is the total integrated luminosity. The bin widths were chosen to be commensurate with the resolution, with typical correct-bin purities above 70%, and the cross section measurements are limited to bins in ξ that contain at least ten entries in the data.
A. Correction factors in the Z/γ * → e + e − channel In the case of the inclusive jet multiplicity, the correction factors vary with the number of jets and are between 1.40 and 1.50. The correction factors for the inclusive jet p T distribution and the p T distribution for the leading jet vary from 1.45 at p T around 30 GeV and 1.50 at p T about 60 GeV to 1.42 at very large p T . The corresponding factors for the p T distribution of the second-leading jet increase from about 1.40 to 1.55 with increasing p T .
The correction factors for the inclusive |y| distribution and the |y| distribution of the leading jet vary from 1.40 for central jets to about 1.60 for very forward jets. The correction factors for the |y| distribution of the second-leading jets are about 1.45 and show a mild rapidity dependence.
The correction factors for the ∆y, ∆φ, and ∆R distributions between the two leading jets increase from 1.30 to 1.50 as the jet separation increases. Finally, the correction factor for the dijet invariant mass distribution varies between 1.40 and 1.55 as m jj increases from 60 GeV to 300 GeV. At very low m jj , the correction factors are about 0.90 and reflect a large sensitivity to the p T thresholds applied in the analysis. Therefore, the cross section as a function of m jj is only reported for m jj > 60 GeV.

B. Correction factors in the
The correction factors for the inclusive jet multiplicity decrease from 1.15 to 1.08 with increasing N jet . The correction factors for the different p T distributions increase from 1.10 to 1.20 as p T increases from 30 GeV to 50 GeV and present a mild p T dependence for p T > 50 GeV. Similarly, the corresponding factors for the different jet |y| distributions vary between 1.15 for central jets and 1.20 for forward jets.
The correction factors for the ∆y, ∆φ, and ∆R distributions, for the two leading jets in events with at least two jets in the final state, vary between 1.10 and 1.20 as the jet separation increases. The correction factors for the m jj distribution vary between 1.10 and 1.20 as m jj increases. As in the electron case, the cross section as a function of m jj is limited to the region m jj > 60 GeV.

X. STUDY OF SYSTEMATIC UNCERTAINTIES
A detailed study of systematic uncertainties is carried out. In the following, a complete description is given for two of the observables: the inclusive cross section as a function of N jet and the inclusive jet cross section as a function of p T , in events with at least one jet in the final state (see Fig. 2). The same sources of systematic uncertainty are considered for the rest of the observables.
• The measured jet energies are increased and decreased by factors between 3% and 10%, depending on p T and η, to account for the absolute jet energy scale (JES) uncertainty, as determined in inclusive jet studies [29]. For a given jet |η|, the jet energy uncertainty tends to decrease with increasing p T , while the uncertainties increase with increasing |η|. An additional 0.1% to 1.5% uncertainty on the jet energy, depending on p T and |η|, is considered for each additional reconstructed primary vertex in the event to account for the uncertainty on the pileup offset subtraction, where the uncertainty decreases (increases) with increasing p T (|η|). Additional uncertainties are included to account for the different quark-and gluon-jet relative population in multi-jets and Z/γ * +jets processes and the presence of close-by jets in the final state, leading to a different average calorimeter response. These effects added in quadrature result in an uncertainty on the measured cross sections that increases from 7% to 22% as N jet increases and from 8% to 12% as p T increases, and constitutes the dominant source of systematic uncertainty for each of the measured distributions. The uncertainty on the jet energy resolution (JER) [29] translates into a 1% uncertainty on the cross section as a function of N jet and into a 1% to 3% uncertainty on the measured cross sections with increasing jet p T and |y|.
• The uncertainty on the estimated multi-jets background in the electron channel translates into an uncertainty on the measured cross sections which rises from 0.6% to 2% as N jet and p T increase. In addition, the background contributions from top quark, W +jets, Z/γ * (→ τ + τ − )+jets, and diboson production processes are varied by +7 −9.6 %, 5%, 5%, and 5%, respectively, to account for the uncertainty on the absolute normalization of the different MC samples. This translates into a less than 1% uncertainty in the measured cross sections. In the Z/γ * (→ µ + µ − )+jets measurements, the impact from the background uncertainties is negligible.
• The correction factors are re-computed using Sherpa instead of ALPGEN to account for possible dependencies on the parton shower, underlying event and fragmentation models, and the PDF sets used in the MC samples. This introduces an uncertainty on the measured cross sections that increases from 0.4% to 4.5% with increasing N jet and p T . In addition, a Bayesian iterative method [31] is used to unfold the data, which accounts for the full migration matrix across bins for a given observable. The ALPGEN MC samples are used to construct the input migration matrices for the different measured distributions and up to three iterations are considered, as optimized separately for each observable using the simulation. The differences with respect to the nominal bin-by-bin correction factors are less than 1% except at very large p T where they vary between 3% and 6%, and are included as an additional source of systematic uncertainty. Altogether, this introduces an uncertainty on the measured cross sections that increases from 0.7% to 7% with increasing N jet and p T .
• The uncertainty on the electron selection is taken into account. It includes uncertainties on the electron absolute energy scale and energy resolution, the uncertainty on the electron identification efficiency, and the uncertainties on the electron reconstruction scale factors applied to the MC simulation. This translates into a 4% uncertainty in the measured Z/γ * (→ e + e − )+jets cross sections, approximately independent of N jet , and jet p T and η.
The uncertainty on the measured cross sections due to the determination of the electron trigger efficiency is negligible.
• The uncertainty on the muon reconstruction efficiency, the muon momentum scale, and the muon momentum resolution translate into a conservative 2% uncertainty in the measured Z/γ * (→ µ + µ − )+jets cross sections, approximately independent of N jet , and jet p T and η. The uncertainty on the muon trigger efficiency introduces a less than 1% uncertainty on the measured cross sections.
For each channel, the different sources of systematic uncertainty are added in quadrature to the statistical uncertainty to obtain the total uncertainty. The total systematic uncertainty increases from 9% to 23% as N jet increases; and from 10% at low p T to 13% at very high p T . Finally, the additional 3.4% uncertainty on the total integrated luminosity [32] is also taken into account.

XI. NEXT-TO-LEADING ORDER PQCD PREDICTIONS
NLO pQCD predictions for Z/γ * (→ e + e − )+jets and Z/γ * (→ µ + µ − )+jets production are computed using the BlackHat program [5]. CTEQ6.6 PDFs [16] are employed and renormalization and factorization scales are set to µ = H T /2, where H T is defined event-by-event as the scalar sum of the p T of all particles and partons in the final state. The anti-k t algorithm with R = 0.4 is used to reconstruct jets at the parton level.
Systematic uncertainties on the predictions related to PDF uncertainties are computed using the Hessian method [33] and are defined as 90% confidence level uncertainties. For the total cross sections, they increase from 2% to 5% with increasing N jet . Additional changes in the PDFs due to the variation of the input value for α s (M Z ) by ±0.002 around its nominal value α s (M Z ) = 0.118 introduce uncertainties on the measured cross sections that increase from 2% to 7% with increasing N jet . These are added in quadrature to the PDF uncertainties. Variations of the renormalization and factorization scales by a factor of two (half) reduce (increase) the predicted cross sections by 4% to 14% as N jet increases.
The theoretical predictions are corrected for QED radiation effects. The correction factors δ QED are determined using ALPGEN MC samples with and without photon radiation in the final state, defined by the lepton four-momentum and photons within a cone of radius 0.1 around the lepton direction. The correction factors are about 2% for the electron and muon channels, and do not present a significant N jet dependence.
The theoretical predictions include parton-to-hadron correction factors δ had that approximately account for nonperturbative contributions from the underlying event and fragmentation into particles. In each measurement, the correction factor is estimated using HERWIG+JIMMY MC samples, as the ratio at the particle level between the nominal distribution and the one obtained by turning off both the interactions between proton remnants and the cluster fragmentation in the MC samples. The non-perturbative correction factors for the inclusive N jet and p T distributions are about 0.99 and exhibit a moderate N jet and p T dependence. However, for very forward jets δ had is about 0.9. The non-perturbative corrections are also computed using PYTHIA-AMBT1 MC samples with different parton shower, fragmentation model, and UE settings. The uncertainty on δ had , defined as the difference between the results obtained with HERWIG/JIMMY-AUET1 and PYTHIA-AMBT1, varies between 2% and 5%.

XII. RESULTS
As mentioned in Section IX, the measured cross sections refer to particle level jets identified using the anti-k t algorithm with R = 0.4, for jets with p T > 30 GeV and |y| < 4.4, and the results are defined in a limited kinematic range for the Z/γ * decay products. The data are compared to the predictions from the different MC event generators implementing Z/γ * (→ e + e − )+jets and Z/γ * (→ µ + µ − )+jets production, as discussed in Section IV, as well as to NLO pQCD predictions, as discussed in Section XI. Tabulated values of the results are available in Tables II to XIII and in the Durham HEP database [34].
A. Inclusive jet multiplicity Figure 3 presents the measured cross sections as functions of the inclusive jet multiplicity (≥ N jet ) for Z/γ * → e + e − and Z/γ * → µ + µ − interactions, in events with up to at least four jets in the final state. The data are well described by the predictions from ALPGEN and Sherpa, and BlackHat NLO pQCD. ALPGEN and Sherpa predictions include a 5% uncertainty from the NNLO pQCD normalization, as discussed in Section IV, and the systematic uncertainty on the BlackHat NLO pQCD predictions is discussed in Section XI. In the case of PYTHIA, the LO pQCD (qq → Z/γ * g and qg → Z/γ * q processes) MC predictions are multiplied by a factor 1.19, as determined from data and extracted from the average of electron and muon results in the ≥ 1 jet bin in Fig. 3. This brings the PYTHIA predictions close to the data. However, for larger N jet , and despite the additional normalization applied, PYTHIA predictions underestimate the measured cross sections.
The measured ratio of cross sections for N jet and N jet − 1 is shown in Fig. 4, compared to the different theoretical predictions. This observable cancels part of the systematic uncertainty and constitutes an improved test of the SM. The ratio is sensitive to the value of the strong coupling, and to the details of the implementation of higher-order matrix elements and soft-gluon radiation contributions in the theoretical predictions. The data indicate that the cross sections decrease by a factor of five with the requirement of each additional jet in the final state. The electron and muon measurements are well described by ALPGEN and Sherpa, and the BlackHat NLO pQCD predictions. PYTHIA predictions underestimate the measured ratios.
B. dσ/dpT and dσ/d|y| The inclusive jet differential cross section dσ/dp T as a function of p T is presented in Fig. 5, for both electron and muon analyses, in events with at least one jet in the final state. The cross sections are divided by the corresponding inclusive Z/γ * cross section times branching ratio σ Z/γ * →ℓ + ℓ − (ℓ = e, µ), separately for Z/γ * → e + e − and Z/γ * → µ + µ − , measured in the same kinematic region for the leptons and consistent with the results in Ref. [30], with the aim of cancelling systematic uncertainties related to lepton identification and the luminosity. The measured differential cross sections decrease by more than two orders of magnitude as p T increases between 30 GeV and 180 GeV. The data are well described by ALPGEN and Sherpa, and the BlackHat NLO pQCD predictions. PYTHIA predictions include the multiplicative factor 1.19 (as described above) and are then divided by the measured σ Z/γ * →ℓ + ℓ − cross sections in this analysis. This results in total normalization factors (×0.0028 pb −1 ) and (×0.0027 pb −1 ) for the electron and muon channels, respectively. PYTHIA shows a slightly softer jet p T spectrum than the data. Similar conclusions are extracted from Fig. 6, where the differential cross sections are presented as a function of the leading-jet p T . Figure 7 shows the measured differential cross sections (1/σ Z/γ * →ℓ + ℓ − )dσ/dp T , for electron and muon channels, as a function of p T of the second leading jet for jets with 30 GeV < p T < 120 GeV, in events with at least two jets in the final state. The measured cross sections decrease with increasing p T , and are again well described by ALPGEN and Sherpa, and the BlackHat NLO pQCD predictions, while PYTHIA does not describe the data. This is expected since PYTHIA only implements pQCD matrix elements for Z/γ * +1 jet production, with the additional parton radiation produced via parton shower.
Inclusive jet differential cross sections (1/σ Z/γ * →ℓ + ℓ − )dσ/d|y| as a function of |y| for jets with p T > 30 GeV are presented in Fig. 8, while Fig. 9 shows the jet measurements as a function of the rapidity of the leading jet. The measured cross sections decrease with increasing |y| and are well described by ALPGEN and the BlackHat NLO pQCD predictions. Sherpa provides a good description of the data in the region |y| < 3.5 but predicts a slightly larger cross section than observed in data for very forward jets. PYTHIA provides a good description of the shape of the measured cross sections in the region |y| < 2.5 but predicts a smaller cross section than the data in the forward region. In Fig. 10, the measured differential cross sections are presented as functions of the |y| of the second leading jet, for events with at least two jets in the final state. The data are described by the predictions from ALPGEN and Sherpa, and BlackHat NLO pQCD, while again PYTHIA does not describe the data.
C. dσ/dm jj The measured differential cross sections (1/σ Z/γ * →ℓ + ℓ − )dσ/dm jj as a function of the invariant mass of the two leading jets in the event for 60 GeV < m jj < 300 GeV are presented in Fig. 11 for both electron and muon channels. The shape of the measured cross section at low m jj is affected by the jet p T threshold in the cross section definition. For m jj > 100 GeV, the measured cross sections decrease with increasing m jj . The measurements are well described by ALPGEN and Sherpa, and the BlackHat NLO pQCD predictions. PYTHIA approximately reproduces the shape of the measured distribution but underestimates the measured cross sections.
D. dσ/d|∆y jj |, dσ/d|∆φ jj |, and dσ/d∆R jj Inclusive dijet cross sections are also measured as a function of the spatial separation of the two leading jets in the final state. Figure 12 shows the measured differential cross section as a function of the rapidity separation of the jets (1/σ Z/γ * →ℓ + ℓ − )dσ/d|∆y jj |, for both the electron and muon analysis, compared to the different predictions. The measured differential cross sections as a function of the azimuthal separation between jets (1/σ Z/γ * →ℓ + ℓ − )dσ/d|∆φ jj | are presented in Fig. 13, and Fig. 14 shows the measured differential cross sections (1/σ Z/γ * →ℓ + ℓ − ) dσ/d∆R jj as a function of the angular separation ∆R jj between the two leading jets in the event. The measurements are well described by ALPGEN and Sherpa, and the BlackHat NLO pQCD predictions, while PYTHIA underestimates the measured cross sections. In particular, PYTHIA underestimates the data for large |∆φ jj | values and for those topologies corresponding to well-separated jets.

E. Combination of electron and muon results
The measured cross section distributions for the Z/γ * (→ e + e − )+jets and Z/γ * (→ µ + µ − )+jets analyses are combined. In this case, the results are not normalized by the inclusive Z/γ * cross section after the combination, with the aim to present also precise absolute jet cross section measurements.
As already discussed, the electron and muon measurements are performed in different fiducial regions for the rapidity of the leptons in the final state. In addition, the QED radiation effects are different in both channels. For each measured distribution, bin-by-bin correction factors, as extracted from ALPGEN Z/γ * (→ e + e − )+jets and Z/γ * (→ µ + µ − )+jets MC samples, are used to extrapolate the measurements to the region p T > 20 GeV and |η| < 2.5 for the leptons, where the lepton kinematics are defined at the decay vertex of the Z boson. The increased acceptance in the lepton rapidities translates into about a 14% and a 5% increase of the measured cross sections in the electron and muon channels, respectively. As already mentioned in Section XI, the correction for QED effects increases the cross sections by about 2%. The uncertainties on the acceptance corrections are at the per mille level, as determined by using Sherpa instead of ALPGEN, and by considering different PDFs among the CTEQ6.6 and MSTW sets. A χ 2 test is performed for each observable to quantify the agreement between the electron and muon results before they are combined, where the statistical and uncorrelated uncertainties are taken into account. The statistical tests lead to probabilities larger than 60% for the electron and muon measurements to be compatible with each other, consistent with slightly conservative systematic uncertainties.
The electron and muon results are combined using the BLUE (Best Linear Unbiased Estimate) [35] method, which considers the correlations between the systematic uncertainties in the two channels. The uncertainties related to the trigger, the lepton reconstruction, and the multi-jets background estimation are considered uncorrelated between the two channels, while the rest of the systematic uncertainties are treated as fully correlated. Figures 15 to 26 show the combined results, and Tables II to XIII collect the final measurements for the electron and muon channels and their combination, together with the multiplicative parton-to-hadron correction factors δ had applied to the BlackHat NLO pQCD predictions (see Section XI). The measurements are well described by the BlackHat NLO pQCD predictions, and by the predictions from ALPGEN and Sherpa. The corresponding χ 2 tests relative to the different predictions, performed separately in each channel and for each observable, lead to χ 2 per degree of freedom values in the range between 0.05 and 2.70. Further details of the combination and the χ 2 tests are presented in the Appendix.

XIII. SUMMARY
In summary, results are reported for inclusive jet production in Z/γ * → e + e − and Z/γ * → µ + µ − events in proton-proton collisions at √ s = 7 TeV. The analysis considers the data collected by the ATLAS detector in 2010 corresponding to a total integrated luminosity of about 36 pb −1 . Jets are defined using the anti-k t algorithm with R = 0.4 and the measurements are performed for jets in the region p T > 30 GeV and |y| < 4.4. Cross sections are measured as a function of the inclusive jet multiplicity, and the transverse momentum and rapidity of the jets in the final state. Measurements are also performed as a function of the dijet invariant mass and the angular separation between the two leading jets in events with at least two jets in the final state. The measured cross sections are well described by NLO pQCD predictions including non-perturbative corrections, as well as by predictions of LO matrix elements of up to 2 → 5 parton scatters, supplemented by parton showers, as implemented in the ALPGEN and Sherpa MC generators.

Appendix -Combined Results
The results for the electron and muon channels are extrapolated to a common acceptance region p T > 20 GeV and |η| < 2.5 for the kinematics of the leptons, defined at the decay vertex of the Z boson before QED radiation. For each bin in a given observable ξ, the measured cross section σ fiducial ξ in each channel is corrected according to where δ QED corrects for QED radiation effects back to the Born level and A extrapolates the result to the new lepton acceptance region. Tables XIV to XVII present the correction factors applied to the measured cross sections, separately for the electron and muon analyses. The results are then combined using the BLUE [35] method that takes into account the correlations between systematic uncertainties in the two channels. The method assumes Gaussian χ 2 distributions and is not directly able to treat the asymmetric systematic uncertainties present in the measured cross sections. Therefore, a modified asymmetric iterative BLUE method is employed.
Three separate BLUE combinations are computed, using as an input the upper, the lower, and the average of the upper and lower uncertainties in the electron and muon channels, leading to three different results here denoted as σ up ξ ± ∆σ up ξ , σ low ξ ± ∆σ low ξ , and σ ave ξ ± ∆σ ave ξ , respectively. The central value for the combined cross section σ ξ , and its upper and lower uncertainties, ∆ + σ ξ and ∆ − σ ξ respectively, are given by the expressions with The BLUE method provides uncertainties on the combined measurement that include both statistical and systematic uncertainties. Finally, χ 2 tests to the data points in each measured cross section before and after extrapolation are performed with respect to the NLO pQCD, ALPGEN, and Sherpa predictions, according to where d j is the measured data point j, t j (s) is the corresponding prediction, ands denotes the vector of standard deviations, s i , for the different independent sources of systematic uncertainty in data and theory, which are considered fully correlated across bins. For each measurement considered, the sums above run over the total number of data points and seven independent sources of systematic uncertainty, and the correlations among systematic uncertainties are taken into account in t j (s). The average of the upper and lower uncertainties in data and theory are employed, and the χ 2 is minimized with respect tos. The results of the χ 2 tests are tabulated in Tables XVIII to XX.    XVI: Multiplicative correction factors, applied to the data in the electron and muon channels, that extrapolate the measured cross sections to the common acceptance region pT > 20 GeV and |η| < 2.5 for the lepton kinematics, defined at the decay vertex of the Z boson before QED radiation.