PolynomialDensity.hh

Go to the documentation of this file.

#ifndef POLYNOMIAL_DENSITY
#define POLYNOMIAL_DENSITY

#include "AbsDensity.hh"
#include "AbsPhaseSpace.hh"
#include "OneDimPhaseSpace.hh"

#include <vector>


class PolynomialDensity : public AbsDensity {

  public: 
  

    PolynomialDensity(const char* pdfName, 
                      OneDimPhaseSpace* thePhaseSpace, 
                      UInt_t maxPower, 
                      TTree* tree, 
                      const char* var, 
                      UInt_t maxEvents = 0);


    PolynomialDensity(const char* pdfName, 
                      AbsPhaseSpace* thePhaseSpace, 
                      UInt_t maxPower, 
                      TTree* tree, 
                      const char* var1, 
                      const char* var2, 
                      UInt_t integEvents, 
                      UInt_t maxEvents = 0);

    virtual ~PolynomialDensity(); 


    Double_t density(std::vector<Double_t> &x);


    AbsPhaseSpace* phaseSpace() { return m_phaseSpace; }
    
  private: 

    void readTuple(TTree* tree, std::vector<TString> &vars, UInt_t maxEvents = 0);
  
    AbsPhaseSpace* m_phaseSpace; 

    UInt_t m_dim; 
    
    Double_t m_par[100]; 
    
    UInt_t m_power; 
}; 

#endif