QBiComponentOptimumFilter.hh

Go to the documentation of this file.

#ifndef _Q_BICOMPONENTOPTIMUMFILTER_HH_
#define _Q_BICOMPONENTOPTIMUMFILTER_HH_
 
#include "QVector.hh"
#include "QVectorC.hh"
#include <vector>
 
namespace Diana { class QRealComplexFFT; };
 
class QBiComponentOptimumFilter {
 
    public:
        
        QBiComponentOptimumFilter(const Diana::QVector& ap1, const Diana::QVector& ap2, const Diana::QVector& an, size_t maxJitter, bool diff=false);
 
        /* @brief destructor */
        virtual ~QBiComponentOptimumFilter();
 
        double GetFilteredNoiseRMS1() const { return fFilteredRMS1; }
 
        double GetFilteredNoiseRMS2() const { return fFilteredRMS2; }
 
        const std::vector<int>& GetJitters() const { return fJitters;}
 
        // Event by event
        QError Filter(const Diana::QVector& p) ;
 
        void GetFiltered(Diana::QVector& f1, Diana::QVector& f2) const;
 
        const std::vector<double>& GetChiSquareWithJitter() const;
 
        int GetJitterAtMinimum() const;
 
        double GetChiSquareAtMinimium() const;
 
        void GetAmplitude(double& a1, double& a2, double& integral) const;
 
        QError GetInterpolated(double& jitter, double& a1, double& a2,double& integral) const;
        
        Diana::QVector GetFitFunction(const double jitter, const double a1, const double a2) const;
 
    private:
 
        QError fCheck;
        size_t fSize;
        bool fIsDiff;
        int fMaxJitter;
        Diana::QRealComplexFFT* fTransformer;
 
        Diana::QVectorC fFilter1;
        Diana::QVectorC fFilter2;
 
        Diana::QVector fAverageNoise;
        std::vector<int> fJitters;
        std::vector<Diana::QVectorC> fAveragePulse1Shifted;
        std::vector<Diana::QVectorC> fAveragePulse2Shifted;     
        Diana::QVector fAveragePulse1, fAveragePulse2;
 
        // event by event;
        Diana::QVector fFiltered1;
        Diana::QVector fFiltered2;
        std::vector<double> fChiSquareWithJitter;
        int fOptimumJitter;
        double fAvgPulseIntegral1;
        double fAvgPulseIntegral2;
        double fAmplitude1;
        double fAmplitude2;
        double fFilteredRMS1;
        double fFilteredRMS2;
 
 
        Diana::QVector NormalizeVector(const Diana::QVector& vec) const;
        void BuildFilter(const Diana::QVector& ap1, const Diana::QVector& ap2);
        QError ComputeAmplitude(const Diana::QVectorC& transformed);
        std::vector<double> ChiSquare(const Diana::QVectorC& p, const Diana::QVector& ap1Filtered, const Diana::QVector& ap2Filtered) const ;
 
        // QRootCintTrigger(QBiComponentOptimumFilter)
 
};
 
#endif