def main(root, visit, ccd, fakes=None, blends=False, listobj=16, filt=None): butler = lsst.daf.persistence.Butler(root) dataId = { 'visit': visit, 'ccd': int(ccd) } if filt is None else { 'tract': visit, 'patch': ccd, 'filter': filt } if fakes is not None: src = matchFakes.getFakeSources(butler, dataId, extraCols=('zeropoint'), radecMatch=fakes) else: src = butler.get('src' if filt is None else 'deepCoadd-src', dataId) if not blends: src = [ s for s in src if (s.get('deblend.nchild') == 0) & (s.get('parent') == 0) ] else: src = [s for s in src if (s.get('deblend.nchild') == 0)] exposure = butler.get('calexp' if filt is None else 'deepCoadd', dataId) if type(listobj) is int: listobj = numpy.random.choice(range(len(src)), listobj, False) srcList = [src[i] for i in listobj] getMosaic(srcList, exposure, 'fakeId' if fakes else 'id')
def main(root, visit, ccds, galType='sersic', output='outputs/'): """ main function for controlling fake source comparison """ butler = dafPersist.Butler(root) fakeTable = None for ccd in ccds: dataId = {'visit': visit, 'ccd': ccd} try: temp = getFakeSources(butler, dataId, tol=1.0, visit=True, ccd=True) except: continue if fakeTable is None: fakeTable = temp.copy(True) else: fakeTable.extend(temp, True) npTable = writeNumpyTable(fakeTable) rerunName = root.split('/')[-2] fitsTable = astropy.table.Table(npTable).write(output + '/' + rerunName + '_galMags.txt', format='ascii')
def main(root, visit, ccd, fakes=None, blends=False, listobj=16, filt=None): butler = lsst.daf.persistence.Butler(root) dataId = {'visit': visit, 'ccd': int(ccd)} if filt is None else {'tract': visit, 'patch': ccd, 'filter': filt} if fakes is not None: src = matchFakes.getFakeSources(butler, dataId, extraCols=('zeropoint'), radecMatch=fakes) else: src = butler.get('src' if filt is None else 'deepCoadd-src', dataId) if not blends: src = [s for s in src if ((s.get('deblend.nchild') == 0) & (s.get('parent') == 0))] else: src = [s for s in src if (s.get('deblend.nchild') == 0)] exposure = butler.get('calexp' if filt is None else 'deepCoadd', dataId) if type(listobj) is int: listobj = numpy.random.choice(range(len(src)), listobj, False) srcList = [src[i] for i in listobj] getMosaic(srcList, exposure, 'fakeId' if fakes else 'id')
def main(root, visit, ccds, galType='sersic', output='outputs/'): """ main function for controlling fake source comparison """ butler = dafPersist.Butler(root) fakeTable = None for ccd in ccds: dataId = {'visit':visit, 'ccd':ccd} try: temp = getFakeSources(butler, dataId, tol=1.0, visit=True, ccd=True) except: continue if fakeTable is None: fakeTable = temp.copy(True) else: fakeTable.extend(temp, True) npTable = writeNumpyTable(fakeTable) rerunName = root.split('/')[-2] fitsTable = astropy.table.Table(npTable).write(output+'/'+rerunName+ '_galMags.txt', format='ascii')