HCalHotTowerFinder/macro/hcal_toweridmacro.C

0001 #pragma once
0002 #if ROOT_VERSION_CODE >= ROOT_VERSION(6,00,0)
0003 #include <fun4all/SubsysReco.h>
0004 #include <fun4all/Fun4AllServer.h>
0005 #include <fun4all/Fun4AllInputManager.h>
0006 #include <fun4all/Fun4AllDstInputManager.h>
0007
0008 #include <fun4all/Fun4AllDstOutputManager.h>
0009 #include <fun4all/Fun4AllOutputManager.h>
0010 #include <fun4all/Fun4AllServer.h>
0011
0012 #include <phool/PHRandomSeed.h>
0013 #include <phool/recoConsts.h>
0014
0015 #include <hcal_towerid/hcal_towerid.h>
0016
0017 R__LOAD_LIBRARY(libfun4all.so)
0018 R__LOAD_LIBRARY(libtowerid.so)
0019 //R__LOAD_LIBRARY(libmakeMBDTrees.so)
0020 #endif
0021 void hcal_toweridmacro(const int nEvents = 5000, const char *listFile = "/sphenix/lustre01/sphnxpro/commissioning/DST/run_22950/dsts_fast/DST-00022950-0018.root", const char *inName = "commissioning.root",const std::string cdbname_i = "cdbname_i.root",const std::string cdbname_o = "cdbname_o.root",  float adccut_i = 150, float adccut_o = 150, float sigmas_lo = 2, float sigmas_hi = 4, float inner_f = 0.05, float outer_f = 0.05)
0022 {
0023   Fun4AllServer *se = Fun4AllServer::instance();
0024   recoConsts *rc = recoConsts::instance();
0025
0026   hcal_towerid *calo = new hcal_towerid(inName,cdbname_i,cdbname_o,adccut_i,adccut_o,sigmas_lo,sigmas_hi,inner_f,outer_f);
0027   se->registerSubsystem(calo);
0028
0029   Fun4AllInputManager *in = new Fun4AllDstInputManager("DSTcalo");
0030   in->AddFile(listFile);
0031   se->registerInputManager(in);
0032
0033   se->run(nEvents);
0034
0035   se->End();
0036   se->PrintTimer();
0037
0038   std::cout << "All done!" << std::endl;
0039
0040   gSystem->Exit(0);
0041 }
0042