sEPDChannelAnalysis/macro/Fun4All_Calo.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 #include <phool/recoConsts.h>
0008 #include <fun4all/Fun4AllNoSyncDstInputManager.h>
0009 #include <fun4all/Fun4AllDstInputManager.h>
0010
0011 #include <fun4all/Fun4AllDstOutputManager.h>
0012 #include <fun4all/Fun4AllOutputManager.h>
0013
0014 // #include <calotowerbuilder/CaloTowerBuilder.h>
0015 #include <caloreco/CaloTowerBuilder.h>
0016 #include <caloreco/CaloWaveformProcessing.h>
0017
0018 #include <ffamodules/FlagHandler.h>
0019 #include <ffamodules/HeadReco.h>
0020 #include <ffamodules/SyncReco.h>
0021 #include <ffamodules/CDBInterface.h>
0022
0023 // #include <runtowerinfo/RunTowerInfo.h>
0024 #include <caloana/CaloAna.h>
0025 #include <fun4all/Fun4AllDstOutputManager.h>
0026 R__LOAD_LIBRARY(libfun4all.so)
0027 R__LOAD_LIBRARY(libfun4allraw.so)
0028 R__LOAD_LIBRARY(libcalo_reco.so)
0029 R__LOAD_LIBRARY(libcaloana.so)
0030 R__LOAD_LIBRARY(libffamodules.so)
0031 #endif
0032
0033 void Fun4All_Calo(const char *filelistcalo = "./test_dst.root", const char *outfile = "/sphenix/user/ecroft/private_prod_52917_test.root")
0034 {
0035
0036
0037
0038
0039   Fun4AllServer *se = Fun4AllServer::instance();
0040   int verbosity = 0;
0041
0042   se->Verbosity(verbosity);
0043   recoConsts *rc = recoConsts::instance();
0044
0045  //===============
0046   // conditions DB flags
0047   //===============
0048
0049   // global tag
0050   rc->set_StringFlag("CDB_GLOBALTAG","MDC2");
0051   // // 64 bit timestamp
0052   rc->set_uint64Flag("TIMESTAMP",6);
0053
0054   Fun4AllInputManager *in2 = new Fun4AllDstInputManager("DSTcalo");
0055   in2->AddListFile(filelistcalo,1);
0056   se->registerInputManager(in2);
0057
0058
0059 CaloAna *ca = new CaloAna("sEPD",outfile,20);
0060 // ca->SetOutputFileName("tim_calo_tree.root");
0061   // // choose CEMC, HCALIN or HCALOUT or whatever you named your
0062   // // calorimeter
0063   se->registerSubsystem(ca);
0064
0065
0066
0067
0068
0069   se->run(-1);
0070   se->End();
0071
0072 }