Measurement of ZZ production in pp collisions at sqrt(s)=7 TeV and limits on anomalous ZZZ and ZZgamma couplings with the ATLAS detector

A measurement of the ZZ production cross section in proton-proton collisions at sqrt(s) = 7 TeV using data recorded by the ATLAS experiment at the Large Hadron Collider is presented. In a data sample corresponding to an integrated luminosity of 4.6 fb-1 collected in 2011, events are selected that are consistent either with two Z bosons decaying to electrons or muons or with one Z boson decaying to electrons or muons and a second Z boson decaying to neutrinos. The ZZ*->llll and ZZ->llnunu cross sections are measured in restricted phase-space regions. These results are then used to derive the total cross section for ZZ events produced with both Z bosons in the mass range 66 to 116 GeV, sigmaZZtot = 6.7 +-0.7 +0.4-0.3 +-0.3 pb, which is consistent with the Standard Model prediction of 5.89+0.22-0.18 pb calculated at next-to-leading order in QCD. The normalized differential cross sections in bins of various kinematic variables are presented. Finally, the differential event yield as a function of the transverse momentum of the leading Z boson is used to set limits on anomalous neutral triple gauge boson couplings in ZZ production.


Introduction
The production of pairs of Z bosons at the Large Hadron Collider (LHC) provides an excellent opportunity to test the predictions of the electroweak sector of the Standard Model (SM) at the TeV energy scale.In the SM, Z boson pairs can be produced via non-resonant processes or in the decay of Higgs bosons.Deviations from SM expectations for the total or differential ZZ production cross sections could be indicative of the production of new resonances decaying to Z bosons or other non-SM contributions.
Non-resonant ZZ production proceeds at leading order (LO) via tand u-channel quark-antiquark interactions, while about 6% of the production proceeds via gluon fusion.The ZZZ and ZZγ neutral triple gauge boson couplings (nTGCs) are absent in the SM, hence there is no contribution from schannel q q annihilation at tree level.These different production processes are shown in figure 1.At the one-loop level, nTGCs generated by fermion triangles have a magnitude of the order of 10 −4 [1].Many models of physics beyond the Standard Model predict values of nTGCs at the level of 10 −4 to 10 −3 [2].The primary signatures of non-zero nTGCs are an increase in the ZZ cross section at high ZZ invariant mass and high transverse momentum of the Z bosons [3].ZZ production has been studied in e + e − collisions at LEP [4][5][6][7][8], in pp collisions at the Tevatron [9-12] and recently in pp collisions at the LHC [13,14].No deviation of the measured total cross section from the SM expectation has been observed, and limits on anomalous nTGCs have been set [8,9,13,14].In searching for the SM Higgs boson, the ATLAS and CMS collaborations observed recently a neutral boson resonance with a mass around 126 GeV [15][16][17].A SM Higgs boson with that mass can decay to two Z bosons only when at least one Z boson is off-shell, and even in this case, the contribution is less than 3%.Searches for high-mass non-SM ZZ resonances have not resulted in any excess above the SM expectations [18].Quark and Gluon Diboson Production Diagrams

t-Channel Diagram
Quark and Gluon Diboson Production Diagrams Figure 1.Leading order Feynman diagrams for ZZ production through the q q and gg initial state at hadron colliders.The s-channel diagram, (c), contains the ZZZ and ZZγ neutral TGC vertices which do not exist in the SM.
in the four-charged-lepton channel: an on-shell ZZ selection denoted by ZZ → + − + − where both Z bosons are required to be within the mass range 66-116 GeV3 and a selection which includes an off-shell Z boson denoted by ZZ * → + − + − where one Z boson is required to be within this mass range and the other can be off-shell and have any mass above 20 GeV.In the + − ν ν channel, the ν ν system is expected to be produced by an off-shell Z boson in 2.6% of the events.Since this fraction is small and only one event selection is used for this channel, it is referred to as ZZ → + − ν ν throughout the paper.The ZZ ( * ) → + − + − channel has an excellent signal-to-background ratio, but it has a branching fraction six times lower than the ZZ → + − ν ν channel; the latter has higher background contributions with an expected signal-to-background ratio around one (after applying the event selections described below).This paper presents the total ZZ production cross section, the fiducial cross section in a restricted phase space for each decay channel (integrated, and as a function of kinematic parameters for the ZZ selections) and limits on anomalous nTGCs using the observed ZZ event yields as a function of the transverse momentum of the leading Z boson 4 .The results presented in this paper supersede the previously published results [13] which were derived with the first 1.02 fb −1 of the dataset used here, only with the ZZ → + − + − decay channel and with the use of the total ZZ event count for the derivation of the limits on anomalous nTGCs.
The total cross section for non-resonant ZZ production is predicted at next-to-leading order (NLO) in QCD to be 6.18 +0.25  −0.18 pb, where the quoted theoretical uncertainties result from varying the factorization and renormalization scales simultaneously by a factor of two whilst using the full CT10 parton distribution function (PDF) error set [19].The cross section is calculated in the on-shell (zero-width) approximation using MCFM [20] with CT10; it includes a 5.8% contribution from gluon fusion.When the natural width of the Z boson is used and both Z bosons are required to be within the Z mass window, the NLO cross section is predicted to be 5.89 +0.22  −0.18 pb.The cross sections given here are calculated at a renormalization and factorization scale equal to half the mass of the diboson system.The total cross section using the zero-width approximation was previously measured to be 8.5 +2.7  −2.3 (stat.)+0.4 −0.3 (syst.)± 0.3 (lumi.)pb [13].This paper is organized as follows: an overview of the ATLAS detector, data, signal and background Monte Carlo (MC) samples used for this analysis is given in section 2; section 3 describes the selection of the physics objects; section 4 describes the fiducial phase space of the measurement, the corresponding ZZ cross section definition and the acceptances of the event selection and fiducial phase space; section 5 explains how the backgrounds to the + − + − and + − ν ν final states are estimated with a combination of simulation and data-driven techniques; section 6 presents the results: cross section, differential cross sections and nTGC limits; finally, a summary of the main results is given in section 7.

The ATLAS detector and data sample
The ATLAS detector [21] is a multipurpose particle detector with a cylindrical geometry.It consists of inner tracking devices surrounded by a superconducting solenoid, electromagnetic and hadronic calorimeters and a muon spectrometer with a toroidal magnetic field.The inner detector, in combination with the 2 T field from the solenoid, provides precision tracking of charged particles in the specific lepton flavours.pseudorapidity range |η| < 2.5 5 .It consists of a silicon pixel detector, a silicon microstrip detector and a straw tube tracker that also provides transition radiation measurements for electron identification in the pseudorapidity range |η| < 2.0.The calorimeter system covers the pseudorapidity range |η| < 4.9.The electromagnetic calorimeter uses liquid argon (LAr) as the active material with lead as an absorber (|η| < 3.2).It identifies electromagnetic showers and measures their energy and position; in the region |η| < 2.5 it is finely segmented and provides electron identification in conjunction with the inner detector which covers the same η region.Hadronic showers are measured in the central rapidity range (|η| < 1.7) by scintillator tiles with iron absorber, while in the end-cap region (1.5 < |η| < 3.2) a LAr calorimeter with a copper absorber is used.In the forward region (3.2 < |η| < 4.9) a LAr calorimeter with a copper absorber for the first layer and tungsten for the last two layers is used for both electromagnetic and hadronic showers.All calorimeters are used to measure jets.The muon spectrometer surrounds the calorimeters; it consists of superconducting air-core toroid magnets, high-precision tracking chambers which provide muon identification and tracking measurement in the pseudorapidity range |η| < 2.7, and separate trigger chambers covering |η| < 2. 4.
A three-level trigger system selects events to be recorded for offline analysis.The events used in this analysis were selected with single-lepton triggers with nominal transverse momentum (p T ) thresholds of 20 or 22 GeV (depending on the instantaneous luminosity of the LHC) for electrons and 18 GeV for muons.The efficiencies of the single-lepton triggers have been determined as a function of lepton pseudorapidity and transverse momentum using large samples of Z → + − events.The trigger efficiencies for events passing the offline selection described below are all greater than 98%.
The measurements presented here uses the full data sample of proton-proton collisions at √ s = 7 TeV recorded in 2011.After data quality requirements, the total integrated luminosity used in the analysis is 4.6 fb −1 with an uncertainty of 3.9% [22].

Simulated data samples
Monte Carlo simulated samples cross-checked with data are used to calculate several quantities used in this measurement, including acceptance, efficiency and some of the background to the ZZ → + − ν ν decay channel.The NLO generator PowhegBox [23,24] with the CT10 PDF set, interfaced to Pythia [25], is used to model the signal for both channels.The LO multi-leg generator Sherpa [26] with the CTEQ6L1 PDF set [27] in comparison with PowhegBox is used to evaluate systematic uncertainties.The contribution from gg → ZZ is modelled by the gg2zz generator [28] interfaced to Herwig [29] to model parton showers and to Jimmy [30] for multiparton interactions.In each case, the simulation includes the interference terms between the Z and γ * diagrams.For both the + − + − and + − ν ν final states, MCFM is used to calculate theoretical uncertainties, and Sherpa is used for the generation of signal samples with neutral triple gauge couplings.
The LO generator Alpgen [31] with CTEQ6L1 PDFs is used to simulate Z+jets, W +jets, Zγ and W γ background events with Jimmy used for multiparton interactions and Herwig for parton showers.The NLO generator MC@NLO [32] with CT10 PDFs is used to model t t background processes as well as W W production.The single-top W t process is modelled with AcerMC [33] with the MSTW2008 PDFs [34].The LO generator Herwig with MSTW2008 PDFs is used to model W Z production.The LO generator Madgraph [35] with CTEQ6L1 PDFs is also used to model Zγ and W γ * events, where Pythia is used for hadronization and showering.
The detector response is simulated [36] with a program based on Geant4 [37].Additional inelastic pp events are included in the simulation, distributed so as to reproduce the number of collisions per bunch-crossing in the data.The detector response to interactions in the out-of-time bunches from pileup is also modelled in the simulation.The results of the simulation are corrected with scale factors determined by comparing efficiencies observed in data to those in the simulated events, and the lepton momentum scale and resolution are finely adjusted to match the observed dilepton spectra in Z → events using a sample of Z bosons.

Event reconstruction and selection
Events are required to contain a primary vertex formed from at least three associated tracks with p T > 400 MeV.

Common lepton selection
Muons are identified by matching tracks (or track segments) reconstructed in the muon spectrometer to tracks reconstructed in the inner detector [38].The momenta of these combined muons are calculated by combining the information from the two systems and correcting for the energy deposited in the calorimeters.The analyses of both decay channels use muons which have full tracks reconstructed in the muon spectrometer with p T > 20 GeV and |η| < 2.5.The ZZ ( * ) → + − + − channel recovers additional ZZ acceptance with minimal additional background using a lower threshold of p T > 7 GeV and by accepting muons with segments reconstructed in the muon spectrometer (in this latter case, the muon spectrometer is used to identify the track as a muon, but its momentum is measured using the inner detector; for the purposes of the discussion below, these muons are also referred to as combined muons).
Electrons are reconstructed from an energy cluster in the electromagnetic calorimeter matched to a track in the inner detector [38]; the transverse momentum is computed from the calorimeter energy and the direction from the track parameters measured in the inner detector.The electron track parameters are corrected for bremsstrahlung energy loss using the Gaussian-sum filter algorithm [39].Electron candidates in the ZZ ( * ) → + − + − (ZZ → + − ν ν) channel are required to have longitudinal and transverse shower profiles consistent with those expected from electromagnetic showers, by satisfying the loose (medium) identification criteria described in ref. [40] reoptimized for the 2011 data-taking conditions.They are also required to have a transverse momentum of at least 7 (20) GeV and a pseudorapidity of |η| < 2.47.
In order to reject non-prompt leptons from the decay of heavy quarks and fake electrons from misidentified jets (charged hadrons or photon conversions), all selected leptons must satisfy isolation requirements based on calorimetric and tracking information and must be consistent with originating from the primary vertex.For the calorimetric isolation the scalar sum of the transverse energies, ΣE T , of calorimeter deposits inside a cone around the lepton, corrected to remove the energy from the lepton and from additional interactions (pile-up), is formed.In the ZZ ( * ) → + − + − (ZZ → + − ν ν) channel, the ΣE T inside a cone of size ∆R = (∆φ) 2 + (∆η) 2 = 0.2 (0.3) around the lepton is required to be no more than 30% (15%) of the lepton p T .For the track isolation, the scalar sum of the transverse momenta, Σp T , of inner detector tracks inside a cone of size ∆R = 0.2 (0.3) around the lepton is required to be no more than 15% of the lepton p T .The wider cone size, in conjunction with the same or tighter requirements on the fraction of extra activity allowed in the cone, corresponds to more stringent isolation requirements applied to the ZZ → + − ν ν channel compared to the ZZ ( * ) → + − + − channel.This reflects the need to reduce the much higher reducible background (predominantly from Z+jets, t t and W W ). To ensure that the lepton originates from the primary vertex, its longitudinal impact parameter |z 0 | is required to be less than 2 mm, and its transverse impact parameter significance (the transverse impact parameter divided by its error), |d 0 /σ d0 |, is required to be less than 3.5 (6) for muons (electrons).Electrons have a worse impact parameter resolution than muons due to bremsstrahlung.Since muons can radiate photons which may then convert to electron-positron pairs, electron candidates within ∆R = 0.1 of any selected muon are not considered.If two electron candidates are within ∆R = 0.1 of each other, the one with the lower p T is removed.

Extended-lepton selection
Two additional categories of muons are considered for the ZZ ( * ) → + − + − channel: forward spectrometer muons with 2.5 < |η| < 2.7 (in a region outside the nominal coverage of the inner detector) and calorimeter-tagged muons with |η| < 0.1 (where there is a limited geometric coverage in the muon spectrometer).Forward spectrometer muons are required to have a full track that is reconstructed in the muon spectrometer; if these muons are also measured in the inner detector, their momentum is measured using the combined information; otherwise, only the muon spectrometer information is used.In either case, such muons are required to have p T > 10 GeV and the ΣE T of calorimeter deposits inside a cone of size ∆R = 0.2 around the muon is required to be no more than 15% of the muon p T , while no requirement is made on Σp T .The same impact parameter requirements as for the combined muons are imposed for the forward muons measured in the inner detector; no such requirement is imposed on those measured in the muon spectrometer only.Calorimeter-tagged muons are reconstructed from calorimeter energy deposits consistent with a muon which are matched to an inner detector track with p T > 20 GeV and are required to satisfy the same impact parameter and isolation criteria as for the combined muons.
The ZZ ( * ) → + − + − channel also uses calorimeter-only electrons with 2.5 < |η| < 3.16 and p T > 20 GeV passing the tight identification requirements [40] for this forward η region, where only the longitudinal and transverse shower profiles in the calorimeters are used for their identification.Their transverse momentum is computed from the calorimeter energy and the electron direction, where the electron direction is computed using the primary vertex position and the shower barycentre position in the calorimeter.Being identified outside the acceptance of the inner detector, no impact parameter requirements can be applied to these calorimeter-only electron candidates, and their charge is not measured.Since only one such electron is allowed in the event, and since all other leptons have their charge measured, the calorimeter-only electron is assigned the charge needed to have two pairs of sameflavour opposite-sign leptons in the event.The requirements described above constrain the additional background introduced by the inclusion of calorimeter-only electrons, and no isolation requirements are imposed on such electrons.
The use of the extended-lepton selection increases the ZZ → + − + − and ZZ * → + − + − acceptance by about 6% from the forward spectrometer muons, 4% from the calorimeter-tagged muons and 6% from the forward electrons.The expected background is kept small by requiring each event to have at most one lepton from each extended-lepton category, and each such lepton to be paired with a non-extended lepton.

Jets and missing transverse momentum
For the ZZ → + − ν ν selection, events which contain at least one well-reconstructed jet are vetoed to reduce background from top-quark production.Jets are reconstructed from topological clusters of energy in the calorimeter [41] using the anti-k t algorithm [42] with radius parameter R = 0.4.The measured jet energy is corrected for detector inhomogeneities and for the non-compensating nature of the calorimeter using p T -and η-dependent correction factors based on Monte Carlo simulations with adjustments from in-situ measurements [43,44].Jets are required to have p T > 25 GeV and |η| < 4.5.In order to minimize the impact of jets from pile-up at high luminosity, the jet vertex fraction is required to be at least 0.75; the jet vertex fraction is defined as the sum of the p T of tracks associated to the jet and originating from the primary vertex, divided by the sum of the p T of all the tracks associated to the jet.If a reconstructed jet and a lepton lie within ∆R = 0.3 of each other, the jet is not considered in the analysis.
The missing transverse momentum E miss T is the imbalance of transverse momentum in the event.A large imbalance in the transverse momentum is a signature of the ZZ → + − ν ν decay channel.The two-dimensional E miss T vector is determined from the negative vectorial sum of reconstructed electron, muon and jet momenta together with calorimeter cells not associated to any object [45].Calorimeter cells are calibrated to the jet energy scale if they are associated with a jet and to the electromagnetic energy scale otherwise.Using calorimeter timing and shower shape information, events that contain jets with p T > 20 GeV and not originating from proton-proton collisions but from e.g.calorimeter signals due to noisy cells are rejected.

ZZ
events are characterized by four high-p T , isolated electrons or muons, in three channels: e + e − e + e − , µ + µ − µ + µ − and e + e − µ + µ − .Selected events are required to have exactly four leptons and to have passed at least a single-muon or single-electron trigger.Each combination of lepton pairs is required to satisfy ∆R( 1 , 2 ) > 0.2, where 1 and 2 are used hereafter to denote a pair of distinct leptons, independent of their Z parent assignment, flavour and charge.To ensure high and well-measured trigger efficiency, at least one lepton must have p T > 20 GeV (25 GeV) for the offline muon (electron) and be matched to a muon (electron) reconstructed online by the trigger system within ∆R = 0.1 (0.15).
Same-flavour, oppositely-charged lepton pairs are combined to form Z candidates.An event must contain two such pairs.In the e + e − e + e − and µ + µ − µ + µ − channels, ambiguities are resolved by choosing the combination which results in the smaller value of the sum of |m + − − m Z | for the two pairs, where m + − is the mass of the dilepton system and m Z is the mass of the Z boson [46].Figure 2 shows the correlation between the invariant mass of the leading (higher p T ) and the sub-leading (lower p T ) lepton pair.The events cluster in the region where both masses are around m Z .At least one lepton pair is required to have invariant mass within the Z mass window, 66 < m + − < 116 GeV.If the second lepton pair satisfies this as well, the event is classified as a ZZ event; if the second pair satisfies m + − > 20 GeV, the event is classified as a ZZ * event.
With the selection described here, 84 ZZ * → + − + − candidates are observed, out of which 66 are classified as ZZ → + − + − candidates.From the 84 (66) ZZ * → + − + − (ZZ → + − + − ) candidates, 8 (7) candidates contain extended leptons.The mass of the leading lepton pair versus the mass of the sub-leading lepton pair.The events observed in the data are shown as solid circles and the ZZ ( * ) → + − + − signal prediction from simulation as boxes.The size of each box is proportional to the number of events in each bin.The region enclosed by the solid (dashed) lines indicates the signal region defined by the requirements on the lepton-pair masses for ZZ (ZZ * ) events, as defined in the text.

ZZ →
ν events are characterized by large missing transverse momentum and two high-p T , isolated electrons or muons.Selected events are required to have exactly two leptons of the same flavour with 76 < m + − < 106 GeV and to have passed at least a single-muon or a single-electron trigger.The mass window is chosen to be tighter than the mass window used for the ZZ ( * ) → + − + − channel in order to reduce the background from t t and W W .The lepton pair is required to have ∆R( + , − ) > 0.3.This requirement reflects the choice of the isolation cone for the leptons.The same trigger matching requirement as in the ZZ ( * ) → + − + − channel is used.
The ZZ → + − ν ν decay channel analysis makes use of several selections to reduce background.The largest background after the mass window requirement consists of Z+jets events, which are associated with non-zero missing transverse momentum when the E miss T is mismeasured or when a b-quark decays to leptons and neutrinos inside of a jet.Since the Z bosons tend to be produced back-to-back, the axial-E miss T (defined as the projection of the E miss and fractional p T difference are shown in figure 3.In order to reduce Z+jets background, the axial-E miss T must be greater than 75 GeV, and the fractional p T difference must be less than 0.4.To reduce background from top-quark production, events which contain at least one reconstructed jet with p T > 25 GeV and |η| < 4.5 are rejected. To reduce background from W Z production, events with a third lepton (electron or muon) with p T greater than 10 GeV are rejected.The shape of the jet multiplicity distribution is well modelled in Monte Carlo simulation as shown in figure 4 for the ZZ → + − + − and ZZ → + − ν ν selections, however, there is an overall excess of about 20% in the ZZ → + − + − selection.With this selection, 87 ZZ → + − ν ν candidates are observed in data.T after all selection requirements, except for the fractional pT difference (the last bin also contains events with fractional pT difference greater than 1).In all plots, the points are data and the stacked histograms show the signal prediction from simulation.The shaded band shows the combined statistical and systematic uncertainties.

Signal acceptance
The Z boson decays to hadrons, neutrinos and charged leptons with branching fractions of 69.9%, 20.0% and 10.1%, respectively [46].The two ZZ decay channels considered in this paper, ZZ → + − + − and ZZ → + − ν ν, have branching fractions of 0.45% and 2.69%, respectively 6 , where decays involving τ leptons are not included in these branching fractions.Some of the ZZ decays produce one or more charged leptons which pass through the uninstrumented regions of the detector, and as such cannot be reconstructed.In order to measure the total ZZ cross section, the measured decays are extrapolated to non-measured parts of the phase-space; this results in the measurement being more dependent on theory predictions.Consequently, two types of cross sections are measured: fiducial and total.The fiducial cross section is the cross section measured within a restricted phase space, and the total cross section is the cross section extrapolated to the total phase space.The total cross section calculation depends on the choice of Z mass range.The cross section is calculated using the Z boson natural width rather than the zero-width approximation, and includes the mass window requirement (66 to 116 GeV) to remove most of the γ * contamination.The ratio of the total cross section calculated with both Z bosons within the mass window to the total cross section calculated using the zero-width approximation is 0.953, as the mass window requirement removes some of the Z bosons in the tails of the mass distribution.

Fiducial region definitions
The fiducial cross section is restricted to a region which is constructed to closely match the instrumented region and the event selection; for simplicity, only the most inclusive requirements on the lepton η and p T are used for the definition of the fiducial phase space.The fiducial cross section σ fid ZZ is calculated as: which depends on a correction factor given by the number of simulated ZZ ( * ) events which satisfy the full event selection divided by the number of ZZ ( * ) events generated in the fiducial region, C ZZ ; the integrated luminosity, L; the number of selected events, N obs ; and the amount of estimated background, N bkg .For the calculation of C ZZ , final states including pairs of oppositely-charged leptons produced from decays of Z → τ + τ − → + − ν νν ν are included in the number of selected events (numerator) since those decays have an identical final state to the signal and are not subtracted as background but are excluded from the fiducial region (denominator) because the fiducial regions are defined only with ZZ ( * ) decays directly to electrons, muons or neutrinos, depending on the channel.The contribution from such τ decays is estimated from Monte Carlo simulation to be < 0.1 % for the ZZ → + − ν ν selection, 0.24±0.01%for the ZZ → + − + − selection and 1.73±0.04%for the ZZ * → + − + − selection.Fiducial requirements are applied at generator level.To reduce the dependence on QED radiation, dressed leptons are always used, for which the lepton four-momentum is summed with the four-momentum of all photons within ∆R = 0.1.
The ZZ → + − + − fiducial region is defined using the following requirements: (i) two pairs of same-flavour opposite-sign electrons or muons, with each lepton satisfying p T > 7 GeV, |η | < 3.16 and at least a distance ∆R = 0.2 from any other selected lepton, i.e., ∆R( 1 , 2 ) > 0.2, and (ii) both dilepton invariant masses within the Z mass window.A ZZ * → + − + − fiducial region is defined with the same criteria as in the ZZ → + − + − case, except that one dilepton invariant mass requirement is relaxed to be greater than 20 GeV.
The ZZ → + − ν ν fiducial region is defined by requiring: (i) two same-flavour opposite-sign electrons or muons, each with p T > 20 GeV, |η | < 2.5, with ∆R( + , − ) > 0.3, (ii) dilepton invariant mass close to the Z boson mass: 76 < m + − < 106 GeV, (iii) dineutrino invariant mass close to the Z boson mass: 66 < m ν ν < 116 GeV, (iv) no jet with p j T > 25 GeV and |η j | < 4.5, and (v) Jets are defined at generator level using the same jet algorithm as used in reconstructed events and including all final state particles after parton showering and hadronization.

Extrapolation to the total phase space
The total ZZ cross section is measured using the ZZ → + − + − and ZZ → + − ν ν selections.The total cross section is calculated using the fiducial acceptance, A ZZ (the fraction of ZZ events with Z bosons in the Z mass window that fall into the fiducial region) and the branching fraction, BF: The fiducial acceptances A ZZ are estimated from Monte Carlo simulation, using PowhegBox for the ZZ → + − ν ν channel and PowhegBox and gg2zz for the ZZ → + − + − channel.The fiducial acceptance of the ZZ → + − ν ν channel is much more constrained than the ZZ → + − + − channel in order to reduce background.Values are given in table 2.

Selection
A ZZ ZZ → + − + − 0.804 ± 0.001 ± 0.010 ZZ → + − ν ν 0.081 ± 0.001 ± 0.004 Table 2. Acceptance AZZ for the two decay channels used for the measurement of the total ZZ production cross section.The first uncertainty is statistical while the second is systematic.

Systematic uncertainties
Table 3 summarizes the systematic uncertainties on C ZZ and A ZZ .For C ZZ in the ZZ ( * ) → + − + − selections, the dominant systematic uncertainties arise from the lepton reconstruction efficiency, the efficiency of the isolation and impact parameter requirements, and the differences in C ZZ estimated by Sherpa and PowhegBox; uncertainties on the trigger efficiency and the lepton energy scale and resolution are small.In the ZZ → + − ν ν channel the dominant C ZZ uncertainties are from uncertainties on the lepton reconstruction efficiency, the lepton energy scale and resolution, and the missing transverse momentum modelling and jet veto uncertainty; uncertainties on the trigger efficiency and due to differences in C ZZ estimated by Sherpa and PowhegBox also contribute.
The uncertainties on C ZZ from the reconstruction efficiency, energy scale and resolution, isolation and impact parameter requirements and trigger efficiency are estimated by varying the data-driven correction factors applied to simulation by their systematic and statistical uncertainties.The systematic uncertainties on events with extended leptons used in the ZZ ( * ) → + − + − channel are slightly higher than in events without them; nevertheless, since their relative contribution is small, the effect on the uncertainty of the combined channels is negligible.The generator systematic uncertainty for C ZZ accounts for the effect of choosing a different renormalization and factorization scale and PDF set.
For A ZZ , the systematic uncertainties are due to theoretical uncertainties which come from the PDFs, the choice of the renormalization and factorization scales, the modelling of the contribution from gg initial states and the parton shower model, as given in table 3.For the ZZ → + − ν ν channel, uncertainties in the efficiency of the jet veto are also taken into account through the calculation of a scale factor; the ratio of the jet veto efficiency in data to that in MC simulation is taken from a sample of single Z events and then applied to ZZ events [47].The systematic uncertainties due to the PDFs and scales are evaluated with MCFM by taking the difference between the A ZZ obtained using the CT10 and MSTW2008 PDF sets, as well as using the 44 CT10 error sets, and by shifting the factorization and renormalization scales up and down by a factor of two.An additional uncertainty is assigned to account for the effect of different modelling at the generator level.Since the ZZ * → + − + − measurement is not used for the total cross section, its A ZZ acceptance is irrelevant and only uncertainty values related to C ZZ are given.
The uncertainty on the integrated luminosity is 3.9% [22].The uncertainty on the background estimates is discussed in the following sections.5 Background estimation Background to the ZZ ( * ) → + − + − signal originates from events with a Z (or W ) boson decaying to leptons accompanied by additional jets or photons (W/Z + X), from top-quark production and from other diboson final states.Such events may contain electrons or muons from the decay of heavyflavoured hadrons, muons from in-flight decay of pions and kaons, or jets and photons misidentified as electrons.The majority of these background leptons are rejected by the isolation requirements.
The background estimate follows a data-driven method in which a sample of events containing three leptons satisfying all selection criteria plus one 'lepton-like jet' is identified; such events are denoted as j.For muons, the lepton-like jets are muon candidates that fail the isolation requirement or fail the impact parameter requirement but not both.For electrons with |η| < 2.47, the lepton-like jets are clusters in the electromagnetic calorimeter matched to inner detector tracks that fail either the full electron selection or the isolation requirement but not both.For electrons with |η| > 2.5, the lepton-like jets are electromagnetic clusters that are reconstructed as electrons but fail the tight identification requirements.The events are otherwise required to satisfy the full event selection, treating the lepton-like jet as if it were a fully identified lepton.The background is then estimated by weighting the j events by a measured factor f , which is the ratio of the probability for a non-lepton to satisfy the full lepton selection criteria to the probability of satisfying the lepton-like jet criteria.The background in which two selected leptons originate from jets is treated similarly, by identifying a data sample with two leptons and two lepton-like jets; such events are denoted as jj.The total number of expected background + − + − events, N (BG), is calculated as: where double counting from j and jj events is accounted for, and the term N (ZZ) is a Monte Carlo estimate correcting for contributions from signal ZZ ( * ) → + − + − events having a real lepton that is classified as a lepton-like jet (the equivalent correction to the term N ( jj) is negligible).
The factor f is measured in a sample of data selected with single-lepton triggers which contain a Z boson candidate: a pair of isolated same-flavour opposite-sign electrons or muons.In these selected events, f is measured, using the lepton and lepton-like jet candidates not assigned to the Z boson, as the ratio of the number of selected leptons to the number of lepton-like jets, after correcting for expected true lepton contributions from W Z and ZZ events using simulation.Independent values as a function of the η and p T of the lepton-like jet are measured, which are then combined assuming they are uncorrelated.The factor f is found to vary from 0.33 ± 0.01 (0.26 ± 0.02) below p T = 10 GeV to 0.09 ± 0.02 (0.46 ± 0.20) above p T = 50 GeV for electrons (muons).The quoted uncertainties are statistical.Then, with the same procedure, a value for f is also derived using the simulated samples of background processes.The difference between the value of f derived in data and in simulation is assigned as a systematic uncertainty on f .The statistical and systematic uncertainties are then added in quadrature to derive a combined uncertainty on f , which varies as a function of p T from 14% (19%) below 10 GeV to 22% (51%) above 50 GeV for electrons (muons).For the muons, the total uncertainty on f is dominated by its statistical uncertainty.The background estimates for the ZZ → + − + − and ZZ * → + − + − selections are 0.9 +1.1 −0.9 (stat.)± 0.7(syst.)and 9.1 ± 2.3(stat.)± 1.3(syst.)events, respectively, as shown in tables 4 and 5.The statistical uncertainty on the background estimate comes from the statistical uncertainty on the numbers of j, jj and ZZ ( * ) → + − + − events used in eq.5.1.The systematic uncertainty results from the combined uncertainty on f .In cases where the overall estimate is negative, the background estimate is described using a truncated Gaussian with mean at zero and standard deviation equal to the estimated statistical and systematic uncertainties added in quadrature.
The extra background induced by the use of the extended leptons in the ZZ ( * ) → + − + − channel is estimated to be negligible in the ZZ → + − + − selection, and about 20% (2 events out of the 9.1 estimated, compared to a signal gain of about 10.6 events out of the 64.4 expected) in the The background is also estimated purely from the simulated samples of background processes, and is predicted to be 1.5 ± 0.4 events for the ZZ → + − + − selection and 8.3 ± 1.3 events for the ZZ * → + − + − selection, with uncertainties being statistical only.These estimates compare well with the data-driven results given in tables 4 and 5.According to the estimate from simulation, the dominant source of background is Z+jets events, with only about a 10% to 20% contribution from other diboson channels (W Z and W W ), and a negligible contribution from events with top quarks.
Differential background distributions are determined by first deriving the shape of the distributions from the background MC samples.This is achieved by selecting events where one Z candidate is ±0.9 (syst.)±1.3(syst.) Table 5. Expected number of background events for the ZZ * → + − + − selection in 4.6 fb −1 of data, for the individual decay modes (columns 2, 3 and 4) and for their combination (last column).If the central value of the estimate is negative, the upper bound on the number of events in that channel is derived as detailed in section 5.1.
required to satisfy the nominal lepton selection, while the other Z candidate is formed by leptons satisfying relaxed criteria for the isolation requirements and transverse impact parameter significance.The shape determined in this way is then scaled such that the total number of events in the distribution is equal to the data-driven background estimate shown in tables 4 and 5.

ZZ → + − ν ν background
There are several sources of background to the ZZ → + − ν ν channel.Processes such as t t, W W , W t or Z → τ + τ − production give two true isolated leptons with missing transverse momentum.Diboson W Z events in which both bosons decay leptonically have three charged leptons, but if one lepton from a W or Z boson decay is not identified, the event has the same signature as the signal.Production of a Z boson in association with jets gives two isolated leptons from the Z boson decay and may have missing transverse momentum if the jet momenta are mismeasured.Finally, production of a W boson in association with jets or photons may satisfy the selection requirements when one of the jets or photons is misidentified as an isolated lepton.All of the backgrounds are measured with data-driven techniques except for W Z and W γ. The total background is estimated to be 46.9 ± 4.8 ± 1.9 events as summarized in table 6.

Backgrounds from t t, W t, W W and Z
The contributions from t t, W t, W W and Z → τ + τ − processes are measured by extrapolating from a control sample of events with one electron and one muon (instead of two electrons or two muons), which otherwise satisfy the full ZZ → + − ν ν selection.This sample is free from signal events.The extrapolation from the eµ channel to the ee or µµ channel uses the relative branching fractions (2 : 1 : 1 for eµ : ee : µµ) as well as the ratio of the efficiencies ee or µµ of the ee or µµ selections to the efficiency eµ of the eµ selection, which differs from unity due to differences in the electron and muon efficiencies.
For the electron channel, this is represented by the equation: where N data eµ is the number of observed eµ events and N sim eµ is the number expected events from processes other than t t, W t, W W and Z → τ + τ − (W Z, ZZ, W +jet, Z+jet and W/γ).Therefore, (N data eµ − N sim eµ ) is the estimate of t t, W t, W W and Z → τ + τ − production in the control sample.The efficiency correction factor, ee / eµ , corrects for the difference between electron and muon efficiency.The efficiency correction factor is measured in data using reconstructed Z → + − events, as where N data ee and N data µµ are the number of observed ee or µµ events in the Z boson mass window, respectively, after all lepton selection requirements and the Z boson mass window requirement are applied.A parallel argument gives N bkg µµ .This procedure is repeated in bins of p Z T in order to obtain the p T distribution of the t t, W t, W W and Z → τ + τ − background.
The dominant uncertainty is statistical (25%), due to the limited number of events in the control samples.Additional uncertainties are due to systematic uncertainties in the normalization of the simulated samples used to correct the eµ contribution (5.5%) and the systematic uncertainty in the efficiency correction factor (4.5%).

Background from W Z production with leptonic decays
Events from leptonic W Z decays may result in an + − E miss T signature when one lepton from the W or Z boson is not reconstructed.The contribution from this process is estimated using the simulated samples described in section 2.1.The estimate is checked using a control region with three highp T isolated leptons.The two dominant processes that contribute to this control region are W Z and Z+jets production, where the W Z boson pair decays to three leptons and a neutrino and the Z+jets contribution has two real leptons from the Z decay and a misidentified lepton from the jet.The technique used to estimate the background in the ZZ ( * ) → + − + − channel is also used to normalize the contribution from Z+jets in the three-lepton control region.The W Z Monte Carlo expectation is consistent with the data.The systematic uncertainties are estimated in the same way as for signal Monte Carlo events.

Background from Z bosons with associated jets
Occasionally events with a Z boson produced in association with jets may have large amounts of missing transverse momentum due to mismeasurement of the momenta of the jets.This background is estimated using events with a high-p T photon and jets as a template, since the mechanism for large missing transverse momentum is the same as in Z+jets events.The events are reweighted such that the photon E T matches the observed Z boson p T and are normalized to the observed Z + jets yield.The procedure is repeated in bins of p Z T in order to obtain the p T distribution of the Z+jets backgrounds.The largest systematic uncertainty is due to the subtraction of W γ, Zγ, t t and W → eν contributions to the γ+jets sample, which is 33% in the ee channel and 37% in the µµ channel.

Background from events with a misidentified lepton
A small contribution to the selected sample is due to events in which one of the two leptons comes from the decay of a W or Z boson (called 'real' below) and the second is a 'fake', corresponding both non-prompt leptons and misidentified π 0 mesons or conversions.
The dominant fake-muon mechanism is the decay of heavy-flavoured hadrons, in which a muon survives the isolation requirements.In the case of electrons, the three mechanisms are heavy-flavour hadron decay, light-flavour jets with a leading π 0 overlapping with a charged particle, and conversion of photons.Processes that contribute are top-quark pair production, production of W bosons in association with jets and multi-jet production.
The 'matrix method' [48] is applied to estimate the fraction of events in the signal regions that contain at least one fake lepton.The method measures the number of fake leptons in backgrounddominated control regions and extrapolates to the ZZ selection region using factors measured in data.The shape of the background is provided by taking the background as uniformly distributed among the bins and treating each bin as statistically uncorrelated.The dominant systematic uncertainty is due to the uncertainty on the extrapolation factors and the limited numbers of events in the control samples, giving a total uncertainty of 63% and 44% in the ee and µµ channels, respectively.

Results
Three types of measurements are presented: • integrated fiducial and total ZZ cross sections; • differential cross sections normalized to the overall measured cross sections for the p Z T and ∆φ( + , − ) of the leading Z boson, and the mass (transverse mass 7 ) of the ZZ system for the ZZ → + − + − (ZZ → + − ν ν) selection; and • limits on the anomalous nTGCs.

Cross section measurements
The expected and observed event yields after applying all selection criteria are shown in table 7 for both channels.Figure 4 shows the jet multiplicity in selected ZZ → + − + − and ZZ → + − ν ν events before the jet veto is applied.Figures 5 and 6 show the transverse momentum and mass of the ZZ system in selected ZZ → + − + − and ZZ * → + − + − events respectively.Figure 7 shows the transverse momentum and mass of the two-charged-lepton system in selected ZZ → + − ν ν events.The shapes of the distributions are consistent with the predictions from the simulation.
The ZZ ( * ) → + − + − and ZZ → + − ν ν fiducial cross sections are determined using a maximum likelihood fitting method, taking into account the integrated luminosity and the C ZZ correction factors discussed in section 4. A Poisson probability function is used to model the number of expected events, multiplied by Gaussian distribution functions which model the nuisance parameters representing systematic uncertainties.The measured fiducial cross sections are: −2.9 (stat.)+1.7 −1.7 (syst.)± 0.5 (lumi.)fb.46.9 ± 4.8 ± 1.9 Table 7. Summary of observed ZZ → + − + − , ZZ * → + − + − and ZZ → + − ν ν candidates in the data, total background estimates and expected signal for the individual decay modes (columns 2 to 4) and for their combination (last column).The quoted uncertainties and limits represent 68% confidence intervals; the first uncertainty is statistical while the second is systematic.The uncertainty on the integrated luminosity (3.9%) is not included.The measured cross sections are compatible with these theoretical values.
The total ZZ cross section is calculated by extrapolating to the full phase space while each Z boson is required to have a mass within the Z mass window.Both ZZ → + − + − and ZZ → + − ν ν events are combined in the maximum likelihood fit, taking into account the known Z branching fractions [46] and the A ZZ kinematic and geometrical acceptances (section 4).Correlated systematic uncertainties between the ZZ → + − + − and ZZ → + − ν ν channels are taken into account in the fit using a single Gaussian for the nuisance parameter for each source of correlated uncertainty.The measured value of the total ZZ cross section is: σ tot ZZ = 6.7 ± 0.7 (stat.)+0.4 −0.3 (syst.)± 0.3 (lumi.)pb.
The result is consistent within errors with the NLO Standard Model total cross section for this process of 5.89 +0.22 −0.18 pb, where the quoted theoretical uncertainties result from varying the factorization and renormalization scales simultaneously by a factor of two and from using the full CT10 PDF error set.

Differential cross sections
The differential cross sections present a more detailed comparison of theory to measurement, allowing a generic comparison of the kinematic distributions to new theories.Variables which are sensitive to new phenomena, such as p Z T , m ZZ and ∆φ( + , − ), are used with bin boundaries chosen to maximize sensitivity to nTGCs.At the same time, the bin widths were chosen to be commensurate with the resolution.
The measured distributions are unfolded back to the underlying distributions, accounting for the effect of detector resolution, efficiency and acceptance, within the fiducial region of each measurement.The unfolding procedure is based on a Bayesian iterative algorithm [49].The algorithm takes as input a prior for the kinematic distribution and iterates using the posterior distribution as prior for the next iteration.The initial prior is taken from the signal Monte Carlo expectation calculated using the PowhegBox generator and three iterations are performed.The uncertainty on the unfolded distributions is dominated by the statistical uncertainty, which is about 30% in most bins.The systematic uncertainty is no more than 5% in any bin.The dependence of the unfolded cross sections on the choice of the initial prior is tested by unfolding the measured distributions using a different generator (Sherpa).The difference between the two is taken as a systematic uncertainty to account for differences in generator modelling (e.g.QCD radiation).The difference in unfolded distributions between three iterations and four iterations is much lower than the statistical uncertainty and it is taken as a further uncertainty on the unfolding procedure.Systematic uncertainties related to detector effects (e.g.lepton reconstruction efficiency) are evaluated using pseudo-experiments.
Figures 8 to 10 show the differential cross sections normalized to the fiducial cross sections for the p Z T and ∆φ( + , − ) of the leading Z boson, and for the mass (transverse mass) of the ZZ system for the ZZ → + − + − (ZZ → + − ν ν) selection.The Standard Model prediction is consistent with the measurement in each case.

Anomalous neutral triple gauge couplings
Anomalous nTGCs for on-shell ZZ production can be parameterized by two CP-violating (f V 4 ) and two CP-conserving (f V 5 ) complex parameters (where V = Z, γ) which are zero in the Standard Model [3].A form-factor parameterization is introduced leading to couplings which vanish at high parton centreof-mass energy √ ŝ: ) n , ensuring partial-wave unitarity.Here, Λ is the energy scale at which physics beyond the Standard Model would be directly observable, f V i0 are the low-energy approximations of the couplings, and n is the form-factor power.Values of n = 3 and Λ = 3 TeV are chosen, so that expected limits are within the values allowed by requiring that unitarity is not violated at LHC energies [3].The results with an energy cutoff Λ = ∞ (i.e.without a form factor) are also presented as a comparison in the unitarity violating scheme.
Limits on anomalous nTGCs are determined using the observed and expected numbers of ZZ → + − + − and ZZ → + − ν ν events binned9 in p Z T , as seen in table 8. Figure 11 shows the observed p Z T distributions, together with the SM expectation and the predicted distributions for nTGC values close to the previous limits obtained by ATLAS [13].Using an increased data sample compared Expected background Expected ZZ signal Observed events ZZ → + − + − 0 < p Z T < 60 GeV 0.6 ± 0.8 ± 0.5 27.9 ± 0. Table 8.Total background, expected signal and observed events as a function of the pT of the leading Z for the ZZ → + − + − and ZZ → + − ν ν selections.For the expected signal and background events, the first uncertainty is statistical and the second is systematic.
with our previous measurement, including the ZZ → + − ν ν channel, and exploiting the differential event yields, the precision is expected to improve by about a factor of five.The dependency of the couplings on the expected number of events in each p Z T bin is parameterized using fully simulated events, generated with Sherpa [26], subsequently reweighted using the Baur-Rainwater [3,50] and BHO [51] MC generators.The next-to-leading-order matrix elements with their nTGC dependence have been extracted from the BHO MC generator for 2 → 5 events and the Baur-Rainwater MC generator for 2 → 4 events and introduced into a framework [52] that enables a calculation of the amplitude given the four vectors and the identity of the incoming and outgoing particles from the hard process.
Confidence intervals for the anomalous triple gauge couplings are determined using the maximum profile likelihood ratio.Limits are set on each coupling, assuming all of the other couplings are zero (as in the Standard Model), and on pairs of couplings assuming the remaining two couplings are zero.Table 9. One-dimensional 95% confidence intervals for anomalous neutral gauge boson couplings, where the limit for each coupling assumes the other couplings are fixed at zero, their SM value.Limits are presented for form factor scales of Λ = 3 TeV and Λ = ∞ and include both statistical and systematic uncertainties; the statistical uncertainties are dominant.
The profile likelihood ratio is calculated for the data, and also for 10000 pseudo-experiments generated using the expected number of events at each point in the one-or two-dimensional nTGC parameter space.A point is rejected if more than 95% of the pseudo-experiments have a larger profile likelihood ratio value than the one observed in data.The systematic errors are included as nuisance parameters.
The resulting limits for each coupling are listed in table 9. Two-dimensional 95% confidence intervals10 are shown in figure 12.The one-dimensional limits are more stringent than those derived from measurements at LEP [8] and the Tevatron [9] and previously by ATLAS [13]; it should be noted that the limits from LEP do not use a form factor, and those from the Tevatron use Λ = 1.2 TeV.A comparison of the LHC limits with those derived from LEP and Tevatron is shown in figure 13.

Conclusions
A measurement of the ZZ ( * ) production cross section in LHC proton-proton collisions at √ s = 7 TeV is presented with data collected by the ATLAS detector, using the ZZ ( * ) → + − + − and ZZ → + − ν ν decay channels.Fiducial cross sections are measured for three production and decay selections, and the results are compatible with the SM expected cross sections.Using the ZZ → + − + − and ZZ → + − ν ν selections, the total ZZ production cross section is determined to be: σ tot ZZ = 6.7 ± 0.7 (stat.)+0.4 −0.3 (syst.)± 0.3 (lumi.)pb.
The result is statistically consistent with the NLO Standard Model prediction of 5.89 +0.22 −0.18 pb, calculated with Z bosons with a mass between 66 and 116 GeV, and supersedes the previous measurements made with part of the same dataset [13].Unfolded distributions of the fiducial cross sections are derived for the p Z T and ∆φ( + , − ) of the leading Z boson and for m ZZ in the ZZ → + − + − selection and the m T in the ZZ → + − ν ν selection.
The event yields as a function of the p T of the leading Z boson for the ZZ → + − + − and ZZ → + − ν ν selections are used to derive 95% confidence intervals for anomalous neutral triple gauge boson couplings.These limits are more stringent than those derived from measurements at LEP [8] and the Tevatron [9].They improve the previous published results from ATLAS [13] by approximately a factor of five and supersede them.Figure 12.Two-dimensional triple gauge coupling limits for form factor scale Λ = ∞.The one-dimensional triple gauge coupling limits are shown as vertical and horizontal lines inside the two-dimensional ellipses, whose shape is determined by the theoretical correlations.For each two-dimensional limit the other TGC parameters are assumed to be zero.Since most of the sensitivity of the measurement is contained in a single bin, the likelihood ratio used to obtain the two-dimensional limits has one effective degree of freedom.[8] The LEP collaborations ALEPH, DELPHI, L3, OPAL, and the LEP Electroweak Working Group, A Combination of preliminary electroweak measurements and constraints on the standard model, hep-ex/0612034.
Figure2.The mass of the leading lepton pair versus the mass of the sub-leading lepton pair.The events observed in the data are shown as solid circles and the ZZ ( * ) → + − + − signal prediction from simulation as boxes.The size of each box is proportional to the number of events in each bin.The region enclosed by the solid (dashed) lines indicates the signal region defined by the requirements on the lepton-pair masses for ZZ (ZZ * ) events, as defined in the text.

T
along the direction opposite to the Z → + − candidate in the transverse plane) is a powerful variable to distinguish ZZ → + − ν ν decays from Z+jets.The axial-E miss T is given by − E miss T • p Z /p Z T , where p Z T is the magnitude of the transverse momentum of the Z candidate.Similarly, the fractional p T difference, |E miss T − p Z T |/p Z T is a good variable to distinguish the two.The axial-E miss T

Figure 3 .
Figure 3.For + − ν ν candidates in all channels figure (a) shows the axial-E miss T after all selection requirements, except for the axial-E miss T , and figure (b) shows the fractional pT difference between E miss T and p ZT after all selection requirements, except for the fractional pT difference (the last bin also contains events with fractional pT difference greater than 1).In all plots, the points are data and the stacked histograms show the signal prediction from simulation.The shaded band shows the combined statistical and systematic uncertainties.

Figure 4 .
Figure 4. (a) Jet multiplicity for the ZZ → + − + − selection and (b) jet multiplicity for the ZZ → + − ν ν selection (with all selections applied but the jet veto).The points represent the observed data.In (a) the ZZ → + − + − background is normalized to the data-driven (dd) estimate, while in (b) the histograms show the prediction from simulation.The shaded band shows the combined statistical and systematic uncertainty on the prediction.

Figure 5 .Figure 6 .Figure 7 .
Figure 5. (a) Transverse momentum p ZZ T and (b) invariant mass m ZZ of the four-lepton system for the ZZ selection.The points represent the observed data and the histograms show the prediction from simulation, where the background is normalized to the data-driven (dd) estimate as described in section 5.1.The shaded band shows the combined statistical and systematic uncertainty on the prediction.

Figure 8 .Figure 9 .
Figure 8. Unfolded ZZ fiducial cross sections in bins of the pT of the leading Z boson for (a) the ZZ → + − + − selection, where a discontinuity is indicated by the parallel pairs of lines, and (b) the ZZ → + − ν ν selection.

Figure 10 .
Figure 10.Unfolded ZZ fiducial cross sections in bins of (a) m ZZ for the ZZ → + − + − selection and (b) m ZZ T for the ZZ → + − ν ν selection.

Figure 11 . 50 3
Figure 11.The leading Z boson transverse momentum distributions for (a) the ZZ → + − + − selection and (b) the ZZ → + − ν ν selection.The observed distributions are shown as filled circles, the SM expected signal and background are shown as filled histograms, and the predicted distributions for four different nTGC samples with form factor scales of Λ = 3 TeV and nTGC coupling values set near the edge of the exclusion set in the 1 fb −1 analysis [13] are shown as dashed lines.

Table 1 .
Correction factors CZZ for each production and decay channel.The first uncertainty is statistical while the second is systematic.

Table 3 .
Summary of systematic uncertainties, as relative percentages of the correction factor CZZ or the acceptance of the fiducial region AZZ .Dashes indicate uncertainties which are not relevant.