compareHists.cc File Reference

#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="")

Function Documentation

getAllKeys()

std::vector<std::string> getAllKeys	(	const TDirectory *	fDir,
		const std::string &	fClassName
	)

Definition at line 16 of file compareHists.cc.

References gather_cfg::cout, submitPVResolutionJobs::key, GetRecoTauVFromDQM_MC_cff::next, mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

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;
} 

getAllObjects()

std::vector<std::string> getAllObjects	(	const TDirectory *	fDir,
		const std::string &	fClassName
	)

Definition at line 30 of file compareHists.cc.

References gather_cfg::cout, GetRecoTauVFromDQM_MC_cff::next, getGTfromDQMFile::obj, mps_fire::result, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                         {
  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;
} 

getObject()

TObject* getObject	(	TDirectory *	fDir,
		const std::vector< std::string > &	fObjectName
	)

Definition at line 44 of file compareHists.cc.

References DMR_cfg::cerr, DeadROC_duringRun::dir, mps_fire::i, and mps_fire::result.

Referenced by 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;
}

main()

int main	(	int	argn,
		char *	argv[]
	)

===============================================================================================================================================================================================

---

variant2: for each run define phi-averaged A for normalization channel (Dref,16) and then, divide Rijk on it, i.e. get RRijk




















































































eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=29 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=28

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:





Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:





Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:





Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

Prepare maps of good/bad channels:

Prepare maps of good/bad channels:

Prepare maps of good/bad channels:

Prepare maps of good/bad channels:

---

CALO JETS

---

PF JETS

Definition at line 104 of file compareHists.cc.

References GCPpyPlots::argv, svgfig::canvas(), DMR_cfg::cerr, gather_cfg::cout, L1TDiffHarvesting_cfi::dir1, L1TDiffHarvesting_cfi::dir2, getAllKeys(), getObject(), compareTotals::hist, makeListRunsInFiles::inputFile, InefficientDoubleROC::inputFileName, makeGifHists(), dataCertificationJetMET_cfi::refFileName, mps_fire::result, AlCaHLTBitMon_QueryRunRegistry::string, and runGCPTkAlMap::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;
}

makeGifHists()

double makeGifHists	(	TH1 *	fHist,
		TH1 *	fRefHist,
		TCanvas *	fCanvas,
		const std::string &	fPrefix = ""
	)

Definition at line 63 of file compareHists.cc.

References visDQMUpload::buf, corrVsCorr::filename, Skims_PA_cff::name, MetAnalyzer::pv(), AlCaHLTBitMon_QueryRunRegistry::string, and runGCPTkAlMap::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;
}