In [6]:

    

from matplotlib import pyplot as plt
import numpy as np
#from numpy.lib.recfunctions import append_fields
from astropy.io import fits
import glob
%matplotlib inline


    

In [21]:

    

class Cluster:
    def __init__(self,prefix):
        self.prefix = prefix
        # read in NSA
        nsafile = prefix+'_NSA.fits'
        self.nsa = fits.getdata(nsafile)
        # read in AGC
        agcfile= prefix+'_matchAGCtoNSA.fits'
        self.agc = fits.getdata(agcfile)
    def plotpositions(self):
        plt.figure()
        plt.plot(self.nsa.RA,self.nsa.DEC,'k.')
    def plotcolormag(self):
        plt.figure()
        self.NUVr=self.nsa.ABSMAG[:,1] - self.nsa.ABSMAG[:,4]
        plt.plot(self.nsa.ABSMAG[:,4],self.NUVr,'bo')


    

In [22]:

    

clusters=['Abell2063']
a2063 = Cluster(clusters[0])


    

In [14]:

    

a2063.plotpositions()


    

In [25]:

    

a2063.plotcolormag()
plt.axis([-23,-16,0,7])


    

    
Out[25]:





[-23, -16, 0, 7]






    

In [ ]: