directory = '/Users/ryanlindley/Research/CIV/' + name
os.chdir(directory)
os.getcwd()
# In[460]:
wl1, f1, gamma1, elem, state = (1548.187, 0.19, 0.00324, 'C', 'IV')
wl2, f2, gamma2, elem, state = (1550.772, 0.0952, 0.00325, 'C', 'IV')
path = '/Users/ryanlindley/Research/CIV/' + name + '/CIV.bin3.linespec'
wl, fl, fe, ctn = auto.np.loadtxt(path).transpose()
lsf = auto.np.loadtxt('/Users/ryanlindley/Research/CIV.old_setup/CIV.lsf')
# In[461]:
new_ctn = auto.find_continuum(name, wl, fl, fe, ctn, [True] * len(wl))
#ctn_flag = auto.mask_wl(wl, [1547.9, 1548.5, 1550.4, 1551.2, 1552.5, 1555]) # add custom ctn if needed
#manual_ctn = auto.find_continuum(name, wl, fl, fe, ctn, ctn_flag)
# In[462]:
def check_larger_continuum(name, wl, ctn):
spec = auto.np.loadtxt('/Users/ryanlindley/Research/CIV/' + name +
'/full.spec')
auto.plt.plot(spec[:, 0], spec[:, 1])
auto.plt.plot(wl, ctn)
auto.plt.xlim(1542, 1556)
auto.plt.ylim(0, 1.5e-14)
auto.plt.show
# In[178]:
wl1, f1, gamma1, elem, state = (1548.187, 0.19, 0.00324, 'C', 'IV')
wl2, f2, gamma2, elem, state = (1550.772, 0.0952, 0.00325, 'C', 'IV')
path = '/Users/ryanlindley/Research/CIV/' + name + '/CIV.bin3.linespec'
wl, fl, fe, ctn = auto.np.loadtxt(path).transpose()
lsf = auto.np.loadtxt('/Users/ryanlindley/Research/CIV.old_setup/CIV.lsf')
# In[180]:
#new_ctn = auto.find_continuum(name, wl, fl, fe, ctn, [True]*len(wl))
ctn_flag = auto.mask_wl(
wl, [1547, 1548.8, 1550, 1551.3]) # add custom ctn if needed
manual_ctn = auto.find_continuum(name, wl, fl, fe, ctn, ctn_flag)
# In[181]:
def check_larger_continuum(name, wl, ctn):
spec = auto.np.loadtxt('/Users/ryanlindley/Research/CIV/' + name +
'/full.spec')
auto.plt.plot(spec[:, 0], spec[:, 1])
auto.plt.plot(wl, ctn)
auto.plt.xlim(1542, 1556)
#auto.plt.ylim(0, 2e-14)
auto.plt.show
check_larger_continuum(name, wl, manual_ctn)
