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)