JCOMPAREHISTOGRAMS::JTestSignificance_t Class Reference

#include <JTestSignificance_t.hh>

Inheritance diagram for JCOMPAREHISTOGRAMS::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

Definition at line 37 of file JTestSignificance_t.hh.

Constructor & Destructor Documentation

JCOMPAREHISTOGRAMS::JTestSignificance_t::JTestSignificance_t ( )

inline

Default constructor.

Definition at line 44 of file JTestSignificance_t.hh.

{}

Member Function Documentation

template<class T >

JTestResult JCOMPAREHISTOGRAMS::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	Significance 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 59 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->GetNbinsZ() ; ++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,
                     JRootObjectID(MAKE_STRING(h1->GetDirectory()->GetPath() << h1->GetName())),
                     JRootObjectID(MAKE_STRING(h2->GetDirectory()->GetPath() << h1->GetName())),
                     parameterName, S , threshold, h3, passed);

      return r;
    };

---

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

• JTestSignificance_t.hh