PatchInfo.h

Go to the documentation of this file.

// PatchInfo.h: Info about a patch
//
// Copyright (C) 2020 ASTRON (Netherlands Institute for Radio Astronomy)
// SPDX-License-Identifier: GPL-3.0-or-later
 
 
#ifndef DP3_SKY_MODEL_PATCHINFO_H
#define DP3_SKY_MODEL_PATCHINFO_H
 
#include <cmath>
#include <string>
 
namespace dp3::sky_model {
 
class PatchInfo {
 public:
  PatchInfo() = default;
 
  PatchInfo(const std::string& name, double ra, double dec, int category,
            double apparentBrightness)
      : itsName(name),
        itsRa(ra),
        itsDec(dec),
        itsCategory(category),
        itsAppBrightness(apparentBrightness) {}
 
  const std::string& getName() const { return itsName; }
 
  double getRa() const { return itsRa; }
 
  double getDec() const { return itsDec; }
 
  int getCategory() const { return itsCategory; }
 
  double apparentBrightness() const { return itsAppBrightness; }
 
  void setRa(double ra) { itsRa = ra; }
 
  void setDec(double dec) { itsDec = dec; }
 
  void setApparentBrightness(double apparentBrightness) {
    itsAppBrightness = apparentBrightness;
  }
 
 private:
  std::string itsName;
  double itsRa = 0.0;
  double itsDec = 0.0;
  int itsCategory = 0;
  double itsAppBrightness = 0.0;
};
 
std::ostream& operator<<(std::ostream& os, const PatchInfo& info);
 
class PatchSumInfo {
 public:
  explicit PatchSumInfo(unsigned int patchId) : itsPatchId(patchId) {}
 
  void add(double ra, double dec, double flux);
 
  double getFlux() const { return itsSumFlux; }
 
  double getRa() const {
    return std::atan2(itsSumY / itsSumFlux, itsSumX / itsSumFlux);
  }
  double getDec() const { return std::asin(itsSumZ / itsSumFlux); }
 
  unsigned int getPatchId() const { return itsPatchId; }
 
 private:
  double itsSumX = 0.0;
  double itsSumY = 0.0;
  double itsSumZ = 0.0;
  double itsSumFlux = 0.0;
  unsigned int itsPatchId = 0;
};
 
void toSkyModel(std::ostream& output, const PatchInfo& patch);
 
}  // namespace dp3::sky_model
 
#endif