LSST Applications 27.0.0,g0265f82a02+469cd937ee,g02d81e74bb+21ad69e7e1,g1470d8bcf6+cbe83ee85a,g2079a07aa2+e67c6346a6,g212a7c68fe+04a9158687,g2305ad1205+94392ce272,g295015adf3+81dd352a9d,g2bbee38e9b+469cd937ee,g337abbeb29+469cd937ee,g3939d97d7f+72a9f7b576,g487adcacf7+71499e7cba,g50ff169b8f+5929b3527e,g52b1c1532d+a6fc98d2e7,g591dd9f2cf+df404f777f,g5a732f18d5+be83d3ecdb,g64a986408d+21ad69e7e1,g858d7b2824+21ad69e7e1,g8a8a8dda67+a6fc98d2e7,g99cad8db69+f62e5b0af5,g9ddcbc5298+d4bad12328,ga1e77700b3+9c366c4306,ga8c6da7877+71e4819109,gb0e22166c9+25ba2f69a1,gb6a65358fc+469cd937ee,gbb8dafda3b+69d3c0e320,gc07e1c2157+a98bf949bb,gc120e1dc64+615ec43309,gc28159a63d+469cd937ee,gcf0d15dbbd+72a9f7b576,gdaeeff99f8+a38ce5ea23,ge6526c86ff+3a7c1ac5f1,ge79ae78c31+469cd937ee,gee10cc3b42+a6fc98d2e7,gf1cff7945b+21ad69e7e1,gfbcc870c63+9a11dc8c8f
LSST Data Management Base Package
Loading...
Searching...
No Matches
Histo4d.h
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#ifndef LSST_JOINTCAL_HISTO4D_H
26#define LSST_JOINTCAL_HISTO4D_H
27
28namespace lsst {
29namespace jointcal {
30
34public:
35 SparseHisto4d() = default;
36 // obvious meanings. NEntries is used as the size of the primary allocation.
37 SparseHisto4d(int n1, double min1, double max1, int n2, double min2, double max2,
38 int n3, double min3, double max3, int n4, double min4, double max4,
39 int nEntries);
41 void fill(const double x[4]);
43 void fill(double x1, double x2, double x3, double x4);
45 int maxBin(double x[4]);
46
48 void zeroBin(double x[4]);
49
51 void binLimits(const double x[4], int idim, double &xMin, double &xMax) const;
52
54 int getNEntries() const { return _ndata; }
55
56 ~SparseHisto4d() = default;
57
58 // private:
59 int code_value(const double x[4]) const;
60 void inverse_code(int code, double x[4]) const;
61 void sort();
62 void print() const;
63
64private:
65 std::vector<int> _data;
66 int _ndata{};
67 int _dataSize{};
68 int _n[4]{};
69 double _minVal[4]{}, _maxVal[4]{};
70 double _scale[4]{};
71 bool _sorted{};
72};
73} // namespace jointcal
74} // namespace lsst
75
76#endif // LSST_JOINTCAL_HISTO4D_H
A class to histogram in 4 dimensions.
Definition Histo4d.h:33
void zeroBin(double x[4])
Definition Histo4d.cc:136
int maxBin(double x[4])
Definition Histo4d.cc:114
void inverse_code(int code, double x[4]) const
Definition Histo4d.cc:77
void fill(const double x[4])
Definition Histo4d.cc:92
void binLimits(const double x[4], int idim, double &xMin, double &xMax) const
return the bin limits of dimension idim (0<=idim<4), around point X.
Definition Histo4d.cc:151
int code_value(const double x[4]) const
Definition Histo4d.cc:67