SavitzkyGolay.h

Go to the documentation of this file.

#ifndef OPALX_SAVITZKY_GOLAY_FILTER_HH
#define OPALX_SAVITZKY_GOLAY_FILTER_HH
 
#include "Filters/Filter.h"
 
class SavitzkyGolayFilter : public Filter {
public:
    SavitzkyGolayFilter(int np, int nl, int nr, int m);
    ~SavitzkyGolayFilter() { ; }
 
    void apply(std::vector<double>& histogram);
    void calc_derivative(std::vector<double>& histogram, const double& h);
 
private:
    int NumberPoints_m;
    int NumberPointsLeft_m;
    int NumberPointsRight_m;
    int PolynomialOrder_m;
    std::vector<double> Coefs_m;
    std::vector<double> CoefsDeriv_m;
};
 
void savgol(
        std::vector<double>& c, const int& np, const int& nl, const int& nr, const int& ld,
        const int& m);
void convlv(
        const std::vector<double>& data, const std::vector<double>& respns, const int& isign,
        std::vector<double>& ans);
void ludcmp(std::vector<double>& a, std::vector<int>& indx, double& d);
void lubksb(std::vector<double>& a, std::vector<int>& indx, std::vector<double>& b);
 
#endif  // OPALX_SAVITZKY_GOLAY_FILTER_HH