Dynamical approach to MPI four-jet production in Pythia

We improve the treatment of Multiple Parton Interactions (MPI) in \textsc{Pythia} by including the \12 mechanism and treating the \22 mechanism in a model-independent way. The \22 mechanism is calculated within the mean field approximation, and its parameters are expressed through Generalized Parton Distributions extracted from HERA data. The parameters related to the transverse parton distribution inside the proton are thus independent of the performed fit. The \12 mechanism is included along the lines of the recently developed perturbative QCD formalism. A unified description of MPI at moderate and hard transverse momenta is obtained within a consistent framework, in good agreement with experimental data measured at 7 TeV. Predictions are shown for the considered observables at 14 TeV. The corresponding code implementing the new MPI approach is available.


I. INTRODUCTION
It is widely realized now that 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-4] until the last decade [5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23], extensive theoretical studies have been carried out. Attempts have been made to incorporate multi-parton collisions in Monte Carlo (MC) event generators [24][25][26][27][28]. 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 have been performed at the Tevatron [29][30][31] and at the LHC [32][33][34][35], showing evidence for MPI at both soft and hard scales. This 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 the effective area Recently, a new approach based in perturbative Quantum Chromodynamics (pQCD) has been developed [17][18][19][20] for describing the MPI and its main ingredients are: • 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 socalled 2 ⊗ 2 mechanism (see Fig. 1 left), an additional 1 ⊗ 2 mechanism ( Fig. 1 right) is included. In this mechanism, which can be described in pQCD, the parton from one of the nucleons splits at some hard scale and creates two hard partons that may 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.
The use of this new formalism at LHC experiments needs its implementation in MC event generators which has not been performed yet. The purpose of the present paper is to make a step ahead towards the implementation of this formalism into MC generators. We use the standard simulation of the MPI implemented in Pythia [26], but with values of σ ef f calculated by using the QCD-based approach of [17][18][19][20], i.e. including 1 ⊗ 2 processes. The current approach used for the description of the MPI in Pythia is based on [25,26].
The Pythia code uses parton distribution functions, dependent on the impact parameter of the collision. From the theoretical point of view these are just one-particle Generalized Parton Distributions GPD 1 (see e.g. [37,38] 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].
Such an approach has, however, a number of difficulties, both conceptual and practical. First of all, a problem arises at the level of mean field approximation. The transverse parton distributions have been extracted from J/Ψ photoproduction measurements at the HERA collider, using QCD factorisation theorems [15,16,37,38]. Hence they can not be treated as free parameters of the model. Secondly, it has been observed that different Pythia parameters are obtained when data sensitive to a different region of the MPI spectrum are used for the fits. For example, it has been shown [39] that different parameters result for fits to UE or hard MPI data. This might be an indication that an additional transverse scale dependence, which is not present in the mean field approach, is needed to describe experimental data on UE and hard MPI simultaneously. Recent improvements in the Pythia MPI model include a dependence of the parton transverse density on the longitudinal momentum fraction (x) [25], but this only accounts for the x values of the hardest dijet. A complete x dependence which considers soft and hard partons may be irrelevant for the UE description where the transverse scales are rather close, but may become important for measurements sensitive to hard MPI.
The approach used in this 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 in Pythia. With these settings, the value of σ (0) ef f would be constant and independent on the scale. In order to implement the x and the scale dependence of σ ef f in collisions where a hard MPI occur, these events are rescaled according to: where σ (0) ef f is the effective cross section in the mean field approach calculated in a model independent way from GPD 1 , parameterized from HERA data [15][16][17], and R corresponds to the correction due to 1 ⊗ 2 mechanism [19,20]. Such an approach is equivalent to using the GPD 1 -based transverse parton densities for parton transverse distributions.
The main result of this paper is that the approach discussed above gives a unified description of both hard MPI and UE experimental data, with good accuracy and few fit parameters. The fit parameters are related to the amount of simulated MPI and of the color string reconnection, and to the separation scale between soft-and hard-scale processes, Q 2 0 whose value is expected to lie in the range 0.5-2 GeV 2 . The transverse scale dependent function R is calculated numerically by solving the nonlinear evolution equation [19,20]. Predictions using this approach are shown later in the paper, and are labeled as "UE Tune Dynamic σ ef f ". Our analysis shows that the values of observable for UE are quite close to the results obtained in a free parton model (mean field approximation), while the inclusion of transverse scale dependent rescaling calculated in pQCD [19] improves the description of hard MPI.
The paper is organized as follows. In section 2, basic theoretical ideas of the used approach are presented, while in section 3 their MC implementation is discussed. In section 4, comparisons for various predictions to observables measured at 7 TeV are shown. In section 5, predictions for these observables are presented for pp collisions at 14 TeV, before drawing the conclusions in Section 6.

II.
A SUMMARY OF THE THEORETICAL BACKGROUND.
The MPI four-jet cross section is characterized by the cross section σ ef f , which corresponds to an effective interaction area [17], and can be written as: where partons 1 and 2 create the first (12), and partons 3 and 4 the second (34) dijet. The pQCD calculation leads to the following expression for σ ef f in terms of two-particle GPD: 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, while the first term corresponds to the conventional case of two partons evolving from a low scale, namely the 2 ⊗ 2 mechanism. This first term 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] G 2 corresponds to the part of double GPD 2 , when both partons are evolved from the initial nonperturbative scale, while [1] G 2 corresponds to the case when one parton evolves up to some hard scale, where it splits to two successive hard partons, each of them in turn participating to the hard dijet event. We refer the reader to [17,18] for the detailed definitions of [1] G 2 and [2] G 2 and their connection to light cone wave functions of the nucleon.
For the two-parton GPD 2 we have: where D(x, Q 2 ) is a conventional parton distribution function (PDF). The use of the mean field approximation results in: and For the two gluon form factor F 2g , we use the exponential parametrization [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 σ ef f 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 the fit of experimental LHC data, but from the fit of single parton GPD 1 . The maximum transversality kinematics i.e. 4Q 2 = x 1 x 2 s for each dijet, have been considered in our approach, being Q the dijet transverse scale, and x 1 , x 2 the Bjorken fractions of the jets.
The second and third terms in Eq. 3 are parameterized as: where 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]. According to the results of [20], the dependence of R on x i in the maximum transversality regime is very weak and can be neglected with high accuracy. The function R also depends on the physical parameter Q 2 0 which corresponds to the separation scale between soft and hard dynamics where the GPD 2 is assumed to factorize.

OBSERVABLES.
In this paper we carry out two types of simulations: one based on the new approach defined in the sections 1 and 2 and one which follows the standard Pythia approach, used for comparison.
Let us recall the standard Pythia approach which is referred as to "UE tune" hereafter. In this study we use the Pythia 8.185 Monte Carlo event generator [24]. It simulates a 2→2 matrix element interfaced to parton shower and Underlying Event (UE). The Pythia 8 event generator uses a simulation of the parton shower ordered in transverse momentum and the Lund string model [42] to implement the hadronization process. The performed study has considered as a starting point the UE simulation implemented in the Pythia 8 tune 4C [26]. This simulation makes use of the CTEQ6L1 [45] parton distribution function and of a simple gaussian as a transverse matter distribution function. A fit to experimental data sensitive to the UE is performed in order to optimize the parameters related to the amount of MPI and colour reconnection in the simulation.
The fit operation has been carried out by using the RIVET [40] software, combined with the PROFESSOR machinery [41]. For the tune, two different observables have been considered at a center-of-mass energy of 7 TeV measured by the ATLAS experiment [36]. They are related to the multiplicity, N chg , and the sum of the transverse momentum, p rmT , of the charged particles in the region transverse to the direction of the leading charged particle in each event. The performed fit has used only the data points corresponding to transverse momenta of the leading charged particle between 2.0 and 15.0 GeV. The exclusion of the very low p T region (≤ 2 GeV) is motivated by the fact that processes at those scales are expected to be dominated by soft physics, including diffractive processes and soft nonperturbative correlations, i.e. along the lines of [18]. The upper cut off is arbitrary, since its variation starting from 5 GeV does not change the values of the observables.
The result of the fit consists of a new set of UE parameters implemented in the "UE tune" hereafter. The values of the Pythia 8 parameters obtained for the "UE tune" after the fit are shown in Table I. The first parameter listed in the table 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 color reconnection among parton strings. The value of σ ef f is found to be around 29.7 mb at 7 TeV; this value is significantly smaller than the one obtained by tuning the correlation observables of the four-jet scenario [35], which is around 19-21 mb. Note that the value of 29.7 mb is quite close to the one determined in mean field approach [17,20].
After fitting the UE observables for the "UE tune" determination, the considered predictions are also tested against measurements sensitive to the hard spectrum of the MPI. Measurements of such type have been conducted by studying correlations between outgoing objects in a proton-proton collision, for instance in four-jet final states measured at 7 TeV by CMS [35]. In this scenario, two dijets have been selected at different transverse momentum; two jets are required to have p T larger than 50 GeV and they are classified as "hard-jet pair", while the so-called "soft-jet pair" is composed by the two other jets selected with p T greater than 20 GeV. Two correlation observables, ∆S and ∆ rel sof t p T that are sensitive to DPS, have been considered. They are, respectively, the azimuthal angle between the two dijet planes and the p T balance between the soft jets and are defined as follows: where pair 1 (pair 2 ) is the hard (soft) jet pair and jet 1 (jet 2 ) is the leading (subleading) soft jet.
Let us now move to the new approach, based on the dynamical pQCD-based formalism, de- We considered both the case of moderate MPI (i.e. MPI at scales of several GeV), relevant for UE, and hard MPI.
For UE we treat separately the events where there is only one hard scattering, which are not rescaled, and the events with additional hard MPI. For the latter events two approaches were checked. First, we rescaled these events according to Equations 9 and 10, taking as Q 1 and Q 2 for the R function the scales of the two hardest scatterings. As shown in Section IV, the influence of this rescaling is very small (less than 5%), with respect to the standard Pythia "UE tune". This may be connected both with the small values of R obtained for UE, and with the fact that the ladder splitting is roughly taken into account for such scales by the large value of the parameter p 0 T ∼ 2 GeV. Subsequently, the second approach tried was to rescale only MPI events starting from the scale of order 4-5 GeV. When we rescale only the MPI starting from this (or a higher) scale, UE observables are not affected at all. At the same time, with this approach we avoid possible double counting effects, since at these scales the regularization formula in Pythia represents an ansatz for higher twist effects, including MPI. Thus, while using Pythia, we can neglect rescaling of MPI in UE, With the current accuracy, any of these two approaches can be used, leading to identical numerical results. This is in agreement with the approach documented in [20], where it was argued that at scales relative to UE the values of σ ef f are close to the ones calculated in the mean field approximation.
We consider now the case of hard MPI, specifically DPS. Two different processes may produce four jets in the final states. The first one is the so-called Single Parton Scattering (SPS) where the four jets are emitted through the same chain while the second one is DPS where the two hard interactions produce one dijet each. A different event topology is expected from these processes: if the four jets are produced through SPS, a high correlation between the objects of the final state is present and this is reflected in their relative configuration in the transverse plane. The direction of the hard jets, for example, is randomized by the emission of the additional two jets within the same chain and their initial p T balance is ruined. Instead, jet pairs coming from DPS events, namely from two independent scatterings, tend to be uncorrelated and their initial back-to-back configuration is less subject to smearing effects coming from additional hard radiation: the jet pairs are expected to exhibit a more balanced configuration in p T and azimuthal angle. In particular, as shown in [35], DPS events add a relevant contribution at low values of ∆S and ∆ rel sof t p T . Here we consider the experimentally relevant example, when the two dijet scales are 50 and 20 GeV.
Similarly as before, the x and scale dependence of σ ef f have been implemented by reweighting on an event-by-event basis the Monte Carlo simulation, as explained above.
In case that only MPI with p T scales smaller than 15 GeV are present in the collision, no x and scale dependence is applied to the σ ef f value of the corresponding event. The choice of 15 GeV is motivated by the fact that we need to treat differently the two contributing processes, SPS and DPS. Events where the two dijets are produced through SPS accompanied by moderate MPI, should not be reweighted [9,14,[17][18][19]. In case a hard MPI occurs in the collision, dynamical σ ef f values are used. In this way, we assume that all collisions with a MPI scale greater than 15 GeV produce the second hard dijet p T > 20 GeV pair selected in the considered four-jet scenario, while MPI at lower scales are below threshold for producing jets with p T > 20 GeV. This approach is generally followed by standard experimental measurements for σ ef f determination, as the ones documented in [32,33]. For our studies, lowering the 15 GeV cut off by 5-10 GeV shows variations of the predictions of DPS-sensitive observables of less than 2%. This is a clear indication of the consistency of our approach.
Various simulation settings have been considered for comparison: • "UE tune": predictions obtained with the parameters listed in Table I and without applying any reweighting of the simulation; this tune uses a constant value of σ ef f , following the standard Pythia approach; • "UE tune Q 2 -dep": predictions obtained with the UE parameters listed in Table I and by applying the scale dependence of σ ef f with Q 2 0 = 1 GeV 2 ; • "UE tune x-dep": predictions obtained with the UE parameters listed in Table I and by applying the x dependence of σ ef f ; • "UE tune Dynamic σ ef f ": predictions obtained with the UE parameters listed in Table I and by applying both the x and the scale dependence with Q 2 0 = 1 GeV 2 .
For the considered "UE tune Dynamic σ ef f ", predictions using Q 2 0 values equal to 0.5, 1 and 2 GeV 2 have been also tested and compared.
A full MC implementation of the presented approach may be different from the one used in this paper, which relies on reweighted events simulated by Pythia. There are at least three reasons for it: • by using the Pythia event generator, all ladders are assumed to evolve independently from the low transverse scale, the initial-state radiation (ISR) being regularized by primordial gluon distribution with a transverse scale equal to p 0 T . No parton ladder splittings are included in this approach; • in Pythia, the geometric picture of the collisions in the impact parameter space corresponds to the 2⊗2 mechanism, while for 1⊗2 mechanism the geometrical picture would be different; • for multi MPI events, namely for events with several MPI within the same collision, we neglect the change of relative weight of 1 ⊗ 2 and 2 ⊗ 2 mechanisms.
In this paper, these effects are neglected. First, the good agreement with experimental data shows that the high regularization scale p 0 T may be a good alternative parametrization of the ladder splittings and of the corresponding changes in ISR at the UE transverse scales. In other words, for UE the high p 0 T , which regularizes the charged particle multiplicity, also approximately fits the change of multiplicity due to ladder splitting. The Pythia regularization formula in this case can be viewed as an ansatz for twist expansion, that may include part of the MPI. Note that the ladder splitting scale is much smaller than the scales of hard dijets created by partons that evolve after the splitting [20,21]. So the effective ladder splitting is partly taken into account for UE by a high p 0 T value. This is the reason why the UE observables change only slightly in the new approach. On the other hand for hard MPI, when the hard splitting scale is much larger than p 0 T , the inaccuracy in accounting for ISR at small p T can be safely neglected.
Second, the direct calculation along the lines of [7] shows that neglecting the change of geometrical picture and of the relative weight between mean field and 1 ⊗ 2 mechanisms, when more than two separate dijet events are present, does not lead to numerical changes.
We conclude that using events simulated with Pythia and reweighted with x-and scaledependent values of σ ef f is a good approximation. In this way, we investigate the influence of changes of σ ef f on MC observables sensitive to UE and DPS.

IV. RESULTS FOR 7 TEV.
In this Section, comparisons between UE-and DPS-sensitive measurements at 7 TeV and various predictions are shown. Figure 2 shows comparisons to ATLAS data [36] on charged particle multiplicity and p T sum in the transverse region as a function of the leading charged particle p T .
Note that these are the observables which have been used in the fitting procedure for the determination of the "UE tune". The measurement is well reproduced by all considered predictions with discrepancies of only up to 10% in the high-p T region (p T > 10 GeV). The intermediate p T region (2 < p T < 10 GeV) is very well reproduced, while all predictions underestimate the first bins at p T > 2 GeV. This effect might be due to a not optimal simulation of diffraction in Pythia 8.
However, no relevant differences are observed for the different σ ef f models. FIG. 2: Charged particle density (left) and p T sum density (right) as a function of the leading charged particle in the transverse regions, measured by the ATLAS experiment at 7 TeV [36]. The data are compared to various predictions: the UE tune with constant σ ef f value (red curve), the UE tune with σ ef f x dependence applied (blue curve), the UE tune with σ ef f scale dependence with Q 2 0 =1.0 GeV 2 applied (black curve) and the UE tune with both σ ef f x and scale dependence with Q 2 0 =1.0 GeV 2 applied (pink curve). The lower panel shows the ratio between the various prediction and the experimental points.
In Figure 3, predictions obtained with different values of the scale Q 2 0 values are shown. All predictions are able to reproduce the measurement at the same good level. From this study, one may conclude that the UE data are not sensitive to the different settings of dynamical dependence applied to σ ef f . Figure 4 shows predictions with the various σ ef f settings considered previously, compared to the normalized cross section distributions as a function of the correlation observables, ∆S and ∆ rel sof t p T , measured in four-jet scenarios [35]. For these variables, the considered models show relevant differences. The static σ ef f dependence ("UE tune") is not able to properly describe the distribution as a function of ∆S; in particular, the region at low values (∆S < 2.5), where a DPS contribution is expected, is underestimated by about 10-18%. By introducing the x dependence for σ ef f ("UE tune x-dep"), the agreement at low values of ∆S does not significantly improve.
When the scale dependence of σ ef f is introduced ("UE tune Q 2 -dep"), the description of the normalized cross section as a function of ∆S gets better with differences not larger than 10%. The   best agreement with the measurement is obtained for predictions where both the x and the Q 2 dependence ("UE tune Dynamic σ ef f ") is included. The normalized cross section as a function of ∆ rel sof t p T is very well reproduced by all considered predictions. However, it has been already observed in [35] that ∆ rel sof t p T is less sensitive to a DPS contribution than ∆S, which uses information from both jet pairs. In Figure 5, predictions obtained with three different values of Q 2 0 (0.5, 1.0 and 2.0 GeV 2 ) are compared to the normalized cross section distributions as a function of ∆S and ∆ rel sof t p T . A considerable level of agreement for the different settings is obtained. Predictions obtained with Q 2 0 = 0.5 GeV 2 are in good agreement with the ∆S measurement but overestimate the first bin of ∆ rel sof t p T . For Q 2 0 = 1 GeV 2 and 2 GeV 2 the agreement tends to improve for ∆ rel sof t p T but is worse for ∆S. However, the measurement of the four-jet correlation observables is not able to discriminate the best choice for the value of Q 2 0 . In order to isolate the DPS contribution from the background produced by 2→4 processes, a dedicated event simulation has been performed with Pythia 8. Events with two hard scatterings within the same pp collision are simulated: the first hard scattering is generated with an exchanged transverse momentum between the outgoing partons,p T , larger than 45 GeV while for the second  ∆ rel sof t p T (right) measured in a four-jet scenario by the CMS experiment at 7 TeV [35]. The data are compared to various predictions: the new UE tune (red curve), the new UE tune with the x dependence applied (blue curve), the new UE tune with only the scale dependence with Q 2 0 =1.0 GeV 2 applied (black curve) and the new UE tune with both x and scale dependence with Q 2 0 =1.0 GeV 2 applied (pink curve). The lower panel shows the ratio between the various prediction and the experimental points. ∆ rel sof t p T (right) measured in a four-jet scenario by the CMS experiment at 7 TeV [35]. The data are compared to various predictions obtained with the new UE tune where both x and scale dependence have been applied with Q 2 0 equal to 1.0 (red curve), 0.5 (blue curve) and 2.0 (black curve) GeV 2 . The lower panel shows the ratio between the various prediction and the experimental points. one,p T is required to be greater than 15 GeV. Figure 6 shows the absolute cross sections predicted by the different settings implemented in the Pythia 8 simulation. The red curve shows the predictions for the UE tune with a static value of σ ef f while the blue, black and pink lines represent the predictions obtained when implementing the dynamical x and Q 2 dependence with y equal to 0.5, 1.0 and 2.0 GeV 2 , respectively. The highest DPS contribution is observed for Q 2 0 = 0.5 GeV 2 and it decreases for increasing Q 2 0 values. The lowest contribution is observed for the static UE tune when no x and Q 2 dependence is applied. The difference between the static and the dynamical σ ef f tune with Q 2 0 = 0.5 GeV 2 is around 80%. The different DPS contributions observed among the considered predictions reflect the decreasing σ ef f values for decreasing Q 2 0 as a function of the scale of the two scatterings (see Appendix of this paper). No significant differences in the shape of these distributions as a function of Q 2 0 are obtained. We observed that predictions of a dynamical σ ef f tune including a x and scale dependence of the transverse parton distribution are fully consistent with experimental data sensitive to moderate and hard MPI. The good agreement obtained for hard MPI is achieved due to contribution of 1 ⊗ 2 mechanism. The contribution from this mechanism is essentially model independent, except for Q 2 0 [20], which is the only new fit parameter, which is expected to lie in the 0.5 − 2 GeV 2 range.
The dynamical σ ef f dependence has been tested for predictions of UE and DPS observables at a center-of-mass energy of 14 TeV. The x and scale dependence of σ ef f follows respectively Equations 9 and 10, similarly as for 7 TeV. Note that the function R in Equation 10 also depends on the center-of-mass energy √ s [18]. Figure 7 shows predictions of charged particle density and the p T sum as a function of the leading charged particle p T , while in Figure  For each plot, the ratio to the predictions obtained with the UE tune with a constant σ ef f value is shown. While for the considered UE observables, a very small change is observed for the various predictions, larger differences are observed when the four-jet correlation observables are investigated. In particular, tunes with a dynamical σ ef f dependence tend to predict a higher  contribution at low ∆S and ∆ rel sof t p T values. These are the regions where a contribution from DPS is expected. The difference between static and dynamical σ ef f dependence is of up to 15% for ∆S < 2.0. Predictions with Q 2 0 equal to 1.0 and 2.0 GeV 2 are very similar to each other, while results obtained with Q 2 0 = 0.5 GeV 2 show a higher contribution at low values of ∆S and ∆ rel sof t p T , where the contribution of hard MPI is expected to be relevant.

VI. CONCLUSIONS
We have developed a new tune "UE tune Dynamic σ ef f " [47], that improves the treatment of hard Multiple Parton Interactions (MPI) in Pythia 8. The tune uses a fit to Underlying Event (UE) data in order to extract the parameters relative to soft MPI and includes values of σ ef f , which contain the 1 ⊗ 2 mechanism. They are calculated directly in the mean field + pQCD approach, as discussed in [20]. The dynamical dependence of σ ef f is not derived from a processdependent fit of the experimental data, but is directly obtained from theoretical calculations [17][18][19][20]. For the parameter Q 2 0 , that separates soft and hard scales, we have considered a range of values 0.5 < Q 2 0 < 2 GeV 2 . At present, the accuracy of the experimental data does not allow to carry a more precise determination, although the central values of the measured observables are better described by 0.5 < Q 2 0 < 1 GeV 2 . We observe that predictions from such tune are in good agreement with experimental measurements at 7 TeV, and for the first time give a consistent description of MPI at both moderate (UE) and hard scales. The results for UE are close to mean field approximation values, as anticipated in [19]. The additional transverse scale dependence of σ ef f , relative to mean field approach, due to 1 ⊗ 2 mechanism, is essential for a unified description of UE and hard MPI.
Predictions, obtained with the new tune for proton-proton collisions at 14 TeV, which are expected to happen within the next LHC phase, are also presented.  [2] N. Paver and D. Treleani, Z. Phys. C 28 187 (1985).