// @(#)root/tmva $Id$
// Author: Tancredi Carli, Dominik Dannheim, Alexander Voigt

/**********************************************************************************
 * Project: TMVA - a Root-integrated toolkit for multivariate data analysis       *
 * Package: TMVA                                                                  *
 * Classes: PDEFoamDensityBase                                                    *
 * Web    : http://tmva.sourceforge.net                                           *
 *                                                                                *
 * Description:                                                                   *
 *      Class PDEFoamDensityBase is an Abstract class representing                *
 *      n-dimensional real positive integrand function                            *
 *      The main function is Density() which provides the event density at a      *
 *      given point during the foam build-up (sampling).                          *
 *                                                                                *
 * Authors (alphabetical):                                                        *
 *      Tancredi Carli   - CERN, Switzerland                                      *
 *      Dominik Dannheim - CERN, Switzerland                                      *
 *      S. Jadach        - Institute of Nuclear Physics, Cracow, Poland           *
 *      Alexander Voigt  - TU Dresden, Germany                                    *
 *      Peter Speckmayer - CERN, Switzerland                                      *
 *                                                                                *
 * Copyright (c) 2008, 2010:                                                      *
 *      CERN, Switzerland                                                         *
 *      MPI-K Heidelberg, Germany                                                 *
 *                                                                                *
 * Redistribution and use in source and binary forms, with or without             *
 * modification, are permitted according to the terms listed in LICENSE           *
 * (http://tmva.sourceforge.net/LICENSE)                                          *
 **********************************************************************************/

#ifndef ROOT_TMVA_PDEFoamDensityBase
#define ROOT_TMVA_PDEFoamDensityBase

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

#ifndef ROOT_TMVA_BinarySearchTree
#include "TMVA/BinarySearchTree.h"
#endif
#ifndef ROOT_TMVA_Event
#include "TMVA/Event.h"
#endif
#ifndef ROOT_TMVA_MsgLogger
#include "TMVA/MsgLogger.h"
#endif

namespace TMVA
{

   // class definition of underlying density
   class PDEFoamDensityBase : public ::TObject
   {
   private:
      std::vector<Double_t> fBox; // range-searching box
      Double_t fBoxVolume;        // volume of range searching box
      Bool_t fBoxHasChanged;      // range searching box has changed

   protected:
      BinarySearchTree *fBst;     // Binary tree to find events within a volume
      mutable MsgLogger *fLogger; //! message logger

      MsgLogger& Log() const { return *fLogger; }

      // calculate volume of fBox
      Double_t GetBoxVolume();

   public:
      PDEFoamDensityBase();
      PDEFoamDensityBase(std::vector<Double_t> box);
      PDEFoamDensityBase(const PDEFoamDensityBase&);
      virtual ~PDEFoamDensityBase();

      // fill event into binary search tree
      void FillBinarySearchTree(const Event* ev);

      // set the range-searching box
      void SetBox(std::vector<Double_t> box) { fBox = box; fBoxHasChanged = kTRUE; }

      // get the range-searching box
      const std::vector<Double_t>& GetBox() const { return fBox; }

      // main function used by PDEFoam
      // returns density at a given point by range searching in BST
      virtual Double_t Density(std::vector<Double_t> &Xarg, Double_t &event_density) = 0;

      ClassDef(PDEFoamDensityBase, 1) // PDEFoam event density interface
   };  //end of PDEFoamDensityBase

}  // namespace TMVA

#endif
 PDEFoamDensityBase.h:1
 PDEFoamDensityBase.h:2
 PDEFoamDensityBase.h:3
 PDEFoamDensityBase.h:4
 PDEFoamDensityBase.h:5
 PDEFoamDensityBase.h:6
 PDEFoamDensityBase.h:7
 PDEFoamDensityBase.h:8
 PDEFoamDensityBase.h:9
 PDEFoamDensityBase.h:10
 PDEFoamDensityBase.h:11
 PDEFoamDensityBase.h:12
 PDEFoamDensityBase.h:13
 PDEFoamDensityBase.h:14
 PDEFoamDensityBase.h:15
 PDEFoamDensityBase.h:16
 PDEFoamDensityBase.h:17
 PDEFoamDensityBase.h:18
 PDEFoamDensityBase.h:19
 PDEFoamDensityBase.h:20
 PDEFoamDensityBase.h:21
 PDEFoamDensityBase.h:22
 PDEFoamDensityBase.h:23
 PDEFoamDensityBase.h:24
 PDEFoamDensityBase.h:25
 PDEFoamDensityBase.h:26
 PDEFoamDensityBase.h:27
 PDEFoamDensityBase.h:28
 PDEFoamDensityBase.h:29
 PDEFoamDensityBase.h:30
 PDEFoamDensityBase.h:31
 PDEFoamDensityBase.h:32
 PDEFoamDensityBase.h:33
 PDEFoamDensityBase.h:34
 PDEFoamDensityBase.h:35
 PDEFoamDensityBase.h:36
 PDEFoamDensityBase.h:37
 PDEFoamDensityBase.h:38
 PDEFoamDensityBase.h:39
 PDEFoamDensityBase.h:40
 PDEFoamDensityBase.h:41
 PDEFoamDensityBase.h:42
 PDEFoamDensityBase.h:43
 PDEFoamDensityBase.h:44
 PDEFoamDensityBase.h:45
 PDEFoamDensityBase.h:46
 PDEFoamDensityBase.h:47
 PDEFoamDensityBase.h:48
 PDEFoamDensityBase.h:49
 PDEFoamDensityBase.h:50
 PDEFoamDensityBase.h:51
 PDEFoamDensityBase.h:52
 PDEFoamDensityBase.h:53
 PDEFoamDensityBase.h:54
 PDEFoamDensityBase.h:55
 PDEFoamDensityBase.h:56
 PDEFoamDensityBase.h:57
 PDEFoamDensityBase.h:58
 PDEFoamDensityBase.h:59
 PDEFoamDensityBase.h:60
 PDEFoamDensityBase.h:61
 PDEFoamDensityBase.h:62
 PDEFoamDensityBase.h:63
 PDEFoamDensityBase.h:64
 PDEFoamDensityBase.h:65
 PDEFoamDensityBase.h:66
 PDEFoamDensityBase.h:67
 PDEFoamDensityBase.h:68
 PDEFoamDensityBase.h:69
 PDEFoamDensityBase.h:70
 PDEFoamDensityBase.h:71
 PDEFoamDensityBase.h:72
 PDEFoamDensityBase.h:73
 PDEFoamDensityBase.h:74
 PDEFoamDensityBase.h:75
 PDEFoamDensityBase.h:76
 PDEFoamDensityBase.h:77
 PDEFoamDensityBase.h:78
 PDEFoamDensityBase.h:79
 PDEFoamDensityBase.h:80
 PDEFoamDensityBase.h:81
 PDEFoamDensityBase.h:82
 PDEFoamDensityBase.h:83
 PDEFoamDensityBase.h:84
 PDEFoamDensityBase.h:85
 PDEFoamDensityBase.h:86
 PDEFoamDensityBase.h:87
 PDEFoamDensityBase.h:88
 PDEFoamDensityBase.h:89
 PDEFoamDensityBase.h:90
 PDEFoamDensityBase.h:91
 PDEFoamDensityBase.h:92
 PDEFoamDensityBase.h:93