// #include "$VMCWORKDIR/KF/KFParticleParformance/KFPartEfficiencies.h"

void CalculateEfficincy(TString workdir="data") {

  TStopwatch timer;
  timer.Start();

 
    // ----  Load libraries   -------------------------------------------------
  gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C");
  basiclibs();
  gSystem->Load("libGeoBase");
  gSystem->Load("libParBase");
  gSystem->Load("libBase");
  gSystem->Load("libCbmBase");
  gSystem->Load("libCbmData");
  gSystem->Load("libField");
  gSystem->Load("libGen");
  gSystem->Load("libPassive");
  gSystem->Load("libMvd");
  gSystem->Load("libSts");
  gSystem->Load("libTrd");
  gSystem->Load("libTof");
  gSystem->Load("libEcal");
  gSystem->Load("libGlobal");
  gSystem->Load("libKF");
  gSystem->Load("libL1");
  gSystem->Load("libMinuit2"); // Needs for rich ellipse fitter
  // ------------------------------------------------------------------------


  KFPartEfficiencies fEfficiency;
  
  const int nDirs = 72;
  
  workdir += "/Signal";
  for(int i=0; i<nDirs; i++)
  {
    TString fileName = workdir;
    fileName += i;
    fileName += "/Efficiency.txt";

    std::fstream file(fileName.Data(),fstream::in);
    KFPartEfficiencies eff;
    eff.AddFromFile(fileName);
    fEfficiency += eff;
  }

  fEfficiency.CalcEff();
  fEfficiency.PrintEff();
}