// $Id: TabletStrokeList.h 2437 2010-08-15 11:06:50Z matevz $

// Copyright (C) 1999-2008, Matevz Tadel. All rights reserved.
// This file is part of GLED, released under GNU General Public License version 2.
// For the licensing terms see $GLEDSYS/LICENSE or http://www.gnu.org/.

#ifndef Tmp1_TabletStrokeList_H
#define Tmp1_TabletStrokeList_H

#include <Glasses/ZNode.h>
#include <Stones/GTSIsoMakerFunctor.h>

#include <Stones/STabletPoint.h>

#include <TKDTree.h>

class TH1I;
class TH2I;

class TabletStrokeList : public ZNode,
			 public GTSIsoMakerFunctor
{
  MAC_RNR_FRIENDS(TabletStrokeList);

public:
  enum PotentialAlgorithm_e
  {
    PA_StrongestPoint,
    PA_TwoStrongestPoints,
    PA_LinearInterpolation,
    PA_CubicInterpolation
  };

private:
  void _init();

protected:
  struct CCoefs { STabletPoint A, B, C, D; };

  // Parameters for iso-triangulation
  PotentialAlgorithm_e mAlgorithm;     // X{GS} 7 PhonyEnum()
  Double_t             mWidth;         // X{GS} 7 Value(-range=>[  0,  1, 1, 1e6])
  Double_t             mPotentialExp;  // X{GS} 7 Value(-range=>[-10,  0, 1, 1e6])
  Double_t             mPressureAlpha; // X{GS} 7 Value(-range=>[  0, 10, 1, 1e6])
  Double_t             mMaxInterPoint; // X{GS} 7 Value(-range=>[  0, 1 , 1, 1e6])
  Bool_t               bMakeCubHistos; // X{GS} 7 Bool();

  // Internal KD-tree and stuff for fast point search.
  TKDTreeIF             *mKDTree;    //!
  vector<STabletPoint*>  mPointRefs; //!
  vector<Float_t>        mArrX;      //!
  vector<Float_t>        mArrY;      //!
  vector<CCoefs>         mCCoefs;    //!
  Double_t               mSearchRad; //!

  void calculate_qubic_coeffs(const vSTabletPoint_t& vec, Int_t i, CCoefs& c);

  void distance_derivatives(Float_t t, const HPointF& P, const CCoefs& c,
			    Float_t& d1, Float_t& d2);

  bool find_closest_points(Double_t x, Double_t y, Double_t z,
			   Int_t id[2], Double_t vm[2]);

  void get_linear_approx(Double_t x, Double_t y, Double_t z,
			 Int_t id[2], Double_t vm[2],
			 Double_t& value, STabletPoint* point=0);

  void get_cubic_approx (Double_t x, Double_t y, Double_t z,
			 Int_t id[2], Double_t vm[2],
			 Double_t& value, STabletPoint* point=0);

  void add_gradient(HPointD& g, Double_t x, Double_t y, Double_t z,
		    Double_t value, STabletPoint* point);

  TH1I *mHN;             //!
  TH2I *mHCubTimevsTime; //!
  TH2I *mHCubTimevsN;    //!
  TH1I *mHCubTime20;     //!
  TH1I *mHDeriv;         //!
  TH1I *mHDeriv20;       //!

public:
  TabletStrokeList(const Text_t* n="TabletStrokeList", const Text_t* t=0);
  virtual ~TabletStrokeList();

  void MakeKDStuff();  // X{E} 7 MButt(-join=>1)
  void ClearKDStuff(); // X{E} 7 MButt()

  void PrintClose(Float_t x, Float_t y, Float_t rad);

  // --- GTSIsoMakerFunctor ---
  virtual void     GTSIsoBegin(GTSIsoMaker* maker, Double_t iso_value);
  virtual Double_t GTSIsoFunc(Double_t x, Double_t y, Double_t z);
  virtual Double_t GTSIsoGradient(Double_t x, Double_t y, Double_t z, HPointD& g);
  virtual void     GTSIsoEnd();

#include "TabletStrokeList.h7"
  ClassDef(TabletStrokeList, 1);
}; // endclass TabletStrokeList

#endif