lsst.jointcal  16.0-26-g3163dd8+1
PhotometryModel.cc
Go to the documentation of this file.
1 // -*- LSST-C++ -*-
2 /*
3  * This file is part of jointcal.
4  *
5  * Developed for the LSST Data Management System.
6  * This product includes software developed by the LSST Project
7  * (https://www.lsst.org).
8  * See the COPYRIGHT file at the top-level directory of this distribution
9  * for details of code ownership.
10  *
11  * This program is free software: you can redistribute it and/or modify
12  * it under the terms of the GNU General Public License as published by
13  * the Free Software Foundation, either version 3 of the License, or
14  * (at your option) any later version.
15  *
16  * This program is distributed in the hope that it will be useful,
17  * but WITHOUT ANY WARRANTY; without even the implied warranty of
18  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
19  * GNU General Public License for more details.
20  *
21  * You should have received a copy of the GNU General Public License
22  * along with this program. If not, see <https://www.gnu.org/licenses/>.
23  */
24 
25 #include "lsst/log/Log.h"
26 
27 #include "lsst/jointcal/CcdImage.h"
29 
30 namespace lsst {
31 namespace jointcal {
32 
33 bool PhotometryModel::validate(CcdImageList const &ccdImageList) const {
34  bool check = true;
35  for (auto const &ccdImage : ccdImageList) {
36  // Don't short circuit so that we log every place the model is negative.
37  check &= checkPositiveOnBBox(*ccdImage);
38  }
39  return check;
40 }
41 
43  bool check = true;
44  auto bbox = ccdImage.getImageFrame();
45  for (auto const &x : {bbox.xMin, bbox.getCenter().x, bbox.xMax})
46  for (auto const &y : {bbox.yMin, bbox.getCenter().y, bbox.yMax}) {
47  // flux, fluxErr, instFlux, instFluxErr all 1
49  star.setInstFluxAndErr(1, 1);
50  double result = transform(ccdImage, star);
51  // Don't short circuit so that we log every place the model is negative.
52  if (result < 0) {
53  LOGLS_WARN(_log, "CcdImage " << ccdImage.getName() << " is negative at (" << x << "," << y
54  << "): " << result);
55  check = false;
56  }
57  }
58  return check;
59 }
60 
61 } // namespace jointcal
62 } // namespace lsst
#define LOGLS_WARN(logger, message)
std::string getName() const
Return the _name that identifies this ccdImage.
Definition: CcdImage.h:79
table::Box2IKey bbox
double x
bool checkPositiveOnBBox(CcdImage const &ccdImage) const
Check that the model is positive on the ccdImage bbox.
The base class for handling stars. Used by all matching routines.
Definition: BaseStar.h:50
Class for a simple mapping implementing a generic AstrometryTransform.
void setInstFluxAndErr(double instFlux, double instFluxErr)
Definition: MeasuredStar.h:85
objects measured on actual images.
Definition: MeasuredStar.h:42
std::unique_ptr< SchemaItem< U > > result
LOG_LOGGER _log
lsst.logging instance, to be created by a subclass so that messages have consistent name...
Handler of an actual image from a single CCD.
Definition: CcdImage.h:64
bool validate(CcdImageList const &ccdImageList) const
Return true if this is a "reasonable" model.
int y
virtual double transform(CcdImage const &ccdImage, MeasuredStar const &measuredStar) const =0
Return the on-sky transformed flux for measuredStar on ccdImage.
Frame const & getImageFrame() const
Frame in pixels.
Definition: CcdImage.h:191