BesselTransformer.hh

Go to the documentation of this file.

#ifndef _BESSEL_TRANSFORMER_HH_
#define _BESSEL_TRANSFORMER_HH_
 
#include "QVector.hh"
#include "QVectorC.hh"
#include "QError.hh"
 
using namespace Diana;
 
class BesselTransformer
{
    public:
        #ifndef __CINT__
        typedef double(*pt2Func)(double*);
        typedef double(*pt2FuncParam)(double* t, double* param);
        #endif
 
        BesselTransformer(double cutFreq, double duration, double samplFreq, double tauRC, bool zeroPad);
 
        const QVectorC& GetResponseFunction() const { return fTransferFunction; }
 
        const QVector& GetDeltaResponse() const { return fDeltaResponse; }
 
        double GetDeltaResponse(double t) const;
 
        double GetExpResponse(double t_0, double tau, double t) const;
 
        QError ConvolveTD(const QVector& pulse, int t_min, QVector& signal) ;
 
        QError ConvolveTD(const QVector& pulse, QVector& signal) {return ConvolveTD(pulse, 0, signal);}
        void ConvolveTD(double (*pulse)(double*)) ;
        #ifndef __CINT__
        void ConvolveTD(pt2FuncParam func,double* param) ;
        void ConvolveFD(pt2FuncParam,double* param) ;
        #endif
        double ConvolvedFuncTD(double* t);
 
        void ConvolveFD(double (*pulse)(double*)) ;
        double ConvolvedFuncFD(double* t);
 
        QError ConvolveFD(const QVector& pulse, QVector& signal) ;
 
    private:
 
        QVector fDeltaResponse;
        QVectorC fTransferFunction;
        QVector fConvolvedVectorTD;
        QVector fConvolvedVectorFD;
        double fCutFreq;
        double fRCFreq;
        double fNormCutFreq;
        double fSamplingFrequency;
        size_t fSize;
        double fDuration;
        // analyitical formula
        size_t fAN;
        double* fAp;
        double* fAq;
        double* fAtheta;
        double* fAphi;
        bool fZeroPad;
        // QRootCintTrigger(BesselTransformer)
 
};
 
#endif