QA_Prod (v1.1)



In [ ]:

    
%matplotlib notebook



In [ ]:

    
# imports
from desispec.qa import qa_prod as dqap

Init

setenv DESI_SPECTRO_DATA /Users/xavier/DESI/DESI_SCRATCH/sim/madrone/



In [ ]:

    
specprod_dir = '/Users/xavier/DESI/DESI_SCRATCH/redux/madrone/'



In [ ]:

    
reload(dqap)
qa_prod = dqap.QA_Prod(specprod_dir)

Remake

In Python



In [ ]:

    
qa_prod.remake_frame_qa(remake_plots=True)

Script

desi_qa_prod --specprod_dir /Users/xavier/DESI/DESI_SCRATCH/redux/madrone --remake_frame 3

Slurp

In Python



In [ ]:

    
qa_prod.slurp(remove=False)

Script

desi_qa_prod --specprod_dir /Users/xavier/DESI/DESI_SCRATCH/redux/madrone/ --slurp --make_frameqa=1

Simple Figs



In [ ]:

    
from desispec.qa import qa_prod as dqap



In [ ]:

    
specprod_dir = '/Users/xavier/DESI/DESI_SCRATCH/redux/madrone/'
reload(dqap)
qa_prod = dqap.QA_Prod(specprod_dir)

Load QA



In [ ]:

    
qa_prod.load_data()



In [ ]:

    
qa_prod.data['20160607'][6]['b0']#['SKYSUB']['QA']

Grab an array of QA values



In [ ]:

    
sky_resid_b, ne_dict = qa_prod.get_qa_array('SKYSUB', 'MED_RESID', channels=['b'])
sky_resid_b



In [ ]:

    
ne_dict

Histogram



In [ ]:

    
plt.clf()
ax = plt.gca()
ax.hist(sky_resid_b)
ax.set_xlim(-1,1)
ax.set_xlabel('MED_SKY_RESID')

Camera Histograms



In [ ]:

    
from desispec.qa import qa_plots as dqqp



In [ ]:

    
reload(dqqp)
dqqp.prod_channel_hist(qa_prod, 'SKYSUB', 'MED_RESID', xlim=(-1,1))

Script

profx.ucolick.org> desi_qa_prod -h                                              
usage: desi_qa_prod [-h] --specprod_dir SPECPROD_DIR
                    [--make_frameqa MAKE_FRAMEQA] [--slurp] [--remove]
                    [--clobber] [--channel_hist CHANNEL_HIST]

Generate Production Level QA

optional arguments:
  -h, --help            show this help message and exit
  --specprod_dir SPECPROD_DIR
                        Path containing the exposures/directory to use
  --make_frameqa MAKE_FRAMEQA
                        Bitwise flag to control remaking the QA files (1) and
                        figures (2) for each frame in the production
  --slurp               slurp production QA files into one?
  --remove              remove frame QA files?
  --clobber             clobber existing QA files?
  --channel_hist CHANNEL_HIST
                        Generate channel histogram(s)



In [ ]: