ManualPolnCalibrator.h

/***************************************************************************
 *
 *   Copyright (C) 2025 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
#ifndef __CalibrationManualPolnCalibrator_H
#define __CalibrationManualPolnCalibrator_H
 
#include "Pulsar/PolnCalibrator.h"
 
namespace Pulsar
{
  class ManualPolnCalibrator : public PolnCalibrator
  {
  public:

    ManualPolnCalibrator (const std::string& filename);

    virtual void load (const std::string& filename);

    void calibrate (Archive* archive);

    void transform (Archive*);

    class Response
    {
    public:


      MJD load (const std::string& str);

      void set_frequency (double f) { frequency = f; }

      double get_frequency () const { return frequency; }

      void set_response (const Jones<float>& J) { response = J; }

      const Jones<float>& get_response () const { return response; }
 
    protected:
 
      double frequency = 0.0; 

      Jones<float> response;
    };

    class FrequencyResponse
    {
    public:

      void set_epoch (const MJD& mjd) { epoch = mjd; }

      MJD get_epoch () const { return epoch; }

      const Response& match (double freq_MHz) const;

      void add (const Response& single) { response.push_back(single); }
 
    protected:

      MJD epoch;

      std::vector<Response> response;
    };

    void add (const MJD& epoch, const Response&);

    const FrequencyResponse& match (const MJD& epoch) const;

    const Jones<float>& get_response (const MJD& epoch, double MHz) const;
 
  protected:

    std::vector<FrequencyResponse> response;

    unsigned current_response = 0;

    std::string ascii_model_filename;

    void transform_work (Archive* arch, const std::string& name, bool invert);
  };
}
 
#endif