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 LOFAR_PARMDB_PATCHINFO_H
#define LOFAR_PARMDB_PATCHINFO_H
 
#include <cmath>
#include <string>
 
namespace dp3 {
 
namespace blob {
class BlobOStream;
class BlobIStream;
}  // namespace blob
 
namespace parmdb {
 
 
class PatchInfo {
 public:
  PatchInfo() {}
 
  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);
 
blob::BlobOStream operator<<(blob::BlobOStream& os, const PatchInfo& info);
 
blob::BlobIStream operator>>(blob::BlobIStream& os, PatchInfo& info);
 
class PatchSumInfo {
 public:
  explicit PatchSumInfo(unsigned int patchId)
      : itsSumX(0),
        itsSumY(0),
        itsSumZ(0),
        itsSumFlux(0),
        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;
  double itsSumY;
  double itsSumZ;
  double itsSumFlux;
  unsigned int itsPatchId;
};
 
void toSkymodel(std::ostream& output, const PatchInfo& patch);
 
 
}  // namespace parmdb
}  // namespace dp3
 
#endif