def plot_source_spectra(name):
plot_source_spectrum(
source=SourceCatalog3FGL()[name], label="Fermi 3FGL", color="r"
)
plot_source_spectrum(
source=SourceCatalog2FHL()[name], label="Fermi 2FHL", color="g"
)
plot_source_spectrum(
source=SourceCatalog1FHL()[name], label="Fermi 1FHL", color="c"
)
plot_source_spectrum(
source=SourceCatalog3FHL()[name], label="Fermi 3FHL", color="b"
)
ax = plt.gca()
ax.set_ylim(1.0e-12, 7.0e-11)
ax.set_xlim(1.0e-4, 2.0)
ax.set_xlabel("Energy (TeV)")
ax.set_ylabel("E^2 dN/dE (erg cm-2 s-1])")
plt.legend(loc=0)
def plot_source_spectra(name):
plot_source_spectrum(source=SourceCatalog3FGL()[name],
label='Fermi 3FGL',
color='r')
plot_source_spectrum(source=SourceCatalog2FHL()[name],
label='Fermi 2FHL',
color='g')
plot_source_spectrum(source=SourceCatalog1FHL()[name],
label='Fermi 1FHL',
color='c')
plot_source_spectrum(source=SourceCatalog3FHL()[name],
label='Fermi 3FHL',
color='b')
ax = plt.gca()
ax.set_ylim(1.e-12, 7.e-11)
ax.set_xlim(1.e-4, 2.)
ax.set_xlabel('Energy (TeV)')
ax.set_ylabel('E^2 dN/dE (erg cm-2 s-1])')
plt.legend(loc=0)
def setup_class(cls):
cls.cat = SourceCatalog1FHL()
def setup_class(cls):
cls.cat = SourceCatalog1FHL()
# Use 1FHL J0534.5+2201 (Crab) as a test source
cls.source_name = "1FHL J0534.5+2201"
cls.source = cls.cat[cls.source_name]
""" Example to show how to plot spectrum of Fermi/LAT sources
"""
import matplotlib.pyplot as plt
from gammapy.catalog import SourceCatalog3FGL, SourceCatalog2FHL, SourceCatalog1FHL, SourceCatalog3FHL
from gammapy.utils.energy import EnergyBounds
# # 2155
name = 'PKS 2155-304'
# load catalogs
fermi_3fgl = SourceCatalog3FGL()
fermi_2fhl = SourceCatalog2FHL()
fermi_1fhl = SourceCatalog1FHL()
fermi_3fhl = SourceCatalog3FHL()
# access crab data by corresponding identifier
src_3fgl = fermi_3fgl[name]
src_2fhl = fermi_2fhl[name]
src_1fhl = fermi_1fhl[name]
src_3fhl = fermi_3fhl[name]
# 3FHL
ax = src_3fgl.spectral_model.plot(src_3fgl.energy_range,
energy_power=2,
label='Fermi 3FGL',
color='r',
flux_unit='erg-1 cm-2 s-1')
# set up an energy array to evaluate the butterfly
emin, emax = src_3fgl.energy_range