//-----------------------------------------------------------
// File and Version Information:
// $Id$
//
// Description:
//      BiCubic spline on 2 B-Spline bases
//
//
// Environment:
//      Software developed for the PANDA Detector at FAIR.
//
// Author List:
//      Felix Boehmer        TUM            (original author)
//
//
//-----------------------------------------------------------

#ifndef BICUBSPLINE_HH
#define BICUBSPLINE_HH

// Base Class Headers ----------------


// Collaborating Class Headers -------
#include "BSpline.h"
#include <vector>


// Collaborating Class Declarations --



class BiCubSpline {

friend class BiCubSplineFitter;

public:

  // Constructors/Destructors ---------
  BiCubSpline(const std::vector<double>* kx,
	      const std::vector<double>* ky,
	      std::vector<std::vector<double>*>* coeffs= NULL);
  ~BiCubSpline();

    
  // Accessors -----------------------

  const std::vector<double>* getKx() {return _kx;}
  const std::vector<double>* getKy() {return _ky;}
  std::vector<std::vector<double>*>* getCoeffs() {return _coeffs;}

  // Modifiers -----------------------
 
  void setKx(const std::vector<double>* kx) {_kx = kx;}
  void setKy(const std::vector<double>* ky) {_ky = ky;}
  void setCoeffs(std::vector<std::vector<double>*>* c) {_coeffs = c;}

  //for ROOT interface:
  //comment: beware there is no possible consistency check since I need
  //to follow TF2's interface restrictions!
  //user is responsible for supplying the right array format
  void setCoeffsByArray(double*);  

  // Operations ----------------------

  double eval(double x, double y);
  

private:

  // Private Data Members ------------
  
  const std::vector<double>* _kx;    //knots in x and y direction
  const std::vector<double>* _ky;
  std::vector<std::vector<double>*>* _coeffs;
  int _xLength, _yLength;
  std::vector<BSpline*> _M;   // M(x)
  std::vector<BSpline*> _N;   // N(y)
  bool _xoutofrange;	          //flag for out of range
  bool _youtofrange;
   

  // Private Methods -----------------

};

#endif

//--------------------------------------------------------------
// $Log$
//--------------------------------------------------------------