Faraday.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2004 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/More/Polarimetry/Pulsar/Faraday.h
 
#ifndef __Calibration_Faraday_H
#define __Calibration_Faraday_H
 
#include "MEAL/Rotation1.h"
 
namespace Calibration {


  class Faraday : public MEAL::Complex2 {
 
  public:
 
    Faraday ();

    std::string get_name () const;

    void set_rotation_measure (const Estimate<double>& rotation_measure);
    Estimate<double> get_rotation_measure () const;

    void set_reference_frequency (double MHz);
    double get_reference_frequency () const;

    void set_frequency (double MHz);
    double get_frequency () const;

    void set_reference_wavelength (double metres);
    double get_reference_wavelength () const;

    void set_wavelength (double metres);
    double get_wavelength () const;

    void set_axis (const Vector<3, double>& axis);

    double get_rotation () const;

    typedef Jones<double> Return;
    void set_measure (double rotation_measure);
    double get_measure () const;
 
  protected:

    MEAL::Rotation1 rotation;

    double reference_wavelength;

    double wavelength;
 
    // ///////////////////////////////////////////////////////////////////
    //
    // MEAL::Complex2 implementation
    //
    // ///////////////////////////////////////////////////////////////////

    void calculate (Jones<double>& result, std::vector<Jones<double> >* grad);
 
  };
 
}
 
#endif