RecoVertex/test/HistoCompare.C

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
0037   Double_t mypv = myoldHisto1->Chi2Test(mynewHisto1,"OU");
0038   std::strstream buf;
0039   std::string value;
0040   buf<<"PV="<<mypv<<std::endl;
0041   buf>>value;
0042
0043   myte->DrawTextNDC(0.2,0.7, value.c_str());
0044
0045   std::cout << "[OVAL] " << myoldHisto1->GetName() << " PV= "<< mypv <<std::endl;
0046   return;
0047
0048 }
0049
0050 HistoCompare::PVCompute(TH2 * oldHisto , TH2 * newHisto , TText * te )
0051 {
0052
0053   myoldHisto2 = oldHisto;
0054   mynewHisto2 = newHisto;
0055   myte = te;
0056
0057   Double_t mypv = myoldHisto2->Chi2Test(mynewHisto2,"OU");
0058   std::strstream buf;
0059   std::string value;
0060   buf<<"PV="<<mypv<<std::endl;
0061   buf>>value;
0062
0063   myte->DrawTextNDC(0.2,0.7, value.c_str());
0064
0065   std::cout << "[OVAL] " << myoldHisto2->GetName() << " PV= " << mypv << std::endl;
0066   return;
0067
0068 }
0069
0070
0071 HistoCompare::PVCompute(TProfile * oldHisto , TProfile * newHisto , TText * te )
0072 {
0073
0074   myoldProfile = oldHisto;
0075   mynewProfile = newHisto;
0076   myte = te;
0077
0078   Double_t mypv = myoldProfile->Chi2Test(mynewProfile,"OU");
0079   std::strstream buf;
0080   std::string value;
0081   buf<<"PV="<<mypv<<std::endl;
0082   buf>>value;
0083
0084   myte->DrawTextNDC(0.2,0.7, value.c_str());
0085
0086   std::cout << "[OVAL] " << myoldProfile->GetName() << " PV = " << mypv << std::endl;
0087   return;
0088
0089 }