#include <Comparator.h>
Public Types | |
enum | Mode { NORMAL, SCALE, RATIO, GRAPH, EFF } |
Public Member Functions | |
Comparator () | |
Comparator (const char *file0, const char *dir0, const char *file1, const char *dir1) | |
TDirectory * | dir0 () |
TDirectory * | dir1 () |
void | Draw (const char *key, Mode mode) |
void | Draw (const char *key0, const char *key1, Mode mode) |
void | DrawGaussSigmaOverMeanSlice (const char *key, const char *key2, const int rebinFactor, Mode mode) |
void | DrawGaussSigmaOverMeanXSlice (const char *key, const int rebinFactor, const int binxmin, const int binxmax, const bool cst_binning, Mode mode) |
void | DrawGaussSigmaSlice (const char *key, const int rebinFactor, Mode mode) |
void | DrawGaussSigmaSlice (const char *key, const int rebinFactor, const int binxmin, const int binxmax, const bool cst_binning, Mode mode) |
void | DrawMeanSlice (const char *key, const int rebinFactor, Mode mode) |
void | DrawSigmaSlice (const char *key, const int rebinFactor, Mode mode) |
void | DrawSlice (const char *key, int binxmin, int binxmax, Mode mode) |
TH1 * | h0 () |
TH1 * | h1 () |
TH1 * | Histo (const char *key, unsigned dirIndex) |
const TLegend & | Legend () |
void | SetAxis (int rebin, float xmin, float xmax) |
void | SetAxis (int rebin) |
void | SetDirs (const char *file0, const char *dir0, const char *file1, const char *dir1) |
set the 2 files, and the directory within each file, in which the histograms will be compared | |
void | SetStyles (Style *s0, Style *s1, const char *leg0, const char *leg1) |
Private Member Functions | |
void | Draw (TH1 *h0, TH1 *h1, Mode mode) |
Private Attributes | |
TDirectory * | dir0_ |
TDirectory * | dir1_ |
TFile * | file0_ |
TFile * | file1_ |
TH1 * | h0_ |
TH1 * | h1_ |
TLegend | legend_ |
int | rebin_ |
bool | resetAxis_ |
Style * | s0_ |
Style * | s1_ |
float | xMax_ |
float | xMin_ |
Definition at line 15 of file Comparator.h.
enum Comparator::Mode |
Comparator::Comparator | ( | ) | [inline] |
Definition at line 27 of file Comparator.h.
Comparator::Comparator | ( | const char * | file0, |
const char * | dir0, | ||
const char * | file1, | ||
const char * | dir1 | ||
) | [inline] |
TDirectory* Comparator::dir0 | ( | ) | [inline] |
TDirectory* Comparator::dir1 | ( | ) | [inline] |
void Comparator::Draw | ( | const char * | key, |
Mode | mode | ||
) |
void Comparator::Draw | ( | TH1 * | h0, |
TH1 * | h1, | ||
Mode | mode | ||
) | [private] |
Definition at line 360 of file Comparator.cc.
References ExpressReco_HICollisions_FallBack::cerr, Styles::FormatHisto(), max(), min, RATIO, and runonSM::text.
{ if( !(h0 && h1) ) { cerr<<"invalid histo"<<endl; return; } TH1::AddDirectory( false ); h0_ = (TH1*) h0->Clone( "h0_"); h1_ = (TH1*) h1->Clone( "h1_"); TH1::AddDirectory( true ); // unsetting the title, since the title of projections // is still the title of the 2d histo // and this is better anyway h0_->SetTitle(""); h1_->SetTitle(""); //h0_->SetStats(1); //h1_->SetStats(1); if(mode!=EFF) { if(rebin_>1) { h0_->Rebin( rebin_); h1_->Rebin( rebin_); } if(resetAxis_) { h0_->GetXaxis()->SetRangeUser( xMin_, xMax_); h1_->GetXaxis()->SetRangeUser( xMin_, xMax_); } } if (mode!=GRAPH) { TPaveStats *ptstats = new TPaveStats(0.7385057,0.720339, 0.9396552,0.8792373,"brNDC"); ptstats->SetName("stats"); ptstats->SetBorderSize(1); ptstats->SetLineColor(2); ptstats->SetFillColor(10); ptstats->SetTextAlign(12); ptstats->SetTextColor(2); ptstats->SetOptStat(1111); ptstats->SetOptFit(0); ptstats->Draw(); h0_->GetListOfFunctions()->Add(ptstats); ptstats->SetParent(h0_->GetListOfFunctions()); //std::cout << "FL: h0_->GetMean() = " << h0_->GetMean() << std::endl; //std::cout << "FL: h0_->GetRMS() = " << h0_->GetRMS() << std::endl; //std::cout << "FL: h1_->GetMean() = " << h1_->GetMean() << std::endl; //std::cout << "FL: h1_->GetRMS() = " << h1_->GetRMS() << std::endl; //std::cout << "FL: test2" << std::endl; TPaveStats *ptstats2 = new TPaveStats(0.7399425,0.529661, 0.941092,0.6885593,"brNDC"); ptstats2->SetName("stats"); ptstats2->SetBorderSize(1); ptstats2->SetLineColor(4); ptstats2->SetFillColor(10); ptstats2->SetTextAlign(12); ptstats2->SetTextColor(4); TText *text = ptstats2->AddText("h1_"); text->SetTextSize(0.03654661); std::ostringstream oss3; oss3 << h1_->GetEntries(); const std::string txt_entries="Entries = "+oss3.str(); text = ptstats2->AddText(txt_entries.c_str()); std::ostringstream oss; oss << h1_->GetMean(); const std::string txt_mean="Mean = "+oss.str(); text = ptstats2->AddText(txt_mean.c_str()); std::ostringstream oss2; oss2 << h1_->GetRMS(); const std::string txt_rms="RMS = "+oss2.str(); text = ptstats2->AddText(txt_rms.c_str()); ptstats2->SetOptStat(1111); ptstats2->SetOptFit(0); ptstats2->Draw(); h1_->GetListOfFunctions()->Add(ptstats2); ptstats2->SetParent(h1_->GetListOfFunctions()); } else { TPaveStats *ptstats = new TPaveStats(0.0,0.0, 0.0,0.0,"brNDC"); ptstats->Draw(); h0_->GetListOfFunctions()->Add(ptstats); ptstats->SetParent(h0_->GetListOfFunctions()); } float min=-999.; float max=+999.; switch(mode) { case SCALE: h1_->Scale( h0_->GetEntries()/h1_->GetEntries() ); case NORMAL: if(s0_) Styles::FormatHisto( h0_ , s0_); if(s1_) Styles::FormatHisto( h1_ , s1_); if( h1_->GetMaximum()>h0_->GetMaximum()) { h0_->SetMaximum( h1_->GetMaximum()*1.15 ); } h0_->Draw(); h1_->Draw("same"); break; case EFF: if(s0_) Styles::FormatHisto( h0_ , s0_); if(s1_) Styles::FormatHisto( h1_ , s1_); // rather arbitrary but useful max= h0_->GetMaximum(); if ( h1_->GetMaximum() > max ) max = h1_->GetMaximum(); if(max > 0.8) max = 1; max *= 1.1 ; min= h0_->GetMinimum(); if ( h1_->GetMinimum() < min ) min = h1_->GetMinimum(); if(min > 0.2) min = 0.; min *= 0.8 ; h0_->SetMaximum( max ); h0_->SetMinimum( min ); h0_->Draw("E"); h1_->Draw("Esame"); break; case GRAPH: if(s0_) Styles::FormatHisto( h0_ , s0_); if(s1_) Styles::FormatHisto( h1_ , s1_); if( h1_->GetMaximum()>h0_->GetMaximum()) { h0_->SetMaximum( h1_->GetMaximum()*1.15 ); } if( h1_->GetMinimum()<h0_->GetMinimum()) { h0_->SetMinimum( h1_->GetMinimum()*1.15 ); } h0_->SetMarkerStyle(21); h0_->SetMarkerColor(2); h1_->SetMarkerStyle(21); h1_->SetMarkerColor(4); h0_->Draw("E1"); h1_->Draw("E1same"); break; case RATIO: h0_->Sumw2(); h1_->Sumw2(); h0_->Divide( h1_ ); if(s0_) Styles::FormatHisto( h0_ , s0_); h0_->Draw(); default: break; } }
void Comparator::Draw | ( | const char * | key0, |
const char * | key1, | ||
Mode | mode | ||
) |
Definition at line 304 of file Comparator.cc.
{ TH1* h0=0; TH1* h1=0; if(mode!=EFF) { h0 = Histo( key0, 0); h1 = Histo( key1, 1); } else { h0 = Histo( key0, 0); TH1 * h0b = Histo( key1, 0); h1 = Histo( key0, 1); TH1 * h1b = Histo( key1, 1); if(rebin_>1) { h0->Rebin( rebin_); h1->Rebin( rebin_); h0b->Rebin( rebin_); h1b->Rebin( rebin_); } if(resetAxis_) { h0->GetXaxis()->SetRangeUser( xMin_, xMax_); h1->GetXaxis()->SetRangeUser( xMin_, xMax_); h0b->GetXaxis()->SetRangeUser( xMin_, xMax_); h1b->GetXaxis()->SetRangeUser( xMin_, xMax_); } h0b->Sumw2(); h0->Sumw2(); h0->Divide(h0,h0b,1.,1.,"B"); h1b->Sumw2(); h1->Sumw2(); h1->Divide(h1,h1b,1.,1.,"B"); } Draw( h0, h1, mode); }
void Comparator::DrawGaussSigmaOverMeanSlice | ( | const char * | key, |
const char * | key2, | ||
const int | rebinFactor, | ||
Mode | mode | ||
) |
Definition at line 261 of file Comparator.cc.
References TH2Analyzer::Average(), dir, and TH2Analyzer::SigmaGauss().
{ TDirectory* dir = dir1_; dir->cd(); TH2D *h2_b = (TH2D*) dir->Get(key2); TH2Analyzer TH2Ana_b(h2_b, rebinFactor); TH1D* meanslice=TH2Ana_b.Average(); TH2D *h2 = (TH2D*) dir->Get(key); TH2Analyzer TH2Ana(h2, rebinFactor); TH1D* ha=TH2Ana.SigmaGauss(); ha->Divide(meanslice); dir = dir0_; dir->cd(); TH2D *h2b = (TH2D*) dir->Get(key); TH2Analyzer TH2Anab(h2b, rebinFactor); TH2D *h2b_b = (TH2D*) dir->Get(key2); TH2Analyzer TH2Anab_b(h2b_b, rebinFactor); TH1D* meansliceb=TH2Anab_b.Average(); TH1D* hb=TH2Anab.SigmaGauss(); hb->Divide(meansliceb); Draw(hb,ha,mode); //Draw(meansliceb,meanslice,mode); }
void Comparator::DrawGaussSigmaOverMeanXSlice | ( | const char * | key, |
const int | rebinFactor, | ||
const int | binxmin, | ||
const int | binxmax, | ||
const bool | cst_binning, | ||
Mode | mode | ||
) |
Definition at line 194 of file Comparator.cc.
References dir, fitFunction_f(), TH2Analyzer::MeanX(), TH2Analyzer::RMS(), and TH2Analyzer::SigmaGauss().
{ TDirectory* dir = dir1_; dir->cd(); TH2D *h2 = (TH2D*) dir->Get(key); TH2Analyzer TH2Ana(h2, binxmin, binxmax, rebinFactor, cst_binning); TH1D* hrms=TH2Ana.RMS(); TH1D* meanXslice = TH2Ana.MeanX(); //for( int i=1; i<=meanXslice->GetNbinsX(); ++i) { // std::cout << "meanXslice->GetBinContent(" << i << ") = " // << meanXslice->GetBinContent(i) << std::endl; // std::cout << "meanXslice->GetBinError(" << i << ") = " // << meanXslice->GetBinError(i) << std::endl; //} //Draw(meanXslice,meanXslice,mode); hrms->Divide(meanXslice); TF1 *fitXfcndgssrms3 = new TF1("fitXfcndgssrms3",fitFunction_f,hrms->GetXaxis()->GetBinLowEdge(1),hrms->GetXaxis()->GetBinUpEdge(hrms->GetNbinsX()),4); fitXfcndgssrms3->SetNpx(500); fitXfcndgssrms3->SetLineWidth(3); fitXfcndgssrms3->SetLineStyle(2); fitXfcndgssrms3->SetLineColor(4); hrms->Fit("fitXfcndgssrms3","0R"); TH1D* ha=TH2Ana.SigmaGauss(); ha->Divide(meanXslice); TF1 *fitXfcndgsse3 = new TF1("fitXfcndgsse3",fitFunction_f,ha->GetXaxis()->GetBinLowEdge(1),ha->GetXaxis()->GetBinUpEdge(ha->GetNbinsX()),4); fitXfcndgsse3->SetNpx(500); fitXfcndgsse3->SetLineWidth(3); fitXfcndgsse3->SetLineStyle(1); fitXfcndgsse3->SetLineColor(4); ha->Fit("fitXfcndgsse3","0R"); dir = dir0_; dir->cd(); TH2D *h2b = (TH2D*) dir->Get(key); TH2Analyzer TH2Anab(h2b, binxmin, binxmax, rebinFactor, cst_binning); TH1D* hrmsb=TH2Anab.RMS(); hrmsb->Divide(meanXslice); TF1 *fitXfcndgssrmsb3 = new TF1("fitXfcndgssrmsb3",fitFunction_f,hrmsb->GetXaxis()->GetBinLowEdge(1),hrmsb->GetXaxis()->GetBinUpEdge(hrmsb->GetNbinsX()),4); fitXfcndgssrmsb3->SetNpx(500); fitXfcndgssrmsb3->SetLineWidth(3); fitXfcndgssrmsb3->SetLineStyle(2); fitXfcndgssrmsb3->SetLineColor(2); hrmsb->Fit("fitXfcndgssrmsb3","0R"); TH1D* hb=TH2Anab.SigmaGauss(); hb->Divide(meanXslice); TF1 *fitXfcndgsseb3 = new TF1("fitXfcndgsseb3",fitFunction_f,hb->GetXaxis()->GetBinLowEdge(1),hb->GetXaxis()->GetBinUpEdge(hb->GetNbinsX()),4); fitXfcndgsseb3->SetNpx(500); fitXfcndgsseb3->SetLineWidth(3); fitXfcndgsseb3->SetLineStyle(1); fitXfcndgsseb3->SetLineColor(2); hb->Fit("fitXfcndgsseb3","0R"); Draw(hb,ha,mode); //Draw(hrms,ha,mode); //Draw(ha,ha,mode); fitXfcndgssrms3->Draw("same"); fitXfcndgsse3->Draw("same"); fitXfcndgssrmsb3->Draw("same"); fitXfcndgsseb3->Draw("same"); }
void Comparator::DrawGaussSigmaSlice | ( | const char * | key, |
const int | rebinFactor, | ||
const int | binxmin, | ||
const int | binxmax, | ||
const bool | cst_binning, | ||
Mode | mode | ||
) |
Definition at line 142 of file Comparator.cc.
References dir, fitFunction_f(), TH2Analyzer::RMS(), and TH2Analyzer::SigmaGauss().
{ TDirectory* dir = dir1_; dir->cd(); TH2D *h2 = (TH2D*) dir->Get(key); TH2Analyzer TH2Ana(h2, binxmin, binxmax, rebinFactor, cst_binning); TH1D* hrms=TH2Ana.RMS(); TF1 *fitfcndgssrms3 = new TF1("fitfcndgssrms3",fitFunction_f,hrms->GetXaxis()->GetBinLowEdge(1),hrms->GetXaxis()->GetBinUpEdge(hrms->GetNbinsX()),4); fitfcndgssrms3->SetNpx(500); fitfcndgssrms3->SetLineWidth(3); fitfcndgssrms3->SetLineStyle(2); fitfcndgssrms3->SetLineColor(4); hrms->Fit("fitfcndgssrms3","0R"); TH1D* ha=TH2Ana.SigmaGauss(); TF1 *fitfcndgsse3 = new TF1("fitfcndgsse3",fitFunction_f,hrms->GetXaxis()->GetBinLowEdge(1),hrms->GetXaxis()->GetBinUpEdge(hrms->GetNbinsX()),4); fitfcndgsse3->SetNpx(500); fitfcndgsse3->SetLineWidth(3); fitfcndgsse3->SetLineStyle(1); fitfcndgsse3->SetLineColor(4); ha->Fit("fitfcndgsse3","0R"); dir = dir0_; dir->cd(); TH2D *h2b = (TH2D*) dir->Get(key); TH2Analyzer TH2Anab(h2b, binxmin, binxmax, rebinFactor, cst_binning); TH1D* hrmsb=TH2Anab.RMS(); TF1 *fitfcndgssrmsb3 = new TF1("fitfcndgssrmsb3",fitFunction_f,hrms->GetXaxis()->GetBinLowEdge(1),hrms->GetXaxis()->GetBinUpEdge(hrms->GetNbinsX()),4); fitfcndgssrmsb3->SetNpx(500); fitfcndgssrmsb3->SetLineWidth(3); fitfcndgssrmsb3->SetLineStyle(2); fitfcndgssrmsb3->SetLineColor(2); hrmsb->Fit("fitfcndgssrmsb3","0R"); TH1D* hb=TH2Anab.SigmaGauss(); TF1 *fitfcndgsseb3 = new TF1("fitfcndgsseb3",fitFunction_f,hrms->GetXaxis()->GetBinLowEdge(1),hrms->GetXaxis()->GetBinUpEdge(hrms->GetNbinsX()),4); fitfcndgsseb3->SetNpx(500); fitfcndgsseb3->SetLineWidth(3); fitfcndgsseb3->SetLineStyle(1); fitfcndgsseb3->SetLineColor(2); hb->Fit("fitfcndgsseb3","0R"); Draw(hb,ha,mode); //Draw(hrms,ha,mode); //Draw(ha,ha,mode); fitfcndgssrms3->Draw("same"); fitfcndgsse3->Draw("same"); fitfcndgssrmsb3->Draw("same"); fitfcndgsseb3->Draw("same"); }
void Comparator::DrawGaussSigmaSlice | ( | const char * | key, |
const int | rebinFactor, | ||
Mode | mode | ||
) |
Definition at line 119 of file Comparator.cc.
References dir, and TH2Analyzer::SigmaGauss().
{ TDirectory* dir = dir1_; dir->cd(); TH2D *h2 = (TH2D*) dir->Get(key); TH2Analyzer TH2Ana(h2,rebinFactor); TH1D* ha=TH2Ana.SigmaGauss(); dir = dir0_; dir->cd(); TH2D *h2b = (TH2D*) dir->Get(key); TH2Analyzer TH2Anab(h2b,rebinFactor); TH1D* hb=TH2Anab.SigmaGauss(); Draw(hb,ha,mode); }
void Comparator::DrawMeanSlice | ( | const char * | key, |
const int | rebinFactor, | ||
Mode | mode | ||
) |
Definition at line 85 of file Comparator.cc.
References TH2Analyzer::Average(), and dir.
{ TDirectory* dir = dir1_; dir->cd(); TH2D *h2 = (TH2D*) dir->Get(key); TH2Analyzer TH2Ana(h2,rebinFactor); TH1D* ha=TH2Ana.Average(); dir = dir0_; dir->cd(); TH2D *h2b = (TH2D*) dir->Get(key); TH2Analyzer TH2Anab(h2b,rebinFactor); TH1D* hb=TH2Anab.Average(); Draw(hb,ha,mode); }
void Comparator::DrawSigmaSlice | ( | const char * | key, |
const int | rebinFactor, | ||
Mode | mode | ||
) |
Definition at line 102 of file Comparator.cc.
References dir, and TH2Analyzer::RMS().
{ TDirectory* dir = dir1_; dir->cd(); TH2D *h2 = (TH2D*) dir->Get(key); TH2Analyzer TH2Ana(h2,rebinFactor); TH1D* ha=TH2Ana.RMS(); dir = dir0_; dir->cd(); TH2D *h2b = (TH2D*) dir->Get(key); TH2Analyzer TH2Anab(h2b,rebinFactor); TH1D* hb=TH2Anab.RMS(); Draw(hb,ha,mode); }
void Comparator::DrawSlice | ( | const char * | key, |
int | binxmin, | ||
int | binxmax, | ||
Mode | mode | ||
) |
Definition at line 47 of file Comparator.cc.
References ExpressReco_HICollisions_FallBack::cerr.
{ static int num = 0; ostringstream out0; out0<<"h0_2d_"<<num; ostringstream out1; out1<<"h1_2d_"<<num; num++; string name0 = out0.str(); string name1 = out1.str(); TH1* h0 = Histo( key, 0); TH1* h1 = Histo( key, 1); TH2* h0_2d = dynamic_cast< TH2* >(h0); TH2* h1_2d = dynamic_cast< TH2* >(h1); if(h0_2d->GetNbinsY() == 1 || h1_2d->GetNbinsY() == 1 ) { cerr<<key<<" is not 2D"<<endl; return; } TH1::AddDirectory( false ); TH1D* h0_slice = h0_2d->ProjectionY(name0.c_str(), binxmin, binxmax, ""); TH1D* h1_slice = h1_2d->ProjectionY(name1.c_str(), binxmin, binxmax, ""); TH1::AddDirectory( true ); Draw( h0_slice, h1_slice, mode); }
TH1* Comparator::h0 | ( | ) | [inline] |
TH1* Comparator::h1 | ( | ) | [inline] |
TH1 * Comparator::Histo | ( | const char * | key, |
unsigned | dirIndex | ||
) |
Definition at line 341 of file Comparator.cc.
References ExpressReco_HICollisions_FallBack::cerr, dir, and h.
{ if(dirIndex>1U) { // dirIndex >= 0, since dirIndex is unsigned cerr<<"bad dir index: "<<dirIndex<<endl; return 0; } TDirectory* dir = 0; if(dirIndex == 0) dir = dir0_; if(dirIndex == 1) dir = dir1_; assert( dir ); dir->cd(); TH1* h = (TH1*) dir->Get(key); if(!h) cerr<<"no key "<<key<<" in directory "<<dir->GetName()<<endl; return h; }
const TLegend& Comparator::Legend | ( | ) | [inline] |
void Comparator::SetAxis | ( | int | rebin, |
float | xmin, | ||
float | xmax | ||
) | [inline] |
Definition at line 47 of file Comparator.h.
References rebin_, resetAxis_, xMax_, and xMin_.
{ rebin_ = rebin; xMin_ = xmin; xMax_ = xmax; resetAxis_ = true; }
void Comparator::SetAxis | ( | int | rebin | ) | [inline] |
Definition at line 57 of file Comparator.h.
References rebin_, and resetAxis_.
{ rebin_ = rebin; resetAxis_ = false; }
void Comparator::SetDirs | ( | const char * | file0, |
const char * | dir0, | ||
const char * | file1, | ||
const char * | dir1 | ||
) |
set the 2 files, and the directory within each file, in which the histograms will be compared
Definition at line 18 of file Comparator.cc.
References cmsRelvalreport::exit.
Referenced by Comparator().
TDirectory* Comparator::dir0_ [private] |
Definition at line 113 of file Comparator.h.
Referenced by dir0().
TDirectory* Comparator::dir1_ [private] |
Definition at line 115 of file Comparator.h.
Referenced by dir1().
TFile* Comparator::file0_ [private] |
Definition at line 112 of file Comparator.h.
TFile* Comparator::file1_ [private] |
Definition at line 114 of file Comparator.h.
TH1* Comparator::h0_ [private] |
Definition at line 117 of file Comparator.h.
Referenced by h0().
TH1* Comparator::h1_ [private] |
Definition at line 118 of file Comparator.h.
Referenced by h1().
TLegend Comparator::legend_ [private] |
Definition at line 123 of file Comparator.h.
Referenced by Legend().
int Comparator::rebin_ [private] |
Definition at line 107 of file Comparator.h.
Referenced by SetAxis().
bool Comparator::resetAxis_ [private] |
Definition at line 110 of file Comparator.h.
Referenced by SetAxis().
Style* Comparator::s0_ [private] |
Definition at line 120 of file Comparator.h.
Style* Comparator::s1_ [private] |
Definition at line 121 of file Comparator.h.
float Comparator::xMax_ [private] |
Definition at line 109 of file Comparator.h.
Referenced by SetAxis().
float Comparator::xMin_ [private] |
Definition at line 108 of file Comparator.h.
Referenced by SetAxis().