GammaLib

// Set input parameters

GRegions      on_regions = load_region("user input");

GRegions      exclusion_regions = load_region("user input");

GObservations obs("obs.xml");

GModels       models("models.xml");

GEbounds      energy_bounds("ebds.fits"');

// Fill spectrum

runs = GCTAOnOffObservations(obs, energy_bounds);

for run in runs.observations:

    // run is a GCTAOnOffSpectrum

    run.set_on_region(on_region)

    run.set_exclusion_regions(exclusion_regions)

    run.compute_off_regions()

    run.fill_spectra() # fills n_on, n_off

run.save("onoff.xml");

// Compute background acceptance 

GCTAOnOffObservations runs("onoff.xml");

for run in runs.observations:

    run.compute_bg_acceptance() // fills bg_acceptance_on,off

run.save("onoff.xml");

// Compute response

GCTAOnOffObservations runs("onoff.xml");

runs.models(models);

for run in runs.observations:

    run.compute_response() # fills gamma_exposure, energy_resolution

run.save("onoff.xml");

// Come back one month later and fit a model

GCTAOnOffObservations runs("onoff.xml");

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:

   GArf gamma_exposure

   GRmf energy_resolution

   GCTABackgroundAcceptanceSpectrum bg_acceptance_on

   GCTABackgroundAcceptanceSpectrum bg_acceptance_off

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 GArf:

   GEbounds energy_bounds

   std::vector<double> effective_area

   double obs_time

   to_arf()

   from_arf()

   fill_from(GResponse response, GModel model)

class GRmf:

   GMatrix energy_resolution_matrix

   to_rmf()

   from_rmf()

class GCTABackgroundAcceptanceSpectrum:

   GEbounds energy_bounds

   std::vector<double> effective_area

   double obs_time

   to_arf()

   from_arf()

   fill_from(GResponse response, GModel model)