def plot_doc(): """ """ img = xFITSImage(os.path.join(XIMPOL_DATA, 'tycho_cmap.fits')) fig = img.plot(show=False) xFITSImage.add_label(fig, 'XIPE 300 ks') plt.show()
def plot_doc(): """ """ img = xFITSImage(os.path.join(XIMPOL_DATA, 'crab_complex_cmap.fits')) fig = img.plot(show=False) xFITSImage.add_label(fig, 'XIPE %d ks' %DURATION/1000.) plt.show()
def plot_doc(): """ """ img = xFITSImage(os.path.join(XIMPOL_DATA, 'casa_cmap.fits')) fig = img.plot(show=False) #Option to draw the psf circle on the count map RAD_PSF = 11/60. fig.show_circles(350.769, 58.873, RAD_PSF/60., lw=2, color='white') fig.add_label(0.73,0.90, 'PSF', relative=True, size='x-large', color='white', horizontalalignment='left') xFITSImage.add_label(fig, 'XIPE 250 ks') plt.show()
def display(): """Display the source model. """ from ximpol.utils.matplotlib_ import pyplot as plt from ximpol.srcmodel.img import xFITSImage print(ROI_MODEL) fig = plt.figure('Energy spectrum') total_spectral_model.plot(logy=True, show=False, label='Total') nonthermal_spectral_model.plot(logy=True, show=False, label='Non-thermal') thermal_spectral_model.plot(logy=True, show=False, label='Thermal') plt.legend(bbox_to_anchor=(0.95, 0.95)) fig = thermal_component.image.plot(show=False) xFITSImage.add_label(fig, 'Chandra 1.5-3.0 keV') fig = nonthermal_component.image.plot(show=False) xFITSImage.add_label(fig, 'Chandra 4.0-6.0 keV') img = xFITSImage(he_img_file_path) fig = img.plot(show=False) polarization_map.build_grid_sample(ROI_MODEL.ra, ROI_MODEL.dec) polarization_map.overlay_arrows(fig) plt.show()