File indexing completed on 2025-08-05 08:13:16
0001 #include <HIJetReco.C>
0002
0003 #include <jetkinematiccheck/JetKinematicCheck.h>
0004
0005 #include <fun4all/Fun4AllDstInputManager.h>
0006 #include <fun4all/Fun4AllInputManager.h>
0007 #include <fun4all/Fun4AllServer.h>
0008 #include <fun4all/SubsysReco.h>
0009
0010 #include <fun4all/Fun4AllDstOutputManager.h>
0011 #include <fun4all/Fun4AllOutputManager.h>
0012 #include <fun4all/Fun4AllServer.h>
0013
0014 #include <phool/PHRandomSeed.h>
0015 #include <phool/recoConsts.h>
0016
0017 R__LOAD_LIBRARY(libfun4all.so)
0018 R__LOAD_LIBRARY(libjetbase.so)
0019 R__LOAD_LIBRARY(libjetbackground.so)
0020 R__LOAD_LIBRARY(libJetKinematicCheck.so)
0021 R__LOAD_LIBRARY(libg4centrality.so)
0022 R__LOAD_LIBRARY(libg4dst.so)
0023
0024 void Fun4All_JetAna(const int nevnt = 0,
0025 const std::string &filelistcalo = "dst_calo_cluster.list",
0026 const std::string &filelistglobal = "dst_global.list",
0027 const std::string &outname = "output.root")
0028 {
0029 Fun4AllServer *se = Fun4AllServer::instance();
0030 int verbosity = 0;
0031
0032 se->Verbosity(verbosity);
0033 recoConsts *rc = recoConsts::instance();
0034
0035 HIJetReco();
0036
0037 JetKinematicCheck *myJetKC = new JetKinematicCheck("AntiKt_Tower_r02_Sub1", "AntiKt_Tower_r03_Sub1", "AntiKt_Tower_r04_Sub1", outname);
0038
0039 myJetKC->setPtRange(5, 100);
0040 myJetKC->setEtaRange(-1.1, 1.1);
0041 se->registerSubsystem(myJetKC);
0042
0043 Fun4AllInputManager *in1 = new Fun4AllDstInputManager("DSTcalo");
0044 in1->AddListFile(filelistcalo, 1);
0045 se->registerInputManager(in1);
0046
0047 Fun4AllInputManager *in2 = new Fun4AllDstInputManager("DSTglobal");
0048 in2->AddListFile(filelistglobal, 1);
0049 se->registerInputManager(in2);
0050
0051 se->run(nevnt);
0052 se->End();
0053
0054 gSystem->Exit(0);
0055 return 0;
0056 }
