def update_lines(self):
""" Update the line indicators with a new redshift.
"""
try:
for l in self.artists['lines']['lines']:
l.remove()
except ValueError:
pass
for art in self.artists['lines']['labels']:
art.remove()
for art in self.artists['lines']['atmos']:
art.remove()
self.artists['lines'] = barak.spec.plotlines(
self.zp1 - 1, self.ax, labels=1, fontsize=10,
lcolor='0.3', lines=self.lines, offsets=False)
if self.DLA is not None:
self.artists['model'].remove()
s = self.spec[self.i]
t,_ = calc_DLA_trans(s.wa, self.zp1-1, 2*self.dvpix,
logN=self.DLA, bHI=15)
self.artists['model'] = self.ax.plot(
s.wa, t*s.co, 'k')[0]
#if not self.showlabels:
# for t in self.artists['lines']['labels']:
# t.set_visible(False)
# for key in self.artists['zlines']:
# for l in self.artists['zlines'][key]:
# try:
# l.remove()
# except ValueError:
# # plot has been removed
# pass
# self.artists['zlines'] = barak.spec.plotlines(
# zp1 - 1, plt.gca(), lines=self.opt.linelist, labels=True)
# plt.draw()
print('z =', self.zp1-1)
fig = plt.figure(figsize=(4.4, 2.5))
fig.subplots_adjust(left=0.13, right=0.97, top=0.95, bottom=0.21)
ax = pl.gca()
sp.fl /= 1000
sp.er /= 1000
sp.co /= 1000
ax.axhline(0, color='0.7')
logN = 16.838, 16.938, 17.038
model = []
for N in logN:
m, ticks = calc_DLA_trans(sp.wa, 2.46408, 7, logN=N, bHI=27)
model.append(m)
plt.fill_between(sp.wa,
sp.co * model[0],
y2=sp.co * model[2],
lw=0,
alpha=0.2,
color='r',
zorder=9)
c0 = between(sp.wa, 3161.52, 3163.82)
ax.plot(sp.wa[c0],
sp.fl[c0],
color='0.7',
fig = plt.figure(figsize=(4.4, 2.5))
fig.subplots_adjust(left=0.13, right=0.97, top=0.95, bottom=0.21)
ax = pl.gca()
sp.fl /= 1000
sp.er /= 1000
sp.co /= 1000
ax.axhline(0, color="0.7")
logN = 16.838, 16.938, 17.038
model = []
for N in logN:
m, ticks = calc_DLA_trans(sp.wa, 2.46408, 7, logN=N, bHI=27)
model.append(m)
plt.fill_between(sp.wa, sp.co * model[0], y2=sp.co * model[2], lw=0, alpha=0.2, color="r", zorder=9)
c0 = between(sp.wa, 3161.52, 3163.82)
ax.plot(sp.wa[c0], sp.fl[c0], color="0.7", lw=0.5, drawstyle="steps-mid", zorder=10)
sp.fl[c0] = np.nan
ax.plot(sp.wa, sp.fl, color="0.4", lw=0.5, drawstyle="steps-mid", zorder=10)
ax.plot(sp.wa, sp.co, "--k", zorder=10, lw=0.5)
plt.plot(sp.wa, sp.co * model[1], "r", lw=1, zorder=10)
ax.set_xlim(3158.1, 3170.9)
