AbsDensity.hh

Go to the documentation of this file.

#ifndef ABS_DENSITY
#define ABS_DENSITY

#include "AbsPhaseSpace.hh"

#include "TMath.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TNtuple.h"
#include "TRandom3.h"

#include <vector>


class AbsDensity {

  public: 


    AbsDensity(const char* pdfName);
    
    virtual ~AbsDensity(); 
    

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


    virtual AbsPhaseSpace* phaseSpace() = 0;


    void slice(std::vector<Double_t> &x, UInt_t num, TH1F* hist);


    void slice(std::vector<Double_t> &x, UInt_t numx, UInt_t numy, TH2F* hist, Bool_t inPhaseSpace = true);

    double transform(TH1F* hist1, TH1F* hist2, double x); 
    

    void project(TH1F* hist);


    void project(TH2F* hist, Bool_t inPhaseSpace = true);


    void setMajorant(Double_t majorant) { m_majorant = majorant; }
    

    void setMaxTries(UInt_t maxTries) { m_maxTries = maxTries; }
    

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


    void generate(TNtuple* tree, UInt_t numEvents); 


    const char* name(void) { return m_name; }


    void setSeed(UInt_t seed = 0) { m_rnd.SetSeed(seed); }

  protected:

    char m_name[256];

    Double_t m_majorant; 

    UInt_t m_maxTries; 

    TRandom3 m_rnd;
};

#endif