Improving on the PYTHIA modelling of equal-scale multi-parton distribution functions

A bstract Multi-parton distribution functions (mPDFs) are non-perturbative objects that are important in the prediction of multiple scattering rates at hadron colliders. In the case where the scales associated with all partons in the mPDF are the same, we have two theoretical constraints on the mPDF. These are symmetry in exchange of the parton indices, and the number and momentum sum rules. In a previous publication JHEP 02 (2023) 090 we found that the equal-scale mPDFs from the P ythia model could not satisfy both of these constraints simultaneously. In this paper we introduce an algorithm for constructing equal-scale mPDFs that is based on the P ythia procedure but has three additional modifications, such that it yields symmetric mPDFs that should satisfy the sum rules to an improved extent. We test the construction for the case of the double and triple parton distribution functions (dPDFs and tPDFs), finding that the sum rules are obeyed to within 10% over the vast majority of the phase space for the scales tested (and deviations only being mildly above this level). We use our dPDFs to compute rapidity asymmetries for same-sign WW and ZZ production via double parton scattering, and compare the results to predictions obtained using P ythia and the GS09 dPDFs of JHEP 03 (2010) 005.