Dynamical approach to MPI four-jet production in Pythia

We modify the treatment of multiple parton interactions (MPI) in Pythia by including the $$1\otimes 2\;$$ mechanism and treating the $$2 \otimes 2\;$$ mechanism in a model-independent way. The $$2 \otimes 2\;$$ 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 $$1\otimes 2\;$$ mechanism is included along the lines of the recently developed formalism in perturbative QCD. 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 also shown for the considered observables at 14 TeV. The corresponding code implementing the new MPI approach is made available.