I think that for this example, it would be really cool to use real data and fit it to model the sky background (with plots to demonstrate). I can provide such an example as a PR once this is merged. SImilarly, we could fit a 1D gaussian above to real spectral data (I can also provide).