2 #ifndef LSST_JOINTCAL_PHOTOMETRY_MAPPING_H 3 #define LSST_JOINTCAL_PHOTOMETRY_MAPPING_H 33 virtual unsigned getNpar()
const = 0;
54 Eigen::Ref<Eigen::VectorXd> derivatives)
const = 0;
62 virtual void offsetParams(Eigen::VectorXd
const &delta) = 0;
105 return _transfo->getNpar();
111 return _transfo->transform(measuredStar.
x, measuredStar.
y, instFlux);
116 Eigen::Ref<Eigen::VectorXd> derivatives)
const override {
120 _transfo->computeParameterDerivatives(measuredStar.
x, measuredStar.
y, instFlux, derivatives);
125 void offsetParams(Eigen::VectorXd
const &delta)
override { _transfo->offsetParams(delta); }
128 Eigen::VectorXd
getParameters()
override {
return _transfo->getParameters(); }
133 for (
unsigned k = 0; k <
getNpar(); ++k) {
134 indices[k] =
index + k;
140 stream <<
"index: " <<
index <<
" params: ";
141 _transfo->dump(stream);
159 _chipMapping(
std::
move(chipMapping)),
160 _visitMapping(
std::
move(visitMapping)) {}
163 unsigned getNpar()
const override {
return _chipMapping->getNpar() + _visitMapping->getNpar(); }
167 double tempFlux = _chipMapping->getTransfo()->transform(measuredStar.
x, measuredStar.
y, instFlux);
168 return _visitMapping->getTransfo()->transform(measuredStar.
getXFocal(), measuredStar.
getYFocal(),
174 Eigen::Ref<Eigen::VectorXd> derivatives)
const override;
178 _chipMapping->offsetParams(delta.segment(0, _chipMapping->getNpar()));
179 _visitMapping->offsetParams(delta.segment(_chipMapping->getNpar(), _visitMapping->getNpar()));
184 Eigen::VectorXd joined(
getNpar());
185 joined << _chipMapping->getParameters(), _visitMapping->getParameters();
194 stream <<
"index: " <<
index <<
" chipMapping: ";
195 _chipMapping->dump(stream);
196 stream <<
"visitMapping: ";
197 _visitMapping->dump(stream);
212 #endif // LSST_JOINTCAL_PHOTOMETRY_MAPPING_H virtual double transform(MeasuredStar const &measuredStar, double instFlux) const =0
Return the on-sky transformed flux at (x,y).
PhotometryMapping(std::shared_ptr< PhotometryTransfo > transfo)
void computeParameterDerivatives(MeasuredStar const &measuredStar, double instFlux, Eigen::Ref< Eigen::VectorXd > derivatives) const override
Compute the derivatives with respect to the parameters (i.e.
Relates transfo(s) to their position in the fitting matrix and allows interaction with the transfo(s)...
virtual Eigen::VectorXd getParameters()=0
std::shared_ptr< PhotometryMapping > getChipMapping() const
void dump(std::ostream &stream=std::cout) const override
Dump the contents of the transfos, for debugging.
void offsetParams(Eigen::VectorXd const &delta) override
Offset the transfo parameters by delta.
Eigen::VectorXd getParameters() override
double transform(MeasuredStar const &measuredStar, double instFlux) const override
Return the on-sky transformed flux at (x,y).
std::shared_ptr< PhotometryMapping > getVisitMapping() const
void setIndex(unsigned i)
Set the index of this mapping in the grand fit.
void getMappingIndices(std::vector< unsigned > &indices) const override
Gets how this set of parameters (of length getNpar()) map into the "grand" fit.
Class for a simple mapping implementing a generic Gtransfo.
void dump(std::ostream &stream=std::cout) const override
Dump the contents of the transfos, for debugging.
void offsetParams(Eigen::VectorXd const &delta) override
Offset the transfo parameters by delta.
objects measured on actual images.
double transform(MeasuredStar const &measuredStar, double instFlux) const override
Return the on-sky transformed flux at (x,y).
virtual ~PhotometryMappingBase()
ChipVisitPhotometryMapping(std::shared_ptr< PhotometryMapping > chipMapping, std::shared_ptr< PhotometryMapping > visitMapping)
unsigned getNpar() const override
Number of total parameters in this mapping.
PhotometryMappingBase & operator=(PhotometryMappingBase const &)=delete
void setFixed(bool _fixed)
Make this mapping's parameters fixed (i.e. not varied during fitting).
unsigned getIndex()
Get the index of this mapping in the grand fit.
std::shared_ptr< PhotometryTransfo > getTransfo() const
virtual void offsetParams(Eigen::VectorXd const &delta)=0
Offset the transfo parameters by delta.
A mapping containing a single photometryTransfo.
virtual void dump(std::ostream &stream=std::cout) const =0
Dump the contents of the transfos, for debugging.
virtual void computeParameterDerivatives(MeasuredStar const &measuredStar, double instFlux, Eigen::Ref< Eigen::VectorXd > derivatives) const =0
Compute the derivatives with respect to the parameters (i.e.
virtual unsigned getNpar() const =0
Number of total parameters in this mapping.
virtual void getMappingIndices(std::vector< unsigned > &indices) const =0
Gets how this set of parameters (of length getNpar()) map into the "grand" fit.
A two-level photometric transform: one for the ccd and one for the visit.
unsigned getNpar() const override
Number of total parameters in this mapping.
Eigen::VectorXd getParameters() override