//#ifndef __CINT__

#include <TROOT.h>
#include <TFile.h>
#include <TH1.h>
#include <TH2.h>
#include <TProfile.h>
#include <TNtuple.h>
#include <TRandom.h>
#include <TChain.h>


#include <iostream>
//#include <iostream.h>
//#include <String.h>


//#include <sstream>
#include <math.h>

//#include <fstream.h>
//#include <stdlib.h>
//#include <string.h>
//#include <iomanip.h>

#include <fstream>

#include <string>
#include <iomanip>
#include <errno.h>
#include <sys/types.h>
#include <sys/stat.h>
using namespace std;
//#endif


#include "anal_rec.h"
//static char psfile[40];

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

  std::string directory, file;
  cout <<"Starting main"<<endl;

  // Check command line arguments
  // cout <<" argc "<<argc<<endl;
  switch (argc) {
  case 1 :
    cout << "Directory : "; cin >> directory;
    cout << "File: "; cin >> file;
    break;
  case 2 :
    directory = argv[1];
    cout << "File: ";cin >> file;
    break;
  case 3 :
    directory = argv[1];
    file = argv[2];
    break;
  default :
    cerr << "tusage: " << argv[0] << " directory file" << endl;
    exit (-1);
  }

#if 0
  if (argc != 3) {
    cerr << "tusage: " << argv[0] << " directory file" << endl;
    exit (-1);
  }
#endif

 // cout << "directory=" << directory << " file=" << file  << endl;
  string path = directory+"/"  + file;
 cout << " file to be analyzed "<< path<<endl;
char file_name[80];
 sprintf(file_name,"%s/%s", directory.c_str(),file.c_str());
// cout <<" file_name "<<file_name<<endl;
  // Check file exists and is readable
/*
  struct stat stat_buffer;
  int result = stat (path.c_str(), &stat_buffer);
  if (result != 0) {
    perror (path.c_str());
    exit (1);
  }
*/

  cout <<" before new tchain "<<endl;
  TChain *fchain=new TChain("Global");
   cout <<" path.c_str() "<<path.c_str()<<endl;
//  fchain->Add(path.c_str());
     fchain->Add(file_name);

  // Substitute file suffix
  file.erase (file.length() - 4, 4);
  file += "ps";
 //sprintf(psfile,"%s",file.c_str());
  std::string root_file;
root_file="hist_"+file;


 cout << " plot are written on  file=" << file  << " root_file "<<root_file<<endl;
///  cout <<" after fchain"<<endl;
  cout <<" in main, nb. of entries "<<fchain->GetEntries()<<endl;
  cout <<" call anal_rec "<<endl;
   anal_rec * es_1 =new anal_rec(fchain);
   es_1->set_psfile_name(file);

//  es_1->Init_hist(file);
  //es_1->book_hist();
  cout <<" avant loop "<<endl;
  es_1->Loop();
  cout <<" in main apres loop"<<endl;
 // return 0;
}