1 #ifndef ROOT_TKDNodeInfo
2 #define ROOT_TKDNodeInfo
8 template <typename Value> class TVectorT;
9 typedef struct TVectorT<Double_t> TVectorD;
10 template <typename Value> class TMatrixT;
11 typedef class TMatrixT<Double_t> TMatrixD;
12 //class TKDInterpolatorBase;
13 class TKDNodeInfo : public TObject
15 //friend class TKDInterpolatorBase;
17 TKDNodeInfo(Int_t ndim = 0);
18 TKDNodeInfo(const TKDNodeInfo & ref);
19 TKDNodeInfo& operator=(const TKDNodeInfo & ref);
20 virtual ~TKDNodeInfo();
21 Double_t CookPDF(const Double_t *point, Double_t &result, Double_t &error);
22 inline Bool_t Has(const Float_t *p) const;
23 void Print(const Option_t * = "") const;
24 void Store(const TVectorD &par, const TMatrixD &cov);
26 Int_t GetSize() const { return fNDim; }
27 Float_t * Data() { return fData; }
28 Float_t * Val() { return fVal; }
29 TMatrixD * Cov() { return fCov; }
30 TVectorD * Par() { return fPar; }
33 void Build(Int_t ndim);
36 Int_t fNDim; // 3 times data dimension
37 Float_t *fData; //[fNDim] node's data
38 Float_t fVal[2]; // measured value for node
39 TMatrixD *fCov; // interpolator covariance matrix
40 TVectorD *fPar; // interpolator parameters
42 ClassDef(TKDNodeInfo, 1) // node info for interpolator
45 //_____________________________________________________________________
46 Bool_t TKDNodeInfo::Has(const Float_t *p) const
50 for(int id=0; id<ndim; id++) if(p[id]>=fData[ndim+2*id] && p[id]<fData[ndim+2*id+1]) n++;
51 if(n==ndim) return kTRUE;