BinnedDensity.hh

Go to the documentation of this file.

#ifndef BINNED_DENSITY
#define BINNED_DENSITY

#include "AbsDensity.hh"

#include "TMath.h"

#include <vector>


class BinnedDensity : public AbsDensity {

  public: 


    BinnedDensity(const char* pdfName, 
                  AbsPhaseSpace* thePhaseSpace, 
                  AbsDensity* d, 
                  UInt_t bins1, 
                  UInt_t bins2 = 0, 
                  UInt_t bins3 = 0, 
                  UInt_t bins4 = 0, 
                  UInt_t bins5 = 0);


    BinnedDensity(const char* pdfName, 
                  AbsPhaseSpace* thePhaseSpace, 
                  std::vector<UInt_t> &binning, 
                  AbsDensity* d);


    BinnedDensity(const char* pdfName, 
                  AbsPhaseSpace* thePhaseSpace, 
                  const char* fileName, 
                  const Double_t cutoff = 1.e50); 

    virtual ~BinnedDensity(); 


    void setCutoff(Double_t cutoff) { m_cutoff = cutoff; }


    void readFromFile(const char* fileName); 


    void readFromTextFile(const char* fileName); 


    void readFromRootFile(const char* fileName); 


    void writeToFile(const char* fileName); 


    void writeToTextFile(const char* fileName); 


    void writeToRootFile(const char* fileName); 


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


    AbsPhaseSpace* phaseSpace() { return m_phaseSpace; }

  private: 


    void init(AbsPhaseSpace* thePhaseSpace, 
              std::vector<UInt_t> &binning, 
              AbsDensity* d);

    std::vector<Double_t> m_map;

    std::vector<UInt_t> m_binning;

    AbsPhaseSpace* m_phaseSpace; 

    AbsDensity* m_density; 

    Double_t m_cutoff; 

};

#endif