MHO_InterpolateSpectralLinePeak.hh

Go to the documentation of this file.

#ifndef MHO_InterpolateSpectralLinePeak_HH__
#define MHO_InterpolateSpectralLinePeak_HH__
 
#include <cmath>
#include <complex>
 
#include "MHO_ContainerDefinitions.hh"
#include "MHO_Operator.hh"
 
namespace hops
{
 
class MHO_InterpolateSpectralLinePeak: public MHO_Operator
{
    public:
        MHO_InterpolateSpectralLinePeak();
        virtual ~MHO_InterpolateSpectralLinePeak(){};
 
        void SetSpecDRData(const visibility_type* spec_dr_data) { fSpecDRData = spec_dr_data; }
 
        void SetWeights(const weight_type* weights) { fWeights = weights; }
 
        void SetMaxBins(int peak_chan, int peak_dr_bin, int peak_freq_bin);
 
        void SetDRAxis(const delay_rate_axis_type* dr_ax) { fDRAxis.Copy(*dr_ax); }
 
        void SetReferenceFrequency(double ref_freq_mhz) { fRefFreqMHz = ref_freq_mhz; }
 
        virtual bool Initialize() override;
        virtual bool Execute() override;
 
        double GetPeakSkyFrequencyMHz() const { return fPeakSkyFreqMHz; }
 
        double GetDelayRate() const { return fDelayRate; }
 
        double GetFringeRate() const { return fFringeRate; }
 
        double GetFringeAmplitude() const { return fFringeAmp; }
 
        double GetFringePhase() const { return fFringePhase; }
 
        double GetPhaseDelay() const { return fPhaseDelay; }
 
    private:
        const visibility_type* fSpecDRData;
        const weight_type* fWeights;
 
        delay_rate_axis_type fDRAxis;
        double fRefFreqMHz;
        double fTotalSummedWeights;
 
        int fPeakChan;
        int fPeakDRBin;
        int fPeakFreqBin;
 
        double fPeakSkyFreqMHz;
        double fDelayRate;
        double fFringeRate;
        double fFringeAmp;
        double fFringePhase;
        double fPhaseDelay;
 
        double parabolic_offset(double fm1, double f0, double fp1) const;
};
 
} // namespace hops
 
#endif