Ugh

Trying to retrace my steps for creating LCS_all_size.fits because there seems to be some discrepancies in how the simard+2011 catalog was matched to my sources.

There are a few galaxies that are in Simard, but their values for B/T where not propagated through to the final fits table - yikes!

So what is the process to make LCS_all_size.fits

LCSmergespiralcats.py input = research/LocalClusters/MasterTables/'+clustername+'mastertable.WithProfileFits.fits'



In [3]:

    
from matplotlib import pyplot as plt
import numpy as np
%matplotlib inline



In [4]:

    
y = np.array([0.,1.,2.,3.,4.,5.],'f')
t = np.array([0.,.54,.72,.91,1.01,1.17],'f')



In [5]:

    
plt.figure()
plt.plot(t,y,'bo')









    Out[5]:





[<matplotlib.lines.Line2D at 0x10fdf2a90>]



In [11]:

    
plt.figure()
plt.plot(t**2,y,'bo')
plt.xlabel('Time Squared')
plt.ylabel('Distance')
c = np.polyfit(t**2,y,1)
tfit = np.arange(0,1.4,.01)
yfit = np.polyval(c,tfit)
plt.plot(tfit,yfit,'r-')
print c









    



[ 3.73263261 -0.00528079]



In [12]:

    
tfit









    Out[12]:





array([ 0.  ,  0.01,  0.02,  0.03,  0.04,  0.05,  0.06,  0.07,  0.08,
        0.09,  0.1 ,  0.11,  0.12,  0.13,  0.14,  0.15,  0.16,  0.17,
        0.18,  0.19,  0.2 ,  0.21,  0.22,  0.23,  0.24,  0.25,  0.26,
        0.27,  0.28,  0.29,  0.3 ,  0.31,  0.32,  0.33,  0.34,  0.35,
        0.36,  0.37,  0.38,  0.39,  0.4 ,  0.41,  0.42,  0.43,  0.44,
        0.45,  0.46,  0.47,  0.48,  0.49,  0.5 ,  0.51,  0.52,  0.53,
        0.54,  0.55,  0.56,  0.57,  0.58,  0.59,  0.6 ,  0.61,  0.62,
        0.63,  0.64,  0.65,  0.66,  0.67,  0.68,  0.69,  0.7 ,  0.71,
        0.72,  0.73,  0.74,  0.75,  0.76,  0.77,  0.78,  0.79,  0.8 ,
        0.81,  0.82,  0.83,  0.84,  0.85,  0.86,  0.87,  0.88,  0.89,
        0.9 ,  0.91,  0.92,  0.93,  0.94,  0.95,  0.96,  0.97,  0.98,
        0.99,  1.  ,  1.01,  1.02,  1.03,  1.04,  1.05,  1.06,  1.07,
        1.08,  1.09,  1.1 ,  1.11,  1.12,  1.13,  1.14,  1.15,  1.16,
        1.17,  1.18,  1.19,  1.2 ,  1.21,  1.22,  1.23,  1.24,  1.25,
        1.26,  1.27,  1.28,  1.29,  1.3 ,  1.31,  1.32,  1.33,  1.34,
        1.35,  1.36,  1.37,  1.38,  1.39])



In [ ]: