lsst.meas.modelfit  21.0.0-2-gecfae73+0341d1150d
UnitTransformedLikelihood.h
Go to the documentation of this file.
1 // -*- lsst-c++ -*-
2 /*
3  * LSST Data Management System
4  * Copyright 2008-2013 LSST Corporation.
5  *
6  * This product includes software developed by the
7  * LSST Project (http://www.lsst.org/).
8  *
9  * This program is free software: you can redistribute it and/or modify
10  * it under the terms of the GNU General Public License as published by
11  * the Free Software Foundation, either version 3 of the License, or
12  * (at your option) any later version.
13  *
14  * This program is distributed in the hope that it will be useful,
15  * but WITHOUT ANY WARRANTY; without even the implied warranty of
16  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17  * GNU General Public License for more details.
18  *
19  * You should have received a copy of the LSST License Statement and
20  * the GNU General Public License along with this program. If not,
21  * see <http://www.lsstcorp.org/LegalNotices/>.
22  */
23 
24 #ifndef LSST_MEAS_MODELFIT_UnitTransformedLikelihood_h_INCLUDED
25 #define LSST_MEAS_MODELFIT_UnitTransformedLikelihood_h_INCLUDED
26 
27 #include <vector>
28 #include <memory>
29 
30 #include "ndarray.h"
31 
32 #include "lsst/pex/config.h"
33 #include "lsst/geom/SpherePoint.h"
37 
42 
43 namespace lsst { namespace meas { namespace modelfit {
44 
52 public:
53 
54  LSST_CONTROL_FIELD(usePixelWeights, bool,
55  "whether to individually weigh pixels using the variance image.");
56 
57  LSST_CONTROL_FIELD(weightsMultiplier, double,
58  "Scaling factor to apply to weights.");
59 
60  explicit UnitTransformedLikelihoodControl(bool usePixelWeights_=false, double weightsMultiplier_=1.0)
61  : usePixelWeights(usePixelWeights_), weightsMultiplier(weightsMultiplier_) {}
62 
63 };
64 
71 public:
79  explicit EpochFootprint(
83  );
84 
88 };
89 
103 public:
104 
107  ndarray::Array<Pixel,2,-1> const & modelMatrix,
108  ndarray::Array<Scalar const,1,1> const & nonlinear,
109  bool doApplyWeights=true
110  ) const override;
111 
123  PTR(Model) model,
124  ndarray::Array<Scalar const,1,1> const & fixed,
125  UnitSystem const & fitSys,
126  geom::SpherePoint const & position,
127  std::vector<PTR(EpochFootprint)> const & epochFootprintList,
129  );
130 
144  PTR(Model) model,
145  ndarray::Array<Scalar const,1,1> const & fixed,
146  UnitSystem const & fitSys,
147  geom::SpherePoint const & position,
148  afw::image::Exposure<Pixel> const & exposure,
149  afw::detection::Footprint const & footprint,
152  );
153 
155 
156 private:
157  class Impl;
158  std::unique_ptr<Impl> _impl;
159 };
160 
161 }}} // namespace lsst::meas::modelfit
162 
163 #endif // !LSST_MEAS_MODELFIT_UnitTransformedLikelihood_h_INCLUDED
#define PTR(...)
An image at one epoch of a galaxy, plus associated info.
afw::image::Exposure< Pixel > const exposure
subregion of exposure that includes footprint
shapelet::MultiShapeletFunction const psf
multi-shapelet model of exposure PSF
afw::detection::Footprint const footprint
footprint of source (galaxy)
EpochFootprint(afw::detection::Footprint const &footprint, afw::image::Exposure< Pixel > const &exposure, shapelet::MultiShapeletFunction const &psf)
Construct a EpochFootprint.
Base class for optimizer/sampler likelihood functions that compute likelihood at a point.
Definition: Likelihood.h:70
Abstract base class and concrete factories that define multi-shapelet galaxy models.
Definition: Model.h:56
Control object used to initialize a UnitTransformedLikelihood.
double weightsMultiplier
"Scaling factor to apply to weights." ;
bool usePixelWeights
"whether to individually weigh pixels using the variance image." ;
UnitTransformedLikelihoodControl(bool usePixelWeights_=false, double weightsMultiplier_=1.0)
A concrete Likelihood class that does not require its parameters and data to be in the same UnitSyste...
void computeModelMatrix(ndarray::Array< Pixel, 2,-1 > const &modelMatrix, ndarray::Array< Scalar const, 1, 1 > const &nonlinear, bool doApplyWeights=true) const override
Evaluate the model for the given vector of nonlinear parameters.
UnitTransformedLikelihood(boost::shared_ptr< Model > model, ndarray::Array< Scalar const, 1, 1 > const &fixed, UnitSystem const &fitSys, geom::SpherePoint const &position, afw::image::Exposure< Pixel > const &exposure, afw::detection::Footprint const &footprint, shapelet::MultiShapeletFunction const &psf, UnitTransformedLikelihoodControl const &ctrl)
Initialize a UnitTransformedLikelihood with data from multiple exposures.
UnitTransformedLikelihood(boost::shared_ptr< Model > model, ndarray::Array< Scalar const, 1, 1 > const &fixed, UnitSystem const &fitSys, geom::SpherePoint const &position, std::vector< boost::shared_ptr< EpochFootprint >> const &epochFootprintList, UnitTransformedLikelihoodControl const &ctrl)
Initialize a UnitTransformedLikelihood with data from multiple exposures.
A simple struct that combines a Wcs and a PhotoCalib.
Definition: UnitSystem.h:42