lsst.synpipe.compareModel Namespace Reference

Functions
def	getParams (record, galType='sersic')
def	getMag (record, fluxType='cmodel')
def	writeNumpyTable (fakeTable)
def	main (root, visit, ccds, galType='sersic', output='outputs/')

Variables
	parser = argparse.ArgumentParser()
	help
	type
	int
	nargs
	str
	dest
	choices
	args = parser.parse_args()
	ccd
	galType
	output

Detailed Description

Compare the input fake galaxy model to the output one

Function Documentation

getMag()

def lsst.synpipe.compareModel.getMag	(		record,
			fluxType = 'cmodel'
	)

return the magnitude and error

Definition at line 33 of file compareModel.py.

def getMag(record, fluxType='cmodel'):
    """
    return the magnitude and error
    """
    flux, fluxerr = record.get(fluxType), record.get(fluxType+'.err')
    mag, magerr = -2.5 * np.log10(flux), 2.5/np.log(10.0)*fluxerr/flux
    mag += record.get('zeropoint')
    return mag, magerr

getParams()

def lsst.synpipe.compareModel.getParams	(		record,
			galType = 'sersic'
	)

return the semi-major axis, axis ratio and position angle (in degrees)

Definition at line 17 of file compareModel.py.

def getParams(record, galType='sersic'):
    """
    return the semi-major axis, axis ratio and position angle (in degrees)
    """
    fluxType = {'dev': 'cmodel.dev',
                'exp': 'cmodel.exp',
                'sersic': 'cmodel',
                'cmodel': 'cmodel'}[galType]

    e = lsst.afw.geom.ellipses.Axes(record.get(fluxType + '.ellipse'))
    q = (e.getB() / e.getA())
    reff = e.getA()
    pa = (e.getTheta() * 180.0 / np.pi)
    return reff, q, pa

main()

def lsst.synpipe.compareModel.main	(		root,
			visit,
			ccds,
			galType = 'sersic',
			output = 'outputs/'
	)

main function for controlling fake source comparison

Examples:

forEachPixel.cc, image1.cc, image2.cc, imageStatistics.cc, maskedImage1.cc, maskedImage2.cc, spatialCellExample.cc, splineInterpolate.cc, and statistics.cc.

Definition at line 83 of file compareModel.py.

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 Exception:
            continue
        if fakeTable is None:
            fakeTable = temp.copy(True)
        else:
            fakeTable.extend(temp, True)

    npTable = writeNumpyTable(fakeTable)
    rerunName = root.split('/')[-2]
    tableFileName = output + '/' + rerunName + '_galMags.txt'
    astropy.table.Table(npTable).write(tableFileName, format='ascii')

writeNumpyTable()

def lsst.synpipe.compareModel.writeNumpyTable

(

fakeTable

)

Writes output to numpy

Definition at line 43 of file compareModel.py.

def writeNumpyTable(fakeTable):
    """Writes output to numpy
    """
    npTable = np.recarray(len(fakeTable),
                          dtype={'names': ['id', 'fakeid', 'visit', 'ccd',
                                           'cmodelMag', 'expMag',
                                           'devMag',
                                           'expReff', 'devReff',
                                           'expBA', 'devBA',
                                           'expPosAng', 'devPosAng',
                                           'cmodelMagErr', 'expMagErr',
                                           'devMagErr', 'nchild', 'parent'],
                                 'formats': [int, int, int, int,
                                             float, float,
                                             float, float, float, float,
                                             float, float, float, float,
                                             float, float, int, int]})
    for indFake, fake in enumerate(fakeTable):
        npTable[indFake]['id'] = fake.get('id')
        npTable[indFake]['fakeid'] = fake.get('fakeId')
        npTable[indFake]['ccd'] = fake.get('ccd')
        npTable[indFake]['visit'] = fake.get('visit')
        nameMatch = {'sersic': 'cmodel',
                     'exp': 'cmodel.exp',
                     'cmodel': 'cmodel',
                     'dev': 'cmodel.dev'}
        for name in ('cmodel', 'exp', 'dev'):
            m1, m2 = getMag(fake, nameMatch[name]+'.flux')
            npTable[indFake][name+'Mag'] = m1
            npTable[indFake][name+'MagErr'] = m2
        for name in ('exp', 'dev'):
            params = getParams(fake, name)
            npTable[indFake][name+'Reff'] = params[0]
            npTable[indFake][name+'BA'] = params[1]
            npTable[indFake][name+'PosAng'] = params[2]
        npTable[indFake]['nchild'] = fake.get('deblend.nchild')
        npTable[indFake]['parent'] = fake.get('parent')
    return npTable

Variable Documentation

args

lsst.synpipe.compareModel.args = parser.parse_args()

Definition at line 116 of file compareModel.py.

ccd

lsst.synpipe.compareModel.ccd

Definition at line 118 of file compareModel.py.

choices

lsst.synpipe.compareModel.choices

Definition at line 113 of file compareModel.py.

dest

lsst.synpipe.compareModel.dest

Definition at line 112 of file compareModel.py.

galType

lsst.synpipe.compareModel.galType

Definition at line 118 of file compareModel.py.

help

lsst.synpipe.compareModel.help

Definition at line 109 of file compareModel.py.

int

lsst.synpipe.compareModel.int

Definition at line 110 of file compareModel.py.

nargs

lsst.synpipe.compareModel.nargs

Definition at line 111 of file compareModel.py.

output

lsst.synpipe.compareModel.output

Definition at line 119 of file compareModel.py.

parser

lsst.synpipe.compareModel.parser = argparse.ArgumentParser()

Definition at line 108 of file compareModel.py.

str

lsst.synpipe.compareModel.str

Definition at line 112 of file compareModel.py.

type

lsst.synpipe.compareModel.type

Definition at line 110 of file compareModel.py.