Measurement of the underlying event in the Drell-Yan process in proton-proton collisions at sqrt(s) = 7 TeV

A measurement of the underlying event (UE) activity in proton-proton collisions at a center-of-mass energy of 7 TeV is performed using Drell--Yan events in a data sample corresponding to an integrated luminosity of 2.2 inverse femtobarns, collected by the CMS experiment at the LHC. The activity measured in the muonic final state (q q-bar to opposite-sign muons) is corrected to the particle level and compared with the predictions of various Monte Carlo generators and hadronization models. The dependence of the UE activity on the dimuon invariant mass is well described by PYTHIA and HERWIG++ tunes derived from the leading jet/track approach, illustrating the universality of the UE activity. The UE activity is observed to be independent of the dimuon invariant mass in the region above 40 GeV, while a slow increase is observed with increasing transverse momentum of the dimuon system. The dependence of the UE activity on the transverse momentum of the dimuon system is accurately described by MADGRAPH, which simulates multiple hard emissions.


Introduction
In hadron-hadron scattering, the "underlying event" (UE) is defined as any hadronic activity that cannot be attributed to the particles originating from the hard scattering, which is characterized by a large momentum transfer, or to the hadronization of initial-and final-state radiation. The UE activity is thus due to the hadronization of partonic constituents, not involved in the hard scattering, that have undergone multiple-parton interactions (MPIs) and to the hadronization of beam remnants that did not participate in other scatterings. These semihard interactions cannot be completely described by perturbative quantum chromodynamics (QCD) and require a phenomenological description involving parameters that must be tuned with the help of data [1].
The experimental study of the UE probes various aspects of hadron production in high energy hadron-hadron collisions. In particular it is sensitive to the interplay of perturbative methods describing the hard process and phenomenological models of the soft interactions that attempt to simultaneously describe MPIs, initial-and final-state radiation, the colour flow between final state partons, and the hadronisation process. Understanding the UE in terms of particle and energy densities will lead to better modelling by Monte Carlo programs that are used in precise measurements of standard model processes and searches for new physics at high energies. The UE affects the estimation of the efficiency of isolation criteria applied to photons and charged leptons, and the energy scale in jet identification. It also affects the reconstruction efficiency for processes like H→ γγ, where the primary vertex is partly determined from the charged particles originating from the UE. Hard MPIs are an important background for new physics searches, e.g. same-sign W production from MPIs [2] is a possible background to the same-sign double lepton SUSY searches [3].
The Compact Muon Solenoid (CMS) [4], ATLAS, and ALICE experiments have carried out UE measurements at centre-of-mass energies ( √ s) of 0.9 TeV and 7 TeV using hadronic events (minimum-bias and single-jet triggered) containing a leading track-jet [5,6] or a leading track [7,8]. The analysis of the central charged particles and forward energy flow correlations in hard processes, e.g. pp → W(Z)X → ν( )X [9], provides supplementary insights into the nature of MPIs. In this paper, we use the Drell-Yan (DY) process [10] with the muonic final state at a centre-of-mass energy of 7 TeV to perform a complementary UE measurement. The DY process with muonic final state is experimentally clean and theoretically well understood, allowing the particles from the UE to be reliably identified. The absence of QCD final-state radiation (FSR) permits a study of different kinematic regions with varying transverse momentum of γ * /Z due to harder or softer initial-state radiation (ISR). The comparison of the UE measurement in DY events with QCD events having a leading track-jet is useful for probing the UE activity in different processes. UE measurements using the DY process have been reported previously in proton-antiproton collisions at √ s = 1.96 TeV [11].
The UE activity at a given centre-of-mass energy is expected to increase with the momentum transfer of the interaction. Events with a harder scale are expected to correspond, on average, to interactions with a smaller impact parameter and, in some models, to more MPIs [12,13]. This increased activity is observed to reach a plateau for high energy scales corresponding to small impact parameter. In this paper we investigate some aspects of the UE modelling in detail by measuring the invariant mass dependence of the UE activity for DY events with small transverse momentum of the DY system. This measurement separates the scale dependence of the UE activity from the ISR effect. The universality of the model parameters, denoted as tunes, implemented in the various MC programs is tested by comparing their predictions with our measurements. The portability of the UE parameters across different event generators, combined in some cases with different parton distribution functions (PDFs), is investigated as well. The modelling of the ISR is studied by measuring the UE activity as a function of the transverse momentum of the DY system. Finally, the dependence of the UE activity on ISR and FSR is determined by comparing the measurements from DY events with previous results from hadronic events containing a leading jet where FSR also plays a role.
The outline of the paper is as follows. Section 2 describes the various observables used in the present study. Section 3 summarizes the different MC models used and corresponding UE parameters. Section 4 presents experimental details: a brief detector description, data samples, event and track selection criteria, correction procedure, and systematic uncertainties. Section 5 presents the results on UE activity measured in DY events and the comparison with the measurements based on a leading track-jet. The main results are summarized in Section 6.

Observables
The UE activity is measured in terms of particle and energy densities. The particle density (1/[∆η∆(∆φ)] N ch ) is computed as the average number of primary charged particles per unit pseudorapidity η and per unit azimuthal separation ∆φ (in radians) between a track and the transverse momentum of the dimuon system. The pseudorapidity is defined as η = − ln(tan(θ/2)), where θ is the polar angle measured with respect to the anticlockwise beam direction. The azimuthal angle φ is measured in the plane perpendicular to the beam axis. The energy density (1/[∆η∆(∆φ)] Σp T ) is expressed in terms of the average of the scalar sum of the transverse momenta of primary charged particles per unit pseudorapidity per unit azimuthal separation. The ratio of the energy and particle densities, as well as the total charged-particle multiplicity N ch and the transverse momentum spectrum are also computed. The chargedparticle multiplicity and transverse momentum distributions are normalized to unit area and to the average number of charged particles per event, respectively. Particles are considered as primary if they originate from the initial proton-proton interaction and are not the decay products of long-lived hadrons with a lifetime exceeding 10 −10 s. Apart from the muons from the DY process, all charged particles in the central region of the detector with pseudorapidity |η| < 2 and with transverse momentum p T > 0.5 GeV/c are considered.
The spatial distribution of the tracks is categorized by the azimuthal separation ∆φ. Particle production in the away region (|∆φ| > 120 • ) is expected to be dominated by the hardest ISR emissions, which balance the dimuon system. The transverse region (60 • < |∆φ| < 120 • ) and towards region (|∆φ| < 60 • ) are more sensitive to soft emissions and, in particular, those due to MPIs. The relevant information about the hard and the soft processes is extracted from the tracking and the muon systems of the CMS detector and thus the derived observables are insensitive to the uncertainties of the calorimetric measurements. The DY events with dimuon mass M µµ around the Z resonance are the least contaminated by background processes (heavyquark, tt, W+jets, and DY → ττ production) [14,15] and best suited for the measurement of the UE activity.
The UE activity is studied as a function of the magnitude of the dimuon transverse momentum (p µµ T = | p µ T,1 + p µ T,2 |) and as a function of M µµ . The dependence of the UE activity on p µµ T for high-mass dimuon pairs effectively probes the ISR spectrum. In order to minimize the background contamination, the p µµ T dependence is studied only in the narrow mass window 81 < M µµ < 101 GeV/c 2 . In contrast to the study of the UE activity in hadronic events using a leading track-jet [5,6], this energy scale is sufficiently large to saturate the MPI contributions. This observation is verified by studying the UE activity as a function of the dimuon mass in a wider mass range, where the total transverse momentum of the dimuon system is kept to a minimum by requiring p µµ T < 5 GeV/c.

Monte Carlo models
The UE dynamics are studied through the comparison of the observables in data with various tunes of PYTHIA6 [16] and its successor PYTHIA8 [17,18]. MADGRAPH (version 5) [19,20], which simulates up to six final-state fermions (including the muons), and POWHEG [21], which includes next-to-leading-order corrections on the hardest emission, are also compared to our measurements. For these two generators, softer emissions are simulated by p T -ordered parton showers using PYTHIA6 tunes and matched with the hard process produced by the generators. Hadronization in PYTHIA6 and PYTHIA8 is based on the Lund string fragmentation model [22]. The measurements are also compared to predictions of the HERWIG++ [23] angular-ordered parton shower and cluster hadronization model [24,25].
The UE contributions from MPIs rely on modelling and tuning of the parameters in the MC generators. The MPI model of PYTHIA relies on two fundamental assumptions [12]: • The ratio of the 2→2 partonic cross section, integrated above a transverse momentum cutoff scale, and the total of the hadronic cross section is a measure of the amount of MPIs. The cutoff scale p 0T is introduced to regularize an otherwise diverging partonic cross section, with Here √ s 0 = 1.8 TeV and is a parameter characterizing the energy dependence of the cutoff scale.
• The number of MPIs in an event has a Poisson distribution with a mean that depends on the overlap of the matter distribution of the hadrons in impact-parameter space.
The MPI model used here [26] includes showering of the MPI process, which is interleaved with the ISR.
The tunes of the models vary mainly in the MPI regularization parameters, p 0T and , in the amount of colour reconnection, and in the PDF used. The Z1 tune [27] of PYTHIA6 adopts the results of a global tuning performed by the ATLAS Collaboration [28] and uses the fragmentation and colour reconnection parameters of the ATLAS AMBT1 tune [29]. The parameters of the Z1 tune related to the MPI regularization cutoff and its energy dependence are adjusted to describe previous CMS measurements of the UE activity in hadronic events [6] and uses the CTEQ5L PDF. The Z2 tune of PYTHIA6 is an update of the Z1 tune using CTEQ6L1 [30], the default used in most CMS generators; the regularization cutoff value at the nominal energy of √ s 0 = 1.8 TeV is optimized to 1.832 GeV/c. The value of the energy evolution parameter for the Z2 tune is 0.275, as for the Z1 tune. The 4C [31] tune of PYTHIA8 follows a similar procedure as the ATLAS AMBT1 tune, but includes ALICE multiplicity data as well. The values of the p 0T ( √ s 0 ) and parameters for the 4C tune are 2.085 GeV/c and 0.19, respectively. The effective value of p 0T at √ s = 7 TeV is about 2.7 GeV/c for both the Z2 and 4C tunes.
A comparison of these models with the measurements is presented in Section 5.

Experimental methods
The present analysis is performed with a sample of proton-proton collisions corresponding to an integrated luminosity of 2.2 fb −1 , collected in March-August 2011 using the CMS detector [4].
Muons are measured in the pseudorapidity range |η| < 2.4 with a detection system consisting of three subsystems: Drift Tubes, Cathode Strip Chambers, and Resistive Plate Chambers. Matching track segments from the muon detector to the tracks measured in the inner tracker results in a transverse momentum resolution between 1% and 5% for p T values up to 1 TeV/c. The tracker subsystem consists of 1440 silicon-pixel and 15 148 silicon-strip detector modules, and it measures charged particle trajectories within the nominal pseudorapidity range |η| < 2.5. The tracker is designed to provide a transverse impact parameter resolution of about 100 µm and a transverse momentum resolution of about 0.7% for 1 GeV/c charged particles at normal incidence (η = 0).
The detector response is simulated in detail using the GEANT4 package [33]. The simulated signal and background events, including heavy-quark, tt, W+jets, and DY → ττ production, are processed and reconstructed in the same manner as collision data.

Event and track selection
The trigger requires the presence of at least two muon candidates. In periods of lower instantaneous luminosity both muons were required to have p T > 7 GeV/c, while in other periods the transverse momentum requirements were 13 GeV/c and 8 GeV/c for the leading and subleading muons, respectively. The trigger efficiency is above 95% for the offline selected DY events with the requirement of 81 < M µµ < 101 GeV/c 2 . The offline selection requires exactly two muons reconstructed in the muon detector and the silicon tracker. Muon candidates are required to satisfy identification criteria based on the number of hits in the muon stations and tracker, transverse impact parameter with respect to the beam axis, and normalized χ 2 of the global fit [15]. The backgrounds from jets misidentified as muons and from semileptonic decays of heavy quarks are suppressed by applying an isolation condition on the muon candidates. The isolation variable I for muons is defined as where the sum is defined in a cone of radius ∆R = (∆φ) 2 + (∆η) 2 = 0.3 around the muon direction; ∆η and ∆φ are the pseudorapidity and azimuthal separation between the muon and tracks or calorimetric towers. Here p T (tracks) is the transverse momentum of tracks, excluding muons, with p T > 1 GeV/c, E T (EM) is the transverse energy deposited in the electromagnetic calorimeter, E T (HAD) is the transverse energy deposited in the hadronic calorimeter, and ρ is the average energy density [34] in the calorimeter and tracker originating from additional inelastic pp interactions (pile-up) in the same bunch crossing as the DY interaction.The calculation of ρ takes into account the number of reconstructed primary vertices in the event; the average value of ρ is 5.6 GeV/c. A muon is considered to be isolated if I < 0.15. Because of the energy density correction, the isolation efficiency is independent of the number of pile-up interactions.
The selected muons are required to have opposite charges, transverse momenta larger than 20 GeV/c, and pseudorapidity |η| < 2.4. Both muons are required to be associated with the same vertex, which is designated as the signal vertex. The selected signal vertex is required to be within ±18 cm of the nominal interaction point as measured along the z direction. At least five tracks are required to be associated with the signal vertex, and the transverse displacement of the signal vertex from the beam axis is required to be less than 2 cm. These criteria select a pure sample of DY events with a total background contribution of less than 0.5% as estimated from simulated events.
Tracks, excluding the selected muons, are considered for the UE measurement if they are well reconstructed in the silicon-pixel and the silicon-strip tracker, have p T > 0.5 GeV/c and |η| < 2, and originate from the signal vertex. To reduce the number of improperly reconstructed tracks, a high purity reconstruction algorithm [35] is used. The high purity algorithm requires stringent cuts on the number of hits, the normalized χ 2 of the track fit, and the consistency of the track originating from a pixel vertex. To reduce the contamination of secondary tracks from decays of long-lived particles and photon conversions, the distances of closest approach between the track and the signal vertex in the transverse plane and in the longitudinal direction are required to be less than 3 times the respective uncertainties. Tracks with poorly measured momenta are removed by requiring σ(p T )/p T < 5%, where σ(p T ) is the uncertainty on the p T measurement.
These selection criteria reject about 10% of primary tracks and 95% of misreconstructed and secondary tracks. The selected tracks have a contribution of about 2% from misreconstructed and secondary tracks.

Corrections and systematic uncertainties
The UE observables, discussed in Section 2, are corrected for detector effects and selection efficiencies. The measured observables are corrected to reflect the activity from all primary charged particles with transverse momentum p T > 0.5 GeV/c and pseudorapidity |η| < 2. The particle and energy densities are corrected using a bin-by-bin technique. In the bin-bybin technique, the correction factor is calculated by taking the bin-by-bin ratio of the particle level and detector level distributions for simulated events and then the measured quantity is multiplied by this correction factor. There is a small growth in the particle and energy densities with increasing p µµ T and M µµ in the towards and transverse regions. Because of this slow growth of densities the bin migration in p µµ T and M µµ has a small effect on the measurements, therefore a bin-by-bin method is considered to be sufficiently precise. There is a fast rise in the energy and particle densities in the away region with the increase of p µµ T , but corrected results using a bin-by-bin method are consistent with correction obtained from a Bayesian [36] technique. The transverse momenta of the charged particles have very good resolution and are corrected using a bin-by-bin method. In this analysis the average of the calculated correction factors from PYTHIA6 Z2, PYTHIA6 D6T, and MADGRAPH Z2 is used to correct the experimental distributions. The maximum deviation from the average correction factor is taken as the model-dependent systematic uncertainty, estimated to be 0.7-1.4% for the particle and energy densities. In the case of charged-particle multiplicity, there is substantial bin migration and the corrected results using the Bayesian [36] and bin-by-bin techniques differ by 10-15%. Therefore the charged-particle multiplicity is corrected using a Bayesian unfolding technique with a response matrix obtained using the PYTHIA6 Z2 tune. The systematic uncertainty related to the correction procedure is calculated by unfolding the data with response matrices obtained using different tunes.
In the analyzed data, there are on average 6-7 collisions in each bunch crossing. Tracks originating from these pile-up interactions cause the UE activity to be overestimated, so the measurements are corrected for the presence of pile-up interactions. The correction factor is calculated as the ratio of the UE activity for simulated events with and without pile-up. The uncertainty in the modelling of the pile-up events is estimated by varying the mean of the expected number of pile-up events by ±1. This uncertainty in pile-up modelling affects the particle and energy densities by 0.3-1.0%. The effect due to pile-up events is small because only the tracks associated with the same vertex as the muon pair are used. The results are also cross-checked with low pile-up 7 TeV data collected during 2010 and the differences are found to be negligible. Table 1: Summary of the systematic uncertainties on the particle and energy densities (in percent). The first three rows show the systematic uncertainties for the particle density in the towards, transverse, and away regions. The last three rows report the systematic uncertainties for the energy density. The numbers outside the parentheses refer to the case where the densities are measured as a function of M µµ and those in the parentheses correspond to the measurements as a function of p µµ T . We also consider possible systematic effects related to trigger requirements, different beam-axis positions in data and simulation, various track selection criteria, muon isolation, and misidentification of tracks. The combined systematic uncertainty related to trigger conditions, the varying beam-axis position, and track selection is less than 0.5%. The systematic uncertainty due to isolation is calculated by removing the isolation condition in the simulated events used for the correction and is found to be 0.8-2.5% for the particle and energy densities.
The yield of secondary tracks originating from the decay of long-lived particles is not correctly predicted by the simulation [37]. To estimate the effect of secondary tracks, a subset of simulated events is created by rejecting tracks that do not have a matching primary charged particle at the generator level. The uncertainty is evaluated by correcting the measurements with this subset of the simulated events, containing fewer secondary tracks, and is found to be 0.7-1.0% for the particle and energy densities.
Though the total contribution of background processes is very small, it affects the measurement at higher p µµ T (50-100 GeV/c) and small M µµ (40-60 GeV/c 2 ) where the contamination from tt and DY→ ττ background processes is 1% and 5%, respectively. The particle and energy densities differ between DY→ ττ and DY→ µµ (the signal process) by 20%. The particle (energy) density for the tt background is two times (four times) that for the signal process. Combination of the differences in the densities for background processes and relative background contributions gives a systematic uncertainty of 0.2-0.9%. Table 1 summarizes the dominant systematic uncertainties on the particle and energy densities. The total systematic uncertainty on the particle and energy densities is in the range 1.5-3.0%, whereas the uncertainties on the track multiplicity and p T spectra reach 10% in the tail (not reported in Table 1). In all figures, inner error bars represent the statistical uncertainty only, while outer error bars account for the quadratic sum of statistical and systematic uncertainties.

Results
The UE activity in DY events, for charged particles with p T > 0.5 GeV/c and |η| < 2.0, is presented as a function of M µµ and p µµ T . The multiplicity and the transverse momentum distributions are also presented for two different sets of events, p µµ T < 5 GeV/c and 81 < M µµ < 101 GeV/c 2 . Finally, the UE activity in the transverse region is compared with that measured in hadronic events using a leading track-jet.

Underlying event in the Drell-Yan process
The energy-scale dependence of the MPI activity is studied by limiting the ISR. To accomplish this we require the muons to be back-to-back in the transverse plane with p µµ T < 5 GeV/c and measure the dependence of the UE activity on the dimuon mass, M µµ . The resulting particle and energy densities are shown in Fig. 1. Because the activity is almost identical in the towards and transverse regions, they are combined as |∆φ| < 120 • . The contribution of ISR to the UE activity is small after requiring p µµ T < 5 GeV/c, as shown by the prediction of HERWIG++ without MPIs. This figure also illustrates the dominant role of MPIs in our current models as they generate more than 80% of the UE activity in these ISR-reduced events. The lack of dependence of the UE activity on M µµ within the range under study (40-140 GeV/c 2 ) indicates that the activity due to MPIs is constant at energy scales down to 40 GeV. The quantitative description by model tunes based on the minimum-bias and UE observables in hadronic events is illustrated by the MC/Data ratios in Fig. 1. In general, PYTHIA6 Z2, PYTHIA8 4C, and HERWIG++ LHC-UE7-2 describe the densities well, whereas the Z2 tune used together with the POWHEG generator underestimates both densities by 5-15%. Both PYTHIA and HERWIG++ model tunes derived from the UE measurement in hadronic events using the leading jet/track approach describe the UE activity in the Drell-Yan events equally well and hence illustrate a certain universality of the underlying event across QCD and electroweak processes in hadronic collisions.
Dependence of the UE activity on the transverse momentum of the dimuon system is shown in Fig. 2 in the towards, transverse, and away regions (top to bottom) for events having M µµ between 81 GeV/c 2 and 101 GeV/c 2 . At this high energy scale, the p µµ T dependence of the UE activity is sensitive to the ISR. The slope in the p µµ T dependence of the UE activity is identical for a model with and without MPIs and is therefore mainly due to ISR. The predictions of HERWIG++ without MPIs underestimate the measurements in the away region as well because the MPIs produce particles uniformly in all directions. The UE activity does not fall to zero when p µµ T → 0 because of the presence of the hard scale set by M µµ . The particle and energy densities in the away region rise sharply with p µµ T and, because of momentum conservation mainly sensitive to the spectrum of the hardest emission, are equally well described by all tunes and generators considered. In the towards and transverse regions there is a slow growth in the particle and energy densities with increasing p µµ T . The energy density increases more than the particle density, implying a continuous increase in the average transverse momentum of the charged particles with p µµ T . This effect is also reflected in the ratio of the energy density to the particle density. The activity in the towards region is qualitatively similar to that in the transverse region. Quantitatively, the activity is higher in the transverse region than the towards region, an effect caused by the spill-over contributions from the recoil activity in the away region, which balances the dimuon system. This observation is visible in particle density energy density ratio of energy and particle densities  The UE activity as a function of the dimuon invariant mass (M µµ ) for events with p µµ T < 5 GeV/c for charged particles having ∆φ < 120 • : (left) particle density; (centre) energy density; (right) ratio of the energy and particle densities. The predictions of PYTHIA6 Z2, POWHEG Z2, PYTHIA8 4C, and HERWIG++ LHC-UE7-2 (with and without MPIs) are also displayed. In the top right plot, the structure around 60-80 GeV/c 2 for HERWIG++ without MPIs reflects the influence of photon radiation by final-state muons, which is enhanced below the Z resonance. Bottom: Ratios of the predictions of various MC models and the measurement. The inner band shows the statistical uncertainity of data whereas the outer band represents the total uncertainty. Fig. 2 at small p µµ T , where the radiation contribution is small and the activity in the transverse region is the same as that in the towards region. Figure 3 presents the ratios of the predictions of various MC models to the measurements for the observables shown in Fig. 2. Statistical fluctuations in the data induce correlated fluctuations for the various MC/data ratios. MADGRAPH in conjunction with PYTHIA6 tune Z2 describes the p µµ T dependence of the UE activity very well, both qualitatively and quantitatively. PYTHIA8 4C and HERWIG++ describe the p µµ T dependence of the particle density within 10-15%, but fail to describe the energy density. PYTHIA8 4C and HERWIG++ agree better with data as p µµ T approaches zero. The combination of the Z2 tune with POWHEG fails to describe the energy density in the towards and transverse regions, but gives a reasonable description of the particle density. This observation, combined with the information in Fig. 1, indicates that the discrepancies are not necessarily due to a flaw in the UE tune, but to an inadequate description of the multiple hard emissions and the different sets of PDFs used with POWHEG. At small p µµ T the comparisons with PYTHIA6 Z2 and POWHEG Z2 are similar to those in Ref. [38], where PYTHIA6 gives a good description of the p µµ T spectrum while POWHEG underestimates the p µµ T . Figure 4 shows the distributions of charged particle multiplicity (top row) and transverse momentum (bottom row). Figure 4 (left) shows a comparison of the normalized distributions in the away, transverse, and towards regions for events satisfying 81 < M µµ < 101 GeV/c 2 . As expected, the transverse and towards regions have fewer charged particles with a softer p T spectrum than the away region. Figure 4 (centre) shows the comparison of the normalized distributions in the transverse region for two different subsets of the selected events, one with 81 < M µµ < 101 GeV/c 2 and one with p µµ T < 5 GeV/c. The charged particle multiplicity is decreased and the p T spectrum is softer when p µµ T < 5 GeV/c is required, because of the reduced contribution of ISR. Figure 4 (right) shows the comparison of the normalized distributions with the predictions of various simulations in the transverse region for events satisfying 81 < M µµ < 101 GeV/c 2 . The charge multiplicity distribution is described well, within 10-15%, by MADGRAPH Z2 and PYTHIA8 4C. The p T spectrum is described within 10-15% by MADGRAPH Z2, whereas PYTHIA8 4C, POWHEG Z2, and HERWIG++ LHC-UE7-2 have softer p T spectra. The various MC programs achieve a similar level of agreement with data in the towards region as in the transverse region.

Comparison with the UE activity in hadronic events
The UE activity was previously measured as a function of leading jet p T in hadronic events for charged particles with pseudorapidity |η| < 2 and with transverse momentum p T > 0.5 GeV/c [6]. Figure 5 shows the comparison of the UE activity measured in the hadronic and the DY events (around the Z peak) in the transverse region as a function of p leading jet T and p µµ T , respectively. For the hadronic events two components are visible: a fast rise for p leading jet T 10 GeV/c due to an increase in the MPI activity, followed by an almost constant particle density and a slow increase in the energy density with p leading jet T . The increase in the particle and energy densities for p leading jet T 10 GeV/c is mainly due to the increase of ISR and FSR. Owing to the presence of a hard energy scale (81 < M µµ < 101 GeV/c 2 ), densities in the DY events do not show a sharply rising part, but only a slow growth with p µµ T due to the ISR contribution.
For p µµ T and p leading jet T > 10 GeV/c, DY events have a smaller particle density with a harder p T spectrum compared to the hadronic events, as can be seen in Fig. 5. This distinction is due to the different nature of radiation in the hadronic and DY events. Drell-Yan events have only initial-particle density energy density ratio of energy and particle densities  T for events satisfying 81 < M µµ < 101 GeV/c 2 : (left) particle density; (centre) energy density; (right) the ratio of the energy density and the particle densities. Predictions of MADGRAPH Z2, POWHEG Z2, PYTHIA8 4C, and HERWIG++ LHC-UE7-2 (with and without MPIs) are superimposed.  Figure 3: Ratios, as functions of p µµ T , of the predictions of various MC models to the measurements in the towards (upper row), transverse (centre row), and away (bottom row) regions for events satisfying 81 < M µµ < 101 GeV/c 2 : (left) particle density; (centre) energy density; (right) the ratio of the energy density and particle densities. The inner band shows the statistical uncertainty on the data whereas the outer band represents the total uncertainty.  state QCD radiation initiated by quarks, which fragment into a smaller number of hadrons carrying a larger fraction of the parent parton energy, whereas the hadronic events have both initial-and final-state QCD radiation predominantly initiated by gluons with a softer fragmentation into hadrons. Similar behavior is observed for the track-jet measurement where the UE activity is higher by 10-20% for gluon-dominated processes, as estimated from simulation. particle density energy density ratio of energy and particle densities  Figure 5: Comparison of the UE activity measured in hadronic and Drell-Yan events (around the Z resonance peak) as a function of p leading jet T and p µµ T , respectively: (left) particle density, (centre) energy density, and (right) ratio of energy and particle densities in the transverse region.

Summary
We have used Drell-Yan events to measure the UE activity in proton-proton collisions at √ s = 7 TeV, which were recorded with the CMS detector at the LHC. The DY process provides a UE measurement where a clean separation of the hard interaction from the soft component is possible. After excluding the muons from the DY process, the towards (|∆φ| < 60 • ) and the transverse (60 • < |∆φ| < 120 • ) regions are both sensitive to initial-state radiation and multiple parton interactions. The DY process provides an effective way to study the dependence of the UE activity on the hard interaction scale, which is related to the invariant mass of the dimuon pair. The influence of the ISR is probed by the dependence on the transverse momentum of the muon pair.
The UE activity is observed to be independent of the dimuon mass above 40 GeV/c 2 , after limiting the recoil activity, which confirms the MPI saturation at this scale. The UE activity in the DY events with no hard ISR is well described by PYTHIA6 and MADGRAPH with the Z2 tune and the CTEQ6L PDF. The Z2 tune does not agree with the data if used with PDFs other than CTEQ6L, as in the case of the POWHEG simulation. The PYTHIA8 4C and HERWIG++ LHC-UE7-2 tunes provide good descriptions of the energy-scale dependence of the UE activity. Thus the dependence of the UE activity on the energy scale is well described by tunes derived from hadronic events, illustrating the universality of MPIs in different processes. This universality is also indicated by the similarity between the UE activity in DY and hadronic events, although these events have different types of radiation. In addition, there is some ambiguity in the definition of the hard scale for both types of events.
The UE activity in the towards and transverse regions shows a slow growth with the transverse momentum of the muon pair and provides an important probe of the ISR. The leading-order matrix element generator MADGRAPH provides a good description of the UE dependence on dimuon transverse momentum. However, PYTHIA, POWHEG, and HERWIG++, which do not simulate the multiple hard emissions with sufficient accuracy, underestimate the energy density, but describe the particle density reasonably well. These measurements provide important input for further tuning or improvements of the Monte Carlo models and also for the understanding of the dynamics of QCD.