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')
