def addReferenceCatalogs(self, rc):
""" addReferenceCatalogs primitive
"""
adOuts = []
try:
log.info(" primitive addReferenceCatalogs")
for ad in rc.get_inputs(style='AD'):
cc = AddReferenceCatalogs(ad, logfile='',
catalogName=rc['catalogName'])
adout = cc.getRefs()
adOuts.append(adout)
rc.report_output(adOuts)
except:
raise FluxCalException("Problems with addReferenceCatalogs")
yield rc
def fluxcal(self, rc):
"""
Fluxcal drives:
- DetectSources
- AddReferenceCatalogs
- CorrelateWithReferenceCatalogs
- CalculateZeropoint
"""
logfile = 'fluxcal.log'
adOuts = []
try:
log.status('*STARTING* Running fluxcal functions')
for ad in rc.get_inputs(style='AD'):
ds = DetectSources (ad, logfile=logfile, sigma=rc['sigma'],
threshold=rc['threshold'], fwhm=rc['fwhm'])
adout = ds.runDS()
cc = AddReferenceCatalogs(adout, logfile=logfile,
catalogName=rc['catalogName'])
adout = cc.getRefs()
corr = CorrelateWithReferenceCatalogs (adout, logfile=logfile,
delta=rc['delta'])
adout = corr.runCorr()
czp = CalculateZeropoint(adout, logfile=logfile,
zplogfile=rc['zplogfile'],
extinction=rc['extinction'])
adout = czp.runZP()
adOuts.append(adout)
rc.report_output(adOuts)
log.status('Fluxcal completed successfully')
except:
raise FluxCalException("Problems with fluxcal")
yield rc
def runFC(self):
log = self.log
adout = self.adout
logfile = self.logfile
adout = self._detsources()
cc = AddReferenceCatalogs(adout, logfile=logfile)
adout = cc.getRefs()
corr = CorrelateWithReferenceCatalogs (adout, logfile=logfile)
adout = corr.runCorr()
czp = CalculateZeropoint(adout, logfile=logfile)
adout = czp.runZP()
self.adout = adout
ofile = 'zp_'+adout.filename
log.info("writing output fits file: "+ofile)
if adout: adout.write(ofile, clobber=True)