MyTUnfoldDensity.h

Go to the documentation of this file.

#ifndef H_MYTUNFOLDDENSITY
#define H_MYTUNFOLDDENSITY
 
#include <TH2.h>
#include <TUnfoldDensity.h>
 
namespace DAS::Unfolding {
 
class MyTUnfoldDensity : public TUnfoldDensity {
 
public:
 
    TH2 * GetA (const char * name)
    {
        auto h = new TH2D(name, "probability matrix", fA->GetNrows(), 0.5, 0.5+fA->GetNrows() ,
                                                      fA->GetNcols(), 0.5, 0.5+fA->GetNcols());
        for (int i = 0; i < fA->GetNrows(); ++i)
        for (int j = 0; j < fA->GetNcols(); ++j) {
            double content = (*fA)[i][j];
            h->SetBinContent(i+1,j+1,content);
        }
        return h;
    }
 
    template<class ...Args> Bool_t MyAddRegularisationCondition (Args ...args)
    {
        return AddRegularisationCondition(args...);
    }
 
    template<class ...Args> MyTUnfoldDensity (Args ...args) :
        TUnfoldDensity(args...)
    {}
};
 
}
#endif