QFFTFir.hh

#ifndef _Q_FFT_FIR_HH_
#define _Q_FFT_FIR_HH_

#include "QVector.hh"
#include "QVectorC.hh"
#include "QError.hh"
#include "QVectorView.hh"
#include "QRealComplexFFT.hh"

using namespace Cuore;

class QFFTFir {

    public:

        enum AlgoType {
            AT_BV,
            AT_BV_Causal,
            AT_SF,
            AT_SF_Infinite,
        };

        QFFTFir(size_t size, double samplFreq, double cutFreq, QFFT::WindowType wt = QFFT::WT_Rect, AlgoType at = AT_BV);

        QError FilterData(const QVector& input, QVector& output, size_t subratio = 1) const; 
        QError FilterDataMinimal(const QVector& input, QVector& output) const; 

        const QVectorC& GetInputFilterFD() const { return fInputFilterFD; }
        const QVector& GetInputFilterTD() const { return fInputFilterTD; }
        double GetInputFilterTD_VPR() const { return fInputFilterTD.GetMin()/fInputFilterTD.GetMax(); }

        const QVectorC& GetFilterFD() const { return fFilterFD; }
        const QVector& GetFilterTD() const { return fFilterTD; }
        
        const QVector&  GetWindow() const { return fWindow; }

        inline const size_t GetFilterSize() const {return fFilterSize;}
        inline const double GetCutFrequency() const {return fCutFreq;}
        inline const double GetSamplingFrequency() const {return fSamplingFreq;}

    private:

        size_t FilterSingleBlock(const QVector& input,QVector& output,size_t subratio, size_t offset) const;
        QVectorC GetInputFilter(size_t size, size_t bin_stop);

        double fCutFreq;
        double fSamplingFreq;
        double fFreqResolution;

        QVectorC fInputFilterFD;
        QVector fInputFilterTD;

        QVector fWindow;

        QVectorC fFilterFD;
        QVector fFilterTD;
        size_t fFilterSize;
        bool fCausal;
        size_t fConvolutionIndex;
        size_t fConvolutionSize;

        mutable QVectorC fFilteredDataFD;
        mutable QVector fFilteredData;
        QRealComplexFFT* fFilterTransformer;

};

#endif

---