Abstract
Beauty and charm quarks are ideal probes of pertubative Quantum Chromodymanics in proton-proton collisions, owing to their large masses. In this paper the role of multi-parton interactions in the production of doubly-heavy hadrons is studied using simulation samples generated with Pythia, a Monte Carlo event generator. Comparisons are made to the stand-alone generators BcVegPy and GenXicc. New methods of speeding up Pythia simulations for events containing heavy quarks are described, enabling the production of large samples with multiple heavy-quark pairs. We show that significantly higher production rates of doubly-heavy hadrons are predicted in models that allow heavy quarks originating from different parton-parton interactions (within the same hadron-hadron collision) to combine to form such hadrons. Quantitative predictions are sensitive to the modelling of colour reconnections. We suggest a set of experimental measurements capable of differentiating these additional contributions.