1 #ifndef __JCOMPAREHISTOGRAMS__JTESTPOISSONLOGLIKELIHOODRATIO__
2 #define __JCOMPAREHISTOGRAMS__JTESTPOISSONLOGLIKELIHOODRATIO__
36 JTest_t(
"Poissonian Log-Likelihood Ratio",
"NLLR")
53 static const double penalty = -1e9;
55 const TH1* h1 =
dynamic_cast<const TH1*
>(o1);
56 const TH1* h2 =
dynamic_cast<const TH1*
>(o2);
58 if (h1 == NULL || h2 == NULL) {
62 if(h1->GetNbinsX() != h2->GetNbinsX() ||
63 h1->GetNbinsY() != h2->GetNbinsY() ||
64 h1->GetNbinsZ() != h2->GetNbinsZ()) {
65 THROW(
JValueOutOfRange,
"JTestPoissonLogLikelihood::test(): Histograms with different binning. The objects: " <<
66 h1->GetName() <<
" and " << h2->GetName() <<
" can not be compared." << endl);
69 TH1* h3 = (TH1*) h1->Clone(h1->GetName() == h2->GetName() ?
75 double LnLratio = 0.0;
77 for (
int i=1 ; i <= h1->GetNbinsX() ; ++i){
78 for (
int j=1 ;
j <= h1->GetNbinsY() ; ++
j){
79 for (
int k=1 ; k <= h1->GetNbinsZ() ; ++k){
83 const double n0 = h1->GetBinContent(i,
j,k);
84 const double n1 = h2->GetBinContent(i,
j,k);
86 if (n0 > 0.0 && n1 > 0.0) {
87 contribution = 2 * (n0 - n1 - n0 * log(n0/n1));
88 }
else if (n0 > 0.0 || n1 > 0.0) {
89 contribution = penalty;
94 h3->SetBinContent(i,
j,k, contribution);
96 LnLratio -= contribution;
101 const bool passed = (LnLratio <
threshold);
105 h3->SetTitle(title.
getTitle().c_str());
123 std::istream&
read(std::istream& in)
override
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
#define MAKE_CSTRING(A)
Make C-string.
#define MAKE_STRING(A)
Make string.
Class dedicated to standardize the title of the graphical objects produced by the JTest_t() derived c...
std::string getTitle() const
Returns a standard string to be used as title of a graphical root object.
Implementation of the Poisson log-likelihood ratio test.
std::istream & read(std::istream &in) override
Read test parameters from input.
JTestPoissonLogLikelihoodRatio()
Default constructor.
double threshold
threshold chi-square to decide if test is passed.
void test(const TObject *o1, const TObject *o2) override
Applies a Poissonian log-likelihood ratio test to the two given histograms.
Interface to read input and write output for TObject tests.
const std::string resultType
test result type
const std::string testName
test name
Auxiliary class to handle file name, ROOT directory and object name.
Exception for accessing a value in a collection that is outside of its range.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
Structure containing the result of the test.