emcalHotTowerFinder/macro/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 #include <cdbobjects/CDBTTree.h>
0015
0016
0017 #include <towerid/towerid.h>
0018
0019 R__LOAD_LIBRARY(libcdbobjects.so)
0020 R__LOAD_LIBRARY(libfun4all.so)
0021 R__LOAD_LIBRARY(libtowerid.so)
0022 //R__LOAD_LIBRARY(libmakeMBDTrees.so)
0023 #endif
0024 void toweridmacro(const int nEvents = 5000, const char *listFile = "/sphenix/lustre01/sphnxpro/commissioning/DST/run_20488/dsts_fast/DST-00020488-0042.root", const char *inName = "commissioning3.root", const std::string cdbname = "test.root", float adccut_sg = 150, float adccut_k = 150,float sigmas_lo = -2, float sigmas_hi = 4, float SG_f = 0.5, float Kur_f = 0.5, float region_f = 1.0)
0025 {
0026   Fun4AllServer *se = Fun4AllServer::instance();
0027   recoConsts *rc = recoConsts::instance();
0028
0029   towerid *calo = new towerid(inName,cdbname,adccut_sg,adccut_k,sigmas_lo,sigmas_hi,nEvents,SG_f,Kur_f);
0030   se->registerSubsystem(calo);
0031
0032   Fun4AllInputManager *in = new Fun4AllDstInputManager("DSTcalo");
0033   in->AddFile(listFile);
0034   se->registerInputManager(in);
0035
0036   se->run(nEvents);
0037
0038   se->End();
0039   se->PrintTimer();
0040
0041   std::cout << "All done!" << std::endl;
0042
0043   gSystem->Exit(0);
0044 }
0045