Use the Lorentz-invariant phase-space measure from the original PDF, including the and factors missing from the local TeX. In the center of mass frame put and , assuming . The spatial Dirac delta distribution sets , leaving
Since , the radial integral gives
The invariant flux factor is , using the Källén function. Dividing the phase space by this flux yields , with the final-state labels retained as in the printed formula.
Let . Of the Mandelstam variables, is
The azimuthal integral contributes , and the angular endpoints are therefore
Changing variables gives the requested expression from the formula supplied in the paper:
There is a normalization qualification: the displayed starting formula integrates over labeled final momenta and contains no . For the physical cross-section of two indistinguishable outgoing quanta of this real scalar field, the identical-particle factor in a final-state phase-space integral divides the full integral by two. With that convention ; equivalently, integrate only one representative of each exchanged pair. The boxed value follows the given formula exactly.