JTestZero Class Reference

#include <JTestZero.hh>

Inheritance diagram for JTestZero:

Public Member Functions
	JTestZero ()
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	JTestZero_2D (TH2 *h1, TH2 *h2, double outliersThreshold, double threshold)

Public Attributes
vector< JTestResult >	results

Private Attributes
double	outliersThreshold
	Fraction of bins allowed to fail. More...
double	threshold
	threshold p-value to decide if test is passed. More...

Detailed Description

Definition at line 22 of file JTestZero.hh.

Constructor & Destructor Documentation

JTestZero::JTestZero ( )

inline

Definition at line 31 of file JTestZero.hh.

{}

Member Function Documentation

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

inlinevirtual

Implements JTest_t.

Definition at line 39 of file JTestZero.hh.

                                   { 
    return in >> outliersThreshold >> threshold;
  };

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

inlinevirtual

Implements JTest_t.

Definition at line 49 of file JTestZero.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 JTestZero::test ( TObject *  o1, TObject *  o2  )

inlinevirtual

Implements JTest_t.

Definition at line 62 of file JTestZero.hh.

                                      {    
    
    if (!(dynamic_cast<TH2*>(o1) == NULL) && !(dynamic_cast<TH2*>(o2) == NULL)) {
                  
      TH2D* h1 = dynamic_cast<TH2D*>(o1);
      TH2D* h2 = dynamic_cast<TH2D*>(o2);
      
      JTestResult r = JTestZero_2D(h1, h2, outliersThreshold, threshold);
      results.push_back(r);
    }
  };

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

inlinevirtual

Implements JTest_t.

Definition at line 80 of file JTestZero.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 JTestZero::clear ( )

inlinevirtual

Implements JTest_t.

Definition at line 93 of file JTestZero.hh.

              {
    results.clear();
  }

JTestResult JTestZero_t::JTestZero_2D ( TH2 *  h1, TH2 *  h2, double  outliersThreshold, double  threshold  )

inlineinherited

Definition at line 37 of file JTestZero_t.hh.

                                                                                         {    
      
    int nx1 = h1->GetNbinsX();
    int nx2 = h2->GetNbinsX();      
    int ny1 = h1->GetNbinsY();
    int ny2 = h2->GetNbinsY();

    if(nx1 != nx2 || ny1 != ny2)
      ERROR("Histograms with different binning. The objects: " << h1->GetName() << " can not be compared." << endl);
        
    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())));
    h3->Reset();

    JTestResult r;

    double failures = 0;
    
    for (int i=1 ; i<nx1+1 ; ++i){
      for (int j=1 ; j<ny1+1 ; ++j){
        
        double m = h1 -> GetBinContent(i,j);
        double n = h2 -> GetBinContent(i,j);
        
        double p1 = 1 - ROOT::Math::poisson_cdf(m,1);
        double p2 = 1 - ROOT::Math::poisson_cdf(n,1);   

        cout << "m " <<  m << " pm " << p1 << " n " << n << " pn " << p2 << endl;

        if ((p1 > threshold && p2 < threshold) || 
            (p1 < threshold && p2 > threshold)){
          h3->Fill(i,j);
          failures+=1./(nx1*ny1);
        }        
      }
    }
      
    (failures > outliersThreshold ? r.passed = false : r.passed = true);
    
    r.message = MAKE_STRING(scientific << setprecision(2)<< "Test: " << "Chi2" <<
                            " HA: " << h1->GetUniqueID() << " " << h1->GetName() <<
                            " HB: " << h2->GetUniqueID() << " " << h2->GetName() <<
                            " Threshold: " << outliersThreshold << " Result: " << failures << " " << (r.passed ? "PASSED" : "FAILED"));

    JResultTitle title("Zero", "Failures[%] ", r.passed , failures * 100);
    
    h3->SetTitle(title.getTitle().c_str());
    
    r.obj = h3;

    return r;      
  };

Member Data Documentation

double JTestZero::outliersThreshold

private

Fraction of bins allowed to fail.

Definition at line 98 of file JTestZero.hh.

double JTestZero::threshold

private

threshold p-value to decide if test is passed.

Definition at line 99 of file JTestZero.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:

• JTestZero.hh