Contents

ZplusJetsPaper.bib \addbibresourcebibtex/bib/ATLAS.bib \addbibresourcebibtex/bib/ConfNotes.bib \AtlasTitleMeasurements of the production cross section of a \Zbosonboson in association with jets in pp collisions at  TeV with the ATLAS detector \AtlasRefCodeSTDM-2016-01 \PreprintIdNumberCERN-EP-2016-297 \AtlasJournalRefEur. Phys. J. C77 (2017) 361 \AtlasDOI10.1140/epjc/s10052-017-4900-z \AtlasAbstract Measurements of the production cross section of a \Zbosonboson in association with jets in proton–proton collisions at  TeV are presented, using data corresponding to an integrated luminosity of  fb collected by the ATLAS experiment at the CERN Large Hadron Collider in 2015. Inclusive and differential cross sections are measured for events containing a \Zbosonboson decaying to electrons or muons and produced in association with up to seven jets with  GeV and . Predictions from different Monte Carlo generators based on leading-order and next-to-leading-order matrix elements for up to two additional partons interfaced with parton shower and fixed-order predictions at next-to-leading order and next-to-next-to-leading order are compared with the measured cross sections. Good agreement within the uncertainties is observed for most of the modelled quantities, in particular with the generators which use next-to-leading-order matrix elements and the more recent next-to-next-to-leading-order fixed-order predictions.
\pdfstringdefDisableCommands\pdfstringdefDisableCommands

## 1 Introduction

The measurement of the production of a \Zbosonboson111Throughout this paper, \Zg-boson production is denoted simply by \Zboson-boson production. in association with jets, \Zj, constitutes a powerful test of perturbative quantum chromodynamics (QCD) [QCD1, QCD2]. The large production cross section and easily identifiable decays of the \Zbosonboson to charged leptonic final states offer clean experimental signatures which can be precisely measured. Such processes also constitute a non-negligible background for studies of the Higgs boson and in searches for new phenomena; typically in these studies, the multiplicity and kinematics of the jets are exploited to achieve a separation of the signal of interest from the Standard Model (SM) \Zjprocess. These quantities are often measured in control regions and subsequently extrapolated to the signal region with the use of Monte Carlo (MC) generators, which are themselves subject to systematic uncertainty and must be tuned and validated using data. Predictions from the most recent generators combine next-to-leading-order (NLO) multi-leg matrix elements with a parton shower (PS) and a hadronisation model. Fixed-order parton-level predictions for \Zjproduction at next-to-next-to-leading order (NNLO) are also available [Boughezal:2015ded, Boughezal:2016isb, Ridder:2015dxa, Ridder:2016nkl].

The \Zjproduction differential cross section was previously measured by the ATLAS [PERF-2007-01], CMS [Chatrchyan:2008aa], and LHCb [Alves:2008zz] collaborations at the CERN Large Hadron Collider (LHC) [1748-0221-3-08-S08001] at centre-of-mass energies of  TeV [Zjets2013, VjetsCms1, VjetsCms2, VjetsCms3, Aaij:2013nxa] and  TeV [Khachatryan:2015ira, Khachatryan:2016crw, AbellanBeteta:2016ugk], and by the CDF and D0 collaborations at the Tevatron collider at  TeV [Aaltonen:2014vma, Abazov:2009pp]. In this paper, proton–proton () collision data corresponding to an integrated luminosity of 3.16 fb, collected at  TeV with the ATLAS detector during 2015, are used for measurements of the \Zboson-boson production cross section in association with up to seven jets within a fiducial region defined by the detector acceptance. The \Zbosonboson is identified using its decays to electron or muon pairs (\Ztoee, \Ztomm). Cross sections are measured separately for these two channels, and for their combination, as a function of the inclusive and exclusive jet multiplicity and the ratio of successive inclusive jet multiplicities , the transverse momentum of the leading jet \pTjfor several jet multiplicities, the jet rapidity \yj, the azimuthal separation between the two leading jets \Dphijj, the two leading jet invariant mass \mjj, and the scalar sum \HTof the transverse momenta of all selected leptons and jets.

The paper is organised as follows. Section 2 contains a brief description of the ATLAS detector. The data and simulated samples as well as the \Zjpredictions used in the analysis are described in Section 3. The event selection and its associated uncertainties are presented in Section 4, while the methods employed to estimate the backgrounds are shown in Section 5. Comparisons between data and Monte Carlo predictions for reconstructed distributions are found in Section 6, while the unfolding procedure is described in Section 7. Section 8 presents the analysis results, the comparisons to predictions, and a discussion of their interpretation. Conclusions are provided in Section 9.

## 2 The ATLAS detector

The ATLAS experiment at the LHC is a multi-purpose particle detector with a forward-backward symmetric cylindrical geometry and nearly coverage in solid angle.222ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the -axis along the beam pipe. The -axis points from the IP to the centre of the LHC ring, and the -axis points upwards. Cylindrical coordinates are used in the transverse plane, being the azimuthal angle around the -axis. The pseudorapidity is defined in terms of the polar angle as . Angular distance is measured in units of . When dealing with massive jets and particles, the rapidity is used, where is the jet/particle energy and is the -component of the jet/particle momentum. It consists of an inner tracking detector, electromagnetic and hadronic calorimeters, and a muon spectrometer. The inner tracker is surrounded by a thin superconducting solenoid magnet and provides precision tracking of charged particles and momentum measurements in the pseudorapidity range . This region is matched to a high-granularity electromagnetic (EM) sampling calorimeter covering the pseudorapidity range , and a coarser granularity calorimeter up to . The hadronic calorimeter system covers the entire pseudorapidity range up to . The muon spectrometer consists of three large superconducting toroids each containing eight coils, a system of trigger chambers, and precision tracking chambers, which provide trigger and tracking capabilities in the range and , respectively. A two-level trigger system [Aaboud:2016leb] is used to select events. The first-level trigger is implemented in hardware and uses a subset of the detector information. This is followed by the software-based high-level trigger system, which runs offline reconstruction, reducing the event rate to approximately .

## 3 Data set, simulated event samples, and predictions

### 3.1 Data set

The data used in this analysis were collected by the ATLAS detector during August to November 2015. During this period, the LHC circulated proton beams with a bunch spacing. The peak delivered instantaneous luminosity was and the mean number of interactions per bunch crossing (hard scattering and pile-up events) was . The data set used in this measurement corresponds to a total integrated luminosity of 3.16 fb.

### 3.2 Simulated event samples

Monte Carlo simulations, normalised to higher-order calculations, are used to estimate most of the contributions from background events, to unfold the data to the particle level, and to compare with the unfolded data distributions. All samples are processed with a Geant4-based simulation [Agostinelli:2002hh] of the ATLAS detector [SOFT-2010-01]. An overview of all signal and background processes considered and of the generators used for the simulation is given in Table 1. Total production cross sections for the samples, their respective uncertainties (mainly coming from parton distribution function (PDF) and factorisation and renormalisation scale variations), and references to higher-order QCD corrections, where available, are also listed in Table 1.

Signal events (i.e. containing a \Zbosonboson with associated jets) are simulated using the \SHERPA v2.2.1 [Gleisberg:2008ta] generator, denoted by \SHERPAtwotwo. Matrix elements (ME) are calculated for up to two additional partons at NLO and up to four partons at leading order (LO) using the Comix [Gleisberg:2008fv] and OpenLoops [Cascioli:2011va] matrix element generators. They are merged with the \SHERPAparton shower [Schumann:2007mg] (with a matching scale of 20 GeV) using the ME+PS@NLO prescription [Hoeche:2012yf]. A five-flavour scheme is used for these predictions. The \NNPDFthreePDF set [Ball:2014uwa] is used in conjunction with a dedicated set of parton-shower-generator parameters (tune) developed by the \SHERPAauthors. This sample is used for the nominal unfolding of the data distributions, to compare to the cross-section measurements, and to estimate the systematic uncertainties.

A simulated sample of \Zjproduction is also produced with the \MADGRAPHaMC(denoted by \MGaMC) v2.2.2 generator [Alwall:2014hca], using matrix elements including up to four partons at leading order and employing the \NNPDFthreePDF set, interfaced to \PYTHIA v8.186 [Sjostrand:2007gs] to model the parton shower, using the CKKWL merging scheme [Lonnblad:2001iq] (with a matching scale of 30 GeV). A five-flavour scheme is used. The A14 [ATL-PHYS-PUB-2014-021] parton-shower tune is used together with the \NNPDFtwothreePDF set [Ball:2012cx]. The sample is denoted by \MGaMCPyCKKWLand is used to provide cross-checks of the systematic uncertainty in the unfolding and to model the small background. In addition, a \MGaMCsample with matrix elements for up to two jets and with parton showers beyond this, employing the \NNPDFthreePDF set and interfaced to \PYTHIA v8.186 to model the parton shower, is generated using the FxFx merging scheme [VjetFxFx] (with a matching scale of 25 GeV [VjetFxFxold]) and is denoted by \MGaMCPyFxFx. This sample also uses a five-flavour scheme and the A14 parton-shower tune with the \NNPDFtwothreePDF set. Both \MGaMCsamples are used for comparison with the unfolded cross-section measurements.

The measured cross sections are also compared to predictions from the leading-order matrix element generator \ALPGEN v2.14 [Mangano:2002ea] interfaced to \PYTHIA v6.426 [Sjostrand:2006za] to model the parton shower, denoted by \ALPGENPysix, using the Perugia2011C [Skands:2010ak] parton-shower tune and the CTEQ6L1 PDF set [Pumplin:2002vw]. A four-flavour scheme is used. Up to five additional partons are modelled by the matrix elements merged with the MLM prescription [Mangano:2002ea] (with a matching scale of 20 GeV). The matrix elements for the production of and events are explicitly included and a heavy-flavour overlap procedure is used to remove the double counting of heavy quarks from gluon splitting in the parton shower.

The \Zboson-boson samples are normalised to the NNLO prediction calculated with the Fewz 3.1 program [Anastasiou:2003ds, FEWZ2, FEWZ3, FEWZ4] with CT10nnlo PDFs [Lai:2010vv].

Contributions from the top-quark, single-boson, and diboson components of the background (described in Section 5) are estimated from the following Monte Carlo samples. Samples of top-quark pair and single top-quark production are generated at NLO with the \POWHEG-Box generator [Nason:2004rx, Frixione:2007vw, Alioli:2010xd] (versions v2 (r3026) for top-quark pairs and v1 for single top quarks (r2556 and r2819 for -and -channels, respectively)) and \PYTHIA v6.428 (Perugia2012 tune [Skands:2010ak]). Samples with enhanced or suppressed additional radiation are generated with the Perugia2012radHi/Lo tunes [Skands:2010ak]. An alternative top-quark pair sample is produced using the \MGaMCgenerator interfaced with \HERWIGpp v2.7.1 [Alwall:2014hca, herwigpp], using the UE-EE-5 tune [Gieseke:2012ft]. The samples are normalised to the cross section calculated at NNLO+NNLL (next-to-next-to-leading log) with the Top++2.0 program [Czakon:2011xx].

The \Wboson-boson backgrounds are modelled using the \MGaMCPyCKKWL v2.2.2 generator, interfaced to \PYTHIA v8.186 and are normalised to the NNLO values given in Table 1. Diboson processes with fully leptonic and semileptonic decays are simulated [MultiBosonPub] using the \SHERPA v2.1.1 generator with the CT10nlo PDF set. The matrix elements contain the doubly resonant , and processes, and all other diagrams with four electroweak vertices. They are calculated for one or zero additional partons at NLO and up to three additional partons at LO and merged with the \SHERPAparton shower using the ME+PS@NLO prescription.

Events involving semileptonic decays of heavy quarks, hadrons misidentified as leptons, and, in the case of the electron channel, electrons from photon conversions are referred to collectively as “multijet events”. The multijet background was estimated using data-driven techniques, as described in Section 5.

Multiple overlaid collisions are simulated with the soft QCD processes of \PYTHIA v.8.186 using the A2 tune [ATLAS-PHYS-PUB-2012-003] and the MSTW2008LO PDF set [Martin:2009iq]. All Monte Carlo samples are reweighted so that the pile-up distribution matches that observed in the data.

### 3.3 Fixed-order predictions

In addition to these Monte Carlo samples, parton-level fixed-order predictions at NLO are calculated by the \BLACKHATSHERPAcollaboration for the production of \Zbosonbosons with up to four partons [bhz3jets, bhz4jets]. The \BLACKHATSHERPApredictions use the CT14 PDF set [Dulat:2015mca] including variations of its eigenvectors at the 68% confidence level, rescaled from 90% confidence level using a factor of . The nominal predictions use a factorisation and renormalisation scale of with uncertainties derived from the envelope of a common variation of the scales by factors of and 2. The effects of PDF and scale uncertainties range from 1% to 4% and from 0.1% to 10%, respectively, for the cross sections of \Zboson-boson production in association with at least one to four partons, and are included in the predictions which are provided by the \BLACKHATSHERPAauthors for the fiducial phase space of this analysis. Since these predictions are defined before the decay leptons emit photons via final-state radiation (Born-level FSR), corrections to the dressed level (where all photons found within a cone of size of the lepton from the decay of the \Zbosonboson are included) are derived from \MGaMCPyCKKWL, separately for each kinematic observable used to measure cross sections, with associated systematic uncertainties obtained by comparing to the \ALPGENPysixgenerator. This correction is needed in order to match the prediction to the lepton definition used in the measurements. The average size of these corrections is approximately %. To bring the prediction from parton to particle level, corrections for the non-perturbative effects of hadronisation and the underlying event are also calculated separately for each observable using the \SHERPA v2.2 generator by turning on and off in the simulation both the fragmentation and the interactions between the proton remnants. The net size of the corrections is up to approximately 10% at small values of \pTjand vanishes for large values of \pTj. An uncertainty of approximately 2% for this correction is included in the total systematic uncertainty of the prediction.

Calculations of cross sections at NNLO QCD have recently become available [Boughezal:2015ded, Boughezal:2016isb, Ridder:2015dxa, Ridder:2016nkl]. In this paper, the results are compared to the calculation, denoted by \Njetti[Boughezal:2015ded, Boughezal:2016isb], which uses a new subtraction technique based on -jettiness [Boughezal:2015dva] and relies on the theoretical formalism provided in soft-collinear effective theory. The predictions, which are provided by the authors of this calculation for the fiducial phase space of this analysis, use a factorisation and renormalisation scale of (where is the invariant mass of the dilepton system) and the CT14 PDF set. The QCD renormalisation and factorisation scales were jointly varied by a common factor of two, and are included in the uncertainties. Non-perturbative and FSR corrections and their associated uncertainties as discussed above are also included in the predictions.

## 4 Event selection

Electron- and muon-candidate events are selected using triggers which require at least one electron or muon with transverse momentum thresholds of or , respectively, with some isolation requirements for the muon trigger. To recover possible efficiency losses at high momenta, additional electron and muon triggers which do not make any isolation requirements are included with thresholds of and , respectively. Candidate events are required to have a primary vertex, defined as the vertex with the highest sum of track , with at least two associated tracks with .

Electron candidates are required to have and to pass “medium” likelihood-based identification requirements [elecperf, ATL-PHYS-PUB-2015-041] optimised for the 2015 operating conditions, within the fiducial region of , excluding candidates in the transition region between the barrel and endcap electromagnetic calorimeters, . Muons are reconstructed for with and must pass “medium” identification requirements [Aad:2016jkr] also optimised for the 2015 operating conditions. At least one of the lepton candidates is required to match the lepton that triggered the event. The electrons and muons must also satisfy \pT-dependent cone-based isolation requirements, using both tracking detector and calorimeter information (described in Refs. [Aad:2015yja] and  [Aad:2014lwa], respectively). The isolation requirements are tuned so that the lepton isolation efficiency is at least 90% for , increasing to 99% at 60 GeV. Both the electron and muon tracks are required to be associated with the primary vertex, using constraints on the transverse impact parameter significance , where is the transverse impact parameter and is its uncertainty, and on the longitudinal impact parameter corrected for the reconstructed position of the primary vertex. The transverse impact parameter significance is required to be less than five for electrons and three for muons, while the absolute value of the corrected multiplied by the sine of the track polar angle is required to be less than 0.5 mm.

Jets of hadrons are reconstructed with the anti- algorithm [Cacciari:2008] with radius parameter using topological clusters of energy deposited in the calorimeters [Aad:2016upy]. Jets are calibrated using a simulation-based calibration scheme, followed by in situ corrections to account for differences between simulation and data [ATLAS-PHYS-PUB-2015-015]. In order to reduce the effects of pile-up contributions, jets with pseudorapidity and are required to have a significant fraction of their tracks with an origin compatible with the primary vertex, as defined by the jet vertex tagger algorithm [ATLAS-CONF-2014-018]. In addition, the expected average energy contribution from pile-up clusters is subtracted according to the catchment area of the jet [Aad:2015ina]. Jets used in the analysis are required to have \pTgreater than  GeV and rapidity .

The overlap between leptons and jets is removed in a two-step process. The first step removes jets closer than to a selected electron, and jets closer than to a selected muon, if they are likely to be reconstructed from photons radiated by the muon. In a second step, electrons and muons are discarded if they are located closer than to a remaining selected jet. This requirement effectively removes events with leptons and jets which are not reliably simulated in the Monte Carlo simulation.

Events containing a \Zboson-boson candidate are selected by requiring exactly two leptons of the same flavour but of opposite charge with dilepton invariant mass in the range . The expected and observed numbers of \Zboson-boson candidates selected for each inclusive jet multiplicity, for , are summarised in Table 2, separately for the \Ztoeeand the \Ztommchannels. The background evaluation that appears in this table is discussed in Section 5. After all requirements, 248,816 and 311,183 \Zonejevents are selected in the electron and muon channels, respectively.

### 4.1 Correction factors and related systematic uncertainties

Some of the object and event selection efficiencies as well as the energy and momentum calibrations modelled by the simulation must be corrected with simulation-to-data correction factors to better match those observed in the data. These corrections and their corresponding uncertainties fall into the following two categories: dependent and not dependent on lepton flavour.

The corrections and uncertainties specific to each leptonic final state (\Ztoeeand \Ztomm) are as follows:

• Trigger: The lepton trigger efficiency is estimated in simulation, with a separate data-driven analysis performed to obtain the simulation-to-data trigger correction factors and their corresponding uncertainties [Aaboud:2016leb].

• Lepton reconstruction, identification, and isolation: The lepton selection efficiencies as determined from simulation are also corrected with simulation-to-data correction factors, with corresponding uncertainties [ATL-PHYS-PUB-2015-041, Aad:2016jkr].

• Energy, momentum scale/resolution: Uncertainties in the lepton calibrations are estimated [Aad:2016jkr] because they can cause a change of acceptance because of migration of events across the \pTthreshold and boundaries.

The corrections and uncertainties common to the electron and muon final states are as follows:

• Jet energy scale and resolution: Uncertainties in the jet energy-scale calibration and resolution have a significant impact on the measurements, especially for the higher jet multiplicities. The jet energy-scale calibration is based on 13 TeV simulation and on in situ corrections obtained from data [ATLAS-PHYS-PUB-2015-015]. The uncertainties are estimated using a decorrelation scheme, resulting in a set of 19 independent parameters which cover all of the relevant calibration uncertainties. The jet energy scale is the dominant systematic uncertainty for all bins with at least one jet. The jet energy-resolution uncertainty is derived by over-smearing the jet energy in the simulation and using the symmetrised variations as the uncertainty.

• Jet vertex tagger: The modelling of the output variable from the jet vertex tagger is validated using data events where the \Zbosonboson recoils against a jet. A percent-level correction is derived and its statistical and systematic uncertainties are used as additional uncertainties in the efficiency to select jets from the primary vertex [ATLAS-CONF-2014-018].

• Pile-up: The imperfect modelling of the effects of pile-up leads to acceptance changes at the percent level for different jet multiplicities. To assess this uncertainty, the average number of interactions per bunch crossing is varied in simulation so that the behaviour of variables sensitive to pile-up matches that observed in data.

• Luminosity: The cross sections have a 2.1% uncertainty from the measurement of the integrated luminosity, which is derived, following a methodology similar to that detailed in Refs. [DAPR-2011-01, Aaboud:2016hhf], from a calibration of the luminosity using beam-separation scans performed in August 2015.

## 5 Background estimation

Contributions from the electroweak (single boson and diboson) and top-quark (single top-quark and top-quark pair) components of the background are estimated using the Monte Carlo samples described in Section 3 with corresponding uncertainties as listed in Table 1. Contributions from multijet events are evaluated with data-driven techniques as described below. A summary of the composition and relative importance of the backgrounds in the candidate \Zjevents is given in Table 2. The overall purity of the \Zj selections (fraction of signal events in the final selection) ranges from 99% in the inclusive sample to 80-85% in the bin.

### 5.1 Top-quark and electroweak backgrounds

The dominant contribution to the background at high jet multiplicities comes from production, with the subsequent leptonic decays of the \Wbosonbosons originating from the top quarks and is evaluated from simulation. An overall uncertainty of 6%, corresponding to the PDF and scale variations on the theoretical predictions of the inclusive cross sections, is assigned (see Table 1). The background estimate is validated through a cross-section measurement of production in the dilepton channel at  TeV [ATLAS-CONF-2015-065] as a function of the jet multiplicity, and the modelling of the additional parton radiation in events by \POWHEGPysixwas found to be in good agreement with this measurement. In addition, a systematic uncertainty in the modelling of the shape of the distributions is derived by modifying the parton-shower intensity in the nominal simulation sample and by comparing to the predictions from the alternative generator \MGaMCHpp(both listed in Table 1). The small contribution from single-top-quark events is also estimated using \POWHEGPysixsamples and assigned a 6% uncertainty.

Diboson production in leptonic and semileptonic final states with at least two leptons of the same flavour constitutes a co-dominant background for high jet multiplicities (see Table 2). The production of bosons in association with jets at  TeV was found to be well modelled by the \SHERPAtwoonegenerator [Aaboud:2016yus]. A 6% uncertainty, again corresponding to PDF and scale variations on the predictions, is assessed. Since in Ref. [Aaboud:2016yus] the measurement is limited by the statistical precision for dibosons  jets (resulting in  hadronic jets for semileptonic diboson decays), an additional systematic uncertainty of 50% in the normalisation of the diboson background is added for \Zsixjs.

Minor background contributions also arise from single-\Wboson-boson production decaying to leptonic final states and from single-\Zboson-boson production in the \Ztotautauprocess, both estimated with simulation and assigned a 5% uncertainty (as given in Table 1).

### 5.2 Multijet background

Background-enriched data control regions are used to estimate the multijet contribution in both the electron and muon channels. They are constructed by loosening the lepton identification and isolation requirements. Templates are built from the dilepton invariant mass distribution, a variable that shows discrimination between multijet background and other processes in regions of its kinematic range, but is largely uncorrelated with the variables used to build the multijet control regions. The templates are subsequently normalised to events passing the \Zboson-boson signal selection.

In the electron channel, the multijet templates are built for each jet multiplicity from events with two same-charge leptons with no isolation requirement, whose identification criteria are looser than those of the signal selection, which the leptons must not satisfy. In the muon channel, the control region is similarly built from events with two leptons which are selected with looser identification requirements than the signal selection and also fail the nominal isolation requirement. In both cases, dedicated triggers better suited to this purpose are used to populate the templates. The small electroweak and top-quark contamination is subtracted using simulated events.

The normalisation of the multijet template is estimated with a log-likelihood fit to the measured dilepton invariant mass distribution for the inclusive \Zbosonselection, using templates for \Ztolland for the electroweak and top-quark background derived from simulation. The fit is performed in the invariant mass windows of  GeV and  GeV for the electron and muon channels, respectively, in order to benefit from the larger multijet contribution in the mass sidebands. The normalisation of the multijet template is allowed to float freely while the remaining non-multijet templates are constrained to be within 6% of the predicted cross sections for these processes as given in Table 1. The multijet fractions are evaluated separately for each jet multiplicity, except for very high jet multiplicities where the templates are statistically limited, and so these fractions are taken from the estimates of the and bins in the electron and muon channels, respectively.

The systematic uncertainties on the multijet background are derived by varying the mass range and bin width of the nominal fit, using the lepton transverse impact parameter as the fitting variable instead of the invariant mass, using alternative simulation samples for the templates, allowing the normalisations of the non-multijet components to vary independently or within a wider range, and varying the lepton resolution and energy/momentum scales. In addition, given the multiple sources of multijet background in the electron channel, an alternative template is constructed by requiring that the electrons fail to meet an isolation criterion instead of failing to meet the nominal signal selection electron identification criterion.

The resulting estimated multijet fractions in each jet multiplicity bin are given in Table 2. Their corresponding total uncertainties are dominated by their systematic components. These systematic components are approximately 70% of the multijet fraction as estimated in the electron and muon channels.

## 6 Kinematic distributions

The level of agreement between data and predictions is evaluated from the comparison of kinematic distributions. Figure 1, which presents the dilepton mass for the \Zonejtopology and the inclusive jet multiplicity, shows how well the \SHERPAtwotwoand \MGaMCPyCKKWLpredictions agree with data. The uncertainty bands shown in these distributions include the statistical uncertainties due to the simulation sample sizes, the event-selection uncertainties described in Section 4.1 (omitting the common 2.1% luminosity uncertainty), and the background normalisation uncertainties described in Section 5.

## 7 Unfolding of detector effects

The cross-section measurements presented in this paper are performed within the fiducial acceptance region defined by the following requirements:

•  GeV,

•  GeV,

•  GeV.

The cross sections are defined at particle (“truth”) level, corresponding to dressed electrons and muons from the \Zbosonbosons. The particle level also includes jets clustered using the anti- algorithm [Cacciari:2008] with radius parameter for final-state particles with decay length  mm, excluding the dressed \Zboson-boson decay products.

The fiducial cross sections are estimated from the reconstructed kinematic observables: jet multiplicity, \pTjfor different jet multiplicities, \yj, \Dphijj, \mjj, and \HT, for events that pass the selection described in Section 4. The expected background components as described in Section 5 are subtracted from the distributions in data. A variable-width binning of these observables is used, such that the purity is at least 50% in each bin and the size of the statistical uncertainty in most of the bins remains below 10%.

An iterative Bayesian unfolding technique [D'Agostini:1994zf], as implemented in the RooUnfold package [RooUnfold], is used to unfold the measurements to the particle level, thereby accounting for detector effects related to inefficiencies, resolution, and systematic biases in the central values of the kinematic variables describing both the leptons and the jets. The iterative unfolding technique updates the initial estimators for the generated (“truth”) distribution in consecutive steps, using Bayes’ theorem in each iteration to derive an unfolding matrix from the initial response matrix (which relates truth and reconstructed distributions of given observables) and the current truth estimator.

The response matrices are constructed using the \SHERPAtwotwo\Ztolljsamples. \SHERPAtwotwois also used to derive the initial truth estimator. In order to enter the response matrix, events must pass the \Zboson-boson selection at generator level and at detector level and contain the number of jets required by the preselection for a given observable at both generator and detector level. Reconstructed jets are required to match the corresponding generator-level jets within a cone of size for all distributions except global quantities such as the jet multiplicity and \HT. A given bin in the response matrix therefore corresponds to the probability that a true jet object in bin is reconstructed in bin of the distribution. Figure 2 illustrates two examples of response matrices. The resulting ratios of detector-level to truth-level event yields are typically 0.65 and 0.8 for the electron and muon channels, respectively.

The background-subtracted data are corrected for the expected fraction of events with reconstructed objects unmatched to any generator object before entering the iterative unfolding. The number of iterations used for the iterative unfolding of each distribution (two) is chosen by unfolding the \SHERPAtwotwosamples reweighted to data and comparing to the generated reweighted distribution. The unfolded event yields are divided by the integrated luminosity of the data sample and the bin width of the distribution in question to provide the final fiducial cross sections. The final result is given by

 σi=1ϵiL∑jUijNdataj(1−funmatchedj), (1)

where is the integrated luminosity, is the reconstruction efficiency for truth bin , corresponds to the number of events observed in data in reconstructed bin and is its fraction of unmatched events calculated from simulation, and is the unfolding matrix calculated after two iterations, using the updated prior from the first iteration and the response matrix.

### 7.1 Systematic uncertainties associated with the unfolding procedure

The limited size of a simulation sample can create biases in the distributions. Systematic uncertainties account for possible residual biases in the unfolding procedure due to, e.g. modelling of the hadronisation in the simulation, migrations into other kinematic distributions not explicitly part of the unfolding, or the finite bin width used in each distribution. The following uncertainties arise from the unfolding procedure.

• The statistical uncertainties of the response matrices derived from \SHERPAtwotwoare propagated to the unfolded cross sections with a toy simulation method. A total of 5000 ensembles (pseudo-experiments) of unfolded samples are generated. For each sample, the number of reconstructed events in each bin is generated randomly according to a Gaussian distribution, where the mean is the nominal number of events before unfolding and the width is its corresponding statistical uncertainty. Unfolding is performed for each ensemble. The widths of resulting distributions are taken as a systematic uncertainty of the unfolding.

• The \SHERPAtwotwosamples are reweighted at generator level, such that the distribution of the leading jet \ptat detector level matches that observed in the data. The modified \SHERPAtwotwosamples are then used to unfold the data again and the variations in the resulting cross sections are used to derive a systematic uncertainty.

• An additional check is performed by unfolding reconstructed \MGaMCPyCKKWLevents using \SHERPAtwotworesponse matrices. The residual non-closure is accounted for by an additional flat uncertainty of 3% for all distributions.

## 8 Results

The measured cross sections, presented in Section 8.1, are calculated in the electron and muon channels separately and the compatibility of the results of the two channels is evaluated. In order to improve the precision of the measurement, these results are then combined, taking into account the correlations of the systematic uncertainties. The comparisons of the combined results to the predictions are presented in Section 8.2.

### 8.1 Results in the individual channels and the combination

The fiducial cross-section measurements in the \Ztoeejand \Ztommjchannels as a function of the inclusive jet multiplicities are presented in Table 3. The data statistical uncertainties are propagated through the unfolding by using pseudo-experiments. As mentioned in Section 7, the systematic uncertainties are propagated through the unfolding via the migration matrices and via the variation of the subtracted background. Table 4 shows the resulting total relative statistical and systematic uncertainties as well as the systematic components (lepton trigger, lepton selection, jet energy scale and resolution, jet vertex tagging, pile-up, luminosity (all described in Section 4.1)), unfolding (described in Section 7), and background (described in Section 5) as a function of the inclusive jet multiplicity, presented separately for the electron and muon channels. The jet energy scale is the dominant systematic uncertainty for all bins with at least one jet.

Figure 3 shows a comparison of the electron and muon channels for the measured fiducial cross section as a function of the inclusive jet multiplicity and of the leading jet \pt for inclusive \Zonejevents. This figure demonstrates that the results in the electron and muon channels are compatible and hence can be combined to improve the precision of the measurement. This figure also shows the result of this combination described below.

The results from the electron and muon channels are combined at dressed level for each distribution separately: inclusive and exclusive jet multiplicities, ratio for successive inclusive jet multiplicities, leading jet \pTfor \Zonefourjevents and jet \pTfor exclusive \Zexcljevents, leading jet rapidity for inclusive \Zonejevents, \HT, \Dphijj, and \mjj. A function whose sum runs over all measurement sets (electrons and muons), measurement points, and some of the uncertainty sources, is used for the combination [Glazov:2005rn, Aaron:2009bp] and distinguishes between bin-to-bin correlated and uncorrelated sources of uncertainties, the latter comprising the statistical uncertainty of the data and the statistical unfolding uncertainty. Uncertainties specific to the lepton flavour and to the background are included in the function, while the remaining, flavour-uncorrelated, systematic uncertainties related to jets, pile-up, luminosity, and unfolding are averaged after the combination.

### 8.2 Comparisons of results to predictions

The cross-section measurement for different inclusive \Zjmultiplicities and their ratios obtained from the combination are found in Tables 5 and 6. Figure 4 shows the comparison of these results with the NLO QCD fixed-order calculations from \BLACKHATSHERPAand with the predictions from \SHERPAtwotwo, \ALPGENPysix, \MGaMCPyCKKWL, and \MGaMCPyFxFx. The plots show the particle-level cross section with the generator predictions normalised to the inclusive NNLO cross sections in the top panel, accompanied by the ratios of the various predictions with respect to the data in the bottom panels. Uncertainties from the parton distribution functions and QCD scale variations are included in the \BLACKHATSHERPApredictions, as described in Section 3.3. A constant 5% theoretical uncertainty is used for \SHERPAtwotwo, \ALPGENPysix, \MGaMCPyCKKWL, and \MGaMCPyFxFx, as described in Table 1. The inclusive jet multiplicity decreases logarithmically while the ratio is flat in the presence of at least one jet. The predictions are in agreement with the observed cross sections and their ratios, except for \SHERPAtwotwo, \ALPGENPysixand \MGaMCPyFxFxfor high jet multiplicity, where a non-negligible fraction of the jets are produced by the parton shower.

The jet transverse momentum is a fundamental observable of the \Zjprocess and probes pQCD over a wide range of scales. Moreover, understanding the kinematics of jets in events with vector bosons associated with several jets is essential for the modelling of backgrounds for other SM processes and searches beyond the SM. The leading jet \pt distribution (which is correlated with the \pTof the \Zbosonboson) in inclusive \Zonefourjevents is shown in Figure 5 and ranges up to 700 GeV. The LO generator \MGaMCPyCKKWLmodels a too-hard jet \pt spectrum. This feature is known from studies of LO generators in collisions at lower centre-of-mass energies [Zjets2013], and can be interpreted as an indication that the dynamic factorisation and renormalisation scale used in the generation is not appropriate for the full jet \pt range. In contrast, the predictions from \BLACKHATSHERPA, \SHERPAtwotwo, and \MGaMCPyFxFx, which are based on NLO matrix elements, are in agreement with the measured cross section within the systematic uncertainties over the full leading jet \pt range. \ALPGENPysixalso shows good agreement with the measured data. The \Njettiprediction models the spectrum for the \Zonejevents well. Uncertainties from the QCD scale variations for the \Njettipredictions are included in the uncertainty band, as described in Section 3.3. For the leading jet rapidity distribution in inclusive \Zonejevents, also shown in this figure, all predictions show good agreement with the measured data within the uncertainties.

The exclusive jet \pt distribution probes the validity of \Zexcljpredictions at increasing QCD scales represented by the jet \pT in the presence of a jet veto at a constant low scale; for a jet \pt range of several hundred GeV, accessible with the current data set, the jet scale is of order ten times larger than the veto scale (30 GeV). Figure 6 demonstrates that all predictions studied are consistent with the data within systematic uncertainties over the full jet \pt range (up to 500 GeV). This figure also shows the measured cross section as a function of the exclusive jet multiplicity, which decreases logarithmically. Similar trends as for the inclusive jet multiplicity (Figure 4) are observed.

Quantities based on inclusive \pTsums of final-state objects, such as \HT, the scalar \pTsum of all visible objects in the final state, are often employed in searches for physics beyond the Standard Model, to enrich final states resulting from the decay of heavy particles. The values \HTor \HT/2 are also commonly used choices for scales for higher-order perturbative QCD calculations. Large values for this quantity can result either from a small number of very energetic particles or from a large number of less energetic particles. Figure 7 shows the measured cross sections as a function of the \HTdistribution (up to 1400 GeV) in inclusive \Zonejevents. The predictions from \SHERPAtwotwo, \ALPGENPysixand \MGaMCPyFxFxdescribe well the \HTdistribution. The prediction from \MGaMCPyCKKWLdescribes well the turn-over in the softer part of the \HTspectrum, but overestimates the contribution at large values of \HT, in line with the overestimate of the cross sections for hard jets. The fixed-order \Zonejprediction from \BLACKHATSHERPAunderestimates the cross section for values of  GeV, as observed in similar measurements at lower centre-of-mass energies [Zjets2013, mySTDM-2012-24], due to the missing contributions from events with higher parton multiplicities, which for large values of \HTconstitute a substantial portion of the data. Agreement is recovered by adding higher orders in perturbative QCD, as demonstrated by the good description of \HTby \Njetti.

Angular relations between the two leading jets and the dijet mass are frequently used to separate either heavier SM particles or beyond-SM physics from the \Zjprocess. Figure 8 shows the differential cross section as a function of azimuthal angular difference between the two leading jets for \Ztwojevents, \Dphijj. The tendency of the two jets to be back-to-back in the transverse plane is well modelled by all predictions. This figure also shows the measured cross sections as a function of the invariant mass \mjjof the two leading jets for \Ztwojevents. The shape of the dijet mass is modelled well by \BLACKHATSHERPA, \SHERPAtwotwo, \ALPGENPysix, and \MGaMCPyFxFx, whereas \MGaMCPyCKKWLshows a harder spectrum.

## 9 Conclusion

Proton–proton collision data at  TeV from the LHC, corresponding to a total integrated luminosity of 3.16 fb, have been analysed by the ATLAS collaboration to study events with \Zbosonbosons decaying to electron or muon pairs, produced in association with one or more jets. The fiducial production cross sections for \Zzerosevenjhave been measured, within the acceptance region defined by  GeV, ,  GeV,  GeV, , and , with a precision ranging from 4% to 30%. Ratios of cross sections for successive jet multiplicities and cross-section measurements as a function of different key variables such as the jet multiplicities, jet \pt for exclusive  jet events, leading jet \pt for \Zonefourjevents, leading jet rapidity for \Zonejevents, , and have also been derived.

The measurements have been compared to fixed-order calculations at NLO from \BLACKHATSHERPAand at NNLO from the \Njetticalculation, and to predictions from the generators \SHERPAtwotwo, \ALPGENPysix, \MGaMCPyCKKWL, and \MGaMCPyFxFx. In general, the predictions are in good agreement with the observed cross sections and cross-section ratios within the uncertainties. Distributions which are dominated by a single jet multiplicity are modelled well by fixed-order NLO calculations, even in the presence of a jet veto at a low scale. The ME+PS generator \MGaMCPyCKKWL, which is based on LO matrix elements, models a too-hard jet spectrum, as observed in  TeV collisions. It however models well the inclusive jet multiplicity distribution over the full multiplicity range. The modelling of the jet \ptand related observables is significantly improved by the ME+PS@NLO generators \SHERPAtwotwoand \MGaMCPyFxFx, which use NLO matrix elements for up to two additional partons. The recent \Njettipredictions describe well key distributions such as the leading jet \pt and \HT.

The results presented in this paper provide essential input for the further optimisation of the Monte Carlo generators of \Zjproduction and constitute a powerful test of perturbative QCD for processes with a higher number of partons in the final state.

## Acknowledgements

We thank CERN for the very successful operation of the LHC, as well as the support staff from our institutions without whom ATLAS could not be operated efficiently.

We acknowledge the support of ANPCyT, Argentina; YerPhI, Armenia; ARC, Australia; BMWFW and FWF, Austria; ANAS, Azerbaijan; SSTC, Belarus; CNPq and FAPESP, Brazil; NSERC, NRC and CFI, Canada; CERN; CONICYT, Chile; CAS, MOST and NSFC, China; COLCIENCIAS, Colombia; MSMT CR, MPO CR and VSC CR, Czech Republic; DNRF and DNSRC, Denmark; IN2P3-CNRS, CEA-DSM/IRFU, France; SRNSF, Georgia; BMBF, HGF, and MPG, Germany; GSRT, Greece; RGC, Hong Kong SAR, China; ISF, I-CORE and Benoziyo Center, Israel; INFN, Italy; MEXT and JSPS, Japan; CNRST, Morocco; NWO, Netherlands; RCN, Norway; MNiSW and NCN, Poland; FCT, Portugal; MNE/IFA, Romania; MES of Russia and NRC KI, Russian Federation; JINR; MESTD, Serbia; MSSR, Slovakia; ARRS and MIZŠ, Slovenia; DST/NRF, South Africa; MINECO, Spain; SRC and Wallenberg Foundation, Sweden; SERI, SNSF and Cantons of Bern and Geneva, Switzerland; MOST, Taiwan; TAEK, Turkey; STFC, United Kingdom; DOE and NSF, United States of America. In addition, individual groups and members have received support from BCKDF, the Canada Council, CANARIE, CRC, Compute Canada, FQRNT, and the Ontario Innovation Trust, Canada; EPLANET, ERC, ERDF, FP7, Horizon 2020 and Marie Skłodowska-Curie Actions, European Union; Investissements d’Avenir Labex and Idex, ANR, Région Auvergne and Fondation Partager le Savoir, France; DFG and AvH Foundation, Germany; Herakleitos, Thales and Aristeia programmes co-financed by EU-ESF and the Greek NSRF; BSF, GIF and Minerva, Israel; BRF, Norway; CERCA Programme Generalitat de Catalunya, Generalitat Valenciana, Spain; the Royal Society and Leverhulme Trust, United Kingdom.

The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN, the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) and BNL (USA), the Tier-2 facilities worldwide and large non-WLCG resource providers. Major contributors of computing resources are listed in Ref. [ATL-GEN-PUB-2016-002].

\printbibliography

The ATLAS Collaboration