def testFitExponRealData(self):
"""
Read timestamps from part of one data file where there was an event.
"""
# We imported FileName at the top of this file. Use this to
# locate the file for this specific observation.
# on turk, the file with the observation is this one:
# /ScienceData/PAL2012/20121211/obs_20121212-1333303.h5
run = 'PAL2012'
sundownDate = '20121211'
obsDate = '20121212'
seq = '133303'
fn = FileName(run, sundownDate, obsDate+"-"+seq)
# The class Cosmic is in ARCONS-pipeline/cosmic/Cosmic.py
cosmic = Cosmic(fn)
t0 = 138595526
t1 = 138595666
dictionary = cosmic.fitExpon(t0,t1)
hist = dictionary['timeHgValues']
bins = np.arange(len(hist))
width = 0.7*(bins[1]-bins[0])
center = (bins[:-1]+bins[1:])/2
#plt.step(center, hist, where = 'post')
plt.plot(bins,hist, label='data')
xFit = np.linspace(0,len(hist),10*len(hist))
pFit = dictionary['pExponFit']
yFit = Cosmic.funcExpon(xFit,*pFit)
plt.plot(xFit,yFit, label='Exponential Fit')
xFit = np.linspace(0,len(hist),10*len(hist))
pFit = dictionary['pGaussFit']
print "pFit=",pFit
yFit = Cosmic.funcGauss(xFit,*pFit)
plt.plot(xFit,yFit, label='Gaussian Fit')
plt.legend()
xLimit = dictionary['xLimit']
plt.axvline(x=xLimit[0], color='magenta', linestyle='-.')
plt.axvline(x=xLimit[1], color='magenta', linestyle='-.')
plt.savefig(inspect.stack()[0][3]+".png")
