LSSTApplications  20.0.0
LSSTDataManagementBasePackage
Chi2.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 <utility>
26 #include <iostream>
27 
28 #include "lsst/jointcal/Chi2.h"
29 
30 namespace lsst {
31 namespace jointcal {
32 
34  double sum = 0;
35  double sum2 = 0;
36  for (auto i : *this) {
37  sum += i.chi2;
38  sum2 += std::pow(i.chi2, 2);
39  }
40  double average = sum / size();
41  double sigma = sqrt(sum2 / size() - std::pow(average, 2));
42  return std::make_pair(average, sigma);
43 }
44 
46  s << "chi2 per star : ";
47  for (auto chi2 : chi2List) {
48  s << *(chi2.star) << " chi2: " << chi2.chi2 << " ; ";
49  }
50  s << std::endl;
51  return s;
52 }
53 
54 } // namespace jointcal
55 } // namespace lsst
lsst::jointcal::Chi2List
Structure to accumulate the chi2 contributions per each star (to help find outliers).
Definition: Chi2.h:100
std::pair< double, double >
std::vector< Chi2Star >::size
Chi2Star size(Chi2Star ... args)
sigma
afw::table::Key< double > sigma
Definition: GaussianPsf.cc:50
std::sqrt
T sqrt(T... args)
Chi2.h
std::ostream
STL class.
lsst::jointcal::operator<<
std::ostream & operator<<(std::ostream &stream, AstrometryMapping const &mapping)
Definition: AstrometryMapping.h:78
lsst::jointcal
Definition: Associations.h:49
lsst
A base class for image defects.
Definition: imageAlgorithm.dox:1
std::endl
T endl(T... args)
lsst::jointcal::Chi2List::computeAverageAndSigma
std::pair< double, double > computeAverageAndSigma()
Compute the average and std-deviation of these chisq values.
Definition: Chi2.cc:33
std::make_pair
T make_pair(T... args)
std::pow
T pow(T... args)