#include <JTestSignificance_t.hh>

Inheritance diagram for JTestSignificance_t:

Public Member Functions
	JTestSignificance_t ()
	Default constructor. More...

template<class T >
JTestResult	JSignificanceTest (T h1, T h2, double threshold, double K, std::string testName, std::string parameterName)
	Significance test for histograms. More...

Detailed Description

Author: bofearraigh; rgruiz

Compares two histograms $ H_a $ and $ H_b $ by calculating the average bin significance

$S = \frac{1}{N} \cdot \sum_{i=0}^{N} \frac{|a_i - K \cdot b_i|}{\sqrt{\sigma ^2(a_i) + K^2 \cdot \sigma ^2(b_i)}}$

where:

is the number of bins in and .
and are the contents of bin in and .
is a normalisation parameter. It should be interpreted as the integral of divided by the integral of
$\sigma$ is calculated assuming that the bin contents follow Poisson distributions.

A threshold value for $ S $ is used to tag the test as passed or failed.

Definition at line 29 of file JTestSignificance_t.hh.

Constructor & Destructor Documentation

JTestSignificance_t::JTestSignificance_t ( )

inline

Default constructor.

Definition at line 36 of file JTestSignificance_t.hh.

36 {}

Member Function Documentation

template<class T >

JTestResult JTestSignificance_t::JSignificanceTest	(	T *	h1,
		T *	h2,
		double	threshold,
		double	K,
		std::string	testName,
		std::string	parameterName
	)

inline

Significance test for histograms.

Parameters

h1	First object
h2	Second object
threshold
K	Normalisation parameter
parameterName	Name of the parameter used to test the histograms
testName	Name of the test used to compare the histograms

Returns: Test result.

Definition at line 51 of file JTestSignificance_t.hh.

                                                                                                                        {    
 
     using namespace std;
     using namespace JPP;
     
     if(h1 -> GetNbinsX() != h2 -> GetNbinsX())
       ERROR("Histograms with different bining. The objects: " << h1 -> GetName() << " can not be compared." << endl);
         
 
     T* h3 = (T*)h1->Clone(h1->GetName()==h2->GetName() ?
                                 MAKE_CSTRING(to_string(h1->GetName())) :
                                 MAKE_CSTRING(to_string(h1->GetName()) + "_VS_" + to_string(h2->GetName())));
     
     h3->Add(h2,-1*K);
     
     double S = 0;
 
     for (int i=0 ; i < h1->GetNbinsX() ; ++i){
       for (int j=0 ; j< h1->GetNbinsY() ; ++j){
         for (int k=0 ; k< h1->GetNbinsY() ; ++k){
           
           double a = h1->GetBinContent(i+1,j+1,k+1);
           double b = h2->GetBinContent(i+1,j+1,k+1);
 
           if(a!=0 || b!=0){  
             S += fabs((a - K*b)/sqrt(a + K*K*b));
           }
         }       
       }
     }
 
     S /= (h1->GetNbinsX()*h1->GetNbinsY()*h1->GetNbinsZ());
 
     bool passed;
     
     (S < threshold ? passed = false : passed = true);
     
     JResultTitle title(testName, parameterName, passed, S);
     
     h3->SetTitle(title.getTitle().c_str());
   
     JTestResult r (testName,
                    string (h1->GetDirectory()->GetPath()).append(h1->GetName()),
                    string (h2->GetDirectory()->GetPath()).append(h2->GetName()),
                    h1->GetDirectory()->GetFile()->GetName(),
                    h2->GetDirectory()->GetFile()->GetName(),
                    parameterName, S , threshold, h3, passed);
 
     return r;
   };

The documentation for this class was generated from the following file:

JTestSignificance_t.hh

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation