PndEmcFadcFilter.h

Go to the documentation of this file.

//-----------------------------------------------------------
//
// Description:
//      FADC FIR (finite impulse response) filter
//  Author Dima Melnychuk
//  (based on original code by A.Mykulyak)
//-----------------------------------------------------------

#ifndef PNDEMCFadcFilter_HXX
#define PNDEMCFadcFilter_HXX

#include <vector>
#include "TObject.h"

class PndEmcAbsPulseshape;

class PndEmcFadcFilter : public TObject {

public:
        enum type {
                arbitrary = 0,
                symmetric = 1,
                antisymmetric = -1
        };

        PndEmcFadcFilter();
        virtual         ~PndEmcFadcFilter();

        virtual void    Filter(const std::vector<Double_t>& in, std::vector<Double_t>& out) const;
        void            SetData(Double_t data[], Int_t i_size, Int_t i_offset);
        Int_t           MinIndex() const;
        Int_t           MaxIndex() const;
        Int_t           Size() const;
        void            Convolute(const PndEmcFadcFilter& filt);
        void            Convolute(const PndEmcFadcFilter& lhs, const PndEmcFadcFilter& rhs);
        void            Clear();
                
        // Different filters 
        //void            SetupIntegrator(Int_t i_width);
        void            SetupMA(Int_t i_width); // Moving average
        void            SetupMWD(Int_t i_width, Double_t tau); // Moving window deconvolution
        void            SetupMatchedFilter(Int_t i_width, PndEmcAbsPulseshape *pulseshape, Double_t sampleRate); // Matched digital filter
        void            SetupBipolarTrapez(Int_t i_rise, Int_t i_flat, Int_t i_width);
        //void            SetupTriangle(Int_t i_rise);
        void            SetupTrapez(Int_t i_rise, Int_t i_flat);
        void            SetupBipolarTriangle(Int_t i_rise);
        void            SetupDifferentiator(Int_t i_lag=0, Int_t i_width=1);
        void            SetupDoubleDifferentiator(Int_t i_npos=1, Int_t i_nneg=1,
                                                                                                                        Int_t i_nzero=0);
        void            SetupPZDifferentiator(Int_t i_lag=0, Double_t d_fac=1.);
        void            SetNormalizeFloating(Double_t d_norm=1.);
        void            SetNormalizeInteger(Int_t i_shift=0);

        std::vector<Double_t> GetWeights() const {return fCoeff;}

private:
        std::vector<Double_t> fCoeff;
        Int_t           fOffset;
        Int_t           fType; // arbitrary, symmetric, antisymmetric
        Bool_t          fIntegerize;
        Double_t        fNormFactor;
        Int_t           fShiftCount;
        
        ClassDef(PndEmcFadcFilter,1)
};

#endif