Wave function Monte Carlo method for polariton condensates