21 unsigned index = firstIndex;
22 if (whatToFit.
find(
"Model") == std::string::npos) {
23 LOGLS_WARN(
_log,
"assignIndices was called and Model is *not* in whatToFit");
28 _fittingChips = (whatToFit.
find(
"ModelChip") != std::string::npos);
29 _fittingVisits = (whatToFit.
find(
"ModelVisit") != std::string::npos);
31 if ((!_fittingChips) && (!_fittingVisits)) {
32 _fittingChips = _fittingVisits =
true;
37 auto mapping = idMapping.second.get();
39 if (mapping->isFixed())
continue;
40 mapping->setIndex(index);
41 index += mapping->getNpar();
46 auto mapping = idMapping.second.get();
47 mapping->setIndex(index);
48 index += mapping->getNpar();
52 idMapping.second->setWhatToFit(_fittingChips, _fittingVisits);
60 auto mapping = idMapping.second.get();
62 if (mapping->isFixed())
continue;
63 mapping->offsetParams(delta.segment(mapping->getIndex(), mapping->getNpar()));
68 auto mapping = idMapping.second.get();
69 mapping->offsetParams(delta.segment(mapping->getIndex(), mapping->getNpar()));
76 idMapping.second.get()->freezeErrorTransform();
79 idMapping.second.get()->freezeErrorTransform();
86 mapping->getMappingIndices(indices);
92 total += idMapping.second->getNpar();
95 total += idMapping.second->getNpar();
102 Eigen::VectorXd &derivatives)
const {
104 mapping->computeParameterDerivatives(measuredStar, measuredStar.
getInstFlux(), derivatives);
110 auto coeffs = transfo->getCoefficients();
112 ndarray::Array<double, 2, 2> chebyCoeffs = allocate(ndarray::makeVector(transfo->getNpar(), 4));
113 Eigen::VectorXd::Index k = 0;
114 auto order = transfo->getOrder();
115 for (ndarray::Size j = 0; j <= order; ++j) {
116 ndarray::Size
const iMax = order - j;
117 for (ndarray::Size i = 0; i <= iMax; ++i, ++k) {
118 chebyCoeffs[k][0] = coeffs[j][i];
119 chebyCoeffs[k][1] = 1;
120 chebyCoeffs[k][2] = i;
121 chebyCoeffs[k][3] = j;
130 idMapping.second->dump(stream);
135 idMapping.second->dump(stream);
144 "ConstrainedPhotometryModel cannot find CcdImage " + ccdImage.
getName());
145 return idMapping->second.get();
148 template <
class ChipTransfo,
class VisitTransfo,
class ChipVisitMapping>
150 afw::geom::Box2D
const &focalPlaneBBox,
int visitOrder) {
156 for (
auto const &
ccdImage : ccdImageList) {
166 if (radius2 < minRadius2) {
167 minRadius2 = radius2;
168 constrainedChip = chip;
170 auto photoCalib =
ccdImage->getPhotoCalib();
177 auto visitTransfo = std::make_unique<VisitTransfo>(visitOrder, focalPlaneBBox);
186 for (
auto const &
ccdImage : ccdImageList) {
193 <<
" visit mappings; holding chip " << constrainedChip <<
" fixed (" 209 return mapping->transform(measuredStar, measuredStar.
getInstFlux());
216 return mapping->transformError(measuredStar, measuredStar.
getInstFlux(), tempErr);
226 assert(mapping !=
nullptr);
231 assert(visitTransfo !=
nullptr);
232 auto focalBBox = visitTransfo->getBBox();
235 auto coeff_f = toChebyMapCoeffs(
236 std::dynamic_pointer_cast<PhotometryTransfoChebyshev>(mapping->
getVisitMapping()->getTransfo()));
248 auto transform = pixToFocal->then(chebyTransform)->then(zoomTransform);
252 auto boundedField = std::make_shared<afw::math::TransformBoundedField>(ccdBBox, *
transform);
253 return std::make_shared<afw::image::PhotoCalib>(oldPhotoCalib->getCalibrationMean(),
254 oldPhotoCalib->getCalibrationErr(), boundedField,
false);
267 return mapping->transform(measuredStar, measuredStar.
getInstMag());
274 return mapping->transformError(measuredStar, measuredStar.
getInstFlux(), tempErr);
285 assert(mapping !=
nullptr);
290 assert(visitTransfo !=
nullptr);
291 auto focalBBox = visitTransfo->getBBox();
294 auto coeff_f = toChebyMapCoeffs(
295 std::dynamic_pointer_cast<PhotometryTransfoChebyshev>(mapping->
getVisitMapping()->getTransfo()));
302 using namespace std::string_literals;
304 ast::MathMap(1, 1, {
"y=pow(10.0,x/-2.5)"s}, {
"x=-2.5*log10(y)"s}));
312 auto transform = pixToFocal->then(chebyTransform)->then(logTransform)->then(zoomTransform);
316 auto boundedField = std::make_shared<afw::math::TransformBoundedField>(ccdBBox, *
transform);
317 return std::make_shared<afw::image::PhotoCalib>(oldPhotoCalib->getCalibrationMean(),
318 oldPhotoCalib->getCalibrationErr(), boundedField,
false);
324 afw::geom::Box2D
const &, int);
327 CcdImageList
const &, afw::geom::Box2D
const &, int);
#define LOGLS_WARN(logger, message)
Photometric offset independent of position, defined as -2.5 * log(flux / fluxMag0).
Relates transfo(s) to their position in the fitting matrix and allows interaction with the transfo(s)...
double transformError(CcdImage const &ccdImage, MeasuredStar const &measuredStar) const override
Return the on-sky transformed flux uncertainty for measuredStar on ccdImage.
std::string getName() const
Return the _name that identifies this ccdImage.
CameraSysPrefix const PIXELS
std::shared_ptr< PhotometryMapping > getChipMapping() const
table::Key< int > detector
double getInstMag() const
T bucket_count(T... args)
nth-order 2d Chebyshev photometry transfo, times the input flux.
Photometric offset independent of position, defined as (fluxMag0)^-1.
table::Key< double > calibrationMean
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...
std::shared_ptr< PhotometryMapping > getVisitMapping() const
std::shared_ptr< afw::image::PhotoCalib > toPhotoCalib(CcdImage const &ccdImage) const override
Return the mapping of ccdImage represented as a PhotoCalib.
int getTotalParameters() const override
Return the total number of parameters in this model.
double computeResidual(CcdImage const &ccdImage, MeasuredStar const &measuredStar) const override
Compute the residual between the model applied to a star and its associated fittedStar.
std::shared_ptr< afw::image::PhotoCalib > getPhotoCalib() const
Return the exposure's photometric calibration.
#define LOGLS_DEBUG(logger, message)
Class for a simple mapping implementing a generic Gtransfo.
double getInstFlux() const
PhotometryMappingBase * findMapping(CcdImage const &ccdImage) const override
Return a pointer to the mapping associated with this ccdImage.
objects measured on actual images.
T dynamic_pointer_cast(T... args)
void initialize(CcdImageList const &ccdImageList, afw::geom::Box2D const &focalPlaneBBox, int visitOrder)
Initialize the chip, visit, and chipVisit mappings by creating appropriate transfos and mappings...
void computeParameterDerivatives(MeasuredStar const &measuredStar, CcdImage const &ccdImage, Eigen::VectorXd &derivatives) const override
Compute the parametric derivatives of this model.
#define LOGLS_INFO(logger, message)
double transform(CcdImage const &ccdImage, MeasuredStar const &measuredStar) const override
Return the on-sky transformed flux for measuredStar on ccdImage.
void dump(std::ostream &stream=std::cout) const override
Dump the contents of the transfos, for debugging.
#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 computeResidual(CcdImage const &ccdImage, MeasuredStar const &measuredStar) const override
Compute the residual between the model applied to a star and its associated fittedStar.
nth-order 2d Chebyshev photometry transfo.
CcdImageKey getHashKey() const
virtual double initialChipCalibration(std::shared_ptr< afw::image::PhotoCalib const > photoCalib)=0
Return the initial calibration to use from this photoCalib.
std::shared_ptr< afw::cameraGeom::Detector > getDetector() const
Handler of an actual image from a single CCD.
std::shared_ptr< afw::image::PhotoCalib > toPhotoCalib(CcdImage const &ccdImage) const override
Return the mapping of ccdImage represented as a PhotoCalib.
double transform(CcdImage const &ccdImage, MeasuredStar const &measuredStar) const override
Return the on-sky transformed flux for measuredStar on ccdImage.
nth-order 2d Chebyshev photometry transfo, plus the input flux.
std::shared_ptr< FittedStar > getFittedStar() const
CameraSys const FOCAL_PLANE
void freezeErrorTransform() override
Once this routine has been called, the error transform is not modified by offsetParams().
virtual double transform(CcdImage const &ccdImage, MeasuredStar const &measuredStar) const =0
Return the on-sky transformed flux for measuredStar on 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.
A two-level photometric transform: one for the ccd and one for the visit.
double transformError(CcdImage const &ccdImage, MeasuredStar const &measuredStar) const override
Return the on-sky transformed flux uncertainty for measuredStar on ccdImage.