dst_eventeval.C

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file dst_eventeval.C

#include <fun4all/Fun4AllInputManager.h>
#include <fun4all/Fun4AllDstInputManager.h>
#include <fun4all/Fun4AllServer.h>
#include <eiceval/EventEvaluatorEIC.h>

R__LOAD_LIBRARY(libfun4all.so)
R__LOAD_LIBRARY(libg4eval.so)
R__LOAD_LIBRARY(libeiceval.so)

namespace EVENT_EVALUATOR
{
  int Verbosity = 0;
  float EnergyThreshold = 0.05;
}  // namespace EVENT_EVALUATOR

void dst_eventeval(const std::string fileList="pb-list.txt", int nEvents=1e5, std::string outputFile="event_eval.root")
{
  Fun4AllServer *se = Fun4AllServer::instance();
  se->Verbosity(0);

  Fun4AllInputManager *hitsin = new Fun4AllDstInputManager("DSTin");
  hitsin->AddListFile(fileList);
  se->registerInputManager(hitsin);

  EventEvaluatorEIC *eval = new EventEvaluatorEIC("EVENTEVALUATOR", outputFile);
  eval->set_reco_tracing_energy_thresholdMC(EVENT_EVALUATOR::EnergyThreshold);
  eval->Verbosity(EVENT_EVALUATOR::Verbosity);

  eval->set_do_TRACKS(true);
  eval->set_do_PROJECTIONS(true);
  eval->set_do_VERTEX(true);

  eval->set_do_HITS(true); // new functionality

  eval->set_do_CEMC(true);
  eval->set_do_EEMC(true);
  eval->set_do_FEMC(true);
  eval->set_do_HCALIN(true);
  eval->set_do_HCALOUT(true);
  eval->set_do_FHCAL(true);
  eval->set_do_CLUSTERS(true);

  eval->set_do_MCPARTICLES(true);
  eval->set_do_HEPMC(true);
  eval->set_do_store_event_level_info(true);

  se->registerSubsystem(eval);

  se->run(nEvents);
  se->End();
  std::cout << "All done" << std::endl;
  delete se;
  gSystem->Exit(0);
}