In [1]:

    

import galaxy_tests
import ResolvedStellarPops as rsp
#import ResolvedStellarPops.PadovaTracks as pc
import os
import numpy as np


    

In [3]:

    

%pylab inline


    

    




Populating the interactive namespace from numpy and matplotlib

In [42]:

    

reload(rsp.Galaxies)


    

    
Out[42]:





<module 'ResolvedStellarPops.Galaxies' from '/home/rosenfield/research/python/ResolvedStellarPops/Galaxies.py'>

In [ ]:

    

fits_src='/home/rosenfield/research/TP-AGBcalib/SNAP/data/angst_no_trim'
fits_files = rsp.fileIO.get_files(fits_src, '*fits')
gal_kw = {'hla': True, 'photsys': 'wfc3snap', 'band': 'opt', 'angst': True,
          'filetype': 'fitstable'}


    

In [33]:

    

ancient_list = ['ddo71', 'kkh37', 'ngc2976', 'ngc404', 'kdg73', 'hs117', 'ddo78']
ancient_fits = np.concatenate([[f for a in ancient_list if a in f] for f in fits_files])


    

In [43]:

    

fig, (axs) = rsp.graphics.GraphicsUtils.setup_multiplot(len(ancient_fits), figsize=(30,30))
axs = np.squeeze(np.concatenate(axs))

for i, fits_file in enumerate(ancient_fits):
    gal = rsp.Galaxies.galaxy(fits_file, **gal_kw)
    gal.plot_cmd(gal.color, gal.mag1, levels=3, threshold=100, ax=axs[i],
                 yfilter=gal.filter1)
    gal.photsys = 'wfc3snap'
    gal.decorate_cmd(ax=axs[i], trgb=True, filter1=gal.filter1)


    

In [ ]: