Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:33:42

0001 #include <iostream.h>
0002 
0003 class HistoCompare {
0004 
0005  public:
0006 
0007   HistoCompare() { std::cout << "Initializing HistoCompare... " << std::endl; } ;
0008 
0009   void PVCompute(TH1 * oldHisto , TH1 * newHisto , TText * te );
0010   void PVCompute(TH2 * oldHisto , TH2 * newHisto , TText * te );
0011   void PVCompute(TProfile * oldHisto , TProfile * newHisto , TText * te );
0012 
0013  private:
0014   
0015   Double_t mypv;
0016 
0017   TH1 * myoldHisto1;
0018   TH1 * mynewHisto1;
0019 
0020   TH2 * myoldHisto2;
0021   TH2 * mynewHisto2;
0022 
0023   TProfile * myoldProfile;
0024   TProfile * mynewProfile;
0025 
0026   TText * myte;
0027 
0028 };
0029 
0030 HistoCompare::PVCompute(TH1 * oldHisto , TH1 * newHisto , TText * te )
0031 {
0032 
0033   myoldHisto1 = oldHisto;
0034   mynewHisto1 = newHisto;
0035   myte = te;
0036  // myoldHisto1->Sumw2();
0037  // mynewHisto1->Sumw2();
0038 
0039   Double_t mypv = myoldHisto1->KolmogorovTest(mynewHisto1);
0040   std::strstream buf;
0041   std::string value;
0042   buf<<"PV="<<mypv<<std::endl;
0043   buf>>value;
0044   
0045   myte->DrawTextNDC(0.2,0.7, value.c_str());
0046 
0047   std::cout << "[OVAL] " << myoldHisto1->GetName() << " PV = " << mypv << std::endl;
0048   return;
0049 
0050 }
0051 
0052 HistoCompare::PVCompute(TH2 * oldHisto , TH2 * newHisto , TText * te )
0053 {
0054 
0055   myoldHisto2 = oldHisto;
0056   mynewHisto2 = newHisto;
0057   myte = te;
0058 
0059   Double_t mypv = myoldHisto2->KolmogorovTest(mynewHisto2);
0060   std::strstream buf;
0061   std::string value;
0062   buf<<"PV="<<mypv<<std::endl;
0063   buf>>value;
0064   
0065   myte->DrawTextNDC(0.2,0.7, value.c_str());
0066 
0067   std::cout << "[OVAL] " << myoldHisto2->GetName() << " PV = " << mypv << std::endl;
0068   return;
0069 
0070 }
0071 
0072 
0073 HistoCompare::PVCompute(TProfile * oldHisto , TProfile * newHisto , TText * te )
0074 {
0075 
0076   myoldProfile = oldHisto;
0077   mynewProfile = newHisto;
0078   myte = te;
0079 
0080   Double_t mypv = myoldProfile->KolmogorovTest(mynewProfile);
0081   std::strstream buf;
0082   std::string value;
0083   buf<<"PV="<<mypv<<std::endl;
0084   buf>>value;
0085   
0086   myte->DrawTextNDC(0.2,0.7, value.c_str());
0087 
0088   std::cout << "[OVAL] " << myoldProfile->GetName() << " PV = " << mypv << std::endl;
0089   return;
0090 
0091 }