AdaptiveSmooth.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2008 by Paul Demorest
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
#ifndef __Pulsar_AdaptiveSmooth_h
#define __Pulsar_AdaptiveSmooth_h
 
#include "Pulsar/Transformation.h"
 
namespace Pulsar {
 
  class Profile;

  class AdaptiveSmooth : public Transformation<Profile> {
 
  public:

    AdaptiveSmooth ();

    virtual ~AdaptiveSmooth ();

    float get_turns () const;

    unsigned get_max_harm() const;

    void set_hold (bool);

    bool get_hold() const;

    void transform(Profile *);

    enum Method { Wiener, Sinc_MSE, Sinc_IC, VonMises };

    void set_method(Method m);
 
  protected:

    std::vector<float> filter;

    int max_harm;

    bool hold;

    float noise_band;

    Method method;
 
  private:

    void compute(const float *fprof, int nbin);

    void compute_wiener(const float *pspec, double sigma2, int nh);

    void compute_sinc_mse(const float *pspec, double sigma2, int nh);

    void compute_sinc_ic(const float *pspec, double sigma2, int nh);

    void compute_lpf(const float *pspec, double sigma2, int nh,
        float (*filter_func)(float,float));

    void compute_vonmises_lpf(const float *pspec, double sigma2, int nh);
 
  };
 
}
 
#endif