LSSTApplications  17.0+124,17.0+14,17.0+73,18.0.0+37,18.0.0+80,18.0.0-4-g68ffd23+4,18.1.0-1-g0001055+12,18.1.0-1-g03d53ef+5,18.1.0-1-g1349e88+55,18.1.0-1-g2505f39+44,18.1.0-1-g5315e5e+4,18.1.0-1-g5e4b7ea+14,18.1.0-1-g7e8fceb+4,18.1.0-1-g85f8cd4+48,18.1.0-1-g8ff0b9f+4,18.1.0-1-ga2c679d+1,18.1.0-1-gd55f500+35,18.1.0-10-gb58edde+2,18.1.0-11-g0997b02+4,18.1.0-13-gfe4edf0b+12,18.1.0-14-g259bd21+21,18.1.0-19-gdb69f3f+2,18.1.0-2-g5f9922c+24,18.1.0-2-gd3b74e5+11,18.1.0-2-gfbf3545+32,18.1.0-26-g728bddb4+5,18.1.0-27-g6ff7ca9+2,18.1.0-3-g52aa583+25,18.1.0-3-g8ea57af+9,18.1.0-3-gb69f684+42,18.1.0-3-gfcaddf3+6,18.1.0-32-gd8786685a,18.1.0-4-gf3f9b77+6,18.1.0-5-g1dd662b+2,18.1.0-5-g6dbcb01+41,18.1.0-6-gae77429+3,18.1.0-7-g9d75d83+9,18.1.0-7-gae09a6d+30,18.1.0-9-gc381ef5+4,w.2019.45
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
friend std::ostream & operator<<(std::ostream &s, Chi2List const &chi2List)
Definition: Chi2.cc:45
T endl(T... args)
afw::table::Key< double > sigma
Definition: GaussianPsf.cc:50
Structure to accumulate the chi2 contributions per each star (to help find outliers).
Definition: Chi2.h:100
A base class for image defects.
T make_pair(T... args)
Chi2Star size(Chi2Star ... args)
T pow(T... args)
T sqrt(T... args)
STL class.
std::pair< double, double > computeAverageAndSigma()
Compute the average and std-deviation of these chisq values.
Definition: Chi2.cc:33