LSST Applications g0f08755f38+c89d42e150,g1635faa6d4+b6cf076a36,g1653933729+a8ce1bb630,g1a0ca8cf93+4c08b13bf7,g28da252d5a+f33f8200ef,g29321ee8c0+0187be18b1,g2bbee38e9b+9634bc57db,g2bc492864f+9634bc57db,g2cdde0e794+c2c89b37c4,g3156d2b45e+41e33cbcdc,g347aa1857d+9634bc57db,g35bb328faa+a8ce1bb630,g3a166c0a6a+9634bc57db,g3e281a1b8c+9f2c4e2fc3,g414038480c+077ccc18e7,g41af890bb2+e740673f1a,g5fbc88fb19+17cd334064,g7642f7d749+c89d42e150,g781aacb6e4+a8ce1bb630,g80478fca09+f8b2ab54e1,g82479be7b0+e2bd23ab8b,g858d7b2824+c89d42e150,g9125e01d80+a8ce1bb630,g9726552aa6+10f999ec6a,ga5288a1d22+065360aec4,gacf8899fa4+9553554aa7,gae0086650b+a8ce1bb630,gb58c049af0+d64f4d3760,gbd46683f8f+ac57cbb13d,gc28159a63d+9634bc57db,gcf0d15dbbd+e37acf7834,gda3e153d99+c89d42e150,gda6a2b7d83+e37acf7834,gdaeeff99f8+1711a396fd,ge2409df99d+cb1e6652d6,ge79ae78c31+9634bc57db,gf0baf85859+147a0692ba,gf3967379c6+02b11634a5,w.2024.45
LSST Data Management Base Package
Loading...
Searching...
No Matches
Public Member Functions | Protected Attributes | List of all members
lsst::jointcal::SimpleAstrometryMapping Class Reference

Class for a simple mapping implementing a generic AstrometryTransform. More...

#include <SimpleAstrometryMapping.h>

Inheritance diagram for lsst::jointcal::SimpleAstrometryMapping:
lsst::jointcal::AstrometryMapping lsst::jointcal::SimplePolyMapping

Public Member Functions

 SimpleAstrometryMapping (AstrometryTransform const &astrometryTransform, bool toBeFit=true)
 
 SimpleAstrometryMapping (SimpleAstrometryMapping const &)=delete
 No copy or move: there is only ever one instance of a given mapping (i.e.. per ccd+visit)
 
 SimpleAstrometryMapping (SimpleAstrometryMapping &&)=delete
 
SimpleAstrometryMappingoperator= (SimpleAstrometryMapping const &)=delete
 
SimpleAstrometryMappingoperator= (SimpleAstrometryMapping &&)=delete
 
virtual void freezeErrorTransform ()
 
std::size_t getNpar () const override
 Number of parameters in total.
 
void getMappingIndices (IndexVector &indices) const override
 Sets how this set of parameters (of length Npar()) map into the "grand" fit Expects that indices has enough space reserved.
 
void transformPosAndErrors (FatPoint const &where, FatPoint &outPoint) const override
 The same as above but without the parameter derivatives (used to evaluate chi^2)
 
void positionDerivative (Point const &where, Eigen::Matrix2d &derivative, double epsilon) const override
 The derivative w.r.t. position.
 
void offsetParams (Eigen::VectorXd const &delta) override
 Remember the error scale and freeze it.
 
Eigen::Index getIndex () const
 position of the parameters within the grand fitting scheme
 
void setIndex (Eigen::Index i)
 Set the index of this mapping in the grand fit.
 
void computeTransformAndDerivatives (FatPoint const &where, FatPoint &outPoint, Eigen::MatrixX2d &H) const override
 Actually applies the AstrometryMapping and evaluates the derivatives w.r.t the fitted parameters.
 
virtual AstrometryTransform const & getTransform () const
 Access to the (fitted) transform.
 
bool getToBeFit () const
 Get whether this mapping is fit as part of a Model.
 
void setToBeFit (bool value)
 Set whether this Mapping is to be fit as part of a Model.
 
void print (std::ostream &out) const override
 Print a string representation of the contents of this mapping, for debugging.
 

Protected Attributes

bool toBeFit
 
Eigen::Index index {}
 
std::shared_ptr< AstrometryTransformtransform
 
std::shared_ptr< AstrometryTransformerrorProp
 
std::unique_ptr< AstrometryTransformLinearlin
 

Detailed Description

Class for a simple mapping implementing a generic AstrometryTransform.

It uses a template rather than a pointer so that the derived classes can use the specifics of the transform. The class simplePolyMapping overloads a few routines.

Definition at line 42 of file SimpleAstrometryMapping.h.

Constructor & Destructor Documentation

◆ SimpleAstrometryMapping() [1/3]

lsst::jointcal::SimpleAstrometryMapping::SimpleAstrometryMapping ( AstrometryTransform const & astrometryTransform,
bool toBeFit = true )
inline

Definition at line 44 of file SimpleAstrometryMapping.h.

46 transform(astrometryTransform.clone()),
47 errorProp(transform),
48 lin(new AstrometryTransformLinear) {}
std::shared_ptr< AstrometryTransform > errorProp
std::shared_ptr< AstrometryTransform > transform
std::unique_ptr< AstrometryTransformLinear > lin

◆ SimpleAstrometryMapping() [2/3]

lsst::jointcal::SimpleAstrometryMapping::SimpleAstrometryMapping ( SimpleAstrometryMapping const & )
delete

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

◆ SimpleAstrometryMapping() [3/3]

lsst::jointcal::SimpleAstrometryMapping::SimpleAstrometryMapping ( SimpleAstrometryMapping && )
delete

Member Function Documentation

◆ computeTransformAndDerivatives()

void lsst::jointcal::SimpleAstrometryMapping::computeTransformAndDerivatives ( FatPoint const & where,
FatPoint & outPoint,
Eigen::MatrixX2d & H ) const
overridevirtual

Actually applies the AstrometryMapping and evaluates the derivatives w.r.t the fitted parameters.

This is grouped into a single call because for most models, evaluating the derivatives w.r.T parameters is not much longer than just transforming

Implements lsst::jointcal::AstrometryMapping.

Reimplemented in lsst::jointcal::SimplePolyMapping.

Definition at line 67 of file SimpleAstrometryMapping.cc.

68 {
69 transformPosAndErrors(where, outPoint);
70 transform->paramDerivatives(where, &H(0, 0), &H(0, 1));
71}
void transformPosAndErrors(FatPoint const &where, FatPoint &outPoint) const override
The same as above but without the parameter derivatives (used to evaluate chi^2)

◆ freezeErrorTransform()

virtual void lsst::jointcal::SimpleAstrometryMapping::freezeErrorTransform ( )
inlinevirtual

Definition at line 56 of file SimpleAstrometryMapping.h.

56 {
57 // from there on, updating the transform does not change the errors.
58 errorProp = transform->clone();
59 }

◆ getIndex()

Eigen::Index lsst::jointcal::SimpleAstrometryMapping::getIndex ( ) const
inline

position of the parameters within the grand fitting scheme

Definition at line 84 of file SimpleAstrometryMapping.h.

◆ getMappingIndices()

void lsst::jointcal::SimpleAstrometryMapping::getMappingIndices ( IndexVector & indices) const
overridevirtual

Sets how this set of parameters (of length Npar()) map into the "grand" fit Expects that indices has enough space reserved.

Implements lsst::jointcal::AstrometryMapping.

Definition at line 34 of file SimpleAstrometryMapping.cc.

34 {
35 if (indices.size() < getNpar()) {
36 indices.resize(getNpar());
37 }
38 for (std::size_t k = 0; k < getNpar(); ++k) {
39 indices[k] = index + k;
40 }
41}
std::size_t getNpar() const override
Number of parameters in total.
T resize(T... args)
T size(T... args)

◆ getNpar()

std::size_t lsst::jointcal::SimpleAstrometryMapping::getNpar ( ) const
inlineoverridevirtual

Number of parameters in total.

Implements lsst::jointcal::AstrometryMapping.

Definition at line 62 of file SimpleAstrometryMapping.h.

62 {
63 if (toBeFit)
64 return transform->getNpar();
65 else
66 return 0;
67 }

◆ getToBeFit()

bool lsst::jointcal::SimpleAstrometryMapping::getToBeFit ( ) const
inline

Get whether this mapping is fit as part of a Model.

Definition at line 97 of file SimpleAstrometryMapping.h.

97{ return toBeFit; }

◆ getTransform()

virtual AstrometryTransform const & lsst::jointcal::SimpleAstrometryMapping::getTransform ( ) const
inlinevirtual

Access to the (fitted) transform.

Reimplemented in lsst::jointcal::SimplePolyMapping.

Definition at line 94 of file SimpleAstrometryMapping.h.

94{ return *transform; }

◆ offsetParams()

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

Remember the error scale and freeze it.

Implements lsst::jointcal::AstrometryMapping.

Definition at line 79 of file SimpleAstrometryMapping.h.

79 {
80 if (toBeFit) transform->offsetParams(delta);
81 }

◆ operator=() [1/2]

SimpleAstrometryMapping & lsst::jointcal::SimpleAstrometryMapping::operator= ( SimpleAstrometryMapping && )
delete

◆ operator=() [2/2]

SimpleAstrometryMapping & lsst::jointcal::SimpleAstrometryMapping::operator= ( SimpleAstrometryMapping const & )
delete

◆ positionDerivative()

void lsst::jointcal::SimpleAstrometryMapping::positionDerivative ( Point const & where,
Eigen::Matrix2d & derivative,
double epsilon ) const
overridevirtual

The derivative w.r.t. position.

Implements lsst::jointcal::AstrometryMapping.

Reimplemented in lsst::jointcal::SimplePolyMapping.

Definition at line 52 of file SimpleAstrometryMapping.cc.

53 {
54 errorProp->computeDerivative(where, *lin, epsilon);
55 derivative(0, 0) = lin->getCoefficient(1, 0, 0);
56 //
57 /* This does not work : it was proved by rotating the frame
58 see the compilation switch ROTATE_T2 in constrainedAstrometryModel.cc
59 derivative(1,0) = lin->getCoefficient(1,0,1);
60 derivative(0,1) = lin->getCoefficient(0,1,0);
61 */
62 derivative(1, 0) = lin->getCoefficient(0, 1, 0);
63 derivative(0, 1) = lin->getCoefficient(1, 0, 1);
64 derivative(1, 1) = lin->getCoefficient(0, 1, 1);
65}

◆ print()

void lsst::jointcal::SimpleAstrometryMapping::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 mapping.

Implements lsst::jointcal::AstrometryMapping.

Definition at line 73 of file SimpleAstrometryMapping.cc.

73{ out << *transform; }

◆ setIndex()

void lsst::jointcal::SimpleAstrometryMapping::setIndex ( Eigen::Index i)
inline

Set the index of this mapping in the grand fit.

Definition at line 87 of file SimpleAstrometryMapping.h.

87{ index = i; }

◆ setToBeFit()

void lsst::jointcal::SimpleAstrometryMapping::setToBeFit ( bool value)
inline

Set whether this Mapping is to be fit as part of a Model.

Definition at line 99 of file SimpleAstrometryMapping.h.

◆ transformPosAndErrors()

void lsst::jointcal::SimpleAstrometryMapping::transformPosAndErrors ( FatPoint const & where,
FatPoint & outPoint ) const
overridevirtual

The same as above but without the parameter derivatives (used to evaluate chi^2)

Implements lsst::jointcal::AstrometryMapping.

Reimplemented in lsst::jointcal::SimplePolyMapping.

Definition at line 43 of file SimpleAstrometryMapping.cc.

43 {
44 transform->transformPosAndErrors(where, outPoint);
45 FatPoint tmp;
46 errorProp->transformPosAndErrors(where, tmp);
47 outPoint.vx = tmp.vx;
48 outPoint.vy = tmp.vy;
49 outPoint.vxy = tmp.vxy;
50}

Member Data Documentation

◆ errorProp

std::shared_ptr<AstrometryTransform> lsst::jointcal::SimpleAstrometryMapping::errorProp
protected

Definition at line 111 of file SimpleAstrometryMapping.h.

◆ index

Eigen::Index lsst::jointcal::SimpleAstrometryMapping::index {}
protected

Definition at line 106 of file SimpleAstrometryMapping.h.

106{};

◆ lin

std::unique_ptr<AstrometryTransformLinear> lsst::jointcal::SimpleAstrometryMapping::lin
protected

Definition at line 114 of file SimpleAstrometryMapping.h.

◆ toBeFit

bool lsst::jointcal::SimpleAstrometryMapping::toBeFit
protected

Definition at line 105 of file SimpleAstrometryMapping.h.

◆ transform

std::shared_ptr<AstrometryTransform> lsst::jointcal::SimpleAstrometryMapping::transform
protected

Definition at line 109 of file SimpleAstrometryMapping.h.


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