JPrintFit.cc

Go to the documentation of this file.

#include <string>
#include <iomanip>
#include <iostream>
#include <set>

#include "TROOT.h"
#include "TFile.h"
#include "TH1.h"
#include "TGraph.h"
#include "TF1.h"

#include "JTools/JRange.hh"
#include "JGizmo/JRootObjectID.hh"
#include "JGizmo/JGizmoToolkit.hh"

#include "Jeep/JPrint.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"


/**
 * \file
 * Auxiliary program to print ROOT fit results.
 * \author mdejong
 */
int main(int argc, char **argv)
{
  using namespace std;
  using namespace JPP;

  JRootObjectID  inputFile;
  string         fcn;
  set<string>    print;
  int            debug;

  try {

    JParser<> zap("Auxiliary program to print ROOT fit results.");

    zap['f'] = make_field(inputFile,  "histogram, e.g: <file name>:<object name>");
    zap['F'] = make_field(fcn,        "function name")           = "user";
    zap['@'] = make_field(print,      "(parameter)+");
    zap['d'] = make_field(debug)        = 1;

    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }


  TObject* p = getObject(inputFile);

  if (p  == NULL) {
    FATAL("No object at " << inputFile << endl);
  };

  TF1* f1 = NULL;

  if (f1 == NULL && dynamic_cast<TF1*>   (p) != NULL) { f1 = dynamic_cast<TF1*>   (p); }
  if (f1 == NULL && dynamic_cast<TH1*>   (p) != NULL) { f1 = dynamic_cast<TH1*>   (p)->GetFunction(fcn.c_str()); };
  if (f1 == NULL && dynamic_cast<TGraph*>(p) != NULL) { f1 = dynamic_cast<TGraph*>(p)->GetFunction(fcn.c_str()); };

  if (f1 == NULL) { 
    FATAL("No function at " << inputFile << " " << fcn << endl);
  };

  for (int i = 0; i != f1->GetNpar(); ++i) {
    if (print.count(f1->GetParName(i)) != 0) {
      cout << FIXED(12,3) << f1->GetParameter(i) << ' '
           << FIXED(12,3) << f1->GetParError (i) << endl;
    }
  }
}