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