JTestRuns_2D Class Reference

#include <JTestRuns_2D.hh>

Inheritance diagram for JTestRuns_2D:

Public Member Functions
	JTestRuns_2D ()
std::istream &	read (std::istream &in)
std::ostream &	write (std::ostream &out)
void	test (TObject *o1, TObject *o2)
virtual void	save (TFile *f, string path)
void	clear ()
JTestResult	JRunsTest (TH1 *h1, TH1 *h2, double threshold, bool logX, bool logY)

Public Attributes
vector< JTestResult >	results

Private Attributes
double	threshold
bool	logX
bool	logY
char	slice

Detailed Description

Definition at line 13 of file JTestRuns_2D.hh.

Constructor & Destructor Documentation

JTestRuns_2D::JTestRuns_2D ( )

inline

Definition at line 19 of file JTestRuns_2D.hh.

{}

Member Function Documentation

std::istream& JTestRuns_2D::read ( std::istream &  in )

inlinevirtual

Implements JTest_t.

Definition at line 27 of file JTestRuns_2D.hh.

                                   { 
    return in >> threshold >> logX >> logY >> slice;
  };

std::ostream& JTestRuns_2D::write ( std::ostream &  out )

inlinevirtual

Implements JTest_t.

Definition at line 37 of file JTestRuns_2D.hh.

                                     {
    for (vector<JTestResult>::const_iterator r = results.begin() ; r != results.end() ; ++r) {
      out << scientific << setprecision(2) << (r->passed ? GREEN : RED) << r->message << endl;
    }
    return out;
  };

void JTestRuns_2D::test ( TObject *  o1, TObject *  o2  )

inlinevirtual

Implements JTest_t.

Definition at line 50 of file JTestRuns_2D.hh.

                                      {    
    
    if (!(dynamic_cast<TH2*>(o1) == NULL) && !(dynamic_cast<TH2*>(o2) == NULL)) {
                  
      TH2D* h1 = dynamic_cast<TH2D*>(o1);
      TH2D* h2 = dynamic_cast<TH2D*>(o2);
      
      if(slice == 'x' || slice == 'X'){
            
        int nSlices1 = h1->GetNbinsX();
        int nSlices2 = h2->GetNbinsX();

        if(nSlices1 != nSlices2)
          ERROR("Histograms with different binning. The objects: " << o1->GetName() << " can not be compared." << endl);

        for (int i=1 ; i<=nSlices1 ; ++i){

          string name = MAKE_STRING(o1->GetName() + to_string("_") + to_string(i));
             
          TH1D* s1 = h1->ProjectionY (name.c_str(),i,i);
          TH1D* s2 = h2->ProjectionY (name.c_str(),i,i);
          
          JTestResult r = JRunsTest(s1, s2, threshold, logX, logY);
          
          results.push_back(r);
        }       
      }else if (slice == 'y' || slice == 'Y'){
           
        int nSlices1 = h1->GetNbinsY();
        int nSlices2 = h2->GetNbinsY();

        if(nSlices1 != nSlices2)
          ERROR("Histograms with different binning. The objects: " << o1->GetName() << " can not be compared." << endl);

        for (int i=1 ; i<=nSlices1 ; ++i){
             
          string name = MAKE_STRING(o1->GetName() + to_string("_") + to_string(i));

          TH1D* s1 = h1->ProjectionX (name.c_str(),i,i);
          TH1D* s2 = h2->ProjectionX (name.c_str(),i,i);
          
          JTestResult r = JRunsTest(s1, s2, threshold, logX, logY);
          results.push_back(r);
        }
      }

    } else if (!(dynamic_cast<TH1*>(o1) == NULL) && !(dynamic_cast<TH1*>(o2) == NULL)) {
      ERROR("For 1D histograms call JChi2_1D: " << o1->GetName() << endl);
    }
  };

virtual void JTestRuns_2D::save ( TFile *  f, string  path  )

inlinevirtual

Implements JTest_t.

Definition at line 106 of file JTestRuns_2D.hh.

                                           {
    
    for (vector<JTestResult>::const_iterator r = results.begin() ; r != results.end() ; ++r){
            
      if (f -> GetDirectory(path.c_str())==0) f->mkdir(path.c_str());
      f->cd(path.c_str());
      r->obj->Write();
    }
  }

void JTestRuns_2D::clear ( )

inlinevirtual

Implements JTest_t.

Definition at line 119 of file JTestRuns_2D.hh.

              {
    results.clear();
  }

JTestResult JTestRuns_t::JRunsTest ( TH1 *  h1, TH1 *  h2, double  threshold, bool  logX, bool  logY  )

inlineinherited

Definition at line 24 of file JTestRuns_t.hh.

                                                                                  {    
    
    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;
  };

Member Data Documentation

double JTestRuns_2D::threshold

private

Definition at line 124 of file JTestRuns_2D.hh.

bool JTestRuns_2D::logX

private

Definition at line 125 of file JTestRuns_2D.hh.

bool JTestRuns_2D::logY

private

Definition at line 126 of file JTestRuns_2D.hh.

char JTestRuns_2D::slice

private

Definition at line 127 of file JTestRuns_2D.hh.

vector<JTestResult> JTest_t::results

inherited

Definition at line 77 of file JTest_t.hh.

---

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

• JTestRuns_2D.hh