JMantisShrimp.cc File Reference

Program to compare root histograms. More...

#include <iostream>
#include <fstream>
#include "TFile.h"
#include "TString.h"
#include "TRegexp.h"
#include "Jeep/JParser.hh"
#include "JSupport/JMeta.hh"
#include "JGizmo/JRootObjectID.hh"
#include "JGizmo/JGizmoToolkit.hh"
#include "JCompareHistograms/JTest_t.hh"
#include "JCompareHistograms/JTestDictionary.hh"

Go to the source code of this file.

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

Detailed Description

Program to compare root histograms.

The input histograms and the test to be applied for each histogram are specified
through the command line.

The histogram names are treated as regular expressions. Therefore, multiple histograms could be tested on one shot.
The command line argument specifies the test to be performed, and the corresponding parameters. See JTestDictionary()

Author

rgruiz

Definition in file JMantisShrimp.cc.

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 78 of file JMantisShrimp.cc.

                                { 

  using namespace JPP;
  using namespace std;

  string steeringFile;
  
  vector<JRootObjectID> input_a;
  vector<JRootObjectID> input_b;
  string test;
  string output;
  string ascii;
  bool onlyFailures;
  
  try {
    
    JParser<> zap("\nProgram to compare root histograms. See the link below this usage for further details.\n");
    
    zap['a'] = make_field(input_a,       "Histogram a");
    zap['b'] = make_field(input_b,       "Histogram b");
    zap['t'] = make_field(test,          "Test parameters");
    zap['o'] = make_field(output,        "output file")      = "out.root";
    zap['t'] = make_field(ascii ,        "output file txt" ) = "";
    zap['w'] = make_field(onlyFailures , "write only failed tests" );
    
    zap(argc,argv);
  }
  catch(const exception &error) {
    ERROR(error.what() << endl);
  }

  JTestDictionary d;

  int testID;
  
  istringstream iss(test);
  iss >> testID;
  d[testID]->read(iss);
  
  vector<TObject*> listOfObjects_a = getListOfObjects(input_a);
  vector<TObject*> listOfObjects_b = getListOfObjects(input_b);

  for(vector<TObject*>::const_iterator object_a = listOfObjects_a.begin() ; object_a != listOfObjects_a.end() ; ++object_a){
    for(vector<TObject*>::const_iterator object_b = listOfObjects_b.begin() ; object_b != listOfObjects_b.end() ; ++object_b){

      d[testID]->test(*object_a,*object_b);
    }
  }
  
  TFile out(output.c_str(),"recreate");
  out.cd();

  d[testID]->write(cout);
  d[testID]->save(&out, "", onlyFailures);

  if(!ascii.empty()){
    
    ofstream results;
    
    results.open (ascii);
    
    d[testID]->write(results, ';', onlyFailures);
    
    results.close();
  }
  
  putObject(&out, JMeta(argc, argv));

  out.Close();
  
  return 0; 
}