EVLAConvFunc.h
Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029 #ifndef SYNTHESIS_TRANSFORM2_EVLACONVFUNC_H
00030 #define SYNTHESIS_TRANSFORM2_EVLACONVFUNC_H
00031
00032 #include <images/Images/ImageInterface.h>
00033 #include <synthesis/TransformMachines2/Utils.h>
00034 #include <synthesis/TransformMachines/BeamCalc.h>
00035 #include <synthesis/TransformMachines2/CFStore.h>
00036 #include <synthesis/TransformMachines2/VLACalcIlluminationConvFunc.h>
00037
00038
00039 #include <synthesis/TransformMachines2/ConvolutionFunction.h>
00040 #include <coordinates/Coordinates/DirectionCoordinate.h>
00041 #include <coordinates/Coordinates/SpectralCoordinate.h>
00042 #include <coordinates/Coordinates/StokesCoordinate.h>
00043 #include <lattices/LatticeMath/LatticeFFT.h>
00044 #include <casa/Logging/LogIO.h>
00045 #include <casa/Logging/LogSink.h>
00046 #include <casa/Logging/LogOrigin.h>
00047
00048 namespace casa {
00049 template<class T> class ImageInterface;
00050 template<class T> class Matrix;
00051
00052 namespace refim{
00053 class EVLAConvFunc : public ConvolutionFunction
00054 {
00055 public:
00056
00057
00058 EVLAConvFunc():
00059 ConvolutionFunction(),bandID_p(-1), polMap_p(), feedStokes_p()
00060 {};
00061 ~EVLAConvFunc() {};
00062 EVLAConvFunc& operator=(const EVLAConvFunc& other);
00063 Int getVLABandID(Double& freq,String&telescopeName);
00064 Bool findSupport(Array<Complex>& func, Float& threshold,Int& origin, Int& R);
00065 void makeConvFunction(const ImageInterface<Complex>& image,
00066 const VisBuffer2& vb,
00067 const Int wConvSize,
00068 const Float pa,
00069 const Float dpa,
00070 CFStore& cfs,
00071 CFStore& cfwts, Bool fillCF=True);
00072 int getVisParams(const VisBuffer2& vb);
00073 Int makePBPolnCoords(const VisBuffer2 &vb,
00074 const Vector<Int>& polMap,
00075 const Int& convSize,
00076 const Int& convSampling,
00077 const CoordinateSystem& skyCoord,
00078 const Int& skyNx, const Int& skyNy,
00079 CoordinateSystem& feedCoord,
00080 Vector<Int>& cfStokes);
00081
00082
00083
00084 void setPolMap(const Vector<Int>& polMap);
00085 void setFeedStokes(const Vector<Int>& feedStokes);
00086 private:
00087 Int bandID_p;
00088 Float Diameter_p, Nant_p, HPBW, sigma;
00089
00090 LogIO& logIO() {return logIO_p;}
00091 LogIO logIO_p;
00092 Vector<Int> polMap_p;
00093 Vector<Int> feedStokes_p;
00094 CountedPtr<IlluminationConvFunc> ATerm_p;
00095 };
00096 };
00097 };
00098 #endif