Photometry model with constraints, \(M(x,y) = M_CCD(x,y)*M_visit(u,v)\). More...

#include <ConstrainedPhotometryModel.h>

Inheritance diagram for lsst::jointcal::ConstrainedPhotometryModel:

Classes
struct	PrepPhotoCalib
	To hold the return of prepPhotoCalib. More...

Public Member Functions
	ConstrainedPhotometryModel (CcdImageList const &ccdImageList, geom::Box2D const &focalPlaneBBox, LOG_LOGGER log, int visitOrder=7, double errorPedestal=0)
	Construct a constrained photometry model. More...

	ConstrainedPhotometryModel (ConstrainedPhotometryModel const &)=delete
	No copy or move: there is only ever one instance of a given model (i.e. per ccd+visit) More...

	ConstrainedPhotometryModel (ConstrainedPhotometryModel &&)=delete

ConstrainedPhotometryModel &	operator= (ConstrainedPhotometryModel const &)=delete

ConstrainedPhotometryModel &	operator= (ConstrainedPhotometryModel &&)=delete

Eigen::Index	assignIndices (std::string const &whatToFit, Eigen::Index firstIndex) override
	Assign indices in the full matrix to the parameters being fit in the mappings, starting at firstIndex. More...

void	offsetParams (Eigen::VectorXd const &delta) override
	Offset the parameters by the provided amounts (by -delta). More...

void	freezeErrorTransform () override
	Once this routine has been called, the error transform is not modified by offsetParams(). More...

void	getMappingIndices (CcdImage const &ccdImage, IndexVector &indices) const override
	Get how this set of parameters (of length Npar()) map into the "grand" fit. More...

std::size_t	getTotalParameters () const override
	Return the total number of parameters in this model. More...

void	computeParameterDerivatives (MeasuredStar const &measuredStar, CcdImage const &ccdImage, Eigen::VectorXd &derivatives) const override
	Compute the parametric derivatives of this model. More...

void	print (std::ostream &out) const override
	Print a string representation of the contents of this mapping, for debugging. More...

	~ConstrainedPhotometryModel ()=default

virtual void	offsetFittedStar (FittedStar &fittedStar, double delta) const =0
	Offset the appropriate flux or magnitude (by -delta). More...

virtual double	computeResidual (CcdImage const &ccdImage, MeasuredStar const &measuredStar) const =0
	Compute the residual between the model applied to a star and its associated fittedStar. More...

virtual double	transform (CcdImage const &ccdImage, MeasuredStar const &measuredStar) const =0
	Return the on-sky transformed flux for measuredStar on ccdImage. More...

virtual double	transformError (CcdImage const &ccdImage, MeasuredStar const &measuredStar) const =0
	Return the on-sky transformed flux uncertainty for measuredStar on ccdImage. More...

virtual double	getRefError (RefStar const &refStar) const =0
	Return the refStar error appropriate for this model (e.g. fluxErr or magErr). More...

virtual double	computeRefResidual (FittedStar const &fittedStar, RefStar const &refStar) const =0
	Return the fittedStar - refStar residual appropriate for this model (e.g. flux - flux or mag - mag). More...

virtual std::shared_ptr< afw::image::PhotoCalib >	toPhotoCalib (CcdImage const &ccdImage) const =0
	Return the mapping of ccdImage represented as a PhotoCalib. More...

std::size_t	getNpar (CcdImage const &ccdImage) const
	Return the number of parameters in the mapping of CcdImage. More...

PhotometryMappingBase const &	getMapping (CcdImage const &ccdImage) const
	Get the mapping associated with ccdImage. More...

bool	validate (CcdImageList const &ccdImageList, int ndof) const
	Return true if this is a "reasonable" model. More...

bool	checkPositiveOnBBox (CcdImage const &ccdImage) const
	Check that the model is positive on the ccdImage bbox. More...

double	getErrorPedestal ()

double	tweakFluxError (jointcal::MeasuredStar const &measuredStar) const
	Add a fraction of the instrumental flux to the instrumental flux error, in quadrature. More...

double	tweakMagnitudeError (jointcal::MeasuredStar const &measuredStar) const
	Add a small magnitude offset to the "instrumental magnitude" error, in quadrature. More...

Protected Types
typedef std::unordered_map< CcdImageKey, std::unique_ptr< ChipVisitPhotometryMapping > >	MapType

typedef std::map< VisitIdType, std::shared_ptr< PhotometryMapping > >	VisitMapType

typedef std::map< CcdIdType, std::shared_ptr< PhotometryMapping > >	ChipMapType

Protected Member Functions
PhotometryMappingBase *	findMapping (CcdImage const &ccdImage) const override
	Return a pointer to the mapping associated with this ccdImage. More...

template<class ChipTransform , class VisitTransform , class ChipVisitMapping >
void	initialize (CcdImageList const &ccdImageList, geom::Box2D const &focalPlaneBBox, int visitOrder)
	Initialize the chip, visit, and chipVisit mappings by creating appropriate transforms and mappings. More...

virtual double	initialChipCalibration (std::shared_ptr< afw::image::PhotoCalib const > photoCalib)=0
	Return the initial calibration to use from this photoCalib. More...

PrepPhotoCalib	prepPhotoCalib (CcdImage const &ccdImage) const
	Helper for preparing toPhotoCalib() More...

Protected Attributes
MapType	_chipVisitMap

VisitMapType	_visitMap

ChipMapType	_chipMap

LOG_LOGGER	_log
	lsst.logging instance, to be created by a subclass so that messages have consistent name. More...

double	errorPedestal

Detailed Description

Photometry model with constraints, \(M(x,y) = M_CCD(x,y)*M_visit(u,v)\).

This model consists of the following components:

A spatially invariant zero point per CCD, constrained across all visits, \(M_CCD\).
A Chebyshev polynomial ( \(a_ij*T_i(x)*T_j(y)\) ) per visit, constrained across all CCDs, \(M_visit\).

Because this model's parameters are degenerate under multiplication by a constant, \(M=(a*M_CCD)*(1/a*M_visit)\), we hold one CCD's zero point fixed to remove that degeneracy.

Definition at line 48 of file ConstrainedPhotometryModel.h.

Member Typedef Documentation

◆ ChipMapType

typedef std::map<CcdIdType, std::shared_ptr<PhotometryMapping> > lsst::jointcal::ConstrainedPhotometryModel::ChipMapType

protected

Definition at line 110 of file ConstrainedPhotometryModel.h.

◆ MapType

typedef std::unordered_map<CcdImageKey, std::unique_ptr<ChipVisitPhotometryMapping> > lsst::jointcal::ConstrainedPhotometryModel::MapType

protected

Definition at line 103 of file ConstrainedPhotometryModel.h.

◆ VisitMapType

typedef std::map<VisitIdType, std::shared_ptr<PhotometryMapping> > lsst::jointcal::ConstrainedPhotometryModel::VisitMapType

protected

Definition at line 107 of file ConstrainedPhotometryModel.h.

Constructor & Destructor Documentation

◆ ConstrainedPhotometryModel() [1/3]

lsst::jointcal::ConstrainedPhotometryModel::ConstrainedPhotometryModel	(	CcdImageList const &	ccdImageList,
		geom::Box2D const &	focalPlaneBBox,
		LOG_LOGGER	log,
		int	visitOrder = `7`,
		double	errorPedestal = `0`
	)

inlineexplicit

Construct a constrained photometry model.

Parameters

	ccdImageList	The list of CCDImages to construct the model for.
	focalPlaneBBox	The bounding box of the camera's focal plane, defining the domain of the visit polynomial.
[in]	log	An lsst::log::Log instance to log messages to.
[in]	visitOrder	The order of the visit polynomial.
[in]	errorPedestal	A pedestal in flux or magnitude to apply to all MeasuredStar flux errors.

Definition at line 60 of file ConstrainedPhotometryModel.h.

             : PhotometryModel(log, errorPedestal), _fittingChips(false), _fittingVisits(false) {
         _chipVisitMap.reserve(ccdImageList.size());
     }

◆ ConstrainedPhotometryModel() [2/3]

lsst::jointcal::ConstrainedPhotometryModel::ConstrainedPhotometryModel ( ConstrainedPhotometryModel const & )

delete

No copy or move: there is only ever one instance of a given model (i.e. per ccd+visit)

◆ ConstrainedPhotometryModel() [3/3]

lsst::jointcal::ConstrainedPhotometryModel::ConstrainedPhotometryModel ( ConstrainedPhotometryModel && )

delete

◆ ~ConstrainedPhotometryModel()

lsst::jointcal::ConstrainedPhotometryModel::~ConstrainedPhotometryModel ( )

default

Member Function Documentation

◆ assignIndices()

Eigen::Index lsst::jointcal::ConstrainedPhotometryModel::assignIndices	(	std::string const &	whatToFit,
		Eigen::Index	firstIndex
	)

overridevirtual

Assign indices in the full matrix to the parameters being fit in the mappings, starting at firstIndex.

Parameters

[in]	whatToFit	String containing parameters to fit.
[in]	firstIndex	Index to start assigning at.

Returns: The highest assigned index.

Implements lsst::jointcal::PhotometryModel.

Definition at line 45 of file ConstrainedPhotometryModel.cc.

                                                                               {
     Eigen::Index index = firstIndex;
     if (whatToFit.find("Model") == std::string::npos) {
         LOGLS_WARN(_log, "assignIndices was called and Model is *not* in whatToFit");
         return index;
     }
  
     // If we got here, "Model" is definitely in whatToFit.
     _fittingChips = (whatToFit.find("ModelChip") != std::string::npos);
     _fittingVisits = (whatToFit.find("ModelVisit") != std::string::npos);
     // If nothing more than "Model" is specified, it means fit everything.
     if ((!_fittingChips) && (!_fittingVisits)) {
         _fittingChips = _fittingVisits = true;
     }
  
     if (_fittingChips) {
         for (auto &idMapping : _chipMap) {
             auto mapping = idMapping.second.get();
             // Don't assign indices for fixed parameters.
             if (mapping->isFixed()) continue;
             mapping->setIndex(index);
             index += mapping->getNpar();
         }
     }
     if (_fittingVisits) {
         for (auto &idMapping : _visitMap) {
             auto mapping = idMapping.second.get();
             mapping->setIndex(index);
             index += mapping->getNpar();
         }
     }
     for (auto &idMapping : _chipVisitMap) {
         idMapping.second->setWhatToFit(_fittingChips, _fittingVisits);
     }
     return index;
 }

◆ checkPositiveOnBBox()

bool lsst::jointcal::PhotometryModel::checkPositiveOnBBox ( CcdImage const & ccdImage ) const

inherited

Check that the model is positive on the ccdImage bbox.

Parameters

ccdImage The ccdImage to test.

Returns: True if the image is positive on a sampling of points of the ccdImage bbox.

Definition at line 49 of file PhotometryModel.cc.

                                                                         {
     bool check = true;
     auto bbox = ccdImage.getImageFrame();
     for (auto const &x : {bbox.xMin, bbox.getCenter().x, bbox.xMax})
         for (auto const &y : {bbox.yMin, bbox.getCenter().y, bbox.yMax}) {
             // flux, fluxErr, instFlux, instFluxErr all 1
             jointcal::MeasuredStar star(jointcal::BaseStar(x, y, 1, 1));
             star.setInstFluxAndErr(1, 1);
             double result = transform(ccdImage, star);
             // Don't short circuit so that we log every place the model is negative.
             if (result < 0) {
                 LOGLS_WARN(_log, "CcdImage " << ccdImage.getName() << " is negative at (" << x << "," << y
                                              << "): " << result);
                 check = false;
             }
         }
     return check;
 }

◆ computeParameterDerivatives()

void lsst::jointcal::ConstrainedPhotometryModel::computeParameterDerivatives	(	MeasuredStar const &	measuredStar,
		CcdImage const &	ccdImage,
		Eigen::VectorXd &	derivatives
	)		const

overridevirtual

Compute the parametric derivatives of this model.

Parameters

[in]	measuredStar	The measured star with the position and flux to compute at.
[in]	ccdImage	The ccdImage containing the measured star, to find the correct mapping.
[out]	derivatives	The computed derivatives. Must be pre-allocated to the correct size.

Implements lsst::jointcal::PhotometryModel.

Definition at line 125 of file ConstrainedPhotometryModel.cc.

                                                                                                {
     auto mapping = findMapping(ccdImage);
     mapping->computeParameterDerivatives(measuredStar, measuredStar.getInstFlux(), derivatives);
 }

◆ computeRefResidual()

virtual double lsst::jointcal::PhotometryModel::computeRefResidual	(	FittedStar const &	fittedStar,
		RefStar const &	refStar
	)		const

pure virtualinherited

Return the fittedStar - refStar residual appropriate for this model (e.g. flux - flux or mag - mag).

Implemented in lsst::jointcal::SimpleMagnitudeModel, lsst::jointcal::SimpleFluxModel, lsst::jointcal::ConstrainedMagnitudeModel, and lsst::jointcal::ConstrainedFluxModel.

◆ computeResidual()

virtual double lsst::jointcal::PhotometryModel::computeResidual	(	CcdImage const &	ccdImage,
		MeasuredStar const &	measuredStar
	)		const

pure virtualinherited

Compute the residual between the model applied to a star and its associated fittedStar.

\[ residual = Model(measuredStar) - fittedStar \]

Parameters

ccdImage	The ccdImage where measuredStar resides.
measuredStar	The measured star position to compute the residual of.

Returns: The residual.

Implemented in lsst::jointcal::SimpleMagnitudeModel, lsst::jointcal::SimpleFluxModel, lsst::jointcal::ConstrainedMagnitudeModel, and lsst::jointcal::ConstrainedFluxModel.

◆ findMapping()

PhotometryMappingBase * lsst::jointcal::ConstrainedPhotometryModel::findMapping ( CcdImage const & ccdImage ) const

overrideprotectedvirtual

Return a pointer to the mapping associated with this ccdImage.

Implements lsst::jointcal::PhotometryModel.

Definition at line 168 of file ConstrainedPhotometryModel.cc.

                                                                                              {
     auto idMapping = _chipVisitMap.find(ccdImage.getHashKey());
     if (idMapping == _chipVisitMap.end())
         throw LSST_EXCEPT(pex::exceptions::InvalidParameterError,
                           "ConstrainedPhotometryModel cannot find CcdImage " + ccdImage.getName());
     return idMapping->second.get();
 }

◆ freezeErrorTransform()

void lsst::jointcal::ConstrainedPhotometryModel::freezeErrorTransform ( )

overridevirtual

Once this routine has been called, the error transform is not modified by offsetParams().

The routine can be called when the mappings are roughly in place. After the call, the transformations used to propagate errors are no longer affected when updating the mappings. This allows an exactly linear fit, which can be necessary for some model+data combinations.

Implements lsst::jointcal::PhotometryModel.

Definition at line 100 of file ConstrainedPhotometryModel.cc.

                                                       {
     for (auto &idMapping : _chipMap) {
         idMapping.second.get()->freezeErrorTransform();
     }
     for (auto &idMapping : _visitMap) {
         idMapping.second.get()->freezeErrorTransform();
     }
 }

◆ getErrorPedestal()

double lsst::jointcal::PhotometryModel::getErrorPedestal ( )

inlineinherited

Definition at line 196 of file PhotometryModel.h.

196 { return errorPedestal; }

◆ getMapping()

PhotometryMappingBase const& lsst::jointcal::PhotometryModel::getMapping ( CcdImage const & ccdImage ) const

inlineinherited

Get the mapping associated with ccdImage.

Definition at line 156 of file PhotometryModel.h.

                                                                             {
         return *(findMapping(ccdImage));
     }

◆ getMappingIndices()

void lsst::jointcal::ConstrainedPhotometryModel::getMappingIndices	(	CcdImage const &	ccdImage,
		IndexVector &	indices
	)		const

overridevirtual

Get how this set of parameters (of length Npar()) map into the "grand" fit.

Parameters

[in]	ccdImage	The ccdImage to look up.
[out]	indices	The indices of the mapping associated with ccdImage.

Implements lsst::jointcal::PhotometryModel.

Definition at line 109 of file ConstrainedPhotometryModel.cc.

                                                                                                        {
     auto mapping = findMapping(ccdImage);
     mapping->getMappingIndices(indices);
 }

◆ getNpar()

std::size_t lsst::jointcal::PhotometryModel::getNpar ( CcdImage const & ccdImage ) const

inlineinherited

Return the number of parameters in the mapping of CcdImage.

Definition at line 153 of file PhotometryModel.h.

153 { return findMapping(ccdImage)->getNpar(); }

lsst::jointcal::PhotometryMappingBase::getNpar

virtual std::size_t getNpar() const =0

Number of total parameters in this mapping.

◆ getRefError()

virtual double lsst::jointcal::PhotometryModel::getRefError ( RefStar const & refStar ) const

pure virtualinherited

Return the refStar error appropriate for this model (e.g. fluxErr or magErr).

Implemented in lsst::jointcal::SimpleMagnitudeModel, lsst::jointcal::SimpleFluxModel, lsst::jointcal::ConstrainedMagnitudeModel, and lsst::jointcal::ConstrainedFluxModel.

◆ getTotalParameters()

std::size_t lsst::jointcal::ConstrainedPhotometryModel::getTotalParameters ( ) const

overridevirtual

Return the total number of parameters in this model.

Implements lsst::jointcal::PhotometryModel.

Definition at line 114 of file ConstrainedPhotometryModel.cc.

                                                                {
     std::size_t total = 0;
     for (auto &idMapping : _chipMap) {
         total += idMapping.second->getNpar();
     }
     for (auto &idMapping : _visitMap) {
         total += idMapping.second->getNpar();
     }
     return total;
 }

◆ initialChipCalibration()

virtual double lsst::jointcal::ConstrainedPhotometryModel::initialChipCalibration ( std::shared_ptr< afw::image::PhotoCalib const > photoCalib )

protectedpure virtual

Return the initial calibration to use from this photoCalib.

Implemented in lsst::jointcal::ConstrainedMagnitudeModel, and lsst::jointcal::ConstrainedFluxModel.

◆ initialize()

template<class ChipTransform , class VisitTransform , class ChipVisitMapping >

template void lsst::jointcal::ConstrainedPhotometryModel::initialize< MagnitudeTransformSpatiallyInvariant, MagnitudeTransformChebyshev, ChipVisitMagnitudeMapping >	(	CcdImageList const &	ccdImageList,
		geom::Box2D const &	focalPlaneBBox,
		int	visitOrder
	)

protected

Initialize the chip, visit, and chipVisit mappings by creating appropriate transforms and mappings.

Definition at line 177 of file ConstrainedPhotometryModel.cc.

                                                                                              {
     // keep track of which chip we want to constrain (the one closest to the middle of the focal plane)
     double minRadius2 = std::numeric_limits<double>::infinity();
     CcdIdType constrainedChip = -1;
  
     // First initialize all visit and ccd transforms, before we make the ccdImage mappings.
     for (auto const &ccdImage : ccdImageList) {
         auto visit = ccdImage->getVisit();
         auto chip = ccdImage->getCcdId();
         auto visitPair = _visitMap.find(visit);
         auto chipPair = _chipMap.find(chip);
  
         // If the chip is not in the map, add it, otherwise continue.
         if (chipPair == _chipMap.end()) {
             auto center = ccdImage->getDetector()->getCenter(afw::cameraGeom::FOCAL_PLANE);
             double radius2 = std::pow(center.getX(), 2) + std::pow(center.getY(), 2);
             if (radius2 < minRadius2) {
                 minRadius2 = radius2;
                 constrainedChip = chip;
             }
             auto photoCalib = ccdImage->getPhotoCalib();
             // Use the single-frame processing calibration from the PhotoCalib as the default.
             auto chipTransform = std::make_unique<ChipTransform>(initialChipCalibration(photoCalib));
             _chipMap[chip] = std::make_shared<PhotometryMapping>(std::move(chipTransform));
         }
         // If the visit is not in the map, add it, otherwise continue.
         if (visitPair == _visitMap.end()) {
             auto visitTransform = std::make_unique<VisitTransform>(visitOrder, focalPlaneBBox);
             _visitMap[visit] = std::make_shared<PhotometryMapping>(std::move(visitTransform));
         }
     }
  
     // Fix one chip mapping, to remove the degeneracy from the system.
     _chipMap.at(constrainedChip)->setFixed(true);
  
     // Now create the ccdImage mappings, which are combinations of the chip/visit mappings above.
     for (auto const &ccdImage : ccdImageList) {
         auto visit = ccdImage->getVisit();
         auto chip = ccdImage->getCcdId();
         _chipVisitMap.emplace(ccdImage->getHashKey(),
                               std::make_unique<ChipVisitMapping>(_chipMap[chip], _visitMap[visit]));
     }
     LOGLS_INFO(_log, "Got " << _chipMap.size() << " chip mappings and " << _visitMap.size()
                             << " visit mappings; holding chip " << constrainedChip << " fixed ("
                             << getTotalParameters() << " total parameters).");
     LOGLS_DEBUG(_log, "CcdImage map has " << _chipVisitMap.size() << " mappings, with "
                                           << _chipVisitMap.bucket_count() << " buckets and a load factor of "
                                           << _chipVisitMap.load_factor());
 }

◆ offsetFittedStar()

virtual void lsst::jointcal::PhotometryModel::offsetFittedStar	(	FittedStar &	fittedStar,
		double	delta
	)		const

pure virtualinherited

Offset the appropriate flux or magnitude (by -delta).

Parameters

fittedStar	The star to update.
delta	The amount to update by.

Implemented in lsst::jointcal::SimpleMagnitudeModel, lsst::jointcal::SimpleFluxModel, lsst::jointcal::ConstrainedMagnitudeModel, and lsst::jointcal::ConstrainedFluxModel.

◆ offsetParams()

void lsst::jointcal::ConstrainedPhotometryModel::offsetParams ( Eigen::VectorXd const & delta )

overridevirtual

Offset the parameters by the provided amounts (by -delta).

The shifts are applied according to the indices given in assignIndices.

Parameters

[in] delta vector of offsets to apply

Implements lsst::jointcal::PhotometryModel.

Definition at line 83 of file ConstrainedPhotometryModel.cc.

                                                                         {
     if (_fittingChips) {
         for (auto &idMapping : _chipMap) {
             auto mapping = idMapping.second.get();
             // Don't offset indices for fixed parameters.
             if (mapping->isFixed()) continue;
             mapping->offsetParams(delta.segment(mapping->getIndex(), mapping->getNpar()));
         }
     }
     if (_fittingVisits) {
         for (auto &idMapping : _visitMap) {
             auto mapping = idMapping.second.get();
             mapping->offsetParams(delta.segment(mapping->getIndex(), mapping->getNpar()));
         }
     }
 }

◆ operator=() [1/2]

ConstrainedPhotometryModel& lsst::jointcal::ConstrainedPhotometryModel::operator= ( ConstrainedPhotometryModel && )

delete

◆ operator=() [2/2]

ConstrainedPhotometryModel& lsst::jointcal::ConstrainedPhotometryModel::operator= ( ConstrainedPhotometryModel const & )

delete

◆ prepPhotoCalib()

ConstrainedPhotometryModel::PrepPhotoCalib lsst::jointcal::ConstrainedPhotometryModel::prepPhotoCalib ( CcdImage const & ccdImage ) const

protected

Helper for preparing toPhotoCalib()

Definition at line 228 of file ConstrainedPhotometryModel.cc.

                                         {
     auto detector = ccdImage.getDetector();
     auto ccdBBox = detector->getBBox();
     ChipVisitPhotometryMapping *mapping = dynamic_cast<ChipVisitPhotometryMapping *>(findMapping(ccdImage));
  
     // There should be no way in which we can get to this point and not have a ChipVisitMapping,
     // so blow up if we don't.
     assert(mapping != nullptr);
     // We know it's a Chebyshev transform because we created it as such, so blow up if it's not.
     auto visitPhotometryTransform = std::dynamic_pointer_cast<PhotometryTransformChebyshev>(
             mapping->getVisitMapping()->getTransform());
     assert(visitPhotometryTransform != nullptr);
     auto focalBBox = visitPhotometryTransform->getBBox();
  
     // Unravel our chebyshev coefficients to build an astshim::ChebyMap.
     auto coeff_f = toChebyMapCoeffs(std::dynamic_pointer_cast<PhotometryTransformChebyshev>(
             mapping->getVisitMapping()->getTransform()));
     // Bounds are the bbox
     std::vector<double> lowerBound = {focalBBox.getMinX(), focalBBox.getMinY()};
     std::vector<double> upperBound = {focalBBox.getMaxX(), focalBBox.getMaxY()};
     afw::geom::TransformPoint2ToGeneric visitTransform(ast::ChebyMap(coeff_f, 1, lowerBound, upperBound));
  
     double chipConstant = mapping->getChipMapping()->getParameters()[0];
  
     // Compute a box that covers the area of the ccd in focal plane coordinates.
     // This is the box over which we want to compute the mean of the visit transform.
     auto pixToFocal = detector->getTransform(afw::cameraGeom::PIXELS, afw::cameraGeom::FOCAL_PLANE);
     geom::Box2D ccdBBoxInFocal;
     for (auto const &point : pixToFocal->applyForward(geom::Box2D(ccdBBox).getCorners())) {
         ccdBBoxInFocal.include(point);
     }
     double visitMean = visitPhotometryTransform->mean(ccdBBoxInFocal);
  
     return {chipConstant, visitTransform, pixToFocal, visitMean};
 }

◆ print()

void lsst::jointcal::ConstrainedPhotometryModel::print ( std::ostream & out ) const

overridevirtual

Print a string representation of the contents of this mapping, for debugging.

This string representation can be very verbose, as it contains all of the parameters of all of the transforms in this model.

Implements lsst::jointcal::PhotometryModel.

Reimplemented in lsst::jointcal::ConstrainedMagnitudeModel, and lsst::jointcal::ConstrainedFluxModel.

Definition at line 154 of file ConstrainedPhotometryModel.cc.

                                                             {
     for (auto &idMapping : _chipMap) {
         out << "Sensor: " << idMapping.first << std::endl;
         idMapping.second->print(out);
         out << std::endl;
     }
     out << std::endl;
     for (auto &idMapping : _visitMap) {
         out << "Visit: " << idMapping.first << std::endl;
         idMapping.second->print(out);
         out << std::endl;
     }
 }

◆ toPhotoCalib()

virtual std::shared_ptr<afw::image::PhotoCalib> lsst::jointcal::PhotometryModel::toPhotoCalib ( CcdImage const & ccdImage ) const

pure virtualinherited

Return the mapping of ccdImage represented as a PhotoCalib.

Implemented in lsst::jointcal::SimpleMagnitudeModel, lsst::jointcal::SimpleFluxModel, lsst::jointcal::ConstrainedMagnitudeModel, and lsst::jointcal::ConstrainedFluxModel.

◆ transform()

virtual double lsst::jointcal::PhotometryModel::transform	(	CcdImage const &	ccdImage,
		MeasuredStar const &	measuredStar
	)		const

pure virtualinherited

Return the on-sky transformed flux for measuredStar on ccdImage.

Parameters

[in]	ccdImage	The ccdImage where measuredStar resides.
	measuredStar	The measured star position to transform.

Returns: The on-sky flux transformed from instFlux at measuredStar's position.

Implemented in lsst::jointcal::SimpleMagnitudeModel, lsst::jointcal::SimpleFluxModel, lsst::jointcal::ConstrainedMagnitudeModel, and lsst::jointcal::ConstrainedFluxModel.

◆ transformError()

virtual double lsst::jointcal::PhotometryModel::transformError	(	CcdImage const &	ccdImage,
		MeasuredStar const &	measuredStar
	)		const

pure virtualinherited

Return the on-sky transformed flux uncertainty for measuredStar on ccdImage.

Identical to transform() until freezeErrorTransform() is called.

Parameters

[in]	ccdImage	The ccdImage where measuredStar resides.
	measuredStar	The measured star position to transform.

Returns: The on-sky flux transformed from instFlux at measuredStar's position.

Implemented in lsst::jointcal::SimpleMagnitudeModel, lsst::jointcal::SimpleFluxModel, lsst::jointcal::ConstrainedMagnitudeModel, and lsst::jointcal::ConstrainedFluxModel.

◆ tweakFluxError()

double lsst::jointcal::PhotometryModel::tweakFluxError ( jointcal::MeasuredStar const & measuredStar ) const

inlineinherited

Add a fraction of the instrumental flux to the instrumental flux error, in quadrature.

Definition at line 199 of file PhotometryModel.h.

                                                                           {
         if (errorPedestal == 0) {
             return measuredStar.getInstFluxErr();
         } else {
             return std::hypot(measuredStar.getInstFluxErr(), measuredStar.getInstFlux() * errorPedestal);
         }
     }

◆ tweakMagnitudeError()

double lsst::jointcal::PhotometryModel::tweakMagnitudeError ( jointcal::MeasuredStar const & measuredStar ) const

inlineinherited

Add a small magnitude offset to the "instrumental magnitude" error, in quadrature.

Definition at line 208 of file PhotometryModel.h.

                                                                                {
         if (errorPedestal == 0) {
             return measuredStar.getInstMagErr();
         } else {
             return std::hypot(measuredStar.getInstMagErr(), errorPedestal);
         }
     }

◆ validate()

bool lsst::jointcal::PhotometryModel::validate	(	CcdImageList const &	ccdImageList,
		int	ndof
	)		const

inherited

Return true if this is a "reasonable" model.

A valid photometry model is positive within each sensor's bounding box.

Parameters

ccdImageList	The ccdImages to test the model validity on.
ndof	The number of degrees of freedom in the fit, e.g. from Fitterbase.computeChi2().

Returns: True if the model is valid on all ccdImages.

Definition at line 33 of file PhotometryModel.cc.

                                                                                {
     bool check = true;
     for (auto const &ccdImage : ccdImageList) {
         // Don't short circuit so that we log every place the model is negative.
         check &= checkPositiveOnBBox(*ccdImage);
     }
     if (ndof < 0) {
         check &= false;
         LOGLS_ERROR(_log, "This model only has "
                                   << ndof << " degrees of freedom, with " << getTotalParameters()
                                   << " total parameters. Reduce the model complexity (e.g. polynomial order)"
                                      " to better match the number of measured sources.");
     }
     return check;
 }

Member Data Documentation

◆ _chipMap

ChipMapType lsst::jointcal::ConstrainedPhotometryModel::_chipMap

protected

Definition at line 111 of file ConstrainedPhotometryModel.h.

◆ _chipVisitMap

MapType lsst::jointcal::ConstrainedPhotometryModel::_chipVisitMap

protected

Definition at line 104 of file ConstrainedPhotometryModel.h.

◆ _log

LOG_LOGGER lsst::jointcal::PhotometryModel::_log

protectedinherited

lsst.logging instance, to be created by a subclass so that messages have consistent name.

Definition at line 223 of file PhotometryModel.h.

◆ _visitMap

VisitMapType lsst::jointcal::ConstrainedPhotometryModel::_visitMap

protected

Definition at line 108 of file ConstrainedPhotometryModel.h.

◆ errorPedestal

double lsst::jointcal::PhotometryModel::errorPedestal

protectedinherited

Definition at line 226 of file PhotometryModel.h.

The documentation for this class was generated from the following files:

/j/snowflake/release/lsstsw/stack/lsst-scipipe-0.7.0/Linux64/jointcal/22.0.1-25-g5bf6245+3ad8ecd50b/include/lsst/jointcal/ConstrainedPhotometryModel.h
/j/snowflake/release/lsstsw/stack/lsst-scipipe-0.7.0/Linux64/jointcal/22.0.1-25-g5bf6245+3ad8ecd50b/src/ConstrainedPhotometryModel.cc

Classes

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Detailed Description

Member Typedef Documentation

◆ ChipMapType

◆ MapType

◆ VisitMapType

Constructor & Destructor Documentation

◆ ConstrainedPhotometryModel() [1/3]

◆ ConstrainedPhotometryModel() [2/3]

◆ ConstrainedPhotometryModel() [3/3]

◆ ~ConstrainedPhotometryModel()

Member Function Documentation

◆ assignIndices()

◆ checkPositiveOnBBox()

◆ computeParameterDerivatives()

◆ computeRefResidual()

◆ computeResidual()

◆ findMapping()

◆ freezeErrorTransform()

◆ getErrorPedestal()

◆ getMapping()

◆ getMappingIndices()

◆ getNpar()

◆ getRefError()

◆ getTotalParameters()

◆ initialChipCalibration()

◆ initialize()

◆ offsetFittedStar()

◆ offsetParams()

◆ operator=() [1/2]

◆ operator=() [2/2]

◆ prepPhotoCalib()

◆ print()

◆ toPhotoCalib()

◆ transform()

◆ transformError()

◆ tweakFluxError()

◆ tweakMagnitudeError()

◆ validate()

Member Data Documentation

◆ _chipMap

◆ _chipVisitMap

◆ _log

◆ _visitMap

◆ errorPedestal