comp_pid.C File Reference

#include "TString.h"
#include "TFile.h"
#include "TCanvas.h"
#include "TKey.h"
#include <iostream>
#include "TH1F.h"
#include "TRandom.h"

Go to the source code of this file.

Functions
bool	comp_pid (TString fn="pid_plot.root", TString fn2="pid_target.root", double minP=0.03, int minev=3, int maxfail=3)

Function Documentation

bool comp_pid	(	TString	fn = "pid_plot.root",
		TString	fn2 = "pid_target.root",
		double	minP = 0.03,
		int	minev = 3,
		int	maxfail = 3
	)

Definition at line 12 of file comp_pid.C.

References c1, f, f2, h, h2, name, next, obj, P, and TString.

{
  TCanvas *c1 = new TCanvas();
  c1->Divide(9,7);
        bool fTest=kFALSE;
        TString templateFile = gSystem->Getenv("VMCWORKDIR");
        templateFile += "/macro/run/";
        templateFile += fn2;
        
        TFile *f=new TFile(fn,"READ");
        if (!f->IsZombie())
        {
                TFile *f2=new TFile(templateFile,"READ");
                        
                TKey *key;
                TIter next(f->GetListOfKeys());
                
                int failcount = 0;
                Int_t yy = 0;
                while ( (key = (TKey*)next()) )
                {
                  c1->cd(yy+1);
                        TObject *obj = key->ReadObj();
                        
                        // only check TH1Fs
                        if (!obj->InheritsFrom("TH1F")) continue;
                        
                        TString name = obj->GetName();
                        TH1F* h  = (TH1F*) obj; 
                        TH1F* h2 = (TH1F*) f2->Get(name);
                        h->SetLineColor(2);
                        h->SetLineWidth(2);
                        h2->SetLineWidth(2);
                        if ( h->GetMaximum()< h2->GetMaximum())
                          {
                            h2->Draw(); h->Draw("same");
                          }
                        else
                          {
                            h->Draw(); h2->Draw("same");
                          }
                        if ( h->GetEntries()<minev ) 
                        {
                                cout << "Histogram (almost) empty : " << name << " \"" << h2->GetTitle() << "\":  N = " <<  h->GetEntries() << endl;
                                failcount++;
                        }
                        else  
                        {
                                double P = h2->KolmogorovTest(h);
                                if ( P<minP )
                                {
                                        cout << "Incompatible distribution: " << name << " \"" << h2->GetTitle() << "\":  P = " << P << endl;
                                        failcount++;
                                }
                        }
                        yy++;
                }
                
                if (failcount<maxfail) fTest = kTRUE;
                
        }
        
        if (fTest){
                cout << " Test passed" << endl;
                cout << " All ok " << endl;  
        }else{
                cout << " Test Failed" << endl;
                cout << " Not Ok " << endl;         
        }
    
   
}