7 LOG_LOGGER _log =
LOG_GET(
"jointcal.PhotometryMapping");
15 Eigen::Ref<Eigen::VectorXd> derivatives)
const {
20 double chipScale = _chipMapping->getTransfo()->transform(measuredStar.
x, measuredStar.
y, 1);
22 _visitMapping->getTransfo()->transform(measuredStar.
getXFocal(), measuredStar.
getYFocal(), 1);
26 if (_nParChips > 0 && !_chipMapping->isFixed()) {
28 Eigen::Ref<Eigen::VectorXd> chipBlock = derivatives.segment(0, _nParChips);
29 _chipMapping->getTransfo()->computeParameterDerivatives(measuredStar.
x, measuredStar.
y, instFlux,
31 chipBlock *= visitScale;
33 if (_nParVisits > 0) {
35 Eigen::Ref<Eigen::VectorXd> visitBlock = derivatives.segment(_nParChips, _nParVisits);
36 _visitMapping->getTransfo()->computeParameterDerivatives(
38 visitBlock *= chipScale;
43 double instFluxErr)
const {
46 _chipMapping->getTransfoErrors()->transform(measuredStar.
x, measuredStar.
y, instFluxErr);
47 return _visitMapping->getTransfoErrors()->transform(measuredStar.
getXFocal(), measuredStar.
getYFocal(),
54 _chipMapping->getMappingIndices(indices);
56 if (_nParVisits > 0) {
60 _visitMapping->getMappingIndices(tempIndices);
61 for (
unsigned k = 0; k < _visitMapping->getNpar(); ++k) {
62 indices.
at(k + _nParChips) = tempIndices.at(k);
69 _nParChips = _chipMapping->getNpar();
74 _nParVisits = _visitMapping->getNpar();
void computeParameterDerivatives(MeasuredStar const &measuredStar, double instFlux, Eigen::Ref< Eigen::VectorXd > derivatives) const override
Compute the derivatives with respect to the parameters (i.e.
void getMappingIndices(std::vector< unsigned > &indices) const override
Gets how this set of parameters (of length getNpar()) map into the "grand" fit.
void setWhatToFit(bool const fittingChips, bool const fittingVisits)
Set whether to fit chips or visits.
double transformError(MeasuredStar const &measuredStar, double instFlux, double instFluxErr) const override
Return the on-sky transformed flux uncertainty for measuredStar on ccdImage.
Class for a simple mapping implementing a generic Gtransfo.
objects measured on actual images.
unsigned getNpar() const override
Number of total parameters in this mapping.