mage_multipanel-redo_as_function.py

''' Same as multipanel-spectrum-ids.py, but plot multiple spectra.  Using for individual galaxies
jrigby, Sept 2016. '''

from builtins import str
import jrr
import matplotlib.pyplot as plt
import numpy as np
mage_mode = "reduction"  
#mage_mode = "released"
methods = ('bystars', 'byneb')  # method of determining systemic redshift

#specs = jrr.mage.wrap_getlist(mage_mode, which_list="all", drop_s2243=True, MWdr=True)
specs = jrr.mage.wrap_getlist(mage_mode, which_list="batch3", drop_s2243=True, MWdr=True)
#specs = jrr.mage.wrap_getlist(mage_mode, which_list="metabatch", drop_s2243=True, MWdr=True)
short_labels = specs['short_label'].unique()
for label in short_labels :
    (sp, resoln, dresoln, LL, z_systemic) = jrr.mage.wrap_open_spectrum(label, mage_mode,  MWdr=True) 
    the_dfs = [sp]
    the_zzs = [z_systemic]  
    the_pdf =  "PDF_Out2/multipanel_" + label + '.pdf'
    zzstr = str(["{0:0.8f}".format(i) for i in the_zzs])
    jrr.plot.echelle_spectrum(the_dfs, the_zzs, LL, outfile=the_pdf, plot_cont=True, norm_by_cont=True, apply_bad=True, title=(label+" "+zzstr), topfid=(1.1,1))#, waverange=(1000,3000))
    plt.close("all")