lsst.pipe.tasks.mocks.simpleMapper Namespace Reference

Classes
class	PersistenceType
class	BypassPersistenceType
class	ExposurePersistenceType
class	SkyMapPersistenceType
class	CatalogPersistenceType
class	SimpleCatalogPersistenceType
class	SourceCatalogPersistenceType
class	ExposureCatalogPersistenceType
class	PeakCatalogPersistenceType
class	SimpleMapping
class	RawMapping
class	SkyMapping
class	TempExpMapping
class	ForcedSrcMapping
class	MapperMeta
class	SimpleMapper

Functions
def	makeSimpleCamera (nX, nY, sizeX, sizeY, gapX, gapY, pixelSize=1.0, plateScale=20.0, radialDistortion=0.925)
def	makeDataRepo (root)

Function Documentation

makeDataRepo()

def lsst.pipe.tasks.mocks.simpleMapper.makeDataRepo

(

root

)

Create a data repository for SimpleMapper and return a butler for it.

Clobbers anything already in the given path.

Definition at line 499 of file simpleMapper.py.

def makeDataRepo(root):
    """
    Create a data repository for SimpleMapper and return a butler for it.
 
    Clobbers anything already in the given path.
    """
    if os.path.exists(root):
        shutil.rmtree(root)
    os.makedirs(root)
    with open(os.path.join(root, "_mapper"), "w") as f:
        f.write("lsst.pipe.tasks.mocks.SimpleMapper\n")
    return lsst.daf.persistence.Butler(root=root)

makeSimpleCamera()

def lsst.pipe.tasks.mocks.simpleMapper.makeSimpleCamera	(		nX,
			nY,
			sizeX,
			sizeY,
			gapX,
			gapY,
			pixelSize = 1.0,
			plateScale = 20.0,
			radialDistortion = 0.925
	)

Create a camera

@param[in] nx: number of detectors in x
@param[in] ny: number of detectors in y
@param[in] sizeX: detector size in x (pixels)
@param[in] sizeY: detector size in y (pixels)
@param[in] gapX: gap between detectors in x (mm)
@param[in] gapY: gap between detectors in y (mm)
@param[in] pixelSize: pixel size (mm) (a float)
@param[in] plateScale: plate scale in arcsec/mm; 20.0 is for LSST
@param[in] radialDistortion: radial distortion, in mm/rad^2
    (the r^3 coefficient of the radial distortion polynomial
    that converts FIELD_ANGLE in radians to FOCAL_PLANE in mm);
    0.925 is the value Dave Monet measured for lsstSim data

Each detector will have one amplifier (with no raw information).

Definition at line 437 of file simpleMapper.py.

):
    """Create a camera
 
    @param[in] nx: number of detectors in x
    @param[in] ny: number of detectors in y
    @param[in] sizeX: detector size in x (pixels)
    @param[in] sizeY: detector size in y (pixels)
    @param[in] gapX: gap between detectors in x (mm)
    @param[in] gapY: gap between detectors in y (mm)
    @param[in] pixelSize: pixel size (mm) (a float)
    @param[in] plateScale: plate scale in arcsec/mm; 20.0 is for LSST
    @param[in] radialDistortion: radial distortion, in mm/rad^2
        (the r^3 coefficient of the radial distortion polynomial
        that converts FIELD_ANGLE in radians to FOCAL_PLANE in mm);
        0.925 is the value Dave Monet measured for lsstSim data
 
    Each detector will have one amplifier (with no raw information).
    """
    pScaleRad = lsst.geom.arcsecToRad(plateScale)
    radialDistortCoeffs = [0.0, 1.0/pScaleRad, 0.0, radialDistortion/pScaleRad]
    focalPlaneToFieldAngle = lsst.afw.geom.makeRadialTransform(radialDistortCoeffs)
 
    ccdBBox = lsst.geom.Box2I(lsst.geom.Point2I(), lsst.geom.Extent2I(sizeX, sizeY))
 
    cameraBuilder = lsst.afw.cameraGeom.Camera.Builder("Simple Camera")
 
    detectorId = 0
    for iY in range(nY):
        cY = (iY - 0.5 * (nY - 1)) * (pixelSize * sizeY + gapY)
        for iX in range(nX):
            cX = (iX - 0.5 * (nX - 1)) * (pixelSize * sizeY + gapX)
            fpPos = lsst.geom.Point2D(cX, cY)
            detectorName = "detector %d,%d" % (iX, iY)
 
            detectorBuilder = cameraBuilder.add(detectorName, detectorId)
            detectorBuilder.setSerial(detectorName + " serial")
            detectorBuilder.setBBox(ccdBBox)
            detectorBuilder.setOrientation(lsst.afw.cameraGeom.Orientation(fpPos))
            detectorBuilder.setPixelSize(lsst.geom.Extent2D(pixelSize, pixelSize))
 
            ampBuilder = lsst.afw.cameraGeom.Amplifier.Builder()
            ampName = "amp"
            ampBuilder.setName(ampName)
            ampBuilder.setBBox(ccdBBox)
            ampBuilder.setGain(1.0)
            ampBuilder.setReadNoise(5.0)
 
            detectorBuilder.append(ampBuilder)
 
            detectorId += 1
 
    cameraBuilder.setTransformFromFocalPlaneTo(lsst.afw.cameraGeom.FIELD_ANGLE, focalPlaneToFieldAngle)
    return cameraBuilder.finish()