#include <HistoCompare.h>

Public Member Functions
TH1 *	Compare (TH1 *h, const TString &hname)

double	GetResult ()

	HistoCompare ()

	HistoCompare (const TString &refFilename)

void	SetChi2Test (bool test=true)

void	SetKGTest (bool test=true)

void	SetReferenceFilename (const TString &filename)

	~HistoCompare ()

Private Attributes
bool	do_nothing_

TFile *	refFile_

TString	refFilename_

TH1 *	refHisto_

TH1 *	resHisto_

double	result_

bool	setChi2Test_

bool	setKGTest_

Detailed Description

class: HistoCompare.h package: Validation/RecoB

author: Victor Bazterra, UIC Francisco Yumiceva, Fermilab (yumic.nosp@m.eva@.nosp@m.fnal..nosp@m.gov)

Definition at line 24 of file HistoCompare.h.

Constructor & Destructor Documentation

HistoCompare::HistoCompare ( )

class: HistoCompare.cc package: Validation/RecoB

author: Victor Bazterra, UIC Francisco Yumiceva, Fermilab (yumic.nosp@m.eva@.nosp@m.fnal..nosp@m.gov)

Definition at line 16 of file HistoCompare.cc.

References do_nothing_, and result_.

                            {
   result_ = 0;
   do_nothing_ = false;
 }

HistoCompare::HistoCompare ( const TString & refFilename )

Definition at line 21 of file HistoCompare.cc.

References gather_cfg::cout, do_nothing_, refFile_, result_, and setChi2Test_.

                                                      {
 
     result_ = 0;
     do_nothing_ = false;
     // open reference file
     refFile_ = new TFile(refFilename);
     if (refFile_->IsZombie()) {
         std::cout << " Error openning file " << refFilename << std::endl;
         std::cout << " we will not compare histograms. " << std::endl;
         do_nothing_ = true;
     } else {
         std::cout << " open file " << refFilename << std::endl;
     }
     
     setChi2Test_ = true;
     
 }

HistoCompare::~HistoCompare ( )

Definition at line 39 of file HistoCompare.cc.

39 {

40 }

Member Function Documentation

TH1 * HistoCompare::Compare	(	TH1 *	h,
		const TString &	hname
	)

Definition at line 43 of file HistoCompare.cc.

References do_nothing_, refFile_, refHisto_, resHisto_, result_, setChi2Test_, and setKGTest_.

                                                        {
 
     if (do_nothing_) return 0;
     
     //std::cout << "histo: " << h->GetName() << " entries: "<< h->GetEntries() << std::endl;
     //std::cout << "hname: " << hname << std::endl;
 
     
     //create a residual histogram
     //if (resHisto_) delete resHisto_;
     resHisto_ = (TH1*) h->Clone(TString(h->GetName())+"_residuals");
     resHisto_->Reset();
 
     // clone input histogram
     TH1 *htemp = (TH1*) h->Clone(TString(h->GetName()));
 
     refFile_->cd();
     //std::cout << "get reference" << std::endl;
     refHisto_ = (TH1*) gDirectory->Get( hname );
     //std::cout << "historef: " << refHisto_->GetName() << " entries: "<< refHisto_->GetEntries() << std::endl;
     //std::cout << "name: " << refHisto_->GetName() << std::endl;
 
 
     // normalize histograms
     htemp->Scale(1./htemp->Integral());
     refHisto_->Scale(1./refHisto_->Integral());
     
     if (setChi2Test_) {
         result_ = refHisto_->Chi2Test(htemp,"UFOF");
     }
     if (setKGTest_) {
         result_ = refHisto_->KolmogorovTest(htemp,"UO");
     }
 
     
     resHisto_->Add( htemp, refHisto_, 1., -1.);
     
     //std::cout << " residual done." << std::endl;
     
     return resHisto_;
     
 }

double HistoCompare::GetResult ( )

inline

Definition at line 54 of file HistoCompare.h.

54 { return result_; };

HistoCompare::result_

double result_

Definition: HistoCompare.h:60

void HistoCompare::SetChi2Test ( bool test = true )

inline

Definition at line 46 of file HistoCompare.h.

References setChi2Test_, and FrontierCondition_GT_autoExpress_cfi::test.

                                      {
         setChi2Test_ = test;
     };

void HistoCompare::SetKGTest ( bool test = true )

inline

Definition at line 50 of file HistoCompare.h.

References setKGTest_, and FrontierCondition_GT_autoExpress_cfi::test.

                                    {
         setKGTest_ = test;
     };

void HistoCompare::SetReferenceFilename ( const TString & filename )

inline

Definition at line 34 of file HistoCompare.h.

References gather_cfg::cout, do_nothing_, lut2db_cfg::filename, refFile_, and refFilename_.

                                                        {
       refFilename_ = filename;
       //if (refFile_) delete refFile_;
       refFile_ = new TFile(refFilename_);
       if (refFile_->IsZombie()) {
           std::cout << " Error openning file " << refFilename_ << std::endl;
           std::cout << " we will not compare histograms. " << std::endl;
           do_nothing_ = true;
       }
       std::cout << " open file" << std::endl;
     };