JPrintFit.cc File Reference

Auxiliary program to print ROOT fit results. More...

#include <string>
#include <iomanip>
#include <iostream>
#include <vector>
#include <set>
#include "TROOT.h"
#include "TFile.h"
#include "TH1.h"
#include "TGraph.h"
#include "TGraph2D.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"

Go to the source code of this file.

Functions
std::istream &	operator>> (std::istream &in, TString &object)
	Read formula from input stream. More...
std::ostream &	operator<< (std::ostream &out, const TString &object)
	Write formula to output stream. More...
int	main (int argc, char **argv)

Detailed Description

Auxiliary program to print ROOT fit results.

Author

mdejong

Definition in file JPrintFit.cc.

Function Documentation

std::istream& operator>> ( std::istream &  in, TString &  object  )

inline

Read formula from input stream.

Parameters

in	input stream
object	formula

Returns

input stream

Definition at line 30 of file JPrintFit.cc.

{
  return object.ReadLine(in);
}

std::ostream& operator<< ( std::ostream &  out, const TString &  object  )

inline

Write formula to output stream.

Parameters

out	output stream
object	formula

Returns

output stream

Definition at line 42 of file JPrintFit.cc.

{
  return out << object.Data();
}

int main	(	int	argc,
		char **	argv
	)

Definition at line 53 of file JPrintFit.cc.

{
  using namespace std;
  using namespace JPP;

  JRootObjectID    inputFile;
  string           fcn;
  vector<Double_t> x;
  set<string>      parameter;
  TString          formula;
  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['x'] = make_field(x,          "(x)+, evaluate function at given point")          = JPARSER::initialised();
    zap['@'] = make_field(parameter,  "(parameter)+, e.g. \"p0\"")                       = JPARSER::initialised();
    zap['F'] = make_field(formula,    "ROOT TFormula (may contain method names of TF1)") = JPARSER::initialised();
    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 && dynamic_cast<TGraph2D*>(p) != NULL) {

    for (TIter i = dynamic_cast<TGraph2D*>(p)->GetListOfFunctions(); (f1 = (TF1*) i()) != NULL; ) {
      if (fcn == f1->GetName()) {
        break;
      }
    }
  }

  if (f1 == NULL || fcn != f1->GetName()) {
    FATAL("No function at " << inputFile << " " << fcn << endl);
  };

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

  if (!x.empty()) {
    cout << f1->EvalPar(x.data()) << endl;
  }

  if (formula != "") {
    cout << SCIENTIFIC(15,5) << getResult(formula, f1) << endl;
  }
}