LSSTApplications
20.0.0
LSSTDataManagementBasePackage
stack
1a1d771
Linux64
jointcal
20.0.0
src
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
33
std::pair<double, double>
Chi2List::computeAverageAndSigma
() {
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
45
std::ostream
&
operator<<
(
std::ostream
& s,
Chi2List
const
& chi2List) {
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)
Generated on Wed Jun 24 2020 18:10:08 for LSSTApplications by
1.8.18