void kratta_calib_reader()
{
  // ----  Load libraries   --------------------------------------------
  gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C");
  basiclibs();
  gSystem->Load("libGeoBase");
  gSystem->Load("libParBase");
  gSystem->Load("libBase");
  
  // -------------------------------------------------------------------
  
  gSystem->Load("libTRootEvent");
  gSystem->Load("libTRootKRATEvent");
  gSystem->Load("libTKratCalibReader");
  gSystem->Load("libTKratCalibAna");
  
  // -------------------------------------------------------------------
  TString InFiles = "runlist_calib.dat" ;
  TString OutFile = "results_kratta_calib_reader.root";

  FairRunAna *fRun= new FairRunAna();
  fRun->SetOutputFile( OutFile );
  
  TKratCalibReader   *pTKratCalibReader  = new TKratCalibReader("TKratCalibReader", InFiles, 0); // 1 - no output, 6 - info about events
  TKratCalibAna      *pTKratCalibAna     = new TKratCalibAna( "TKratCalibAna", 2 );

  // ------ Task configuration -----------------------------------------
  pTKratCalibAna->SetModuleToAnalyse(7);

  //pTKratCalibReader->SetVerboseLevel(6);
  
  pTKratCalibReader->SetVerboseLevel(10);
  //pTKratCalibReader->SetKratContInName ("KRATTAEVENT"); 
  pTKratCalibReader->SetKratContOutName ("KRATTAEVENT_FF"); 
  pTKratCalibReader->SaveOutputIntoTree();
    
  // -------------------------------------------------------------------
  fRun->AddTask( pTKratCalibReader );
  //fRun->AddTask( pTKratCalibAna );
  
  // -------------------------------------------------------------------
  fRun->Init();
  
  Long64_t events = pTKratCalibReader->GetNumberOfEventsInInputFiles();
  cout << "All events in all files = " << events << " able to be analysed" << endl;
  cout << "Press enter to start... "  << flush;
  getchar();
  
  // -------------------------------------------------------------------
  
  fRun->Run(0, events ); /// second argument can't be larger than number of events, otherwise it crashes.
  
  cout << "TEST_PASSED" << endl; /// for testing only  
  cout << "Done" << endl;
  
}