GammaLib

#! /usr/bin/env python

import gammalib

# Load OG only

obs = gammalib.GObservations('obs.xml')

print('=== Observation without response ===')

print(obs[0])

obs.save('write_obs.xml')

# Load OG with response group

obs = gammalib.GObservations('obs_rspgrp.xml')

print('')

print('=== Observation with response group ===')

print(obs[0])

obs.save('write_obs_rspgrp.xml')

# Save IRF FITS file (we need this later)

obs[0].response().save('irf.fits', True)

# Load OG with response file

obs = gammalib.GObservations('obs_rspfile.xml')

print('')

print('=== Observation with response file for continuum IRF ===')

print(obs[0])

obs.save('write_obs_rspfile.xml')

# Save IRF FITS file (we need this later)

obs[0].response().save('irf.fits', True)

# Load OG with response group for line response

obs = gammalib.GObservations('obs_rspgrp_line.xml')

print('')

print('=== Observation with response file for line IRF ===')

print(obs[0])

obs.save('write_obs_rspgrp_line.xml')