QTriggerOptimumFilter.hh

Go to the documentation of this file.

#ifndef _Q_TRIGGEROPTIMUMFILTER_HH_
#define _Q_TRIGGEROPTIMUMFILTER_HH_
 
#include "QTrigger.hh"
#include <vector>
#include "QVector.hh"
#include "QVectorC.hh"
#include "QRealComplexFFT.hh"
 
 
Q_BEGIN_NAMESPACE
 
class QTriggerOptimumFilter : public QTrigger{
    public:
        QTriggerOptimumFilter(const QBaseModule* Reader,const string triglab) : QTrigger(Reader,triglab){optimum_filter_built=false;read_parameters=false;triggerAlgo="OPTIMUMFILTER";}
        QTriggerOptimumFilter(const QBaseModule* Reader) : QTrigger(Reader){optimum_filter_built=false;read_parameters=false;triggerAlgo="OPTIMUMFILTER";}
        
            QTriggerOptimumFilter():QTrigger(){optimum_filter_built=false;read_parameters=false;triggerAlgo="OPTIMUMFILTER";}
        
        ~QTriggerOptimumFilter(){}
 
        void TriggerStream(const vector<double>& stream);
        
        void Reset(){};
 
            bool BuildOptimumFilter();
            
            void SetCFGParameters(string apinput="",string aninput="",int ds=1,double trigLevl=5.,double receed=2.,string APOwner="AveragePulses",string ANOwner="NoiseAvgPowerSpectrum"){
            average_noise=aninput;
                average_pulse=apinput;
                dataset = ds;   
                triggerLevel=trigLevl;
                receed_lvl=receed;
                apowner = APOwner;
            anowner = ANOwner;
                read_parameters=true;
                
        }
        
        void SetCFGParameters(vector<TString> params){
            SetCFGParameters(params[0].Data(),params[1].Data(),params[2].Atoi(),params[3].Atof(),params[4].Atof(),params[5].Data(),params[6].Data());
        }
       
        std::string GetParameters(const bool useDesc=false)const{
            std::ostringstream out;
            out.str("");
            if(useDesc){
                out<<"APFile :      "<<average_pulse<<std::endl;
                out<<"ANFile :      "<<average_noise<<std::endl;
                out<<"Dataset:      "<<dataset<<std::endl;
                out<<"TriggerLevel: "<<triggerLevel<<std::endl;
                out<<"Receed:       "<<receed_lvl<<std::endl;
                out<<"APOwner:      "<<apowner<<std::endl;
                out<<"ANOwner:      "<<anowner<<std::endl;
            }else{
                out<<"[ "<<average_pulse<<" "
                         <<average_noise<<" "
                         <<dataset<<" "
                         <<triggerLevel<<" "
                         <<receed_lvl<<" "
                         <<apowner<<" "
                         <<anowner<<" ]";
            }
            return out.str();
        }
 
 
        void GetParametersFromCFG();
        
            const Diana::QVector& Filter(const QVector& window_stream);
            
            const Diana::QVector& GetAveragePulse(){return ap;}
            const Diana::QVector& GetAveragePulseDoubled(){return ap_d;}
            const Diana::QVector& GetAveragePulseDoubledSmoothed(){return smoothed_ap;}
            const Diana::QVector& GetAveragePulseSmoothing(){return smoothing_ap;}
            const Diana::QVector& GetAveragePulseFiltered(){return apfilt;}
            
            const Diana::QVector& GetFilterTD(){return fFilterTD;}
            const Diana::QVector& GetFilterTDDoubled(){return fFilterTD_d;}
            
            const Diana::QVector& GetAverageNoise(){return an;}
            const Diana::QVector& GetFilteredWindow(){return filt_stream_window;}
            
            const Diana::QVectorC& GetFilterFD(){return fFilter1;}
            const Diana::QVectorC& GetFilterFDDoubled(){return fFilterFD_d;}
            const Diana::QVectorC& GetAveragePulseFD(){return ap_FD;}
            
            double GetTriggerLevel(){return triggerLevel*rms;}
            
            int GetIndexDifference(){return IndexDifference;}
            
            double GetAPMax(){return ap_max_norm;}
 
        ClassDef(QTriggerOptimumFilter,0);
        
        private:
            
            string average_noise;
            string average_pulse;
            //string extraLabel;
        string apowner;
        string anowner;
            double triggerLevel;
            bool read_parameters;
            
            int dataset;
            bool optimum_filter_built;
            
            Diana::QVector filt_stream_window;//
            Diana::QVector fFilterTD_d;//
            Diana::QVector fFilterTD;//
        
        Diana::QVector an;//
 
        Diana::QVector ap;//
        Diana::QVector ap_d;//
        Diana::QVector smoothed_ap;//
        Diana::QVector apfilt;//
        Diana::QVector smoothing_ap;//
 
            Diana::QVectorC fFilterFD_d;
            Diana::QVectorC fFilter1;
            Diana::QVectorC ap_FD;
        
            int IndexDifference;
            double receed_lvl;
 
            Diana::QRealComplexFFT* fTransformer;
            Diana::QRealComplexFFT* fTransformer_d;
            double ap_max_norm;
 
};
Q_END_NAMESPACE
 
#endif