lsst.jointcal  22.0.1-22-g41e933f+587b875011
StarList.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 #ifndef STARLIST__CC
26 #define STARLIST__CC
27 
28 #include "lsst/pex/exceptions.h"
29 
30 #include "lsst/jointcal/Frame.h"
31 #include "lsst/jointcal/StarList.h"
32 #include "lsst/jointcal/BaseStar.h"
35 
37 
38 namespace lsst {
39 namespace jointcal {
40 
41 template <class Star>
43  typedef StarList<Star>::Element E;
44  this->sort([](const E &e1, const E &e2) { return (e1->getFlux() > e2->getFlux()); });
45 }
46 
47 template <class Star>
48 void StarList<Star>::cutTail(const int nKeep) {
49  int count = 0;
50  auto si = this->begin();
51  for (; si != this->end() && count < nKeep; ++count, ++si)
52  ;
53  while (si != this->end()) {
54  si = this->erase(si);
55  }
56 }
57 
58 template <class Star>
59 void StarList<Star>::extractInFrame(StarList<Star> &out, const Frame &frame) const {
60  for (auto const &star : *this) {
61  if (frame.inFrame(*star)) {
62  out.push_back(std::make_shared<Star>(*star));
63  }
64  }
65 }
66 
67 template <class Star>
69  copy.clearList();
70  for (auto const &si : *this) copy.push_back(std::make_shared<Star>(*si));
71 }
72 
73 // Explicit instantiations
74 template class StarList<BaseStar>;
75 template class StarList<FittedStar>;
76 template class StarList<MeasuredStar>;
77 } // namespace jointcal
78 } // namespace lsst
79 
80 #endif /* STARLIST__CC */
int end
T begin(T... args)
rectangle with sides parallel to axes.
Definition: Frame.h:38
bool inFrame(double x, double y) const
inside?
Definition: Frame.cc:120
std::lists of Stars.
Definition: StarList.h:58
void extractInFrame(StarList< Star > &out, const Frame &frame) const
copy the part of the std::list which is included in the frame at the end of another std::list
Definition: StarList.cc:59
void cutTail(const int nKeep)
cuts the end of the std::list
Definition: StarList.cc:48
void copyTo(StarList< Star > &copy) const
clears copy and makes a copy of the std::list to copy
Definition: StarList.cc:68
void fluxSort()
a model routine to sort the std::list
Definition: StarList.cc:42
T copy(T... args)
T count(T... args)
def erase(frame=None)
Class for a simple mapping implementing a generic AstrometryTransform.
T push_back(T... args)
T sort(T... args)