Measurement of distributions sensitive to the underlying event in inclusive Z-boson production in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$pp$$\end{document}pp collisions at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=7$$\end{document}s=7 TeV with the ATLAS detector

A measurement of charged-particle distributions sensitive to the properties of the underlying event is presented for an inclusive sample of events containing a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z$$\end{document}Z-boson, decaying to an electron or muon pair. The measurement is based on data collected using the ATLAS detector at the LHC in proton–proton collisions at a centre-of-mass energy of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$7$$\end{document}7 TeV with an integrated luminosity of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4.6$$\end{document}4.6 fb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-1}$$\end{document}-1. Distributions of the charged particle multiplicity and of the charged particle transverse momentum are measured in regions of azimuthal angle defined with respect to the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z$$\end{document}Z-boson direction. The measured distributions are compared to similar distributions measured in jet events, and to the predictions of various Monte Carlo generators implementing different underlying event models.


Introduction
In order to perform precise Standard Model measurements or to search for new physics phenomena at hadron colliders, it is important to have a good understanding of not only the shortdistance hard scattering process, but also of the accompanying activity -collectively termed the underlying event (UE). This includes partons not participating in the hard-scattering process (beam remnants), and additional hard scatters in the same proton-proton collision, termed multiple parton interactions (MPI). Initial and final state gluon radiation (ISR, FSR) also contribute to the UE activity. It is impossible to unambiguously separate the UE from the hard scattering process on an event-by-event basis. However, distributions can be measured that are sensitive to the properties of the UE.
The soft interactions contributing to the UE cannot be calculated reliably using perturbative quantum chromodynamics (pQCD) methods, and are generally described using different phenomenological models, usually implemented in Monte Carlo (MC) event generators. These models contain e-mail: atlas.publications@cern.ch many parameters whose values and energy dependences are not known a priori. Therefore, the model parameters must be tuned to experimental data to obtain insight into the nature of soft QCD processes and to optimise the description of UE contributions for studies of hard-process physics.
Measurements of distributions sensitive to the properties of the UE have been performed in proton-proton ( pp) collisions at √ s = 900 GeV and 7 TeV in ATLAS [1-5], ALICE [6] and CMS [7,8]. They have also been performed in pp collisions in events with jets and in Drell-Yan events at CDF [9,10] at centre-of-mass energies of √ s = 1.8 TeV and 1.96 TeV .
This paper reports a measurement of distributions sensitive to the UE, performed with the ATLAS detector [11] at the LHC in pp collisions at a centre-of-mass energy of 7 TeV. The full dataset acquired during 2011 is used, corresponding to an integrated luminosity of 4.64 ± 0.08 fb −1 . Events with a Z -boson candidate decaying into an electron or muon pair were selected, and observables constructed from the final state charged particles (after excluding the lepton pair) were studied as a function of the transverse momentum 1 of the Z -boson candidate, p Z T . This paper is organised as follows: the definitions of the underlying event observables are given in Sect. 2. The ATLAS detector is described in Sect. 3. In Sect. 4, the MC models used in this analysis are discussed. Sections 5 and 6 describe the event selection, and the correction for the effect of multiple proton-proton interactions in the same bunch crossing (termed pile-up). The correction of the data to the 1 The ATLAS reference system is a Cartesian right-handed coordinate system, with the nominal collision point at the origin. The anticlockwise beam direction defines the positive z-axis, while the positive x-axis is defined as pointing from the collision point to the center of the LHC ring and the positive y-axis points upwards. The azimuthal angle φ is measured around the beam axis, and the polar angle θ is measured with respect to the z-axis. The pseudorapidity is given by η = − ln tan(θ/2). Transverse momentum is defined relative to the beam axis. particle level, and the combination of the electron and muon channel results are described in Sect. 7. Section 8 contains the estimation of the systematic uncertainties. The results are discussed in Sect. 9 and finally the conclusions are presented in Sect. 10.

Underlying event observables
Since there is no final-state gluon radiation associated with a Z -boson, lepton-pair production consistent with Z -boson decays provides a cleaner final-state environment than jet production for measuring the characteristics of the underlying event in certain regions of phase space. The direction of the Z -boson candidate is used to define regions in the azimuthal plane that have different sensitivity to the UE, a concept first used in [12]. As illustrated in Fig. 1, the azimuthal angular difference between charged tracks and the Z -boson, | φ| = |φ − φ Z -boson |, is used to define the following three azimuthal UE regions: -| φ| < 60 • , the toward region, -60 • < | φ| < 120 • , the transverse region, and -| φ| > 120 • , the away region.
These regions are well defined only when the measured p Z T is large enough that, taking into account detector reso-  These are defined for each azimuthal region under consideration except for p Z T lution, it can be used to define a direction. The away region is dominated by particles balancing the momentum of the Z -boson except at low values of p Z T . The transverse region is sensitive to the underlying event, since it is by construction perpendicular to the direction of the Z -boson and hence it is expected to have a lower level of activity from the hard scattering process compared to the away region. The two opposite transverse regions may be distinguished on an event-by-event basis through their amount of activity, as measured by the sum of the charged-particle transverse momenta in each of them. The more or less-active transverse regions are then referred to as trans-max and trans-min, respectively, with the difference between them on an event-by-event basis for a given observable defined as trans-diff [13,14]. The activity in the toward region, which is similarly unaffected by additional activity from the hard scatter, is measured in this analysis, in contrast to the underlying event analysis in dijet events [5].
The observables measured in this analysis are derived from the number, N ch , and transverse momenta, p T , of stable charged particles in each event. They have been studied both as one-dimensional distributions, inclusive in the properties of the hard process, and as profile histograms which present the dependence of the mean value of each observable (and its uncertainty) on p Z T . The observables are summarised in Table  1. The mean charged-particle transverse momentum is constructed on an event-by-event basis and is then averaged over all events to calculate the observable mean p T .

The ATLAS detector
The ATLAS detector [11] covers almost the full solid angle around the collision point. The components that are relevant for this analysis are the tracking detectors, the liquid-argon (LAr) electromagnetic sampling calorimeters and the muon spectrometer.
The inner tracking detector (ID) has full coverage in azimuthal angle φ and covers the pseudorapidity range |η| < 2.5. It consists of a silicon pixel detector (pixel), a semiconductor tracker (SCT) and a straw-tube transition radiation tracker (TRT). These detectors are located at a radial distance from the beam line of 50.5-150, 299-560 a nd 563-1,066 mm , respectively, and are contained within a 2 T axial magnetic field. The inner detector barrel (end-cap) consists of 3 (2 × 3) pixel layers, 4 (2 × 9) layers of double-sided silicon strip modules, and 73 (2 × 160) layers of TRT strawtubes. These detectors have position resolutions typically of 10, 17 a nd 130 µm for the r -φ coordinates (only for TRT barrel), respectively. The pixel and SCT detectors provide measurements of the r -z coordinates with typical resolutions of 115 a nd 580 µm, respectively. The TRT acceptance is |η| < 2.0. A track traversing the barrel typically has 11 silicon hits (3 pixel clusters and 8 strip clusters) and more than 30 straw-tube hits.
A high-granularity lead, liquid-argon electromagnetic sampling calorimeter [15] covers the pseudorapidity range |η| < 3.2. Hadronic calorimetry in the range |η| < 1.7 is provided by an iron scintillator-tile calorimeter, consisting of a central barrel and two smaller extended barrel cylinders, one on either side of the central barrel. In the end-caps (|η| > 1.5), the acceptance of the LAr hadronic calorimeters matches the outer |η| limits of the end-cap electromagnetic calorimeters. The LAr forward calorimeters provide both electromagnetic and hadronic energy measurements, and extend the coverage to |η| < 4.9.
The muon spectrometer (MS) measures the deflection of muon tracks in the large superconducting air-core toroid magnets in the pseudorapidity range |η| < 2.7. It is instrumented with separate trigger and high-precision tracking chambers. Over most of the η-range, a precision measurement of the track coordinates in the principal bending direction of the magnetic field is provided by monitored drift tubes. At large pseudorapidities, cathode strip chambers with higher granularity are used in the innermost plane over the range 2.0 < |η| < 2.7. The ATLAS trigger system consists of a hardware-based Level-1 (L1) trigger and a software-based High Level Trigger, subdivided into the Level-2 (L2) and Event-Filter (EF) [16] stages. In L1, electrons are selected by requiring adjacent electromagnetic (EM) trigger towers exceed a certain E T threshold, depending on the detector η. The EF uses the offline reconstruction and identification algorithms to apply the final electron selection in the trigger. The Z → e + e − events are selected in this analysis by using a dielectron trigger in the region |η| < 2.5 with an electron transverse energy, E T , threshold of 12 GeV. The muon trigger system, which covers the pseudorapidity range |η| < 2.4, consists of resistive plate chambers in the barrel (|η| < 1.05) and thin gap chambers in the end cap regions (1.05 < |η| < 2.4). Muons are reconstructed in the EF combining L1 and L2 information. The Z → μ + μ − events in this analysis are selected with a first-level trigger that requires the presence of a muon candidate reconstructed in the muon spectrometer with transverse momentum of at least 18 GeV. The trigger efficiency for the events selected as described in Sect. 5 is very close to 100 %.

Monte Carlo simulations
Monte Carlo event samples including a simulation of the ATLAS detector response are used to correct the measurements for detector effects, and to estimate systematic uncertainties. In addition, predictions of different phenomenological models implemented in the MC generators are compared to the data corrected to the particle level. Samples of inclusive Z → e + e − and Z → μ + μ − events were produced using the leading order (LO) Pythia 6 [17], Pythia 8 [18], Herwig++ [19,20], Sherpa [21], Alpgen [22] and next to leading order (NLO) Powheg [23] event generators, including various parton density function (PDF) parametrisations. The Alpgen and Sherpa matrix elements are generated for up to five additional partons, thereby filling the phase space with sufficient statistics for the full set of measured observables. It should be noted, that since the measurements are all reported in bins of p Z T , the results presented in this paper are not sensitive to the predicted shape of the p Z T spectrum, even though they are sensitive to jet activity in the event. Table 2 lists the different MC models used in this paper.
Pythia 6, Pythia 8 and Herwig++ are all leadinglogarithmic parton shower (PS) models matched to leadingorder matrix element (ME) calculations, but with different ordering algorithms for parton showering, and different hadronization models. In scattering processes modelled by lowest-order perturbative QCD two-to-two parton scatters, with a sufficiently low p T threshold, the partonic jet cross-section exceeds that of the total hadronic cross-section. This can be interpreted in terms of MPI. In this picture, the ratio of the partonic jet cross-section to the total crosssection is interpreted as the mean number of parton interactions per event. This is implemented using phenomenological models [24], which include (non-exhaustively) further low-p T screening of the partonic differential cross-section, and use of phenomenological transverse matter-density profiles inside the hadrons. The connection of colour lines between partons, and the rearrangement of the colour structure of an event by reconnection of the colour strings, are implemented in different ways in these phenomenological models.
The Pythia 6 and Pythia 8 generators both use p Tordered parton showers, and a hadronisation model based on the fragmentation of colour strings. The Pythia 8 generator adds to the Pythia 6 MPI model by interleaving not only the ISR emission sequence with the MPI scatters, but also the FSR emissions. The Herwig++ generator implements a cluster hadronization scheme with parton showering ordered  [22] between the matrix element and parton shower to avoid double-counting of jet production mechanisms. A related matching process is used to interface Pythia 6 to the next-to-leading-order (NLO) Powheg generator, where the matching scheme avoids both doublecounting and NLO subtraction singularities [27,28]. Different settings of model parameters, tuned to reproduce existing experimental data, have been used for the MC generators. The Pythia 6, Pythia 8, Herwig + Jimmy, Her-wig++ and Sherpa tunes have been performed using mostly Tevatron and early LHC data. The parton shower generators used with Alpgen and Powheg do not use optimised tunes specific to their respective parton shower matching schemes.
For the purpose of correcting the data for detector effects, samples generated with Sherpa (with the CTEQ6L1 PDF and the corresponding UE tune), and Pythia 8 tune 4C [36] were passed through ATLFAST2 [37], a fast detector simulation software package, which used full simulation in the ID and MS and a fast simulation of the calorimeters. Comparisons between MC events at the reconstructed and particle level are then used to correct the data for detector effects. Since the effect of multiple proton-proton interactions is corrected using a data-driven technique (as described in Sect. 6), only single proton-proton interactions are simulated in these MC samples.

Event selection
The event sample was collected during stable beam conditions, with all detector subsystems operational. To reject contributions from cosmic-ray muons and other non-collision backgrounds, events are required to have a primary vertex (PV). The PV is defined as the reconstructed vertex in the event with the highest p 2 T of the associated tracks, consistent with the beam-spot position (spatial region inside the detector where collisions take place) and with at least two associated tracks with p T > 400 MeV.
Electrons are reconstructed from energy deposits measured in the EM calorimeter and associated to ID tracks. They are required to satisfy p T > 20 GeV and |η| < 2.4, excluding the transition region 1.37 < |η| < 1.52 between the barrel and end-cap electromagnetic calorimeter sections. Electron identification uses shower shape, track-cluster association and TRT criteria [38]. Muons are reconstructed from track segments in the MS associated to ID tracks [39]. They are required to have p T > 20 GeV and |η| < 2.4. Both electrons and muons are required to have longitudinal impact parameter multiplied by sin θ of the ID track, |z 0 | sin θ < 10 mm with respect to the PV. The dilepton invariant mass of oppositely charged leptons, m ll , is required to be in the region 66 < m ll < 116 GeV at this stage. No explicit isolation requirement is applied to the muons, but in the case of electrons, some isolation is implied by the identification algorithm. The correction for this effect is discussed in Sect. 7.3.
The tracks in the calculation of UE observables satisfy the following criteria [40]: -p T > 0.5 GeV and |η| < 2.5; -a minimum of one pixel and six SCT hits; -a hit in the innermost pixel layer, if the corresponding pixel module was active; -transverse and longitudinal impact parameters with respect to the PV, |d 0 | < 1.5 mm and |z 0 | sin θ < 1.5 mm, respectively; -for tracks with p T > 10 GeV, a goodness of fit probability greater than 0.01 in order to remove mis-measured tracks.
The tracks corresponding to the leptons forming the Zboson candidate are excluded.

Correction for pile-up
The average expected number of pile-up events per hardscattering interaction (μ) was typically in the range 3−12 in the 2011 dataset. Of the tracks selected by the procedure described above and compatible with the PV of the hard-scattering event, up to 15 % originate from pile-up, as described below. Due to the difficulty in modelling accurately the soft interactions in pp collisions and the fact that pile-up conditions vary significantly over the data-taking period, a data-driven procedure has been derived to correct the measured observables for the pile-up contribution.
The measured distribution of any track-based observable can be expressed as the convolution of the distribution of this variable for the tracks originating from the Zboson production vertex, with the distribution resulting from the superimposed pile-up interactions. The pile-up contribution is estimated from data by sampling tracks originating from a vertex well separated from the hard-scattering PV. In each event, the pile-up contribution to a given observable is derived from tracks selected with the same longitudinal and transverse impact parameter requirements as the PV tracks, but with respect to two points located at z distances of +2 cm and −2 cm from the hard-scattering PV. The shift of 2 cm relative to the PV introduces a bias in the density of the pile-up interactions. This is corrected on the basis of the shape of the distribution of the z distance between pairs of interactions in the same bunch crossing. This distribution is well approximated by a Gaussian with variance σ = √ 2σ BS , where σ BS ≈ 6 cm is the effective longitudinal variance of the interaction region averaged over all events. Pile-up distributions are thus obtained for each observable and are deconvoluted from the corresponding measured distributions at the hard-scattering PV.
The stability of the pile-up correction for different beam conditions is demonstrated in Fig. 2. The figure compares the distributions of the average charged particle multiplicity density, N ch /δη δφ as a function of p Z T , before and after pile-up correction, for two sub-samples with an average of 3.6 and 6 interactions per bunch crossing ( μ ), respectively. Each distribution is normalised to that obtained for the full sample after pile-up correction. The dependence of the normalised charged multiplicity distributions on p Z T which can be seen before correction in Fig. 2 reflects the fact that actual   contributions to this observable depend on p Z T , while the pileup contribution is independent of p Z T . The pile-up corrected results agree to better than 2 %, a value much smaller than the size of the correction, which may be as large as 20 % for this observable in low p Z T bins for the data-taking periods with the highest values of μ . The systematic uncertainty arising from this procedure is discussed in Sect. 8.

Unfolding to particle level, background corrections and channel combination
After correcting for pile-up, an iterative Bayesian unfolding [41] of all the measured observables to the particle level is performed. This is followed by a correction of the unfolded distributions for the small amount of background from other physics processes. At this point, the electron and muon measurements are combined to produce the final results.

Unfolding
The measurements are presented in the fiducial region defined by the Z -boson reconstructed from a pair of oppositely charged electrons or muons each with p T > 20 GeV and |η| < 2.4 and with a lepton pair invariant mass in the range 66 < m ll < 116 GeV.
The results in Sect. 9 are presented in the Born approximation, using the leptons before QED FSR to reconstruct the Z -boson. These results are also provided in HEPDATA [42] using dressed leptons. These are defined by adding vectorially to the 4-momentum of each lepton after QED FSR the 4momenta of any photons not produced in hadronic decays and found within a cone of R = 0.1 around the lepton, where the angular separation R is given by The UE observables are constructed from stable charged particles with p T > 0.5 GeV and |η| < 2.5, excluding Zboson decay products. Stable charged particles are defined as those with a proper lifetime τ > 0.3 × 10 −10 s, either directly produced in pp interactions or from the subsequent decay of particles with a shorter lifetime.
Bayesian iterative unfolding was used to correct for residual detector resolution effects. This method requires two inputs: an input distribution of the observable (the MC generator-level distribution is used for this), and a detector response matrix which relates the uncorrected measured distribution in this observable to that defined at the event generator level, also termed the particle level. The detector response matrix element, S i j is the probability that a particular event from bin i of the particle-level distribution is found in bin j of the corresponding reconstructed distribution, and is obtained using simulation. For the profile histogram observables in this paper, a two-dimensional (2D) histogram was created with a fine binning for the observable of interest, such that each unfolding bin corresponds to a region in the 2D space.
The unfolding process is iterated to avoid dependence on the input distribution: the corrected data distribution produced in each iteration is used as the input for the next. In this analysis, four iterations were performed since this resulted only in a small residual bias when tested on MC samples while keeping the statistical uncertainties small. The unfolding uses the Sherpa simulation for the input distributions and unfolding matrix. In the muon channel, the MC events are reweighted at the particle level in terms of a multi-variable distribution constructed for each distribution of interest using the ratio of data to detector-level MC, so that the detectorlevel MC closely matches the data. This additional step is omitted in the electron channel for the reasons discussed in Sect. 7.3.
The dominant correction to the data is that related to track reconstruction and selection efficiencies, in particular at lowp T . After the selection described in Sect. 5, the rate of fake tracks (those constructed from tracker noise and/or hits which were not produced by a single particle) is found to be very small. This, as well as a small contribution of secondaries (i.e. tracks arising from hadronic interactions, photon conversions to electron-positron pairs, and decays of long-lived particles) is corrected for by the unfolding procedure.

Backgrounds
The background to the Z -boson signal decaying into a lepton pair consists of a dominant component from multijet production, smaller components from other physics sources, and a very small component from non-collision backgrounds. A fully data-driven correction procedure has been developed and applied directly to the unfolded distributions to take into account the influence of the backgrounds.
The primary vertex requirement removes almost all of the beam-induced non-collision background events. Similarly, the impact parameter requirements on the leptons reduce the cosmic-ray background to a level below 0.1 % of the signal. These residual backgrounds were considered as negligible in the analysis.
The pp collision backgrounds to Z → e + e − or Z → μ + μ − decays were found to be of the order of a few percent of the signal in the mass window [43]. The resonant backgrounds from W Z, Z Z and Z γ pair production with a Z boson decaying into leptons were estimated from simulated samples and found to amount to less than 0.2 % of the selected events. Their impact on the underlying event observables is negligible and they were not considered further here.
The contribution from the non-resonant backgrounds (i.e. from all other pp collision processes) is larger, typically between 1 and 2 % of the signal, depending on the p Z T range considered, and is dominated by multijet production with a combination of light-flavour jets misidentified as electrons and heavy-flavour jets with a subsequent semileptonic decay of a charm or beauty hadron. This contribution is estimated to correspond to 0.5 % of the signal for Z → e + e − decays and to 1-2 % of the signal for Z → μ + μ − decays. The background in the electron channel is somewhat lower because of the implicit isolation requirement imposed on the electrons through the electron identification requirements. Smaller contributions to the non-resonant background arise from diboson, tt and single top production and amount to less than 0.3 % of the signal, increasing to 1 % at p Z T > 50 GeV. The still smaller contributions from processes such as W or Z production with jets, where a jet is misidentified as a lepton, are treated in the same way as the multijet background. These contributions amount to less than 0.1 % of the signal sample.
The non-resonant background is corrected for by studying the UE observables as a function of m ll , the half-width of the mass window around the Z -boson signal peak. Since the distributions of UE observables in non-resonant background processes are found to be approximately constant as a function of the dilepton mass and the background shape under the Z -boson mass peak is approximately linear, the background contribution to any UE observable is approximately proportional to m ll . Thus, the background contribution can be corrected for by calculating the UE observables for different values of m ll , chosen here to be between 10 and 25 GeV, and extracting the results which could be measured for a pure signal with m ll → 0. This procedure is performed separately for each bin of the distributions of interest.
The validity of the linear approximation for the m ll dependence of the background contribution was checked for  all observables studied in this analysis. An example is presented in Fig. 3, where the m ll dependence is shown for one bin of the p T differential distribution, as obtained in the toward region for 30 < p Z T < 35 GeV and shown separately for the electron and muon channels. The values plotted in Fig. 3 are normalised to the corrected combined value. The values of the observables in the muon channel increase linearly with m ll . The difference in the slope observed between the muon and the electron samples is due to the larger background in the muon channel, as discussed above. A straight line is fitted through the points obtained for the various m ll values shown in Fig. 3 for each channel. For each bin in the observable and p Z T , the muon and electron channels values agree with each other after extrapolating to m ll = 0 within the uncertainties of the fit procedure, which are represented by the shaded areas and include the statistical and systematic uncertainties from the fit itself (as discussed in Sect. 8, as well as the relevant correlations.
The effect of the background on the unfolded distributions can be summarised as follows: in the case of the electron channel, which has less background than the muons, the background in the average values of p T and N ch is below 1 %. The absence of any isolation requirement applied to the muons leads to significantly higher background levels in certain regions, with corrections ranging from as high as 6-8 % for the average values of p T in the toward region at high p Z T , to about 1 % for the average values of N ch . The back-ground correction is done after unfolding to avoid resolution issues present at the detector level.

Combination of the electron and muon channels
Before combining the electron and muon channels, the analysis must correct for a bias over a limited region of the phase space which affects the measurements in the electron channel when one of the electrons is close to a jet produced in association with the Z boson. This bias is observed at high p Z T , mostly in the toward region and to a lesser extent in the transverse region, and affects the p T distribution for high values of p T , typically p T > 30 GeV. It arises from the imperfect modelling of the electron shower shape variables in the simulation, which leads to an underestimate of the electron identification efficiency for electrons close to jets. The bias on the observable can be as large as 50 % for p T = 100 GeV. Since it is not reproduced precisely enough by the simulation of the electron shower, in the relevant narrow regions of phase space a tightened isolationcriterion was applied to electrons to exclude the mismodelled event configurations and the proper geometric correction was deduced from the muon channel unaffected by jet overlap. The combined results for electrons and muons in the affected bins are assigned a larger uncertainty, since the contribution of events from the electron-decay channel is significantly reduced leading to a larger overall uncertainty. The most significant effect is observed for the p T > 100 GeV in the toward and transverse region.
As discussed in Sect. 2 and in Sect. 7.1, the electron and muon results are unfolded and then combined, both as Bornlevel lepton pairs and as dressed lepton pairs, and accounting for the uncorrelated and correlated terms in the systematic uncertainties between the channels (as described in Sect. 8).
Combining the dressed electron and muon pairs induces < 0.1 % additional systematic uncertainty on the UE observables compared to the Born level results. Figure 4 illustrates the excellent agreement between the fully unfolded and corrected UE observables for the electron and muon channels, once the specific correction procedure described above has been applied to the electron channel in the limited phase space regions where significant hadronic activity occurs close to one of the electrons. As shown for the specific region 20 < p Z T < 50 GeV in Fig. 4, the differential distributions for p T and N ch agree within statistical uncertainties over most of the range of relevance, except for high values of p T , where the electron bias has been corrected as described above, and where the total uncertainty on the combined measurement has been enlarged as shown by the shaded error band in the ratio plot. The shape of the p T distribution in the region around 1 GeV reflects the p T threshold of 0.5 GeV applied in the track selection. The following sources of uncertainty have been assessed for the measured distributions after all corrections and unfolding. Table 3 summarises the typical sizes of the systematic uncertainties for the UE observables as a function of p Z T .
Lepton selection: systematic uncertainties due to the lepton selection efficiencies have been assessed using MC simulation. The data are first unfolded using the nominal MC samples, then with samples corresponding to a ±1σ variation of the efficiencies [43]. These uncertain-ties are assumed to be uncorrelated between the electron and muon channels. The resulting uncertainty is less than 1 % for all observables over most of the kinematic range.
Track reconstruction: the systematic uncertainty on the track reconstruction efficiency originating from uncertainties on the detector material description is estimated as in Ref.
[40] for |η| > 2.1. The typical value for |η| < 2.1 is ±1 % while it is approximately 5 % for |η| > 2.1. The effect of this uncertainty on the final results is less than 2 %. This uncertainty is fully correlated between the electron and muon channels. Impact parameter requirement: the fraction of secondary particles (i.e. those originating from decays and interactions in the inner detector material) in data is reproduced by the MC simulation to an accuracy of ∼ 10-20 %, obtained by comparing d 0 distributions in MC and in the data corrected for pile-up. To assess the corresponding systematic uncertainty, the track impact parameter requirements on |d 0 | and |z 0 |sinθ are varied from the nominal values of 1.5 to 1.0 and 2.5 mm, resulting in fractions of secondaries varying between 0.5 to 4.0 %, and the resulting distributions are unfolded using MC samples selected with the same impact parameter requirements. The maximum residual difference of 2 % or less between these unfolded distributions and the nominal unfolded distribution is taken as the uncertainty arising from this requirement. This uncertainty is also fully correlated between the electron and muon channels.
Pile-up correction: the pile-up correction uncertainty originates from the uncertainty in the pile-up density fitted along with the spatial distribution of tracks originating from pile-up, and the difference between the pile-up densities measured for Z -boson and for randomly triggered events. In addition to these, the stability of the correction method with respect to the instantaneous luminosity was estimated by performing the correction procedure independently on datasets with different average numbers of reconstructed vertices, as shown in Fig. 2. The total uncertainty due to the pile-up correction is taken to be the quadratic combination of the uncertainties from these sources, and it is at most 2 % for the average underlying event observables. The overall uncertainty is fully correlated between the electron and muon channels.
Background correction: the uncertainty is evaluated by comparing the results of the linear fit to those obtained using a second-order polynomial. This uncertainty is at most 2 % for the maximum background uncertainty on p T , which is the most strongly affected variable, and is assumed to be uncorrelated between the electron and muon channels. Any potential correlation arising from the common tt and diboson backgrounds is neglected  because they become sizable only for p Z T > 100 GeV, where the total uncertainty is dominated by the statistical uncertainity on the background.
Unfolding: the uncertainty due to the model-dependence of the unfolding procedure is taken from the degree of non-closure between the Pythia 8 initial particlelevel distributions and the corresponding detector-level Pythia 8 distributions unfolded and corrected using the Sherpa sample, which was reweighted to agree with Pythia 8 at the detector level. This uncertainty varies between 0.5 and 3 % for the profile distributions, and is assumed to be uncorrelated between the electron and muon channels.
Bias due to implicit isolation: this uncertainty is estimated by varying the electron isolation requirement used to derive the correction discussed in Sect. 7.3. The uncertainty is assigned to the electron channel and does not exceed ∼1 % for the profile distributions.
Other potential sources of systematic uncertainty have been found to be negligible. The total uncertainty in each measured bin is obtained by propagating the systematic component of the error matrix through the channel combination. For the differential distributions in Sect. 9.2, the unfolding model dependent uncertainty increases to about 5 %, resulting in slightly larger overall systematic uncertainties.

Overview of the results
The results are shown in Sect. 9.2, first for the differential distributions of charged particle p T and N ch in intervals of p Z T , and then for the same distributions for a representative p Z T range compared to MC model predictions. The normalised quantities, N ch /δη δφ and p T /δη δφ, are obtained The angular areas for the transverse, toward, and away region observables are δφ δη = (2 × π/3) × (2 × 2.5) = 10π/3, while for trans-max/min/diff, δφ δη = 5π/3.
Since the away region is dominated by the jets balancing the p Z T [43], the focus will be on the toward, transverse, trans-max and trans-min regions. In the transverse region, the extra jet activity is more likely to be assigned to the transmax region. Assuming the same flat UE activity in trans-min and trans-max regions, the trans-diff region, the difference between the observables measured in trans-max and transmin regions, is expected to be dominated by the hard scattering component. In Sect. 9.3 profile histograms are shown. Finally, in Sect. 9.4, the results are compared to previous measurements from ATLAS where distributions sensitive to the underlying event were measured as a function of the kinematics of either the leading charged particle [1], or the leading jet [5].  Fig. 6. At values below p T /δη δφ of 0.1 GeV, the distributions exhibit a decrease, which is independent of p Z T . This is followed by a sharp increase at higher p T /δη δφ, which is an artifact of requiring at least two tracks with p T of at least 0.5 GeV in every event. Then a broad distribution can be seen extending to p T /δη δφ of about 1 GeV, followed  Table 2). For p T /δη δφ < 0.1 GeV, there is a large spread in the predictions of the MC models relative to the data, with Powheg providing the best description. The intermediate region with 0.1 < p T /δη δφ < 1 GeV, is well reproduced by most of the MC models. For the higher p T /δη δφ ranges, most of the MC models underestimate the number of events, with the exception of Sherpa and Alpgen, which have previously been shown to provide good models of multi-jet produced in association with a Zboson [43]. This observation may indicate that even the transmin region is not free of additional jets coming from the hard scatter.
The distributions of the charged particle multiplicity density in the four UE regions are shown in Figs. 9 and 10 for the same p Z T intervals used in Figs. 5 and 6, respectively. The distributions in the transverse, toward and trans-max regions exhibit similar features, with the exception of the largest multiplicities, which are suppressed in the trans-min region, com- pared to the trans-max one. In the trans-min region, as for the p T /δη δφ distribution, limited dependence on p Z T is observed at low multiplicity. The suppression of large multiplicities in the trans-min region is more pronounced in the lower p Z T intervals. The comparison of these multiplicity distributions to various MC models, in the same p Z T interval, between 20-50 GeV, is shown in Figs. 11 and 12 for all the UE regions. In contrast to the p T /δη δφ distributions, none of the MC models, except Pythia 8, describes the data distributions, in particular for N ch /δη δφ > 2.

Average distributions
The evolution of the event activity in the four UE regions with the hard scale can be conveniently summarised by the average value of the UE observables as a function of p Z T . In Fig. 13 Fig. 11 Comparisons of data and MC predictions for charged particle multiplicity density, N ch /δη δφ, for Z -boson transverse momentum, p Z T , in the interval 20-50 GeV, in the toward (a) and transverse (b) regions. The bottom panels in each plot show the ratio of MC predictions to data. The shaded bands represent the combined statistical and systematic uncertainties, while the error bars show the statistical uncertainties in the away region. This difference increases with increasing p Z T . The away region density is large due to the presence in most cases of a jet balancing the Z -boson in p T . The density in the transverse region is seen to be systematically higher than that in the toward region, which can be explained by the fact that for high p Z T , additional radiated jets balancing p Z T affect the transverse region more than the toward region [43]. The difference between the three regions disappears at low p Z T due to the fact that the UE regions are not well defined with respect to the actual Z -boson direction.
In Fig. 13, p T /δη δφ is seen to rise much faster as a function of p Z T in the trans-max region than in the trans-min region. The slowing down of the rise of p T /δη δφ at high The comparison of the p T /δη δφ distribution as a function of p Z T with the predictions of various MC models is shown in Figs. 14 and 15 in the UE regions sensitive to the underlying event characteristics. For clarity of comparison, the statistically least significant p Z T > 210 GeV bin is omitted. The variation in the range of predictions is quite wide, although less so than for the differential p T distributions. The best description of the transverse  13 The average values of charged particle scalar p T density, p T /δη δφ , as a function of Z -boson transverse momentum, p Z T , in the transverse, toward and away regions (a), and in the trans-max, transmin and trans-diff regions (b). The results are plotted at the center of each p Z T bin. The error bars depict combined statistical and systematic uncertainties and trans-max regions is given by Sherpa, followed by Pythia 8, Alpgen and Powheg. The observation that the multi-leg and NLO generator predictions are closer to the data than most of the pure parton shower generators suggests that these regions are affected by the additional jets coming from the hard interaction. Jet multiplicities in events with a Z -boson have been studied by the LHC experiments [43], and they are well described by Sherpa and Alpgen.
The discrepancy between the Pythia 8 AU2 tune and the Pythia 6 Perugia tune possibly indicates the effect of using LHC UE data for the former in addition to the shower model improvement. In the trans-min region, which is the most sensitive to the UE, none of the models fully describe the data. Apart from Herwig++, and Sherpa, which predicts a faster rise of p T than observed in data, the other generators model the data better in the trans-min region than they do in the transverse or trans-max regions. This possibly indicates that in the LO shower generators the underlying event is well modelled but perturbative jet activity is not.
In Fig. 16, N ch /δη δφ is shown as a function of p Z T in the different UE regions. The profiles behave in a similar way to p T /δη δφ . However, the trans-diff N ch /δη δφ activity is lower than that for trans-min, while for p T /δη δφ , it is the other way around. This indicates that the trans-diff region, which is a measure of extra activity in the trans-max region over the trans-min region, is populated by a few particles In Figs. 17 and 18, in which various MC model predictions are compared to N ch /δη δφ as a function of p Z T , a different pattern from that of p T /δη δφ is observed. The Pythia 6 Perugia 2011C tune and Alpgen provide the closest predictions in all three regions. Sherpa, Pythia 8 and Powheg predict higher average multiplicities, with Sherpa being the farthest from the data. On the other hand, Herwig++ mostly underestimates the data. The p T /δη δφ and N ch /δη δφ distributions as functions of p Z T in the trans-diff region are compared with the MC model predictions in Fig. 19. While all MC models, except for Herwig++ predict the multiplicity fairly well, only Sherpa and Alpgen predict the p T average values well in certain ranges. The better modelling of this region by MC models with additional jets coming from matrix element rather than from parton shower again confirms that the trans-diff region is most sensitive to the additional radiated jets.
The difficulty of describing the p T /δη δφ and N ch /δη δφ average values simultaneously in MC models is reflected in the comparison of data and MC model predictions for p T in Fig. 20. The p T as a function of p Z T is reasonably described by Alpgen and Sherpa for high p Z T , while all the other models predict softer spectra. The correlation of p T with N ch , shown in Fig. 21, follows the pattern Overall the Pythia 8 prediction is the closest to the data, followed by Pythia 6 and Powheg, although for N ch < 5, all three have much softer distributions than the data. The other models do well in this low N ch region, but are then much lower than the data for high N ch .
From all the distributions considered, it can be inferred that the jets radiated from the hard scatter will affect the  (b) Fig. 18 Comparison of data and MC predictions for charged particle multiplicity density average values, N ch /δη δφ , as a function of Zboson transverse momentum, p Z T , in the trans-max (a) and trans-min (b) regions. The bottom panels in each plot show the ratio of MC predictions to data. The shaded bands represent the combined statistical and systematic uncertainties, while the error bars show the statistical uncertainties underlying event observables and therefore these must be properly reproduced in order to obtain an accurate MC description of the UE. The UE region least affected by the presence of extra jets is the trans-min region.

Comparison with other ATLAS measurements
The results from this analysis are compared to the results obtained when the leading object is either a charged particle [1] or a hadronic jet [5]. The underlying event analysis with a leading charged particle was performed with the early The differential N ch /δη δφ and p T /δη δφ distributions for leading jet and Z -boson events are compared in Figs. 22 and 23 for the trans-max and trans-min regions. While the N ch /δη δφ distributions are similar, a clear difference is observed in the high tails of the p T /δη δφ distribution, which are more populated in Z -boson events than in jet events. This difference was traced to the definition of the leading object. In the case of jets, the accompanying activity can never contain jets with a p T higher than that of the leading jet, whereas there is no such restriction for Z -boson   20 Comparison of data and MC predictions for charged particle mean p T as a function of Z -boson transverse momentum, p Z T , in the toward (a) and transverse (b) regions. The bottom panels in each plot show the ratio of MC predictions to data. The shaded bands represent the combined statistical and systematic uncertainties, while the error bars show the statistical uncertainties events. As a test, the average p T was determined for Zboson events after rejecting all events in which at the detector level there was a jet with p T higher than the p Z T , with jets selected as in [5]. The average was found to be about 20-30 % lower than for the standard selection, and the average values in jet and Z -boson events are in close agreement in this case.
The hard scales used for the analyses are different and the choice of the main observable used to assess the evolution of the underlying event reflects this to a certain extent in the figures. Nevertheless, certain common qualita-   21 Comparison of data and MC predictions for charged particle mean p T as a function of charged particle multiplicity, N ch , in the toward (a) and transverse (b) regions. The bottom panel in each plot shows the ratio of MC predictions to data. The shaded bands represent the combined statistical and systematic uncertainties, while the error bars show the statistical uncertainties tive features can be observed by comparing p T /δη δφ and N ch /δη δφ as functions of the leading object p T in the transverse region, and also separated into the transmax/min regions as shown in Figs. 24 and 25. The measurements with a leading jet are complementary to the measurements with a leading track, and a smooth continuation at 20 GeV is observed (in Fig. 24), corresponding to the lowest jet p T for which the jet measurement could be performed and the highest leading track momentum included in the leading track analysis. Where the p T of the leading object is less than 50 GeV, a large difference is observed both for the N ch and p T average values between the jet and Z -boson measurements in Fig. 24; the increase Although the N ch density is similar in the underlying event associated with a jet to that with a Z -boson for higher values of the hard scale (≥50 GeV), there are residual differences in the average p T densities. The activity in events with a Z -boson is systematically higher than that in events with jets. From the behaviour of the underlying event properties in the trans-max/min regions in Fig. 25, this difference originates mostly from the trans-max region, due to selection bias discussed previously in this section. The trans-min region is very similar between the two measurements, despite the different hard scales, indicating again that this region is least sensitive to the hard interaction and most sensitive to the MPI component.

Conclusion
Measurements sensitive to the underlying event have been presented, using an inclusive sample of Z -boson decays, obtained from a data set collected in proton-proton collisions at the LHC corresponding to an integrated luminosity of 4.6 fb −1 . The transverse and toward regions with respect to the reconstructed Z -boson are most sensitive to the underlying event. The transverse region was further subdivided into trans-max and trans-min regions on an event-by-event basis depending on which one had a higher p T ; this subdivision provides additional power to discriminate between the different processes contributing to the underlying event models.
The results show the presence of a hard component in the p T distribution of particles, presumably originating from extra jet activity associated with the Z -boson production. It and Zboson transverse momentum, p Z T , in the transverse region. The error bars in each case show the combined statistical and systematic uncertainties. The insets show the region of transition between the leading charged particle and leading jet results in more detail is observed in all the investigated regions, with the transmin region least affected by it. The average underlying event activity increases with p Z T , until it reaches a plateau, which is again most prominent in the trans-min region. The results have been compared to a number of MC models, using several tunes of commonly used underlying event models. MC model predictions qualitatively describe the data well, but with some significant discrepancies, providing precise information sensitive to the choices of parameters used in the various underlying-event models. Careful tuning of these parameters in the future may improve the description of the data by the different models in future LHC measurements and studies. The study of such variables in Z -boson events provides a probe of the underlying event which is complementary to that from purely hadronic events. A comparison between them shows similar underlying event activity for the transmin region.