PolynomialFit.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2008 by Paul Demorest
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
#ifndef __Pulsar_PolynomialFit_h
#define __Pulsar_PolynomialFit_h
 
#include "Estimate.h"
#include "Reference.h"
 

class PolynomialFit : public Reference::Able {
 
public:

  PolynomialFit ();

  virtual ~PolynomialFit ();

  void reset();

  void set_degree(int n) { ndeg=n; calcd=false; }

  int get_degree() { return ndeg; }

  void add_data(double x, Estimate<double> y);

  void compute();

  double evaluate(double x);

  double evaluate_deriv(double x);

  double evaluate_integ(double x0, double x1);

  double evaluate_moment(double x0, double x1, int n);

  double get_rchi2();

  Estimate<double> get_coeff(int n);
 
protected:

  std::vector<double> x;

  std::vector< Estimate<double> > y;

  int ndeg;

  double x_avg;

  double chi2;

  bool calcd;

  std::vector< Estimate<double> > coeffs;

  bool check_range(double x);
 
private:
 
};
 
#endif