GammaLib

def hap:

    # HAP

    on_region = load_region('user input')

    exclusion_regions = load_region('user input')

    GObservations obs = load_obs('user input filename')

    models = GModels('models.xml')

    energy_bounds = GEbounds('user input filename')

    runs = GCTAOnOffObservations(obs, energy_bounds, models)

    runs.set_debug_level(3)

    for run in runs.spectra:

        # run is a GCTAOnOffSpectrum

        run.set_on_region(on_region)

        run.set_exclusion_regions(exclusion_regions)

        #maker = ReflectedRegionComputer()

        #maker.compute_off_regions(run)

        run.compute_off_regions()

        run.fill_spectra() # fills n_on, n_off

    # if you like make a separate loop / ctool

    for run in runs.observations:

        run.compute_bg_acceptance() # fills bg_acceptance_on,off

    for run in runs.observations:

        run.compute_response() # fills gamma_exposure, energy_resolution

    runs.save()

    # come back one month later and fit a different model

    # FitSpectrum

    runs.load()

    runs.fit_parametric_model(models)

    runs.compute_flux_points(models)

class GCTAOnOffObservations:

    """Container for spectra; one spectrum per run"""

    std::vector<GCTAOnOffObservation> observations

    def compute_total_observations():

        """Sum run-spectra into total spectrum.

        Store again in spectra member"""

class GCTAOnOffObservation:

    """Represents all data / response / analysis parameters

    for one run"""

    Regions on_regions

    Regions off_regions

    Regions on_exclusion_regions

    Regions off_exclusion_regions

    GCTAOnOffSpectrum spectrum

    GCTAOnOffResponse response

    compute_significance()

    compute_flux()

    load()

    save()

class GCTAOnOffResponse:

    GCTABackgroundAcceptanceSpectrum bg_acceptance_on

    GCTABackgroundAcceptanceSpectrum bg_acceptance_off

    ARF gamma_exposure

    GEnergyResolution energy_resolution

class GCTAOnOffSpectrum:

    GCountsSpectrum n_on

    GCountsSpectrum n_off

class GCountsSpectrum:

    GEbounds energy_bounds

    std::vector<double> counts

    to_pha()

    from_pha()

    fill_from(GEventList events)

class ARF:

    GEbounds energy_bounds

    std::vector<double> effective_area

    double obs_time

    to_arf()

    from_arf()

    fill_from(GResponse response, GModel model)

class GCTABackgroundAcceptanceSpectrum:

    GEbounds energy_bounds

    std::vector<double> effective_area

    double obs_time

    to_arf()

    from_arf()

    fill_from(GResponse response, GModel model)

class GEnergyResolution:

    GMatrix energy_resolution_matrix

    to_rmf()

    from_rmf()