EvtPolInt.hh

Go to the documentation of this file.

//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of models developed at BES collaboration
//      based on the EvtGen framework.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/BesCopyright
//      Copyright (A) 2006      Ping Rong-Gang, Pang Cai-Ying@IHEP
//
// Module: EvtEvtPloInt.hh
//
// Description: Routine to deal with polynomial interpolation
//
// Modification history:
//
//    Ping R.-G.    December, 2010   Module created
//
//------------------------------------------------------------------------
#ifndef EVTPOLINT_HH
#define EVTPOLINT_HH
 
#include <cmath>
#include <iostream>
#include <string>
#include <vector>
using std::endl;
using std::vector;
 
class EvtPolInt {
public:
  EvtPolInt( vector<double> xi, vector<double> yi, double x ) {
    xx   = x;
    size = xi.size();
    vx.clear();
    vy.clear();
    for ( int i = 0; i < size; i++ )
    {
      vx.push_back( xi[i] );
      vy.push_back( yi[i] );
    }
  }
  virtual ~EvtPolInt() {}
 
  void polynomial(); // polynomial interpolation see: Numerical recipes in C++ page 82
  void ratint();
  vector<double> spline();
  void           splint();
  double         getvalue();
  double         geterror();
 
private:
  vector<double> vx, vy;
  double         xx, value, error;
  int            size;
};
 
#endif