We present a study on using Markov Chain Monte Carlo (MCMC) techniques to explore the high-dimensional and multi-modal phase space of scattering events at high-energy particle colliders. To this end, we combine the BAT.jl package that provides implementations of a variety of MCMC algorithms with the Sherpa event generator framework. We discuss technical aspects of the implementation and the resulting algorithm and present first results for the process of $Z+3$ jets production at the LHC.