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 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 to
the case for single π0 production where significant deviations between
theory and experiment have been observed.Comment: 22 pages, 15 figures; 3 references added, one figure modified for
clarit