In [1]:
import numpy as np
import matplotlib.pyplot as plt
import os
In [2]:
jpas_filters_path = '/home/mldantas/Dropbox/DoutoradoIAG/Challenge/Filters/JPAS_filters'
jplus_filters_path = '/home/mldantas/Dropbox/DoutoradoIAG/Challenge/Filters/JPLUS_SDSS_filters'
jpas_filters_list = np.loadtxt('/home/mldantas/Dropbox/DoutoradoIAG/Challenge/Filters/jpas_filters_list.txt', dtype=str)
jplus_filters_list = np.loadtxt('/home/mldantas/Dropbox/DoutoradoIAG/Challenge/Filters/jplus_filters_list.txt', dtype=str)
In [7]:
for each_filter in jpas_filters_list:
wavelength = np.loadtxt(os.path.join(jpas_filters_path, each_filter), dtype=float, usecols=[0])
flux = np.loadtxt(os.path.join(jpas_filters_path, each_filter), dtype=float, usecols=[1])
plt.plot(wavelength, flux, '-', linewidth=0.5)
# plt.fill_between(wavelength, flux, interpolate=True)
plt.xlabel('Wavelength ($\AA$)', fontsize=15)
plt.ylabel('Throughput', fontsize=15)
plt.gcf()
plt.savefig('/home/mldantas/Dropbox/DoutoradoIAG/Challenge/Filters/jpas_filters.pdf', dpi = 100)
plt.show()
plt.clf()
In [6]:
for each_filter in jplus_filters_list:
wavelength = np.loadtxt(os.path.join(jplus_filters_path, each_filter), dtype=float, usecols=[0])
flux = np.loadtxt(os.path.join(jplus_filters_path, each_filter), dtype=float, usecols=[1])
plt.plot(wavelength, flux, '-')
plt.xlabel('Wavelength ($\AA$)', fontsize=15)
plt.ylabel('Throughput', fontsize=15)
plt.gcf()
plt.savefig('/home/mldantas/Dropbox/DoutoradoIAG/Challenge/Filters/jplus_filters.pdf', dpi = 100)
plt.show()
In [ ]: