QOptimumFilter.hh

Go to the documentation of this file.

#ifndef _Q_OPTIMUMFILTER_HH_
#define _Q_OPTIMUMFILTER_HH_
 
#include "QVector.hh"
#include "QVectorC.hh"
#include <vector>
 
namespace Diana { class QRealComplexFFT; };
 
class QOptimumFilter {
 
    public:
        enum JitterMode {
            J_ABSMAX, 
            J_LOCMAX,
            J_FIXED,
            J_NOJITTER,
        };
 
        QOptimumFilter(const Diana::QVector& ap1, const Diana::QVector& an, int maxJitter, bool useDiff = false, bool debugOn = false);
 
        /* @brief destructor */
        virtual ~QOptimumFilter();
 
 
        int GetMaxJitter() const { return fMaxJitter; }
 
        double GetCutOffFrequency() const { return fCutOffFrequency; }
 
        Diana::QVector GetAveragePulse() const { return fAveragePulse1; }
 
        Diana::QVector GetAveragePulseEnergySpectrum() const { return fAveragePulseEnergySpectrum; }
 
        Diana::QVector GetNoisePowerSpectrum() const { return fAverageNoise; }
 
        Diana::QVectorC GetFilterFD() const { return fFilter1; }
 
        Diana::QVector GetFilterTD() const { return fFilterTD; }
 
        Diana::QVector GetFilteredNoise() const { return fFilteredAverageNoise; }
 
        double GetFilteredNoiseRMS() const { return fFilteredRMS1; }
 
        Diana::QVector GetSNRFD() const { return fSNRFreq; }
 
        // Event by event
        QError Filter(const Diana::QVector& p) ;
 
        const Diana::QVector&  GetFiltered() const ;
 
        Diana::QVector  GetFilteredShifted() const ;
 
 
        int  GetShift() const {return fMaxPos;}
 
        QError SetJitter(JitterMode mode, int start = 0);
 
        QError Get(double& jitter, double& chi2, double& amplitude,double& integral, double& left, double& right) const;
        QError GetHighSNRChi2(const double& threshold, double& chi2) const;
        QError GetInterpolated(double& jitter, double& chi2, double& amplitude,double& integral, double& left, double& right) const;
 
        std::vector<double> Get()const;
 
        Diana::QVector GetFitFunction(const double jitter, const double a1) const;
 
        std::vector<Diana::QVector> fDebugVectors,fDebugSpectra;
 
    protected:
        QOptimumFilter(const Diana::QVector& ap1, const Diana::QVector& an);
        Diana::QVector NormalizeVector(const Diana::QVector& vec) const;
        virtual QError ManipulateInputs();
        QError BuildFilter();
        QError fCheckForFilteredSamples;
        QError fCheckForBuiltFilter;
        double fFilteredRMS1;
        double fFilteredAutoCorrelation;
        Diana::QVector fAveragePulse1;
        size_t fSize;
        int fMaxPos;
        int fMaxJitter;
        bool fDebugOn;
        bool fUseDiff;
        Diana::QVectorC fFilter1;
 
 
    private:
 
        QError Parabola(double xa,double xb, double xc, double ya,double yb, double yc,double &min, double& minpos) const;
       Diana::QRealComplexFFT* fTransformer;
        double fCutOffFrequency;
        int fAvgPulseHalfTimeToMax;
        int fFilteredAveragePulseHalfWidth;
        double fFilteredAveragePulseHalfWidthValueLeft;
        double fFilteredAveragePulseHalfWidthValueRight;
 
 
 
        Diana::QVector fAverageNoise;
        Diana::QVector fAverageNoiseForChi2;
        Diana::QVector fAveragePulseEnergySpectrum;
        double fAvgPulseIntegral1;
        std::vector<Diana::QVectorC> fAveragePulse1Shifted;
        Diana::QVector fSNRFreq;
        Diana::QVector fSNRNorm;
 
 
        size_t fEffectiveFilterLength;
        Diana::QVector fFilteredAverageNoise;
        Diana::QVector fFilterTD;
        std::vector<bool> fValidSample;
 
        // event by event;
        Diana::QVector fFiltered1;
        Diana::QVectorC fTransformedWaveform;
        int fOptimumJitter;
 
 
 
 
        // QRootCintTrigger(QOptimumFilter)
 
};
 
#endif