#include <PhotometryMapping.h>

Inheritance diagram for lsst::jointcal::ChipVisitMagnitudeMapping:

Public Member Functions
	ChipVisitMagnitudeMapping (std::shared_ptr< PhotometryMapping > chipMapping, std::shared_ptr< PhotometryMapping > visitMapping)

double	transformError (MeasuredStar const &measuredStar, double value, double valueErr) const override
	Return the on-sky transformed flux uncertainty for measuredStar on ccdImage. More...

void	computeParameterDerivatives (MeasuredStar const &measuredStar, double value, Eigen::Ref< Eigen::VectorXd > derivatives) const override
	Compute the derivatives with respect to the parameters (i.e. More...

std::size_t	getNpar () const override
	Number of total parameters in this mapping. More...

double	transform (MeasuredStar const &measuredStar, double value) const override
	Return the on-sky transformed flux for measuredStar on ccdImage. More...

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

Eigen::VectorXd	getParameters () override

void	getMappingIndices (IndexVector &indices) const override
	Gets how this set of parameters (of length getNpar()) map into the "grand" fit. More...

void	setWhatToFit (bool const fittingChips, bool const fittingVisits)
	Set whether to fit chips or visits. More...

void	dump (std::ostream &stream=std::cout) const override
	Dump the contents of the transforms, for debugging. More...

std::shared_ptr< PhotometryMapping >	getChipMapping () const

std::shared_ptr< PhotometryMapping >	getVisitMapping () const

std::size_t	getNParChip () const

std::size_t	getNParVisit () const

void	setFixed (bool _fixed)
	Make this mapping's parameters fixed (i.e. not varied during fitting). More...

bool	isFixed ()

Eigen::Index	getIndex ()
	Get the index of this mapping in the grand fit. More...

void	setIndex (Eigen::Index i)
	Set the index of this mapping in the grand fit. More...

Protected Attributes
std::size_t	_nParChip

std::size_t	_nParVisit

std::shared_ptr< PhotometryMapping >	_chipMapping

std::shared_ptr< PhotometryMapping >	_visitMapping

Eigen::Index	index

bool	fixed

Detailed Description

Definition at line 297 of file PhotometryMapping.h.

Constructor & Destructor Documentation

◆ ChipVisitMagnitudeMapping()

lsst::jointcal::ChipVisitMagnitudeMapping::ChipVisitMagnitudeMapping	(	std::shared_ptr< PhotometryMapping >	chipMapping,
		std::shared_ptr< PhotometryMapping >	visitMapping
	)

inline

Definition at line 299 of file PhotometryMapping.h.

301 : ChipVisitPhotometryMapping(chipMapping, visitMapping) {}

lsst::jointcal::ChipVisitPhotometryMapping::ChipVisitPhotometryMapping

ChipVisitPhotometryMapping(std::shared_ptr< PhotometryMapping > chipMapping, std::shared_ptr< PhotometryMapping > visitMapping)

Definition: PhotometryMapping.h:215

Member Function Documentation

◆ computeParameterDerivatives()

void lsst::jointcal::ChipVisitMagnitudeMapping::computeParameterDerivatives	(	MeasuredStar const &	measuredStar,
		double	value,
		Eigen::Ref< Eigen::VectorXd >	derivatives
	)		const

overridevirtual

Compute the derivatives with respect to the parameters (i.e.

the coefficients).

Parameters

[in]	measuredStar	The measured star position to transform.
[in]	value	The instrument flux or magnitude to compute the derivative at.
[out]	derivatives	The computed derivatives, in the same order as the deltas in offsetParams.

Implements lsst::jointcal::PhotometryMappingBase.

Definition at line 119 of file PhotometryMapping.cc.

                                                                                                        {
     // TODO DM-12161: possible optimization is to merge transform and computeDerivatives,
     // and/or save these intermediate calculations when transforming flux to use in derivatives.
     // Like what AstrometryMappings do with `computeTransformAndDerivatives` vs. `transformPosAndErrors`.
 
     // NOTE: See DMTN-036 for the math behind this.
     if (getNParChip() > 0 && !_chipMapping->isFixed()) {
         // The chip derivatives start at 0, independent of the full-fit indices.
         Eigen::Ref<Eigen::VectorXd> chipBlock = derivatives.segment(0, getNParChip());
         _chipMapping->getTransform()->computeParameterDerivatives(measuredStar.x, measuredStar.y, instFlux,
                                                                   chipBlock);
     }
     if (getNParVisit() > 0) {
         // The visit derivatives start at the last chip derivative, independent of the full-fit indices.
         Eigen::Ref<Eigen::VectorXd> visitBlock = derivatives.segment(getNParChip(), getNParVisit());
         _visitMapping->getTransform()->computeParameterDerivatives(
                 measuredStar.getXFocal(), measuredStar.getYFocal(), instFlux, visitBlock);
     }
 }

◆ dump()

void lsst::jointcal::ChipVisitPhotometryMapping::dump ( std::ostream & stream = std::cout ) const

inlineoverridevirtualinherited

Dump the contents of the transforms, for debugging.

Implements lsst::jointcal::PhotometryMappingBase.

Definition at line 261 of file PhotometryMapping.h.

                                                            {
         stream << "index: " << index << " chipMapping: ";
         _chipMapping->dump(stream);
         stream << "visitMapping: ";
         _visitMapping->dump(stream);
     }

◆ freezeErrorTransform()

void lsst::jointcal::ChipVisitPhotometryMapping::freezeErrorTransform ( )

inlineoverridevirtualinherited

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::PhotometryMappingBase.

Definition at line 234 of file PhotometryMapping.h.

                                          {
         _chipMapping->freezeErrorTransform();
         _visitMapping->freezeErrorTransform();
     }

◆ getChipMapping()

std::shared_ptr<PhotometryMapping> lsst::jointcal::ChipVisitPhotometryMapping::getChipMapping ( ) const

inlineinherited

Definition at line 268 of file PhotometryMapping.h.

268 { return _chipMapping; }

lsst::jointcal::ChipVisitPhotometryMapping::_chipMapping

std::shared_ptr< PhotometryMapping > _chipMapping

Definition: PhotometryMapping.h:279

◆ getIndex()

Eigen::Index lsst::jointcal::PhotometryMappingBase::getIndex ( )

inlineinherited

Get the index of this mapping in the grand fit.

Definition at line 115 of file PhotometryMapping.h.

115 { return index; }

lsst::jointcal::PhotometryMappingBase::index

Eigen::Index index

Definition: PhotometryMapping.h:122

◆ getMappingIndices()

void lsst::jointcal::ChipVisitPhotometryMapping::getMappingIndices ( IndexVector & indices ) const

overridevirtualinherited

Gets how this set of parameters (of length getNpar()) map into the "grand" fit.

Expects that indices has enough space reserved.

Implements lsst::jointcal::PhotometryMappingBase.

Definition at line 39 of file PhotometryMapping.cc.

                                                                              {
     if (indices.size() < getNpar()) indices.resize(getNpar());
     // If we're fitting the chip mapping, fill those indices.
     if (_nParChip > 0) {
         _chipMapping->getMappingIndices(indices);
     }
     // If we're fitting the visit mapping, fill those indices.
     if (_nParVisit > 0) {
         // TODO DM-12169: there is probably a better way to feed a subpart of a std::vector
         // (maybe a view or iterators?)
         IndexVector tempIndices(_visitMapping->getNpar());
         _visitMapping->getMappingIndices(tempIndices);
         // We have to insert the visit indices starting after the chip indices.
         for (std::size_t k = 0; k < _visitMapping->getNpar(); ++k) {
             indices.at(k + _nParChip) = tempIndices.at(k);
         }
     }
 }

◆ getNpar()

std::size_t lsst::jointcal::ChipVisitPhotometryMapping::getNpar ( ) const

inlineoverridevirtualinherited

Number of total parameters in this mapping.

Implements lsst::jointcal::PhotometryMappingBase.

Definition at line 224 of file PhotometryMapping.h.

224 { return _nParChip + _nParVisit; }

lsst::jointcal::ChipVisitPhotometryMapping::_nParVisit

std::size_t _nParVisit

Definition: PhotometryMapping.h:276

lsst::jointcal::ChipVisitPhotometryMapping::_nParChip

std::size_t _nParChip

Definition: PhotometryMapping.h:276

◆ getNParChip()

std::size_t lsst::jointcal::ChipVisitPhotometryMapping::getNParChip ( ) const

inlineinherited

Definition at line 271 of file PhotometryMapping.h.

271 { return _nParChip; }

lsst::jointcal::ChipVisitPhotometryMapping::_nParChip

std::size_t _nParChip

Definition: PhotometryMapping.h:276

◆ getNParVisit()

std::size_t lsst::jointcal::ChipVisitPhotometryMapping::getNParVisit ( ) const

inlineinherited

Definition at line 272 of file PhotometryMapping.h.

272 { return _nParVisit; }

lsst::jointcal::ChipVisitPhotometryMapping::_nParVisit

std::size_t _nParVisit

Definition: PhotometryMapping.h:276

◆ getParameters()

Eigen::VectorXd lsst::jointcal::ChipVisitPhotometryMapping::getParameters ( )

inlineoverridevirtualinherited

Implements lsst::jointcal::PhotometryMappingBase.

Definition at line 240 of file PhotometryMapping.h.

                                            {
         Eigen::VectorXd joined(getNpar());
         joined << _chipMapping->getParameters(), _visitMapping->getParameters();
         return joined;
     }

◆ getVisitMapping()

std::shared_ptr<PhotometryMapping> lsst::jointcal::ChipVisitPhotometryMapping::getVisitMapping ( ) const

inlineinherited

Definition at line 269 of file PhotometryMapping.h.

269 { return _visitMapping; }

lsst::jointcal::ChipVisitPhotometryMapping::_visitMapping

std::shared_ptr< PhotometryMapping > _visitMapping

Definition: PhotometryMapping.h:280

◆ isFixed()

bool lsst::jointcal::PhotometryMappingBase::isFixed ( )

inlineinherited

Definition at line 101 of file PhotometryMapping.h.

101 { return fixed; }

lsst::jointcal::PhotometryMappingBase::fixed

bool fixed

Definition: PhotometryMapping.h:124

◆ setFixed()

void lsst::jointcal::PhotometryMappingBase::setFixed ( bool _fixed )

inlineinherited

Make this mapping's parameters fixed (i.e. not varied during fitting).

Definition at line 100 of file PhotometryMapping.h.

100 { fixed = _fixed; }

lsst::jointcal::PhotometryMappingBase::fixed

bool fixed

Definition: PhotometryMapping.h:124

◆ setIndex()

void lsst::jointcal::PhotometryMappingBase::setIndex ( Eigen::Index i )

inlineinherited

Set the index of this mapping in the grand fit.

Definition at line 118 of file PhotometryMapping.h.

118 { index = i; }

lsst::jointcal::PhotometryMappingBase::index

Eigen::Index index

Definition: PhotometryMapping.h:122

◆ setWhatToFit()

void lsst::jointcal::ChipVisitPhotometryMapping::setWhatToFit	(	bool const	fittingChips,
		bool const	fittingVisits
	)

inherited

Set whether to fit chips or visits.

This must be called before anything that depends on knowing the number of parameters in the fit, such as offsetParams(), getParameters(), or computeParameterDerivatives().

Parameters

fittingChips	Fit the chip transform.
fittingVisits	Fit the visit transform.

Definition at line 58 of file PhotometryMapping.cc.

                                                                                                {
     if (fittingChips) {
         _nParChip = _chipMapping->getNpar();
     } else {
         _nParChip = 0;
     }
     if (fittingVisits) {
         _nParVisit = _visitMapping->getNpar();
     } else {
         _nParVisit = 0;
     }
 }

◆ transform()

double lsst::jointcal::ChipVisitPhotometryMapping::transform	(	MeasuredStar const &	measuredStar,
		double	value
	)		const

inlineoverridevirtualinherited

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

Parameters

[in]	measuredStar	The measured star position to transform.
[in]	value	The instrument flux or magnitude to transform.

Returns: The on-sky value () transformed from value at measuredStar's position.

Implements lsst::jointcal::PhotometryMappingBase.

Definition at line 227 of file PhotometryMapping.h.

                                                                                     {
         double temp = _chipMapping->getTransform()->transform(measuredStar.x, measuredStar.y, value);
         return _visitMapping->getTransform()->transform(measuredStar.getXFocal(), measuredStar.getYFocal(),
                                                         temp);
     }

◆ transformError()

double lsst::jointcal::ChipVisitMagnitudeMapping::transformError	(	MeasuredStar const &	measuredStar,
		double	value,
		double	valueErr
	)		const

overridevirtual

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

Matches the underlying PhotometryTransform's transformError() until freezeErrorTransform() is called.

Parameters

[in]	measuredStar	The measured star position to transform.
[in]	value	The flux or magnitude to transform.
[in]	valueErr	The flux or magnitude uncertainty to transform.

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

Note: This method takes instFlux and instFluxErr: the error calculation has to use fluxes to get the math right.

Implements lsst::jointcal::PhotometryMappingBase.

Definition at line 112 of file PhotometryMapping.cc.

                                                                            {
     // The transformed error is s_mout = 2.5/ln(10) * instFluxErr / instFlux
     // because the other components of the mapping (f0, the polynomials) disappear in the partial derivative.
     return 2.5 / std::log(10.0) * instFluxErr / instFlux;
 }

Member Data Documentation

◆ _chipMapping

std::shared_ptr<PhotometryMapping> lsst::jointcal::ChipVisitPhotometryMapping::_chipMapping

protectedinherited

Definition at line 279 of file PhotometryMapping.h.

◆ _nParChip

std::size_t lsst::jointcal::ChipVisitPhotometryMapping::_nParChip

protectedinherited

Definition at line 276 of file PhotometryMapping.h.

◆ _nParVisit

std::size_t lsst::jointcal::ChipVisitPhotometryMapping::_nParVisit

protectedinherited

Definition at line 276 of file PhotometryMapping.h.

◆ _visitMapping

std::shared_ptr<PhotometryMapping> lsst::jointcal::ChipVisitPhotometryMapping::_visitMapping

protectedinherited

Definition at line 280 of file PhotometryMapping.h.

◆ fixed

bool lsst::jointcal::PhotometryMappingBase::fixed

protectedinherited

Definition at line 124 of file PhotometryMapping.h.

◆ index

Eigen::Index lsst::jointcal::PhotometryMappingBase::index

protectedinherited

Definition at line 122 of file PhotometryMapping.h.

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

/j/snowflake/release/lsstsw/stack/Linux64/jointcal/19.0.0-2-gb96a1c4+3/include/lsst/jointcal/PhotometryMapping.h
/j/snowflake/release/lsstsw/stack/Linux64/jointcal/19.0.0-2-gb96a1c4+3/src/PhotometryMapping.cc

Public Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ChipVisitMagnitudeMapping()

Member Function Documentation

◆ computeParameterDerivatives()

◆ dump()

◆ freezeErrorTransform()

◆ getChipMapping()

◆ getIndex()

◆ getMappingIndices()

◆ getNpar()

◆ getNParChip()

◆ getNParVisit()

◆ getParameters()

◆ getVisitMapping()

◆ isFixed()

◆ setFixed()

◆ setIndex()

◆ setWhatToFit()

◆ transform()

◆ transformError()

Member Data Documentation

◆ _chipMapping

◆ _nParChip

◆ _nParVisit

◆ _visitMapping

◆ fixed

◆ index