EJones.h

Go to the documentation of this file.

//# EJones.h: Antenna-pattern related calibration term: E Jones
//# Copyright (C) 1996,1997,2000,2001,2002,2003
//# Associated Universities, Inc. Washington DC, USA.
//#
//# This library is free software; you can redistribute it and/or modify it
//# under the terms of the GNU Library General Public License as published by
//# the Free Software Foundation; either version 2 of the License, or (at your
//# option) any later version.
//#
//# This library is distributed in the hope that it will be useful, but WITHOUT
//# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
//# FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Library General Public
//# License for more details.
//#
//# You should have received a copy of the GNU Library General Public License
//# along with this library; if not, write to the Free Software Foundation,
//# Inc., 675 Massachusetts Ave, Cambridge, MA 02139, USA.
//#
//# Correspondence concerning AIPS++ should be adressed as follows:
//#        Internet email: aips2-request@nrao.edu.
//#        Postal address: AIPS++ Project Office
//#                        National Radio Astronomy Observatory
//#                        520 Edgemont Road
//#                        Charlottesville, VA 22903-2475 USA
//#
//#

#ifndef SYNTHESIS_EJONES_H
#define SYNTHESIS_EJONES_H

#include <casa/aips.h>
#include <casa/Containers/Record.h>
#include <casa/BasicSL/Complex.h>
#include <synthesis/MeasurementComponents/SolvableVisCal.h>

namespace casa { //# NAMESPACE CASA - BEGIN

// Forward declaration
class VisEquation;

// **********************************************************
//  EGainCurve
//

class EGainCurve : public SolvableVisJones {
public:

  // Constructor
  EGainCurve(VisSet& vs);
  EGainCurve(casacore::String msname,casacore::Int MSnAnt,casacore::Int MSnSpw);
  EGainCurve(const MSMetaInfoForCal& msmc);

  virtual ~EGainCurve();

  // Return the type enum
  virtual Type type() { return VisCal::E; };

  // Return type name as string
  virtual casacore::String typeName()     { return "EGainCurve"; };
  virtual casacore::String longTypeName() { return "EGainCurve (Gain(elev) corrections)"; };

  // Gain curve coeffs are casacore::Float parameters
  virtual VisCalEnum::VCParType parType() { return VisCalEnum::REAL; };

  // Type of Jones matrix (
  Jones::JonesType jonesType() { return Jones::Diagonal; };

  // Local specialization of setApply to extract gain curves from table
  void setApply(const casacore::Record& applypar);

  // Cal Library version
  void setCallib(const casacore::Record& applypar,const casacore::MeasurementSet& selms);

  // Generate gain curves caltable via specify
  void setSpecify(const casacore::Record& specify);
  void specify(const casacore::Record& specify);

  // Guess (throws error because we don't yet solve for this)
  virtual void guessPar(VisBuffer& vb);

protected:

  // EGainCurve has eight casacore::Float pars per ant (4 per pol)
  virtual casacore::Int nPar() { return 8; };

  // Jones matrix elements are NOT trivial
  virtual casacore::Bool trivialJonesElem() { return false; };

  // Access to z.a. data
  casacore::Vector<casacore::Double>& za() { return za_; };

  // Detect zenith angle for this vb
  virtual void syncMeta(const VisBuffer& vb);
  virtual void syncMeta2(const vi::VisBuffer2& vb);

  // Calculate parameters (in this case, the z.a.)
  virtual void calcPar();

  // Calculate the EGainCurve matrix for all ants
  virtual void calcAllJones();

private:

  // avoid compiler warnings
  using SolvableVisJones::setApply;  

  // Name of the gain curve source table
  casacore::String gainCurveSrc_;

  // The observation start time (for selecting gain curves)
  casacore::Double obstime_;

  // The ant name list (for selecting gain curves)
  casacore::Vector<casacore::String> antnames_;

  // The spw ref frequencies (for selecting gain curves)
  casacore::Vector<casacore::Double> spwfreqs_;

  // The spw band names (for selecting gain curves)
  casacore::Vector<casacore::String> spwbands_;

  // Zenith angles
  casacore::Vector<casacore::Double> za_;

  // Effeciency samples
  casacore::Vector<casacore::Double> eff_;

};


} //# NAMESPACE CASA - END

#endif