python-programs-1/superposition-spetrum-curve-JPAS2017.py at master · AngelGSoto/python-programs-1 · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
'''
Superponient the JPAS transmission curve and spectrum of a symbiotic
'''
from __future__ import print_function
import numpy as np
import glob
import json
import matplotlib.pyplot as plt
import seaborn as sns
import os

def load_file(filename):
    wll, ress = [], []
    data = np.loadtxt(filename, delimiter=None, converters=None, skiprows=0,
                                       usecols=None, unpack=False, ndmin=0)
    for i in data:
        wl = str(i[0])
        res = str(i[1])
        wll.append(wl)
        ress.append(res)
    return wll, ress

pattern = "JPAS_J*.tab"
#pattern = "Jv0915_8100.res"
file_list = glob.glob(pattern)

pattern1 = "JPAS_*SDSS.tab"
file_list1 = glob.glob(pattern1)

fig = plt.figure(figsize=(12, 7))
ax1 = fig.add_subplot(111)
for f in file_list:
    x, y = load_file(f)
    #plt.plot(x,y, c="gray")#, label=f.split('_t')[0])
    #plt.xticks(x, file_list, rotation=90)

for f in file_list1:
    x, y = load_file(f)
    ax1.plot(x,y, linewidth=2.3)#, label=f.split('_t')[0])

# Spectrum
def sys(spectra):
    #datadir = "../../Halo-PNe-spectros/"
    datadir = "../../PN-Tesis-Didactica/"
    file_ = spectra
    x = np.loadtxt(os.path.join(datadir, file_), delimiter = None, skiprows = 0, usecols = None,
                   unpack = False, dtype = np.dtype([('Wl', '|f8'), ('Flux', 'f8')]))
    return x['Wl'], x['Flux']

x, y = sys("PN-M157-medium.dat")
#x1, y1 = sys("LHa_115_N_60.txt")


y /=80.0e-15
#y1 /=1.5e-14

#print(x, y)
lgd_kws = {'frameon': True, 'fancybox': True, 'shadow': True}
sns.set(style="dark")
#plt.xlim(6e3, 10500)
# plt.xlim(3e3, 10500)
# plt.ylim(0.0, 0.7)
plt.tick_params(axis='x', labelsize=15)
plt.tick_params(axis='y', labelsize=15)
#plt.plot(x, y, linewidth=1.3, color="k")#, label="Lin 358")
#plt.plot(x1, y1, linewidth=1.3, color="green", label="LHa 115")
#plt.legend(fontsize='x-small')
plt.xlabel("Wavelength($\AA$)", fontsize= 18)
plt.ylabel("Transmission", fontsize= 18)
#plt.legend()

plt.savefig('JPAS2017-filter-M157-medium.jpg')