The Metropolis algorithm for on-shell four-momentum phase space

We present several implementations of the Metropolis method, an adaptive Monte Carlo algorithm, which allow for the calculation of multi-dimensional integrals over arbitrary on-shell four-momentum phase space. The Metropolis technique reveals itself very suitable for the treatment of high energy processes in particle physics, particularly when the number of final state objects and of kinematic constraints on the latter gets larger. We compare the performances of the Metropolis algorithm with those of other programs widely used in numerical simulations.