Created a load_scaled_data function to handle datasets with non-unity

AREASCAL values (i.e., XMM-RGS data).  ISIS got that "right" only for
unbinned data.  This reimplements the implementation so that normal
binning/statistics can be applied.