Perform a Monte Carlo calculation of the goodness-of-fit for the C-statistic
XSPEC> goodness [<no. of realizations>]
This command simulates <no. of realizations> spectra based on the model and
writes out the percentage of these simulations with C less than that for the
data. If the observed spectrum was produced by the model then this number
should be around 50%.