/*
 * Macro:   unpack_kratta_simReader
 * Created: 2013-04-30
 * Autor:   sebastian.kupny@uj.edu.pl
 * 
 * Description:
 *   Unpack_kratta_simReader analyse MC Simulation data (simulations results)
 * 
 * 1. Input:
 * Example of runlist file structure:
3
/home/sebastian/Datasrc/MonteCarloSimulations/mcs_csi_p_130204_153746/r3bsim_p_Ekin0100MeV_p0445MeV.root
/home/sebastian/Datasrc/MonteCarloSimulations/mcs_csi_p_130204_153746/r3bsim_p_Ekin0030MeV_p0239MeV.root
/home/sebastian/Datasrc/MonteCarloSimulations/mcs_csi_p_130204_153746/r3bsim_p_Ekin0110MeV_p0467MeV.root

 */ 

void unpack_kratta_simReader( TString A_runlist = "", TString A_descr = "" )
{
  // ----  Load libraries   -------------------------------------------------
  gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C");
  basiclibs();
  //gSystem->Load("libGeoBase");
  //gSystem->Load("libParBase");
  gSystem->Load("libBase");

  // -----Load specific libs for kratta analysis--------------------------
  
  ///Load libs for kratta analysis
  gSystem->Load("libASYEOSData");
  gSystem->Load("libTKratProgressBar");
  gSystem->Load("libTKratSimKapReader");

  gSystem->Load("libTKratSimKapAna1");
  gSystem->Load("libTKratSimKapAna2");
    
  // ----------------------------------------------------------------------
  // Style setup (can be moved to rootlogon)
  gStyle->SetPalette(1);
  
  // ------------------------------------------------------------------------  
  TString InFiles           = "runlist_sim-local_p.dat" ;
  TString outFile           = "simana_0509.root";
  TString outFileWithHist   = "simana_hist_0509.root";
  
  
  if( A_runlist != "" ){
    InFiles = A_runlist;
  }
  if( A_descr != "" ){
    outFile         = Form( "simana_%s.root", A_descr.Data() );
    outFileWithHist = Form( "simana_%s_hist.root", A_descr.Data() );
  } 
  

  FairRunAna *fRun= new FairRunAna();
  fRun->SetOutputFile( outFile );
  
  TKratSimKapReader        *pTKratSimKapReader = new TKratSimKapReader("TKratSimKapReader", InFiles, 2);
  TKratSimKapAna1          *pTKratSimKapAna1   = new TKratSimKapAna1("TKratSimKapAna1", 1 );
  TKratSimKapAna2          *pTKratSimKapAna2   = new TKratSimKapAna2("TKratSimKapAna2", 1 );
  //--------------------------------------------------------------------  
  pTKratSimKapReader->DisableProgressBar();
  pTKratSimKapAna1->SetFileWithHistogramsName( outFileWithHist );
  cout <<  pTKratSimKapAna1->GetFileWithHistogramsName() << endl;
  cout << "InFiles="<<  InFiles << endl;
  cout << "outFile="<< outFile << endl;
  cout << "outFileWithHist="<< outFileWithHist << endl;
  
  //--------------------------------------------------------------------  
  fRun->AddTask( pTKratSimKapReader );
  fRun->AddTask( pTKratSimKapAna1 );
  fRun->AddTask( pTKratSimKapAna2 );
  
  //--------------------------------------------------------------------  
  fRun->Init();
  Long64_t eventss = 15000;
  eventss = pTKratSimKapReader->GetNumberOfEventsInInputFiles();
  cout << "Events in files = " << eventss << endl;

  //--------------------------------------------------------------------  
  cout << "Press any key and enter to start..."  << endl;
  getchar();
  
  fRun->Run(0, eventss );
  
  cout << "Done" << endl;
  
}