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;
 
    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