/*
Program reading data used in KNOSSOS
*/


#include <iostream>
#include <iomanip>

// For ROOT

#include "TFile.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TF1.h"

void ResidualsInfo(std::string inputlog="/data/FAIRsorf/pandaroot/karavdina/myKnossos/logForRES"){

//  Initialize

  // gROOT->Reset();
  TString outfile= "ResidualsINFO.root";
  TFile* fi = new TFile(outfile,"RECREATE");  
  TH1* hres = new TH1F("hres","residuals",1e3,-0.1,0.1);
//  Load files
 std::ifstream input;
 std::string tmptxt;
 cout<<"File with residuals: "<<inputlog<<endl;
 input.open(inputlog.c_str());
 double resval;
 while(!input.eof()){
   input>>tmptxt;
   //  cout<<"tmptxt = "<<tmptxt<<endl;
   if(tmptxt[0]=='!' && tmptxt[1]!='!'){
     cout<<"tmptxt = "<<tmptxt<<endl;
     input>>resval;
     cout<<"resval = "<<resval<<endl;
     hres->Fill(resval);
   }
 }
 // input.close();
 hres->Write();
 fi->Close();
}