sending trigger decisions (true and false) for every event

[u/mrichter/AliRoot.git] / STAT / TKDNodeInfo.h

#ifndef ROOT_TKDNodeInfo
#define ROOT_TKDNodeInfo

#ifndef ROOT_TObject
#include "TObject.h"
#endif

#ifndef ROOT_TBox
#include "TBox.h"
#endif

#ifndef ROOT_TMarker
#include "TMarker.h"
#endif

template <typename Value> class TVectorT;
typedef struct TVectorT<Double_t> TVectorD;
template <typename Value> class TMatrixT;
typedef class TMatrixT<Double_t> TMatrixD;
class TKDNodeInfo : public TObject
{
  //friend class TKDInterpolatorBase;
public:
  class TKDNodeDraw : public TBox {
  public:
    TKDNodeDraw();
    ~TKDNodeDraw() {;}
    void  Draw(Option_t* option = "");
    void  SetNode(TKDNodeInfo*, UChar_t s, UChar_t ax1=0, UChar_t ax2=1);
  private:
    TKDNodeDraw(const TKDNodeDraw & ref);
    TKDNodeDraw& operator=(const TKDNodeDraw & ref);

    TMarker  fCOG;      // COG of the node

    ClassDef(TKDNodeDraw, 1)   // graphical representation of TKDNodeInfo
  };

  TKDNodeInfo(Int_t ndim = 0);
  TKDNodeInfo(const TKDNodeInfo & ref);
  TKDNodeInfo& operator=(const TKDNodeInfo & ref);
  virtual  ~TKDNodeInfo();
  Double_t  CookPDF(const Double_t *point, Double_t &result, Double_t &error);
  inline Bool_t    Has(const Float_t *p) const;
  void      Print(const Option_t * = "") const;
  void      Store(const TVectorD &par, const TMatrixD &cov);

  Int_t GetSize() const { return fNDim; }
  Float_t *  Data() { return fData; } 
  Float_t *  Val() { return fVal; } 
  TMatrixD * Cov() { return fCov; }
  TVectorD * Par() { return fPar; }

protected:
  void      Build(Int_t ndim);

private:
  Int_t     fNDim;          // 3 times data dimension
  Float_t   *fData;         //[fNDim] node's data
  Float_t   fVal[2];        // measured value for node
  TMatrixD  *fCov;          // interpolator covariance matrix
  TVectorD  *fPar;          // interpolator parameters

  ClassDef(TKDNodeInfo, 1)  // node info for interpolator
};

//_____________________________________________________________________
Bool_t TKDNodeInfo::Has(const Float_t *p) const
{
  Int_t n = 0;
  Int_t ndim = fNDim/3;
  for(int id=0; id<ndim; id++) if(p[id]>=fData[ndim+2*id] && p[id]<fData[ndim+2*id+1]) n++;
  if(n==ndim) return kTRUE;
  return kFALSE;
}


#endif