bestChiSq > 100 and ngoodobs > 100
In [22]:
import matplotlib.pyplot as plt
%matplotlib notebook
import numpy as np
from astropy.table import Table as tbl
In [12]:
source = tbl.read('allsky4.tbl', format = 'ipac')
In [24]:
ret = plt.hist(source['bestmeanmag'], range = (12,27), bins = 50)
plt.xlabel('Mean Mag')
plt.ylabel('Counts')
Out[24]:
This is showing the best mean magnitude observed for all sources. There appears to be a linear trend for the counts between ~14 and ~18.5 magnitude stars. Then, the distrubtion turns into a more typical Gaussian-like shape. However, this may be a exponential trend upward that just falls off at higher magnitudes because of the telescope being used and the time it takes to observe fainter stars. The latter seems like a more probable explanation as the number of stars increases at higher magnitudes. Additionally, I am somewhat concerned about the area surrounding the 26th magnitude as this seems fainter than what could be measured, and there is a broad area in which very few observations were made. These seem to be extraneous points due to some unknown source.
Below are the counts and cent
In [25]:
print(ret)
In [ ]: