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_EVLACONVFUNC_H
00030 #define SYNTHESIS_EVLACONVFUNC_H
00031
00032 #include <images/Images/ImageInterface.h>
00033 #include <synthesis/TransformMachines/Utils.h>
00034 #include <synthesis/TransformMachines/BeamCalc.h>
00035 #include <synthesis/TransformMachines/CFStore.h>
00036 #include <synthesis/TransformMachines/VLACalcIlluminationConvFunc.h>
00037
00038
00039 #include <synthesis/TransformMachines/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 class VisBuffer;
00052 class EVLAConvFunc : public ConvolutionFunction
00053 {
00054 public:
00055
00056
00057 EVLAConvFunc():
00058 ConvolutionFunction(),bandID_p(-1), polMap_p(), feedStokes_p()
00059 {};
00060 ~EVLAConvFunc() {};
00061 EVLAConvFunc& operator=(const EVLAConvFunc& other);
00062 Int getVLABandID(Double& freq,String&telescopeName);
00063 Bool findSupport(Array<Complex>& func, Float& threshold,Int& origin, Int& R);
00064 void makeConvFunction(const ImageInterface<Complex>& image,
00065 const VisBuffer& vb,
00066 const Int wConvSize,
00067 const Float pa,
00068 const Float dpa,
00069 CFStore& cfs,
00070 CFStore& cfwts, Bool fillCF=True);
00071 int getVisParams(const VisBuffer& vb);
00072 Int makePBPolnCoords(const VisBuffer&vb,
00073 const Vector<Int>& polMap,
00074 const Int& convSize,
00075 const Int& convSampling,
00076 const CoordinateSystem& skyCoord,
00077 const Int& skyNx, const Int& skyNy,
00078 CoordinateSystem& feedCoord,
00079 Vector<Int>& cfStokes);
00080
00081
00082
00083 void setPolMap(const Vector<Int>& polMap);
00084 void setFeedStokes(const Vector<Int>& feedStokes);
00085 private:
00086 Int bandID_p;
00087 Float Diameter_p, Nant_p, HPBW, sigma;
00088
00089 LogIO& logIO() {return logIO_p;}
00090 LogIO logIO_p;
00091 Vector<Int> polMap_p;
00092 Vector<Int> feedStokes_p;
00093 CountedPtr<IlluminationConvFunc> ATerm_p;
00094 };
00095 };
00096 #endif