Double parton scattering is usually suppressed by single parton scattering. However, in some specific regions of phase-space, the differential cross-sections are comparable. Also, for a given final state, it might happen that the double parton scattering is the dominant contribution if the single parton scattering is suppressed by a higher multiplicity of couplings. For these reasons, it turns out to be necessary to include double parton scattering in event generators in order to give a better description of the data at high energy scales such as at the LHC.
Double parton scattering requires the use of double parton distributions and many efforts have been made during the last decade to produce realistic sets of those distributions, despite the lack of experimental data. In this work, we explore some directions to include those double parton distributions within the conventional partons showers.