|
|
|
#include <iostream>#include <string>#include <vector>#include "TROOT.h"#include "TFile.h"#include "TKey.h"#include "TH1F.h"#include "TH2F.h"#include "TCanvas.h"#include "TPostScript.h"#include "TPaveText.h"#include "TStyle.h"#include "TText.h"Go to the source code of this file.
Functions | |
| std::vector< std::string > | getAllKeys (const TDirectory *fDir, const std::string &fClassName) |
| std::vector< std::string > | getAllObjects (const TDirectory *fDir, const std::string &fClassName) |
| TObject * | getObject (TDirectory *fDir, const std::vector< std::string > &fObjectName) |
| int | main (int argn, char *argv[]) |
| double | makeGifHists (TH1 *fHist, TH1 *fRefHist, TCanvas *fCanvas, const std::string &fPrefix="") |
| std::vector<std::string> getAllKeys | ( | const TDirectory * | fDir, |
| const std::string & | fClassName | ||
| ) |
Definition at line 16 of file compareHists.cc.
References gather_cfg::cout, combine::key, and query::result.
Referenced by main().
{
std::cout << "getAllKeys-> " << fDir->GetName() << ", " << fClassName << std::endl;
// fDir->ls();
std::vector<std::string> result;
TIter next (fDir->GetListOfKeys ());
for (TKey* key = 0; (key = (TKey *) next());) {
std::cout << "key from list: " << key->GetName() << '/' << key->GetClassName () << std::endl;
if (fClassName == key->GetClassName ()) {
result.push_back (std::string (key->GetName ()));
}
}
return result;
}
| std::vector<std::string> getAllObjects | ( | const TDirectory * | fDir, |
| const std::string & | fClassName | ||
| ) |
Definition at line 30 of file compareHists.cc.
References gather_cfg::cout, VarParsing::obj, and query::result.
{
std::cout << "getAllObjects-> " << fDir->GetName() << ", " << fClassName << std::endl;
// fDir->ls();
std::vector<std::string> result;
TIter next (fDir->GetList ());
for (TObject* obj = 0; (obj = (TObject *) next());) {
std::cout << "name from list: " << obj->GetName() << '/' << obj->ClassName () << std::endl;
if (fClassName == obj->ClassName ()) {
result.push_back (std::string (obj->GetName ()));
}
}
return result;
}
| TObject* getObject | ( | TDirectory * | fDir, |
| const std::vector< std::string > & | fObjectName | ||
| ) |
Definition at line 44 of file compareHists.cc.
References ExpressReco_HICollisions_FallBack::cerr, dir, i, and query::result.
Referenced by HcalLedAnalyzer::beginJob(), HcalPedestalAnalyzer::endJob(), and main().
{
TObject* result = 0; // nothing so far
TDirectory* dir = fDir;
for (unsigned i = 0; i < fObjectName.size (); ++i) {
dir->GetObject (fObjectName[i].c_str(), result);
if (result) {
if (i < fObjectName.size () - 1) {
dir = (TDirectory*) result;
result = 0;
}
}
else {
std::cerr << "getObject-> Can not find (sub)dir/object " << fObjectName[i] << " in directory " << dir->GetName () << std::endl;
return 0;
}
}
return result;
}
| int main | ( | int | argn, |
| char * | argv[] | ||
| ) |
------------------ CALO JETS ------------------
------------------ PF JETS ------------------
Definition at line 104 of file compareHists.cc.
References svgfig::canvas(), ExpressReco_HICollisions_FallBack::cerr, gather_cfg::cout, getAllKeys(), getObject(), estimatePileup::hist, analyzePatCleaning_cfg::inputFile, makeGifHists(), query::result, and indexGen::title.
{
int result = 0; // OK
std::string inputFileName (argv[1]);
std::string refFileName (argv[2]);
std::string globalTitle = argn > 2 ? argv[3] : "";
std::cout << "Processing file " << inputFileName << std::endl;
TFile* inputFile = TFile::Open (inputFileName.c_str());
TFile* refFile = TFile::Open (refFileName.c_str());
if (inputFile) {
std::cout << "ls for the file:" << std::endl;
inputFile->ls ();
std::vector<std::string> dirName1 = getAllKeys (inputFile, "TDirectory");
for (unsigned idir = 0; idir < dirName1.size(); idir++) {
TDirectory* dir1 = 0;
inputFile->GetObject (dirName1[idir].c_str(), dir1);
if (dir1) {
std::vector<std::string> dirName2 = getAllKeys (dir1, "TDirectory");
for (unsigned idir2 = 0; idir2 < dirName1.size(); ++idir2) {
TDirectory* dir2 = 0;
dir1->GetObject (dirName2[idir2].c_str(), dir2);
if (dir2) {
std::vector<std::string> histKeys = getAllKeys (dir2, "TH1F");
// output
gStyle->SetOptStat (kFALSE);
TCanvas canvas ("Jets","Jets",800,600);
TPostScript ps ((dirName2[idir2]+std::string(".ps")).c_str(), -112);
ps.Range(29.7 , 21.0);
for (unsigned ihist = 0; ihist < histKeys.size (); ++ihist) {
TH1* hist = 0;
dir2->GetObject (histKeys[ihist].c_str(), hist);
if (hist) {
std::vector<std::string> histPathName;
histPathName.push_back (dirName1[idir]);
histPathName.push_back (dirName2[idir2]);
histPathName.push_back (histKeys[ihist]);
TH1* refhist = (TH1*) getObject (refFile, histPathName);
if (refhist) {
std::string title = globalTitle.empty () ? dirName2[idir2] : globalTitle;
double pv = makeGifHists (hist, refhist, &canvas, title);
std::cout << "pv for hist " << dirName1[idir] << '/' << dirName2[idir2] << '/' << histKeys[ihist] << " is " << pv << std::endl;
ps.NewPage();
}
}
else {
std::cerr << "Can not get histogram " << histKeys[ihist] << std::endl;
}
}
}
}
}
else {
std::cerr << "Can not find dir1: " << dirName1[idir] << std::endl;
}
}
}
else {
std::cerr << " Can not open input file " << inputFileName << std::endl;
result = 1;
}
return result;
}
| double makeGifHists | ( | TH1 * | fHist, |
| TH1 * | fRefHist, | ||
| TCanvas * | fCanvas, | ||
| const std::string & | fPrefix = "" |
||
| ) |
Definition at line 63 of file compareHists.cc.
References lut2db_cfg::filename, AlCaRecoCosmics_cfg::name, and indexGen::title.
Referenced by main().
{
double pv = fHist->KolmogorovTest (fRefHist, "OU");
// set style
TPad pad ("pad", "pad", 0, 0, 1, 0.9, 0);
pad.SetLogy ();
pad.Draw();
char buf [1024];
sprintf (buf, "%s: Kolmogorov Test PV = %5.3f", fPrefix.c_str(), pv);
TPaveText title (0.3,0.85,0.95, 0.99, buf);
title.SetFillColor(pv > 0.01 ? 3 : 2);
// TText* t1 = title.AddText (fPrefix.c_str());
sprintf (buf, "Kolmogorov Test PV = %6.4f", pv);
// TText* t2 = title.AddText (buf);
// t2->SetTextSize(0.3);
title.Draw();
pad.cd();
fHist->Sumw2 ();
fHist->Scale (fRefHist->GetSumOfWeights () / fHist->GetSumOfWeights ());
fHist->SetMarkerStyle (21);
fHist->SetMarkerSize (0.7);
fRefHist->SetLineColor (2);
fRefHist->SetFillColor (42);
std::string name = fRefHist->GetTitle ();
int blank = name.rfind (' ');
if (blank >= 0) name.erase (0, blank+1);
fHist->SetXTitle (name.c_str());
fHist->SetTitle ("");
fRefHist->Draw ();
fHist->Draw ("e1p,same");
std::string filename = name + ".gif";
fCanvas->Print (filename.c_str());
fCanvas->Update ();
return pv;
}
1.7.3