SingleAxis.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2004 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// psrchive/More/Polarimetry/Pulsar/SingleAxis.h
 
#ifndef __Calibration_SingleAxis_H
#define __Calibration_SingleAxis_H
 
#include "MEAL/ProductRule.h"
#include "MEAL/Complex2.h"
 
#include "Estimate.h"
#include "Stokes.h"
 
namespace MEAL
{
  template<typename> class Gain;
  class Boost1;
  class Rotation1;
}
 
namespace Calibration {


  class SingleAxis : public MEAL::ProductRule<MEAL::Complex2> {
 
  public:

    SingleAxis ();

    SingleAxis (const SingleAxis& s);

    const SingleAxis& operator = (const SingleAxis& s);

    ~SingleAxis ();

    SingleAxis* clone () const override;

    void set_cyclic (bool flag = true);


    Estimate<double> get_gain () const;

    Estimate<double> get_diff_gain () const;
    
    Estimate<double> get_diff_phase () const;

    void set_gain (const Estimate<double>& gain);

    void set_diff_gain (const Estimate<double>& gamma);
    
    void set_diff_phase (const Estimate<double>& phi);
 
 
    // Non-free parameters

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

    Vector<3, double> get_axis () const;

    void solve (const std::vector< Estimate<double> >& output);

    void invert ();

    const SingleAxis& operator *= (const SingleAxis& other);
 
    // ///////////////////////////////////////////////////////////////////
    //
    // Model implementation
    //
    // ///////////////////////////////////////////////////////////////////

    std::string get_name () const;

    MEAL::Gain<Complex2>* get_gain_transformation ();

    MEAL::Boost1* get_boost_transformation ();

    MEAL::Rotation1* get_rotation_transformation ();
 
  protected:

    Reference::To< MEAL::Gain<Complex2> > gain;

    Reference::To<MEAL::Boost1> boost;

    Reference::To<MEAL::Rotation1> rotation;
 
  private:

    void init ();
 
  };
 
}
 
#endif