SumDensity.hh

Go to the documentation of this file.

#ifndef SUM_DENSITY
#define SUM_DENSITY

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

#include "TMath.h"

#include <vector>


class SumDensity : public AbsDensity {

  public: 
  

    SumDensity(const char* pdfName, 
                      AbsPhaseSpace* thePhaseSpace, 
                      std::vector<AbsDensity*> &densityComponents, 
                      std::vector<Double_t> &weights);


    SumDensity(const char* pdfName, 
                      AbsPhaseSpace* thePhaseSpace, 
                      AbsDensity* d1, 
                      AbsDensity* d2, 
                      Double_t w1 = 1., 
                      Double_t w2 = 1.);


    SumDensity(const char* pdfName, 
                      AbsPhaseSpace* thePhaseSpace, 
                      AbsDensity* d1, 
                      AbsDensity* d2, 
                      AbsDensity* d3, 
                      Double_t w1 = 1., 
                      Double_t w2 = 1., 
                      Double_t w3 = 1.);


    SumDensity(const char* pdfName, 
                      AbsPhaseSpace* thePhaseSpace, 
                      AbsDensity* d1, 
                      AbsDensity* d2, 
                      AbsDensity* d3, 
                      AbsDensity* d4, 
                      Double_t w1 = 1., 
                      Double_t w2 = 1., 
                      Double_t w3 = 1., 
                      Double_t w4 = 1.);

    virtual ~SumDensity(); 


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


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

    void init(AbsPhaseSpace* thePhaseSpace, 
              std::vector<AbsDensity*> &densityComponents, 
              std::vector<Double_t> &weights); 

    AbsPhaseSpace* m_phaseSpace; 
    
    std::vector<AbsDensity*> m_densityComponents; 
    
    std::vector<Double_t> m_weights; 
    
    UInt_t m_dim; 

}; 

#endif