12 double instFluxFromMag(
double mag) {
return std::pow(10, mag / -2.5); }
19 unsigned ipar = firstIndex;
20 for (
auto const &i :
_myMap) {
21 auto mapping = i.second.get();
22 mapping->setIndex(ipar);
23 ipar += mapping->getNpar();
30 auto mapping = i.second.get();
31 mapping->offsetParams(delta.segment(mapping->getIndex(), mapping->getNpar()));
37 i.second->freezeErrorTransform();
44 if (indices.
size() < mapping->getNpar()) indices.
resize(mapping->getNpar());
45 indices[0] = mapping->getIndex();
51 total += i.second->getNpar();
58 Eigen::VectorXd &derivatives)
const {
60 mapping->computeParameterDerivatives(measuredStar, measuredStar.
getInstFlux(), derivatives);
65 stream << i.first <<
": ";
66 i.second->dump(stream);
75 "SimplePhotometryModel cannot find CcdImage " + ccdImage.
getName());
76 return i->second.get();
81 for (
auto const &
ccdImage : ccdImageList) {
82 auto photoCalib =
ccdImage->getPhotoCalib();
84 auto transfo = std::make_shared<FluxTransfoSpatiallyInvariant>(photoCalib->getCalibrationMean());
96 return mapping->transform(star, star.
getInstFlux());
102 return mapping->transformError(star, star.
getInstFlux(), tempErr);
108 return std::make_unique<afw::image::PhotoCalib>(calibration, oldPhotoCalib->getCalibrationErr());
113 for (
auto const &
ccdImage : ccdImageList) {
114 auto photoCalib =
ccdImage->getPhotoCalib();
116 double calib = magFromFlux(photoCalib->getCalibrationMean());
117 auto transfo = std::make_shared<MagnitudeTransfoSpatiallyInvariant>(calib);
130 return mapping->transform(star, star.
getInstMag());
136 return mapping->transformError(star, star.
getInstMag(), tempErr);
141 double calibration = instFluxFromMag(
findMapping(ccdImage)->getParameters()[0]);
143 return std::make_unique<afw::image::PhotoCalib>(calibration, oldPhotoCalib->getCalibrationErr());
double transform(CcdImage const &ccdImage, MeasuredStar const &measuredStar) const override
Return the on-sky transformed flux for measuredStar on ccdImage.
std::shared_ptr< afw::image::PhotoCalib > toPhotoCalib(CcdImage const &ccdImage) const override
Return the mapping of ccdImage represented as a PhotoCalib.
Relates transfo(s) to their position in the fitting matrix and allows interaction with the transfo(s)...
virtual Eigen::VectorXd getParameters()=0
std::string getName() const
Return the _name that identifies this ccdImage.
int getTotalParameters() const override
Return the total number of parameters in this model.
SimpleMagnitudeModel(CcdImageList const &ccdImageList, double errorPedestal_=0)
std::shared_ptr< afw::image::PhotoCalib > toPhotoCalib(CcdImage const &ccdImage) const override
Return the mapping of ccdImage represented as a PhotoCalib.
double getInstMag() const
void computeParameterDerivatives(MeasuredStar const &measuredStar, CcdImage const &ccdImage, Eigen::VectorXd &derivatives) const override
Compute the parametric derivatives of this model.
PhotometryMappingBase * findMapping(CcdImage const &ccdImage) const override
Return the mapping associated with this ccdImage.
void getMappingIndices(CcdImage const &ccdImage, std::vector< unsigned > &indices) const override
Get how this set of parameters (of length Npar()) map into the "grand" fit.
std::shared_ptr< afw::image::PhotoCalib > getPhotoCalib() const
Return the exposure's photometric calibration.
Class for a simple mapping implementing a generic Gtransfo.
SimpleFluxModel(CcdImageList const &ccdImageList, double errorPedestal_=0)
double getInstFlux() const
double computeResidual(CcdImage const &ccdImage, MeasuredStar const &measuredStar) const override
Compute the residual between the model applied to a star and its associated fittedStar.
objects measured on actual images.
double transform(CcdImage const &ccdImage, MeasuredStar const &measuredStar) const override
Return the on-sky transformed flux for measuredStar on ccdImage.
#define LOGLS_INFO(logger, message)
unsigned assignIndices(std::string const &whatToFit, unsigned firstIndex) override
Assign indices in the full matrix to the parameters being fit in the mappings, starting at firstIndex...
#define LSST_EXCEPT(type,...)
double tweakFluxError(jointcal::MeasuredStar const &measuredStar) const
Add a fraction of the instrumental flux to the instrumental flux error, in quadrature.
LOG_LOGGER _log
lsst.logging instance, to be created by a subclass so that messages have consistent name...
void offsetParams(Eigen::VectorXd const &delta) override
Offset the parameters by the provided amounts (by -delta).
double tweakMagnitudeError(jointcal::MeasuredStar const &measuredStar) const
Add a small magnitude offset to the "instrumental magnitude" error, in quadrature.
double computeResidual(CcdImage const &ccdImage, MeasuredStar const &measuredStar) const override
Compute the residual between the model applied to a star and its associated fittedStar.
CcdImageKey getHashKey() const
void dump(std::ostream &stream=std::cout) const override
Dump the contents of the transfos, for debugging.
Handler of an actual image from a single CCD.
std::shared_ptr< FittedStar > getFittedStar() const
double transformError(CcdImage const &ccdImage, MeasuredStar const &measuredStar) const override
Return the on-sky transformed flux uncertainty for measuredStar on ccdImage.
Photometric response model which has a single photometric factor per CcdImage.
double transformError(CcdImage const &ccdImage, MeasuredStar const &measuredStar) const override
Return the on-sky transformed flux uncertainty for measuredStar on ccdImage.
void freezeErrorTransform() override
Once this routine has been called, the error transform is not modified by offsetParams().