FaradayConstraint.h

Go to the documentation of this file.

#ifndef DP3_FARADAY_ROTATION_CONSTRAINT_H_
#define DP3_FARADAY_ROTATION_CONSTRAINT_H_
 
#include "Constraint.h"
 
#include "base/CalType.h"
#include "common/PhaseLineFitter.h"
 
#include <vector>
#include <ostream>
 
namespace dp3::ddecal {
 
class FaradayConstraint final : public Constraint {
 public:
  FaradayConstraint(base::CalType diagonal_solution_type,
                    std::optional<double> max_rotation_value)
      : diagonal_solution_type_(diagonal_solution_type),
        max_rotation_value_(max_rotation_value) {}
 
  void Apply(SolutionSpan& solutions, double time) override;
 
  void Initialize(size_t n_antennas,
                  const std::vector<uint32_t>& solutions_per_direction,
                  const std::vector<double>& frequencies) override;
 
  void SetWeights(const std::vector<double>& weights) override;
 
  void SetSubSolutionWeights(
      const std::vector<std::vector<double>>& sub_solution_weights) override;
 
  std::vector<ConstraintResult> GetResult() const override { return results_; }
 
 private:
  void PerformFit(SolutionSpan& solutions, size_t sub_solution, size_t antenna,
                  std::vector<common::phase_fitting::FitSample>& scratch_space);
 
  std::vector<ConstraintResult> results_;
  std::vector<std::vector<double>> sub_solution_weights_;
  std::vector<double> frequencies_;
  common::phase_fitting::SlopeFitRange fit_range_;
  base::CalType diagonal_solution_type_;
  std::optional<double> max_rotation_value_;
};
 
}  // namespace dp3::ddecal
 
#endif