mode.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2016 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
// epsic/src/mode.h
 
#ifndef __epsic_mode_h
#define __epsic_mode_h
 
#include "Stokes.h"
#include "Spinor.h"
#include "Matrix.h"
#include "BoxMuller.h"
#include "Minkowski.h"
 
#include <iostream>
 
namespace epsic
{
 
  /***************************************************************************
   *
   *  a single source of electromagnetic radiation
   *
   ***************************************************************************/
 
  class mode
  {
  public:
    mode ();
    virtual ~mode () { }
 
    virtual void set_Stokes (const Stokes<double>& mean);
    virtual Stokes<double> get_Stokes () { return mean; }

    virtual Stokes<double> get_mean () const { return mean; }

    virtual Matrix<4,4, double> get_covariance () const
    { return Minkowski::outer (mean, mean); }

    virtual Matrix<4,4, double> get_crosscovariance (unsigned ilag) const
    { if (ilag>0) return 0; else return get_covariance(); }

    virtual Spinor<double> get_field ();

    virtual BoxMuller* get_normal () { return normal; }
    virtual void set_normal (BoxMuller* n) { normal = n; }

    const Jones<double>& get_polarizer () const { return polarizer; }
      
  private:
    Stokes<double> mean;
    Jones<double> polarizer;
 
    BoxMuller* normal;
    double rms;
  };
 
  class mode_decorator : public mode
  {
  protected:
    mode* source;
 
  public:
    mode_decorator (mode* s) { source = s; }
    mode* get_source () { return source; }
 
    void set_Stokes (const Stokes<double>& mean) { source->set_Stokes(mean); }
    Stokes<double> get_Stokes () { return source->get_Stokes(); }
 
    Matrix<4,4,double> get_covariance() const { return source->get_covariance(); }
    Stokes<double> get_mean () const { return source->get_mean(); }
 
    Spinor<double> get_field () { return source->get_field(); }
    BoxMuller* get_normal () { return source->get_normal(); }
    void set_normal (BoxMuller* n) { source->set_normal(n); }
  };
 
  class field_transformer : public mode_decorator
  {
  public:
    field_transformer (mode* s) : mode_decorator (s) { }
 
    Spinor<double> get_field () { return transform( source->get_field() ); }
 
    virtual Spinor<double> transform ( const Spinor<double>& ) = 0;
  };
 
} // namespace epsic
 
#endif // ! defined __epsic_mode_h