Test full model fits/errors, refactor code

Run ipynb with all regions/mu values for SN 1006, checking free fits
based on best grid points (2, 3 bands).  Roughly compute error from
grids and perform fits for more accurate error bounds.