Dynamical approach to MPI in W+dijet and Z+dijet production within the PYTHIA event generator

The new numerical approach that includes 1⊗2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\otimes 2$$\end{document} mechanisms is applied to double parton scattering (DPS) in W+dijet and Z+dijet final-state production in proton–proton collisions at the LHC. By using the underlying event (UE) simulation from a pythia 8.205 tune extracted in hadronic events, we show that, like in the case of a four-jet final state, the inclusion of 1⊗2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\otimes 2$$\end{document} mechanisms improves the description of experimental data measured at 7 TeV. The analysis is based on applying an event-by-event reweighting factor to a standard pythia 8.205 sample, by using the theoretical value of σeff\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _\mathrm{eff}$$\end{document}, which includes corrections due to 1⊗2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\otimes 2$$\end{document} mechanisms. In addition, predictions for proton–proton collisions at a center-of-mass energy of 13 TeV are shown for DPS-sensitive observables. The relevant code, used for this analysis, is publicly available at the following link: http://desy.de/~gunnep/SigmaEffectiveDependence/.


I. INTRODUCTION
Hard Multiple Parton Interactions (MPI) play an important role in the description of inelastic proton-proton (pp) collisions at high center-of-mass energies.Starting from the eighties [1][2][3][4] until the last decade , extensive theoretical studies have been performed.Significant progress was made on the simulation of multi-parton collisions in Monte Carlo (MC) event generators [23][24][25][28][29][30][31][32].Multiple parton interactions can serve as a probe for non-perturbative correlations between partons in the nucleon wave function and are crucial for determining the structure of the Underlying Event (UE) at Large Hadron Collider (LHC) energies.Moreover, they constitute an important background for new physics searches at the LHC.A large number of experimental measurements has been released at the Tevatron [33][34][35] and at the LHC [36][37][38][39][40][41][42], showing a clear evidence of MPI at both soft and hard scales.The latter case is usually referred to as "Double Parton Scattering" (DPS), which involves two hard scatterings within the same hadronic collision.
The cross section of such an event is generally expressed in terms of σ eff 26].In the socalled "mean-field approximation", the cross section σ eff is the effective area which measures the transverse distribution of partons inside the colliding hadrons and their overlap in a collision.
Recently, a new approach based in perturbative Quantum Chromodynamics (pQCD) has been developed [17][18][19][20] for describing MPI.Its main ingredients are listed in the following: • the MPI cross sections are expressed through new objects, namely double Generalized Parton distributions (GPD 2 ); • besides the conventional mean-field parton model approach to MPI, represented by the so- included.In this mechanism, which can be calculated in pQCD, a parton from one of the nucleons splits at some hard scale and creates two hard partons that participate in MPI.
This mechanism leads to a significant transverse-scale dependence of MPI cross sections; • the contribution of the 2⊗2 mechanism to GPD 2 is calculated in a mean-field approximation with model-independent parameters.By using the pythia 8 event generator, this new formalism has been implemented in the simulation of MPI in pp collisions at LHC energies for four-jet final states [43].The new approach has been developed by including a dynamic σ eff calculation and has two basic differences from the conventional pythia MPI simulation [28].On the one hand, the MPI cross section is calculated in mean-field approach by using the factorization property of GPD 2 , which connects them to the more conventional one-parton GPD, which have been measured at HERA and parameterized in [15][16][17].On the other hand, corrections due to 1 ⊗ 2 mechanisms are implemented by solving nonlinear pQCD evolution equations [19,20].The value of σ eff is obtained according to: where σ (0) eff is the effective cross section in the mean-field approach calculated in a modelindependent way, and R(Q 2  1 , Q 2 2 , Q 2 0 ) is calculated by solving iteratively the nonlinear evolution equation, as explained in detail in [19,20].The two scales Q 1 , Q 2 are the transverse scales of the two hardest dijet systems produced within the same pp collision, while Q 2 0 = 0.5 − 1 GeV 2 is the scale which divides soft and hard processes in the MPI formalism [20].In such an approach, results obtained for MPI cross sections are model independent, and do not need additional fit parameters for characterizing the MPI.It has been shown [43] that the new formalism implemented in pythia 8 -called in [43] and hereafter "UE Tune Dynamic σ eff " -achieves a consistent description of both soft and hard MPI in four-jet final states.
The aim of this paper is to extend the MPI approach described in [43] for W+dijet and Z+dijet final states.Various predictions are compared to experimental data on W+dijet (hereafter referred to as Wjj) [36,37]  Unlike a four-jet final state, it has been shown [36,37] that higher-order contributions from single parton scattering (SPS) are crucial for a consistent description of DPS-sensitive observables in the Wjj channel.These are not included in the matrix element simulated by the pythia event generator.Their neglection will lead to incorrect reults, and strong disagreement with experimental data.Other event generators, like madgraph [44] and powheg [45], which simulate multileg and higher-order matrix elements, are more suitable for studies of the Wjj and Zjj final states.The madgraph event generator includes both NLO and NNLO real corrections to inclusive W and Z processes but no virtual ones, while powheg generates a full NLO calculation of Wjj and Zjj cross sections with both real and virtual corrections.The matrix elements generated by madgraph and powheg at the parton level are then interfaced to the UE simulation provided by pythia 8.In particular, the applied UE simulation uses the parameters of the so-called "UE Tune" [43], extracted from UE data in hadronic events at transverse scales between 2 and 5 GeV.Even though the hard scales are much higher than 5 GeV in inclusive W and Z boson events, our model-independent calculation of σ eff in Wjj and Zjj final states shows values which are very similar to σ eff obtained in the four-jet channel.This looks however a pure numerical coincidence but motivates the choice of using the "UE Tune" also for W and Z boson events.However, it is important to note that the value of the rescaling function R in Eq. ( 1) calculated for Wjj and Zjj is very different from the corresponding rescaling function R for four-jet final states.It is only the ratio between σ eff and R that is numerically similar between Wjj, Zjj, and four-jet processes.
The paper is organized in the following way.In chapter II, the basic formalism and the numerical calculation of σ eff is described, along with its MC implementation.In chapter III, predictions of observables on Wjj and Zjj final states are compared for different settings of the UE simulation, while chapter IV considers variables on inclusive W and inclusive Z processes.Chapter V shows predictions of variables in the same previous final states at 14 TeV, while summary and conclusions are given in chapter VI.Theoretical dependence of σ eff in Wjj and Zjj as a function of the dijet transverse scale is investigated in appendix A for 7 and 14 TeV center-of-mass energies, while in appendix B, results obtained with the madgraph MC event generator are considered.

A. Theoretical tools
The approach of the paper is based on the calculation of the MPI cross section by means of the effective cross section σ eff .In the case of Wjj or Zjj channels, the MPI cross section can be written as: where partons 1 and 2 create the gauge boson, and partons 3 and 4 the dijet system.The pQCD calculation leads to the following expression for σ eff in terms of two-particle GPD: The scale Q 1 of first hard collision is kept fixed to Wjj and Zjj production.The second and third terms in Eq. (3) correspond to the 1 ⊗ 2 mechanism, when two partons are generated from the splitting of a parton at a hard scale after evolution.The first term corresponds to the conventional case of two partons evolving from a low scale, namely the 2 ⊗ 2 mechanism and can be calculated in the mean field approximation [15][16][17].The momentum ∆ is conjugated to the relative distance between the two participating partons.The full double GPD is a sum of two terms: Here [2] GP D 2 corresponds to the part of the GPD 2 , referring to the occurrence when both partons are evolved from an initial nonperturbative scale.The [1] GP D 2 function corresponds to the case when one parton evolves up to some hard scale, at which it then splits into two successive hard partons, each of them participating in turn to the hard dijet event [55].
The difference with respect to [43] is that the partons 1 and 2 are quarks in the case of W and Z boson events (u and d for W and uū or d d for Z production), instead of gluons in the case of four-jet final states.Hence, the GPD in Wjj and Zjj is defined as where D(x, Q 2 ) is a conventional parton distribution function (PDF).The use of the mean-field approximation results into: and The numerical analysis of HERA data shows that the gluonic and quark radii of the nucleon are of similar size [46].Hence, here we neglect the difference between the two-gluon form factors and its quark analogues.By assuming no difference between the initial partons, the structure functions cancel out in the mean-field calculation of σ eff .Consequently, the 2 ⊗ 2 part of the calculation can be done using the two-gluon form factors only.For the two-gluon form factor F 2g , we use the exponential parametrization described in [16].In fact, it leads to the same numerical results as the dipole form [15], but it is more convenient for calculations.This parametrization is unambiguously fixed by J/Ψ diffractive charmonium photo/electro production at HERA.The functions D are the conventional nucleon structure functions and F 2g can be parameterized as: where B g (x)= B 0 + 2K Q •log(x 0 /x), with x 0 ∼ 0.0012, B 0 = 4.1 GeV −2 and K Q = 0.14 GeV −2 .In our implementation the central values of the parameters B 0 and K Q [16] have been used, which are known with an accuracy of ∼ 8%.Integrating over ∆ 2 , we obtain for the part of σ eff corresponding to the first term in Eq. ( 3): where x 1−4 are the longitudinal momentum fractions of the partons participating in the 2 ⊗ 2 mechanism.This cross section corresponds to the free parton model and is model independent in the sense that its parameters are determined not from a fit to experimental LHC data, but from a fit to single parton GPD.The maximum transversality kinematics for the dijet system, i.e.
4Q 2 = x 3 x 4 s, have been considered in our approach, being Q the dijet transverse scale, and x 3 , x 4 the Bjorken fractions of the two jets.Concerning the dependence of σ eff on the parton scales, the results documented in [20] show that rescaling factors R in Eq. ( 9) are different for different final states, e.g.four-jet, Wjj and Zjj.Hence, we calculated the rescaling factors separately for the two considered channels, Wjj and Zjj.More details on the obtained values of σ eff as a function of the dijet scale are described in Appendix B.

B. Monte Carlo implementation and definition of experimental observables
The standard simulation of MPI implemented in pythia 8 [30] is considered, but with values of σ eff calculated by using the QCD-based approach of [17][18][19][20], i.e. including 1 ⊗ 2 processes.
The simulation of the MPI in pythia is based on [29,30].The pythia code uses single parton distribution functions, dependent on the impact parameter between the two colliding partons.From a theoretical point of view, these are just GPD 1 (see e.g.[47,48] for a review).The parameters set in the pythia simulation relative to the transverse parton density are extracted from fits to experimental data on UE, sensitive to the contribution of the MPI.This procedure is closely related to mean-field-based schemes, see e.g.[17].
The approach developed in [43] and used in the present paper combines the standard pythia MPI model with the one of [17][18][19][20].We use a single gaussian to model the matter distribution function of the protons.With these settings, the value of σ eff would be constant and independent of the scale.In order to implement the σ eff dependence as a function of the parton momentum fraction x and of the scale, the events where a hard MPI occur, are rescaled according to Eq. ( 1).
Two types of simulations are considered: one based on the new approach defined in Section 1 and called "Dynamic σ eff ", and one which follows the standard pythia 8 approach without any rescaling, called "UE tune".
The UE tune [43] has been extracted from fits to UE data in hadronic final states for scales of the leading charged particle in the range 2-5 GeV 2 and its parameters are listed in Table I The first parameter, MultipartonInteractions:pT0Ref, refers to the value of transverse momentum, p 0 T , defined at √ s = 7 TeV, used for the regularization of the cross section in the infrared limit, according to the formula 1/p 4 T → 1/(p 2 T + p 0 2 T ) 2 .The second parameter is the probability of colour reconnection among parton strings.The value of σ eff is found to be around 29.7 mb at 7 TeV and is quite close to the one determined in the mean-field approach [17,20].
Two different sets of observables are studied in final states with W or Z bosons: the first one consists of variables sensitive to hard MPI, the second one includes observables which are mostly influenced by the UE, namely by MPI at moderate scales.Final states with a Z or a W boson are separately investigated at the stable-particle level by using the rivet framework [49].For the UE study, an inclusive Z or W boson production is required.The Z boson is reconstructed through its muonic decay: two muons with p T > 20 GeV in |η| < 2.4 are required with an invariant mass between 81 and 101 GeV 2 .For the W-boson selection, a final state with one muon with p T > 30 GeV and |η| < 2.1 and a missing transverse energy of 30 GeV is required.In the Wjj and Zjj channel, two jets clustered with the anti-k T algorithm with p T > 20 GeV in |η| < 2.0 are added to the selection of, respectively, the W and the Z boson.The various final states for which the considered predictions are tested are summarized in Table II.The following observables are investigated for the study of DPS in Wjj and Zjj final states: where boson may be the W or Z boson, jet 1,2 is the jet pair and jet 1 (jet 2 ) is the leading (subleading) jet.
The study of the UE contribution has been performed through the usual a-la-Rick-Field strategy [50].The direction of the reconstructed boson identifies the direction of the hard scattering and defines different regions in the plane transverse to the beam direction: the "toward" region (|∆φ| < 60 • ), two transverse regions (60 < |∆φ| < 120 • ) and the "away" region (|∆φ| > 120 • ).
Only observables measured in the transverse regions have been considered in this study, since they are the ones which are most affected by the UE contribution.The observables refer to the amount of number of charged particles and of their transverse momentum and are: • charged particle multiplicity density (N ch ); • transverse momentum sum density (Σp T ).
Charged particles which contribute to these quantities are selected in each event within a region of pseudorapidity |η| < 2.5 with a lower p T cut of 500 MeV.
The x and scale dependence of σ eff has been implemented in the considered predictions by reweighting on an event-by-event basis the MC simulation in presence of a hard MPI (p T > 15 GeV).The x dependence is given by Eq. ( 9), where x 1,2 are taken as the longitudinal momentum fractions of the partons participating in the hardest scattering (W-or Z-boson production), while x 3,4 refer to the longitudinal momentum fractions of the partons participating in the hardest MPI.
The scale dependence is expressed by Eq. ( 9), where R takes for Q 2 the scale of the hardest MPI.
Different values of Q 2 0 have been considered in the range between 0.5 and 1 GeV 2 .Predictions with the following simulation settings are considered for comparison: • "UE Tune" [43]: predictions obtained without applying any reweighting of the simulation; this tune uses a constant value of σ eff , following the standard Pythia approach, and its parameters have been extracted by fits to the UE measurement in hadronic final states; • "UE Tune x-dep": predictions obtained with the parameters of the UE tune and by applying the x dependence of σ eff ; • "UE Tune Dynamic σ eff ": predictions obtained with the parameters of the UE tune and by applying both x and scale dependence for σ eff values; two different tunes are shown, corresponding to values of Q 2 0 equal to 0.5 and 1 GeV 2 .As it was done in four-jet final states [43], the dependence of the cross section on the p T of the outgoing partons is assumed to be the same in 2 ⊗ 2 and 1 ⊗ 2 production mechanisms.This might affect differential distributions as a function of the jet balance, i.e. ∆ rel p T .This approach is implemented for various Monte Carlo event generators which use different matrix-element calculations: • pythia 8 [28], which implements a 2→2 LO ME (q q → Z and q q → W for, respectively, Z-and W-boson production), where additional hard partons in the final state are generated through the parton-shower simulation in a leading-log approximation; • powheg [45] interfaced to pythia 8, which implements a 2→4 NLO ME; • madgraph [44] interfaced to pythia 8, which implements a 2→4 LO ME, where up to four partons in addition to the Z-or the W-boson are simulated within the ME calculation.
The pythia 8 sample uses the CTEQ6L1 [53] PDF set, while the powheg sample has been generated with the CT10NLO [54] PDF set.For the madgraph sample, the CTEQ6L1 [53] PDF set has been used and the matching and merging scale between matrix element (ME) and parton shower (PS) have been set to, respectively, 10 and 20 GeV in the MLM formalism [51].
Predictions obtained with the considered event generators have been compared to the ∆S and ∆ rel p T observables, measured at 7 TeV by the CMS experiment in the Wjj channel [37] and they are shown in Fig. 2. Predictions obtained with powheg and madgraph interfaced to pythia 8 UE Tune are able to follow the shape of the measured points, while pythia 8 does not describe at the same level of the agreement.In particular, higher-order contributions fill the region of the phase space which is most sensitive to the DPS signal.This effect, already observed in [37], is a clear indication of the need of higher-order matrix elements to give a reasonable description of DPSsensitive observables.Predictions obtained with powheg and madgraph without the simulation of MPI, also shown in Figure 2, are not able to follow the trend of the measured ∆S and ∆ rel p T .
In particular, they underestimate the region of ∆S < 2 by about 50-70% and the region of ∆ rel p T < 0.15 by about 10-20%.These are the regions of the phase space where a signal from hard MPI is expected to contribute.The large discrepancy observed between data and predictions without the simulation of MPI clearly indicates the need of MPI contributions in the current models for a good description of DPS-sensitive observables.
In the following sections, results are shown by using the powheg event generator interfaced to pythia 8, which consistently includes both real and virtual NLO corrections for hard Wjj and Zjj processes.Comparisons with simulations obtained with madgraph, which was used for experimental extractions of σ eff by the CMS collaboration [37], are documented in Appendix A, while predictions with pythia 8 standalone are dropped from the discussion.

III. DPS-SENSITIVE OBSERVABLES IN W+DIJET AND Z+DIJET FINAL STATES
In this section, comparisons of various predictions for DPS-sensitive observables in Wjj and Zjj are shown at 7 TeV.In Fig. 3, comparisons to data measured by the CMS experiment [37][38][39] in the Wjj channel at 7 TeV are considered.They refer to the normalized distributions of the correlation observables ∆S (left) and ∆ rel p T (right).Predictions of powheg interfaced to pythia 8 UE Tune are considered with different σ eff dependence applied: no reweighting, only x-dependent σ eff values calculated in mean-field approach, x-and scale-dependent σ eff values with The inclusion of contributions from 1⊗2 production mechanisms in the predictions with dynamic σ eff values improves the agreement with the measurement for the ∆S observable.They follow the decreasing shape of the observable better than the predictions obtained with the UE Tune without any rescaling and the UE Tune with only x-dependence applied.The ∆ rel p T observable is in good agreement with every prediction, except at values ∆ rel p T < 0.15, where the curves are above the data by about 10-25%.This might be due to the fact that in our approach we assume the ∆ rel p T dependence of the cross section to be the same in 2 ⊗ 2 and 1 ⊗ 2 production mechanisms.This Tune with different σ eff dependence applied: no reweighting applied (red line), x-dependent σ eff values (blue line), x-and scale-dependent σ eff values with Q 2 0 = 0.5 GeV 2 (green line) and x-and scale-dependent σ eff values with Q 2 0 = 1 GeV 2 (pink line).Also shown are the ratios of these tunes to the data.
does not need to be true [17][18][19][20].It would be interesting to study if indeed different assumptions on the differential 1 ⊗ 2 and 2 ⊗ 2 cross sections bring to significant difference and play a role at low values of ∆ rel p T .However, this demands additional analytical and numerical work, and will be done elsewhere [52].Also, it would be interesting to study the agreement of the considered predictions for a different scale of the dijet system, different from 20 GeV as in this case.

IV. UE OBSERVABLES IN INCLUSIVE W AND Z BOSON EVENTS
Predictions obtained with the considered tunes are also tested for UE observables in inclusive W and Z boson events.This kind of events are sensitive to MPI at moderate scales.In Fig. 5, predictions on charged-particle multiplicity and p T sum densities are shown for inclusive W events in the transverse region as a function of the W-boson p T .Curves obtained with powheg interfaced to pythia 8 UE Tune and implementing different σ eff dependence differ less than 2% from each other.This effect is very similar to the one observed in hadronic events, documented in [43].
In Fig. 6, various predictions obtained with powheg interfaced to pythia 8 UE Tune are shown of the two UE observables in the transverse region as a function of the Z-boson p T and compared to the measurement performed by the CMS experiment [50].As seen in inclusive W events, the difference among the considered curves is of the order of 2%.All predictions are able to follow the data points reasonably well at all scales with differences up to 10%.change is observed for variables sensitive to the contribution of moderate MPI and predictions with or without dynamic σ eff values are able to reproduce the data at the same good level of agreement.

V. PREDICTIONS OF DPS-SENSITIVE OBSERVABLES AT 14 TEV
In this Section, predictions of DPS-sensitive observables at 14 TeV are shown for Wjj and Zjj final states.Only the powheg event generator is considered with the same settings used for comparisons at 7 TeV.The energy extrapolation of the p 0 T value at 14 TeV is applied through the parameter of the pythia 8 tune 4C [30].In 7, predictions at 14 TeV are shown for the normalized distributions of the correlation observables ∆S and ∆ rel p T in the Wjj and Zjj channel, obtained with powheg interfaced to pythia 8 UE Tune with different σ eff dependence applied: no σ eff reweighting, x-dependent σ eff values, and x-and scale-dependent σ eff values with Q 2 0 = 0.5 GeV 2 and Q 2 0 = 1 GeV 2 .Very similar behaviour is observed for the two channels.Predictions obtained without any rescaling differ of about 10-15% from the curves which include the x and scale dependence of σ eff in the regions of phase space where a DPS signal is expected to contribute, namely ∆S < 2 and ∆ rel p T < 0.2.No relevant difference is observed in case a value of Q 2 0 equal to 0.5 or 1.0 GeV 2 is used.A higher DPS sensitivity might result for a different jet selection.
For instance, bigger differences by about 20-25% are observed between predictions obtained with powheg interfaced topythia 8 with and without event reweighting, in case the two jets are selected with a rapidity separation ∆η > 6, in association with a W or a Z boson.A requirement of a large |∆η| indeed suppresses the contribution of SPS processes and increases the sensitivity to DPS contributions.

VI. CONCLUSIONS
The new tune, developed in [43] for the analysis of multiple parton interactions (MPI) in four-jet final states and including contributions from 1 ⊗ 2 mechanisms, is compatible with measurements sensitive to double parton scattering (DPS) and MPI at moderate scales in W+dijet and Z+dijet channels.In order to properly treat events with a W or a Z boson with associated jets, it is necessary to include higher-order contributions within the matrix element calculation.In this paper, simulations using the powheg event generator interfaced to the underlying event (UE) simulation provided by pythia 8 are considered.Predictions using dynamic σ eff values dependent on the longitudinal momentum fractions and on the scale of the process improve the description of correlation observables measured in W+dijet final states.The experimental accuracy achieved so far does not allow to make conclusion on the best value of the scale separating soft and hard MPI, Values of Q 2 0 between 0.5 and 1 are able to describe the measurement at the same level of agreement.Differences by about 15% are observed between these predictions and measured data if the jet balance in transverse momentum p T (∆ rel p T ) is considered.This might be due to the assumption made in our simulation that differential distributions in 1⊗2 mechanisms have the same p T dependence as the conventional 2 ⊗ 2 diagram, as a simple generalization of the formula given in [18][19][20].Although explicit formulae for differential distributions of 1 ⊗ 2 production mechanisms are known [19,20], their actual implementation still demands additional work, both analytical and numerical.describe better the measurement, especially at values of ∆S < 2 and at ∆ rel p T < 0.2.A slight better agreement than predictions obtained with the powheg matrix element is observed for the normalized cross section as a function of ∆ rel p T .As for the powheg case, it is not possible to discriminate the best value of Q 2 0 , due to the large experimental uncertainty.In Fig. 8(bottom), predictions with the same settings are tested on ∆S and ∆ rel p T in the Zjj channel.As the results for powheg, differences between predictions with dynamic σ eff values are 10-15% above the curve using no σ eff rescaling at ∆S < 2 and at ∆ rel p T < 0.2.
In Fig. 9, the observables sensitive to the UE are investigated by comparing different madgraph predictions.In Fig. 9(top), predictions on the charged particle multiplicity and p T sum densities as a function of the W-boson p T are shown in the inclusive W channel, while in Fig. 9(bottom) the same predictions are compared to the CMS measurement in the inclusive Z channel.Very similar conclusions as drawn when considering powheg can be extracted.Differences among the various predictions are observed only of the order of 2% and, in the case of inclusive Z production, all of them are able to reproduce the trend of the measured points.
In conclusion, DPS-sensitive observables are better described by predictions using dynamic σ eff values, while variables sensitive to MPI at moderate scales are not strongly affected by σ eff variation within our approach.This is also the case for powheg.Considering the fact that same conclusions hold for madgraph and powheg, real NLO corrections are dominant in Wjj and Zjj final states, while virtual ones have a low impact.

FIG. 2 :
FIG.2: CMS data[37] at 7 TeV for the normalized distributions of the correlation observables ∆S (left) and ∆ rel p T (right) in the W+dijet channel, compared to predictions generated with pythia 8 UE Tune, madgraph and powheg interfaced to pythia 8 UE Tune.Predictions obtained with madgraph and powheg interfaced to pythia 8 UE Tune, without the simulation of the MPI are also compared to the measurement.The ratios of these predictions to the data are shown in the lower panels.

2 10 − 2 10 − 1 ∆ 4 ∆
FIG. 4: Predictions at 7 TeV for the normalized distributions of the correlation observables ∆S (left) and ∆ rel p T (right) in the Z+dijet channel, of simulations performed with powheg interfaced to pythia 8 UE Tune with different σ eff dependence applied: no σ eff reweighting applied (red line), x-dependent σ eff values (blue line), x-and scale-dependent σ eff values with Q 2 0 = 0.5 GeV 2 (green line) and x-and scale-dependent σ eff values with Q 2 0 = 1 GeV 2 (pink line).Also shown are the ratios of each curve to the predictions of the UE Tune.

2 N 2 NFIG. 5 : 2 N 2 NFIG. 6 :
FIG.5: Predictions for the (left) charged-particle and (right) p T sum densities in the transverse regions as defined by the W-boson in proton-proton collisions at 7 TeV.Simulations obtained with powheg interfaced to pythia 8 UE Tune are considered with different σ eff dependence applied: no reweighting applied (red line), x-dependent σ eff values (blue line), x-and scale-dependent σ eff values with Q 2 0 = 0.5 GeV 2 (green line) and x-and scale-dependent σ eff values with Q 2 0 = 1 GeV 2 (pink line).Also shown are the ratios of these tunes to the predictions of the UE Tune.

FIG. 8 : 2 N 2 N 2 N 2 NFIG. 9 :FIG. 10 :
FIG. 8: top CMS data at 7 TeV for the normalized distributions of the correlation observables ∆S (left) and ∆ rel p T (right) in the W+dijet channel, compared to predictions of madgraph interfaced to pythia 8 UE Tune with different σ eff dependence applied: reweighting applied (red line), x-dependent σ eff values (blue line), x-and scale-dependent σ eff values with Q 2 0 = 0.5 GeV 2 (green line) and x-and scale-dependent σ eff values with Q 2 0 = 1 GeV 2 (pink line).Also shown are the ratios of these predictions to the data.bottom Predictions at 7 TeV for the normalized distributions of the correlation observables ∆S (left) and ∆ rel p T (right) in the Z+dijet channel, of simulations performed with madgraph interfaced to pythia 8 UE Tune with different σ eff dependence applied: no σ eff reweighting applied (red line), x-dependent σ eff values (blue line), x-and scale-dependent σ eff values with Q 2 0 = 0.5 GeV 2 (green line) and x-and scale-dependent σ eff values with Q 2 0 = 1 GeV 2 (pink line).Also shown are the ratios of these tunes to the predictions of the UE Tune.

TABLE I :
. Pythia 8 parameters obtained after the fit to the UE observables.The value of pT0Ref is given at a reference energy of 7 TeV.Values of σ eff at 7 and 14 TeV are also shown in the table.