// Macro to reconstruct the detector output in eicroot. 

eic_reco(Int_t nEvents = 10){
    
    // ========================================================================
    // Verbosity level (0=quiet, 1=event level, 2=track level, 3=debug)
    Int_t iVerbose = 0;
    TString MCFile = "data/eic_tau_p_without_selection.mc.root";
    TString parFile = "data/params_tau_p_without_selection.mc.root";
    // ----  Load libraries   -------------------------------------------------
    gROOT->LoadMacro("$VMCWORKDIR/example/gconfig/basiclibs.C");
    basiclibs();
    gSystem->Load("libFairTools");
    gSystem->Load("libFairDB");
    gSystem->Load("libGeoBase");
    gSystem->Load("libParBase");
    gSystem->Load("libTrkBase");
    gSystem->Load("libProof");
    gSystem->Load("libProofPlayer");
    gSystem->Load("libBase");
    gSystem->Load("libMCStack");
    gSystem->Load("libGen");
    gSystem->Load("libfield");
    gSystem->Load("libRho");
    gSystem->Load("libEicTools");
    gSystem->Load("libEICPASSIVE");
    gSystem->Load("libFairBatr");
    gSystem->Load("libFairEMCA");
    gSystem->Load("libFairDirc");
    gSystem->Load("libFairEmlf");
    gSystem->Load("libFairHACA");
    gSystem->Load("libFairHtck");
    gSystem->Load("libFairRich");
    gSystem->Load("libFairSilm");
    gSystem->Load("libFairTrac");

    // ------------------------------------------------------------------------
    // Output file
    FairFileNameCreator creator(MCFile.Data());
    TString DigiFile = creator.GetDigiFileName(false).c_str();
    TString outFile = creator.GetRecoFileName(false).c_str();
    
    std::cout << "MCFile  : " << MCFile.Data()<< std::endl;
    std::cout << "DigiFile: " << DigiFile.Data()<< std::endl;
    std::cout << "RecoFile: " << outFile.Data()<< std::endl;
    
    // -----   Timer   --------------------------------------------------------
    TStopwatch timer;
    timer.Start();
    // ------------------------------------------------------------------------
    
    // -----   Reconstruction run   -------------------------------------------
    FairRunAna *fRun= new FairRunAna();
    fRun->SetInputFile(MCFile);
    fRun->AddFriend(DigiFile);
    fRun->SetOutputFile(outFile);
    // ------------------------------------------------------------------------
    
    
    
    // -----  Parameter database   --------------------------------------------
    FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
    FairParRootFileIo* parInput1 = new FairParRootFileIo(kTRUE);
    parInput1->open(parFile.Data(),"UPDATE");
    rtdb->setFirstInput(parInput1);
    
    // -----    Default MVD hit producer   --------------------------------------------
    
 //   PndMvdClusterTask* mvdmccls = new PndMvdClusterTask();
 //   mvdmccls->SetVerbose(iVerbose);
 //   fRun->AddTask(mvdmccls);
    
    // -----    OR separate tasks   ---------------------------------------------------
    //Double_t chargecut = 1.e5;
    //PndMvdStripClusterTask* mvdmccls = new PndMvdStripClusterTask();
    //mvdmccls->SetVerbose(iVerbose);
    //fRun->AddTask(mvdmccls);
    //PndMvdPixelClusterTask* mvdClusterizer = new PndMvdPixelClusterTask();
    //mvdClusterizer->SetVerbose(iVerbose);
    //fRun->AddTask(mvdClusterizer);
    
    rtdb->setOutput(parInput1);
    rtdb->print();
    // =====                 End of HitProducers                           =====
    // =========================================================================
    //   PndMvdGeoPar* geoPar  = (PndMvdGeoPar*)(rtdb->getContainer("PndMvdGeoPar")); 
    
    // -----   Intialise and run   --------------------------------------------
    fRun->Init();
    
    fRun->Run(0,nEvents);
    // ------------------------------------------------------------------------
    
    rtdb->saveOutput();
    rtdb->print();
    // -----   Finish   -------------------------------------------------------
    timer.Stop();
    Double_t rtime = timer.RealTime();
    Double_t ctime = timer.CpuTime();
    cout << endl << endl;
    cout << "Macro finished succesfully." << endl;
    cout << "Output file is "    << outFile << endl;
    cout << "Parameter file is " << parFile << endl;
    cout << "Real time " << rtime << " s, CPU time " << ctime << " s" << endl;
    cout << endl;
    // ------------------------------------------------------------------------

    
    
}