Implements the Monte Carlo Wave Function method for quantum optics