startWl = 4000
endWl = 5000
hp=True
flat = True
outputFN = 'ImageStack_%s_%s-%s'%(intTime,startWl, endWl)
if flat==True:
outputFN+='_flat'
if hp==True:
outputFN+='_hp'
#outputFN+='.h5'
identifier=outputFN
print outputFN
LC=LightCurve(fileID=identifier,path=path,targetName=None,run=None,verbose=True,showPlot=False)
LC.makeImageStack(imageStackFilename='',dt=intTime,wvlStart=startWl,wvlStop=endWl,
weighted=True, fluxWeighted=False, getRawCount=False,
scaleByEffInt=True, deadTime=100.e-6)
LC.makeAllCentroidFiles(centroidFilenames=['',''])
print LC.centroids
print LC.flags
LC.makeLightCurve(photometryType='aper')
'''
LC = LightCurve(path,verbose=verbose,showPlot=showPlot)
obsFNs = LC.getFNs()
#print 'obs[10]: ',obsFNs[10:11][0].obs()
obsFiles = generateObsObjectList(obsFNs, wvlLowerLimit=startWl, wvlUpperLimit=endWl, loadHotPix=hp, loadWvlCal=True, loadFlatCal=flat, loadSpectralCal=False)
print 'numObs: ',len(obsFiles)
#outputFN = '%ss_%s-%s'%(intTime,startWl, endWl)
outputFN = '%ss_%s'%(intTime,filt)
if flat==True:
outputFN+='_flat'
if flux==True:
outputFN+='_flux'
if hp==True:
outputFN+='_hp'
identifier=outputFN
print outputFN
#replace everything below here with lightCurve.py example under __main()__
=======
LC=LightCurve(fileID=identifier,path=path,targetName=None,run=None,verbose=True,showPlot=False)
LC.makeImageStack(imageStackFilename='',dt=intTime,wvlStart=startWl,wvlStop=endWl,
weighted=True, fluxWeighted=True, getRawCount=False,
scaleByEffInt=True, deadTime=100.e-6, filterName=filt)
print "\n\n###########################\nFINISHED MAKING IMAGE STACKS\n#############################\n\n"
LC.makeAllCentroidFiles(centroidFilenames=['',''])
'''
>>>>>>> 7edaa38518bd8e99dd76c9a379439a900670e4b2
LC = LightCurve(path,verbose=verbose,showPlot=showPlot)
obsFNs = LC.getFNs()
#print 'obs[10]: ',obsFNs[10:11][0].obs()
obsFiles = generateObsObjectList(obsFNs, wvlLowerLimit=startWl, wvlUpperLimit=endWl, loadHotPix=hp, loadWvlCal=True, loadFlatCal=flat, loadSpectralCal=False)
#outputFN = '%ss_%s-%s'%(intTime,startWl, endWl)
outputFN = '%ss_%s'%(intTime,filt)
if flat==True:
outputFN+='_flat'
if flux==True:
outputFN+='_flux'
if hp==True:
outputFN+='_hp'
#outputFN+='.h5'
identifier=outputFN
print outputFN
LC=LightCurve(fileID=identifier,path=path,targetName=None,run=None,verbose=True,showPlot=False)
'''
LC.loadImageStack()
LC.loadAllCentroidFiles()
print LC.centroids
print LC.flags
#print "Finished loading centroid files"
LC.makeLightCurve(photometryType='aper', interpolation='cubic')
'''
LC.loadLightCurve(photometryFilename='',photometryType='aper')
photometryDict=LC.photometry_dict
showPlot=True
startWl = 4000
endWl = 11000
hp=True
flat = True
outputFN = 'ImageStack_%s-%s'%(startWl, endWl)
if flat==True:
outputFN+='_flat'
if hp==True:
outputFN+='_hp'
outputFN+='.h5'
print outputFN
LC = LightCurve(path,verbose=verbose,showPlot=showPlot)
obsFNs = LC.getFNs(fromObs = True, fromImageStack=False)
#print 'obs[10]: ',obsFNs[10:11][0].obs()
obsFiles = generateObsObjectList(obsFNs, wvlLowerLimit=startWl, wvlUpperLimit=endWl, loadHotPix=hp, loadWvlCal=True, loadFlatCal=flat, loadSpectralCal=False)
print 'numObs: ',len(obsFiles)
im_dict = getImages(fromObsFile=True, fromPhotonList=False, fromImageStack=False, obsFiles=obsFiles, integrationTime=10, weighted=flat, fluxWeighted=False, getRawCount=False)
images=im_dict['images']
pixIntTimes=im_dict['pixIntTimes']
print 'numImages: ',len(images)
writeImageStack(images=images, pixIntTimes=pixIntTimes, startTimes=im_dict['startTimes'], intTimes=im_dict['intTimes'], path=path, outputFilename=outputFN)
