LSST Applications g2079a07aa2+86d27d4dc4,g2305ad1205+a659bff248,g2bbee38e9b+3c60f8fe34,g337abbeb29+3c60f8fe34,g33d1c0ed96+3c60f8fe34,g3502564af9+d77d6d1350,g3a166c0a6a+3c60f8fe34,g487adcacf7+25d9892218,g4be5004598+d77d6d1350,g50ff169b8f+96c6868917,g52b1c1532d+585e252eca,g591dd9f2cf+4d81263f9a,g5cd07815a0+980d2b1c3b,g607f77f49a+d77d6d1350,g858d7b2824+d77d6d1350,g88963caddf+83e433e629,g99cad8db69+a4d3c48eeb,g9ddcbc5298+9a081db1e4,ga1e77700b3+bcf1af89ad,ga57fefb910+9a39d7b2d7,gae0086650b+585e252eca,gb065fddaf9+4f9fd82a2c,gb0e22166c9+60f28cb32d,gb363559e06+d84b1d3d07,gb3b7280ab2+4563d032e1,gb4b16eec92+babe958938,gba4ed39666+c2a2e4ac27,gbb8dafda3b+ed6854b564,gc120e1dc64+b72d212f87,gc28159a63d+3c60f8fe34,gc3e9b769f7+921dbcd359,gcf0d15dbbd+9a39d7b2d7,gdaeeff99f8+f9a426f77a,gddc38dedce+585e252eca,ge79ae78c31+3c60f8fe34,w.2024.21
LSST Data Management Base Package
Loading...
Searching...
No Matches
Classes | Variables
lsst.pipe.tasks.mergeMeasurements Namespace Reference

Classes

class  MergeMeasurementsConnections
 

Variables

 schema : `lsst.afw.table.Schema`, optional
 
 initInputs : `dict`, optional
 
 catalogs : `lsst.afw.table.SourceCatalog`
 

Variable Documentation

◆ catalogs

lsst.pipe.tasks.mergeMeasurements.catalogs : `lsst.afw.table.SourceCatalog`
_DefaultName = "mergeCoaddMeasurements"
ConfigClass = MergeMeasurementsConfig

inputDataset = "meas"
outputDataset = "ref"

def __init__(self, schema=None, initInputs=None, **kwargs):
    super().__init__(**kwargs)

    if initInputs is not None:
        schema = initInputs['inputSchema'].schema

    if schema is None:
        raise ValueError("No input schema or initInputs['inputSchema'] provided.")

    inputSchema = schema

    self.schemaMapper = afwTable.SchemaMapper(inputSchema, True)
    self.schemaMapper.addMinimalSchema(inputSchema, True)
    self.instFluxKey = inputSchema.find(self.config.snName + "_instFlux").getKey()
    self.instFluxErrKey = inputSchema.find(self.config.snName + "_instFluxErr").getKey()
    self.fluxFlagKey = inputSchema.find(self.config.snName + "_flag").getKey()

    self.flagKeys = {}
    for band in self.config.priorityList:
        outputKey = self.schemaMapper.editOutputSchema().addField(
            "merge_measurement_%s" % band,
            type="Flag",
            doc="Flag field set if the measurements here are from the %s filter" % band
        )
        peakKey = inputSchema.find("merge_peak_%s" % band).key
        footprintKey = inputSchema.find("merge_footprint_%s" % band).key
        self.flagKeys[band] = pipeBase.Struct(peak=peakKey, footprint=footprintKey, output=outputKey)
    self.schema = self.schemaMapper.getOutputSchema()

    self.pseudoFilterKeys = []
    for filt in self.config.pseudoFilterList:
        try:
            self.pseudoFilterKeys.append(self.schema.find("merge_peak_%s" % filt).getKey())
        except Exception as e:
            self.log.warning("merge_peak is not set for pseudo-filter %s: %s", filt, e)

    self.badFlags = {}
    for flag in self.config.flags:
        try:
            self.badFlags[flag] = self.schema.find(flag).getKey()
        except KeyError as exc:
            self.log.warning("Can't find flag %s in schema: %s", flag, exc)
    self.outputSchema = afwTable.SourceCatalog(self.schema)

def runQuantum(self, butlerQC, inputRefs, outputRefs):
    inputs = butlerQC.get(inputRefs)
    dataIds = (ref.dataId for ref in inputRefs.catalogs)
    catalogDict = {dataId['band']: cat for dataId, cat in zip(dataIds, inputs['catalogs'])}
    inputs['catalogs'] = catalogDict
    outputs = self.run(**inputs)
    butlerQC.put(outputs, outputRefs)

def run(self, catalogs):

Definition at line 201 of file mergeMeasurements.py.

◆ initInputs

lsst.pipe.tasks.mergeMeasurements.initInputs : `dict`, optional

Definition at line 131 of file mergeMeasurements.py.

◆ schema

lsst.pipe.tasks.mergeMeasurements.schema : `lsst.afw.table.Schema`, optional
pseudoFilterList = pexConfig.ListField(
    dtype=str,
    default=["sky"],
    doc="Names of filters which may have no associated detection\n"
    "(N.b. should include MergeDetectionsConfig.skyFilterName)"
)
snName = pexConfig.Field(
    dtype=str,
    default="base_PsfFlux",
    doc="Name of flux measurement for calculating the S/N when choosing the reference band."
)
minSN = pexConfig.Field(
    dtype=float,
    default=10.,
    doc="If the S/N from the priority band is below this value (and the S/N "
    "is larger than minSNDiff compared to the priority band), use the band with "
    "the largest S/N as the reference band."
)
minSNDiff = pexConfig.Field(
    dtype=float,
    default=3.,
    doc="If the difference in S/N between another band and the priority band is larger "
    "than this value (and the S/N in the priority band is less than minSN) "
    "use the band with the largest S/N as the reference band"
)
flags = pexConfig.ListField(
    dtype=str,
    doc="Require that these flags, if available, are not set",
    default=["base_PixelFlags_flag_interpolatedCenter", "base_PsfFlux_flag",
             "ext_photometryKron_KronFlux_flag", "modelfit_CModel_flag", ]
)
priorityList = pexConfig.ListField(
    dtype=str,
    default=[],
    doc="Priority-ordered list of filter bands for the merge."
)
coaddName = pexConfig.Field(
    dtype=str,
    default="deep",
    doc="Name of coadd"
)

def validate(self):
    super().validate()
    if len(self.priorityList) == 0:
        raise RuntimeError("No priority list provided")


class MergeMeasurementsTask(pipeBase.PipelineTask):

Definition at line 129 of file mergeMeasurements.py.