The production of pairs of hadrons in hadronic collisions is studied using a next-to-leading-order Monte Carlo program based on the phase space slicing technique. Up-to-date fragmentation functions based on fits to CERN LEP data are employed, together with several versions of current parton distribution functions. Good agreement is found with data for the dihadron mass distribution. A comparison is also made with data for the dihadron angular distribution. The scale dependence of the predictions and the dependence on the choices made for the fragmentation and parton distribution functions are also presented. The good agreement between theory and experiment is contrasted with the case for single ${\pi }^0$ production where significant deviations between theory and experiment have been observed.

• Received 2 October 2001

DOI:https://doi.org/10.1103/PhysRevD.65.034011