QMultiChannelDecorrelator.hh

Go to the documentation of this file.

#ifndef _Q_MULTI_CHANNEL_DECORRELATOR_
#define _Q_MULTI_CHANNEL_DECORRELATOR_
 
#include "QChannelCovariance.hh"
#include "QVectorI.hh"
#include "QRealComplexFFT.hh"
 
class QMultiChannelDecorrelator {
 
    public:
        QMultiChannelDecorrelator(const int chan, const QChannelCovariance* cov);
        virtual ~QMultiChannelDecorrelator();
        
        void SetSideChannels(const std::vector<int>& channels);
 
        const std::vector<int>& GetSideChannels() const { return fSideChannels; }
        const Diana::QVector& GetOriginalPowerSpectrum() const {return fOriginalPowerSpectrum;} ;
 
        const Diana::QVector& GetFilteredPowerSpectrum() const {return fDecorrelatedPowerSpectrum;}  ;
 
        virtual void Filter(const Diana::QVector& input, const std::vector<const Diana::QVector*>& sideInputs, Diana::QVector& output) const;
        virtual void Filter(const std::vector<const Diana::QVector*>& inputs, const std::vector<const Diana::QVector*>& sideInputs, Diana::QVector& output) const;
 
        const std::vector<int>& GetChannels() const { return fChannels; }
 
        void SetScale(double scale) { fScale = scale; } 
 
        double GetOriginalResolution() const { return fOriginalResolution; }
 
        double GetFilteredResolution() const { return fFilteredResolution; }
 
    protected:
        QMultiChannelDecorrelator(const int chan, const QChannelCovariance* cov, bool init);
        QMultiChannelDecorrelator(const std::vector<int>& channels, const QChannelCovariance* cov, bool init);
        virtual void Init(const std::vector<int>& allChannels);
        std::vector<int>::iterator FindChannel(std::vector<int>& channels,const int chan);
        double ComputeResolution(const Diana::QVector& powerSpectrum) const;
 
        std::vector<int> fChannels;
        const QChannelCovariance* fCovarianceMatrix;
        std::vector<int> fAvailableSideChannels;
        Diana::QRealComplexFFT* fTransformer;
 
        std::vector<int> fSideChannels;
        std::vector<Diana::QVectorC> fWeights;
        Diana::QVector fOriginalPowerSpectrum;
        Diana::QVector fDecorrelatedPowerSpectrum;
        double fScale;
        double fOriginalResolution;
        double fFilteredResolution;
    private:
        void Constructor(const std::vector<int>& channels, const QChannelCovariance* cov, bool init);
 
        // QRootCintTrigger(QMultiChannelDecorrelator)
};
 
 
#endif