void hitcount(TString FileName = "test", Int_t nEvents = 0) {
    // ========================================================================

//    gROOT->LoadMacro("$VMCWORKDIR/gconfig/rootlogon.C");
//    rootlogon();

    Int_t iVerbose = 0;

    // PndFileNameCreator creator(FileName.Data());
    // std::cout <<    creator.GetParFileName()    << std::endl;
    // std::cout <<    creator.GetSimFileName()        << std::endl;
    // std::cout <<    creator.GetDigiFileName()       << std::endl;
    // std::cout <<    creator.GetRecoFileName()       << std::endl;
    // std::cout <<    creator.GetCustomFileName("blubb") << std::endl;


    // // TString simFile  = creator.GetSimFileName();
    // TString simFile = "sim_complete.root";
    // // TString parFile     = creator.GetParFileName() ;
    // TString parFile = "simparams.root";
    // // ----  Load libraries   -------------------------------------------------

    // // gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
    // // ------------------------------------------------------------------------
    // // Output file
    // //    PndFileNameCreator creator(simFile.Data());
    // // TString DigiFile = creator.GetDigiFileName();
    // TString digiFile = "digi_complete.root";
    // // TString RecoFile = creator.GetRecoFileName();
    // TString RecoFile = "reco_complete.root";
    // //    TString FitFile = creator.GetTrackFindingFileName(false).c_str();

    // // TString outFile = creator.GetCustomFileName("ana");
    // TString outFile = "mc_hitcount.root"

    // std::cout << "simFile  : " << simFile.Data() << std::endl;
    // std::cout << "DigiFile: " << DigiFile.Data() << std::endl;
    // std::cout << "RecoFile: " << RecoFile.Data() << std::endl;
    //    std::cout << "FitFile:  " << FitFile.Data() << std::endl;

    // gSystem->Load("libMCMatch");
    // gSystem->Load("libMCMatchExamples");

    TString parFile = "simparams.root";
    TString simFile = "sim_complete.root";
    TString digiFile = "digi_complete.root";
    TString recoFile = "reco_complete.root";
    // TString anaParFile = "";
    TString outputFile = "mc_hitcount.root";

    // -----   Timer   --------------------------------------------------------
    TStopwatch timer;
    timer.Start();
    // ------------------------------------------------------------------------



    // -----   Reconstruction run   -------------------------------------------
    FairRunAna *fRun = new FairRunAna();

    fRun->SetInputFile(simFile);
    fRun->AddFriend(digiFile);
    fRun->AddFriend(recoFile);
    //  fRun->AddFriend(FitFile);
    //  fRun->AddFriend(PidFile);

    fRun->SetOutputFile(outputFile);
    // ------------------------------------------------------------------------

    TString anaParFile = gSystem->Getenv("VMCWORKDIR");
    anaParFile +=  "/macro/params/all.par";
    FairParAsciiFileIo *parIo1 = new FairParAsciiFileIo();
    parIo1->open(anaParFile.Data(), "in");

    // -----  Parameter database   --------------------------------------------
    FairRuntimeDb *rtdb = fRun->GetRuntimeDb();

    FairParRootFileIo *parInput1 = new FairParRootFileIo(kTRUE);
    parInput1->open(parFile.Data(), "UPDATE");

    rtdb->setFirstInput(parInput1);
    rtdb->setSecondInput(parIo1);
    Bool_t kParameterMerged = kTRUE;
    rtdb->setOutput(parInput1);
    rtdb->print();


    // -----------------------------------------------------------------------



    // =========================================================================
    // ======                       Hit Producers                         ======
    // =========================================================================



    PndMCMatchCreatorTask *mcMatch = new PndMCMatchCreatorTask();
    mcMatch->SetVerbose(1);
    fRun->AddTask(mcMatch);

    PndMCIdealTrackFinder *idealtrackfinder = new PndMCIdealTrackFinder();
    fRun->AddTask(idealtrackfinder);

    PndMcHitCount * hitCountTask = new PndMcHitCount();
    fRun->AddTask(hitCountTask);

    fRun->Init();
    fRun->Run(0, nEvents);

    // ------------------------------------------------------------------------

    // rtdb->saveOutput();
    // -----   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 "    << outputFile << endl;
    cout << "Parameter file is " << parFile << endl;
    cout << "Real time " << rtime << " s, CPU time " << ctime << " s" << endl;
    cout << endl;
    // ------------------------------------------------------------------------


}