JTestRuns_t.hh

Go to the documentation of this file.

#ifndef __JTESTRUNS_T__
#define __JTESTRUNS_T__

#include <iostream>
#include "JTest_t.hh"

using namespace std;

/*
 * Runs test.
 */
class JTestRuns_t 
{ 
public:  
  
  JTestRuns_t(){}

  /*
   * Tests the statistical compatibility of two ROOT TObjects
   * 
   * \param h1 First object
   * \param h2 Second object
   */
  JTestResult JRunsTest(TH1* h1, TH1* h2, double threshold, bool logX, bool logY) {    
    
    JTestResult r;
    
    if(h1 -> GetNbinsX() != h2 -> GetNbinsX())
      ERROR("Histograms with different bining. The objects: " << h1 -> GetName() << " can not be compared." << endl);

    double R = h1->Integral();
    double T = h2->Integral();
    
    int    n = 1;
    double p = 0;
    double q = 0;
    
    bool a = ((T/R)*h1->GetBinContent(1) - h2->GetBinContent(1)) < 0;

    (a ? p++ : q++);

    for (int i = 2 ; i<h1->GetNbinsX() ; ++i){
      
      bool b = ((T/R)*h1->GetBinContent(i) - h2->GetBinContent(i)) < 0;

      (b ? p++ : q++);

      if (b != a){
        n++;
        a=b;
      }
    }

    double N = 1 + 2*p*q/(p+q) ;
    double s = sqrt( 2*p*q*(2*p*q-p-q)/(p+q)/(p+q)/(p+q-1) );
    double d = (n-N)/s;

    (fabs(d) > threshold ? r.passed = false : r.passed = true);
    
    r.message = MAKE_STRING(scientific << setprecision(2)<< "Test: " << "Runs" <<
                            " HA: " << h1->GetUniqueID() << " " << h1->GetName() <<
                            " HB: " << h2->GetUniqueID() << " " << h2->GetName() <<
                            " Threshold: " << threshold << " Result: " << fabs(d) << " " << (r.passed ? "PASSED" : "FAILED"));

    TH2D* h3 = (TH2D*)h1->Clone(h1->GetName()==h2->GetName() ?
                                MAKE_CSTRING(to_string(h1->GetName())) :
                                MAKE_CSTRING(to_string(h1->GetName()) + "_VS_" + to_string(h2->GetName())));
    
    for (int i=1 ; i<h1->GetNbinsX() ; ++i) {
      h3->SetBinContent(i , (T/R)*h1->GetBinContent(i) - h2->GetBinContent(i)); 
    }
    
    JResultTitle title("Runs", "Deviation[#sigma] ", r.passed , fabs(d));
    
    h3->SetTitle(title.getTitle().c_str());    
    
    r.obj = h3;
    
    return r;
  };
}; 
  
#endif