PointSource.h

Go to the documentation of this file.

// PointSource.h: Point source model component with optional spectral index and
// rotation measure.
//
// Copyright (C) 2020 ASTRON (Netherlands Institute for Radio Astronomy)
// SPDX-License-Identifier: GPL-3.0-or-later
 
#ifndef DPPP_POINTSOURCE_H
#define DPPP_POINTSOURCE_H
 
#include <vector>
 
#include "ModelComponent.h"
#include "Stokes.h"
#include "base/Direction.h"
 
#include <memory>
 
namespace dp3 {
namespace base {
 
 
 
class PointSource : public ModelComponent {
 public:
  typedef std::shared_ptr<PointSource> Ptr;
  typedef std::shared_ptr<const PointSource> ConstPtr;
 
  PointSource(const Direction &position);
  PointSource(const Direction &position, const Stokes &stokes);
 
  const Direction &direction() const override { return itsDirection; }
  void setDirection(const Direction &position);
 
  void setStokes(const Stokes &stokes);
 
  template <typename T>
  void setSpectralTerms(double refFreq, bool isLogarithmic, T first, T last);
 
  void setRotationMeasure(double fraction, double angle, double rm);
 
  Stokes stokes(double freq) const;
  Stokes stokes() const;
  const std::vector<double> &spectrum() const { return itsSpectralTerms; }
  double referenceFreq() const { return itsRefFreq; }
 
  void accept(ModelComponentVisitor &visitor) const override;
 
  bool hasLogarithmicSI() const { return itsHasLogarithmicSI; }
 
  bool hasSpectralTerms() const;
  bool hasRotationMeasure() const;
 
  double polarizationAngle() const { return itsPolarizationAngle; }
  double polarizedFraction() const { return itsPolarizedFraction; }
  double rotationMeasure() const { return itsRotationMeasure; }
 
 private:
  Direction itsDirection;
  Stokes itsStokes;
  double itsRefFreq;
  std::vector<double> itsSpectralTerms;
  double itsPolarizedFraction;
  double itsPolarizationAngle;
  double itsRotationMeasure;
  bool itsHasRotationMeasure;
  bool itsHasLogarithmicSI;
};
 
 
// -------------------------------------------------------------------------- //
// - Implementation: PointSource                                            - //
// -------------------------------------------------------------------------- //
 
template <typename T>
void PointSource::setSpectralTerms(double refFreq, bool isLogarithmic, T first,
                                   T last) {
  itsRefFreq = refFreq;
  itsHasLogarithmicSI = isLogarithmic;
  itsSpectralTerms.clear();
  itsSpectralTerms.insert(itsSpectralTerms.begin(), first, last);
}
 
}  // namespace base
}  // namespace dp3
 
#endif