To Do

  • Use the Python logging module to produce output and allow for a verbosity setting. [Easy]

  • Move from using obsid to a more-generic label for the multi-data-set case (i.e. when there are multiple data sets in a given annulus)

  • Allow the theta value to vary between data set in each annulus (for the multi-data-set case)

  • Support elliptical annuli

  • Create and use more generalized ModelStack and DataStack classes to allow for general mixing models. [Hard]

  • Add summary methods (similar to Sherpa’s show series of commands)

  • Add a specialized version of plot_source_component which takes a general model expression (i.e. without the _annulus suffix) and displays all combos per annulus. A bit tricky to get right. There’s also a plot_model_component.