run.C

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

#if ROOT_VERSION_CODE >= ROOT_VERSION(6,00,0)

#include <GlobalVariables.C>

#include <fun4all/SubsysReco.h>
#include <fun4all/Fun4AllServer.h>
#include <fun4all/Fun4AllInputManager.h>
#include <fun4all/Fun4AllDstInputManager.h>
#include <fun4all/Fun4AllDstOutputManager.h>

#include </sphenix/u/weihuma/install/include/trackpidassoc/ElectronPid.h>
#include </sphenix/u/weihuma/install/include/trackpidassoc/TrackPidAssoc.h>


R__LOAD_LIBRARY(libfun4all.so)
R__LOAD_LIBRARY(libtrackpid.so)

#endif

void run(
  const char *inputFile = "/sphenix/sim/sim01/sphnxpro/MDC1/embed/embedDST_sHijing_0_20fm_50kHz_bkg_0_20fm-0000000001-04594.root",
  const string &outputroot = "/sphenix/u/weihuma/RunOutput/embedDST_sHijing_upsilon_0_20fm",
  bool output_ntuple = true
)
{
  gSystem->Load("libg4dst");
  gSystem->Load("libtrackpid");

  Fun4AllServer *se = Fun4AllServer::instance();
  se->Verbosity(1);

  ElectronPid *ePid = new ElectronPid("ElectronPid",outputroot+"_ElectronPid.root");
  ePid->set_output_ntuple(output_ntuple);
  ePid->Verbosity(0);
  ePid->setTrackcutlimits(1,1,20,10);//nmvtx, nintt, ntpc, quality
  ePid->setEMOPcutlimits(0.7,1.5);
  ePid->setHinOEMcutlimit(0.2);
  ePid->setPtcutlimit(2.0,30.0);
  ePid->setHOPcutlimit(0.3);
  se->registerSubsystem(ePid);

  Fun4AllInputManager *in = new Fun4AllDstInputManager("DST_TRACKS");
  in->Verbosity(1);
  //in->fileopen(inputFile);
 //in->AddListFile("filelist_0_4d88fm.txt"); //sHijing with Upsilon embeded;
 // in->AddListFile("filelist_0_12fm.txt"); //sHijing with Upsilon embeded;
  in->AddListFile("filelist_0_20fm.txt"); //sHijing with Upsilon embeded;
  
  se->registerInputManager(in);

  if(output_ntuple) {
        Fun4AllOutputManager *outePid = new Fun4AllDstOutputManager("outePid",outputroot+"_ElectronPid_DST.root");
        outePid->AddNode("TrackPidAssoc");
        se->registerOutputManager(outePid);
        outePid->Verbosity(1);
        outePid->Print();
  }

  se->run();
  se->End();

}