JConvertToPDF1D.cc File Reference

Auxiliary program to convert 1D histograms to PDFs. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <vector>
#include "TROOT.h"
#include "TFile.h"
#include "TKey.h"
#include "TH1.h"
#include "TProfile.h"
#include "TString.h"
#include "TRegexp.h"
#include "JGizmo/JRootObjectID.hh"
#include "JGizmo/JGizmoToolkit.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Detailed Description

Auxiliary program to convert 1D histograms to PDFs.

The option -f corresponds to <file name>:<object name>.

Possible operations (option -O <operation>:

• N
  make PDF;
• W
  divide bin contents by width of bin;
• E
  apply operation also to bin errors;

  Author

  mdejong

Definition in file JConvertToPDF1D.cc.

Function Documentation

main()

int main	(	int	argc,
		char **	argv )

Definition at line 36 of file JConvertToPDF1D.cc.

{
  using namespace std;
  using namespace JPP;
  
  vector<JRootObjectID> inputFile;
  string                outputFile;
  string                option;
  int                   debug;
 
  try {
 
    JParser<> zap("Auxiliary program to convert 1D histograms to PDFs.");
 
    zap['f'] = make_field(inputFile,  "<input file>:<object name>");
    zap['O'] = make_field(option);
    zap['o'] = make_field(outputFile, "ROOT file with histogram(s)") = "pdf.root";
    zap['d'] = make_field(debug)      = 1;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
 
  vector<TH1*> listOfObjects;
 
  for (vector<JRootObjectID>::const_iterator input = inputFile.begin(); input != inputFile.end(); ++input) {
 
    DEBUG("Input: " << *input << endl);
 
    TDirectory* dir = getDirectory(*input);
 
    if (dir == NULL) {
      ERROR("File: " << input->getFullFilename() << " not opened." << endl);
      continue;
    }
 
    const TRegexp regexp(input->getObjectName());
 
    TIter iter(dir->GetListOfKeys());
 
    for (TKey* key; (key = (TKey*) iter.Next()) != NULL; ) {
 
      const TString tag(key->GetName());
 
      DEBUG("Key: " << tag << " match = " << tag.Contains(regexp) << endl);
 
      // option match
 
      if (tag.Contains(regexp) && isTObject(key)) {
 
        TObject* object = key->ReadObj();
        
        TH1*  h1 = dynamic_cast<TH1*>(object);
 
        try { h1 = dynamic_cast<TProfile&>(*object).ProjectionX(); } catch(exception&) {}
        
        if (h1 != NULL) {
          listOfObjects.push_back(h1);
        } else {
          ERROR("Incompatible object " << object->GetName() << endl);
        }
      }
    }
  }
 
  TFile out(outputFile.c_str(), "recreate");
 
  for (vector<TH1*>::iterator h1 = listOfObjects.begin(); h1 != listOfObjects.end(); ++h1) {
 
    convertToPDF(**h1, option);
    
    (*h1)->Write();
  }
  
  out.Write();
  out.Close();
}