JetNconstituents/macro/Fun4All_JetNconst.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 <g4centrality/PHG4CentralityReco.h>
0016
0017
0018 #include <HIJetReco.C>
0019
0020
0021 #include <jetnconstituents/JetNconstituents.h>
0022
0023 // #include <jetvalidation/JetValidation.h>
0024
0025 R__LOAD_LIBRARY(libfun4all.so)
0026 R__LOAD_LIBRARY(libg4jets.so)
0027 R__LOAD_LIBRARY(libjetbackground.so)
0028 R__LOAD_LIBRARY(libJetNconstituents.so)
0029 R__LOAD_LIBRARY(libg4centrality.so)
0030 R__LOAD_LIBRARY(libg4dst.so)
0031
0032
0033 #endif
0034
0035 void Fun4All_JetNconst(const char *filelistcalo = "/sphenix/u/tmengel/dst_calo_cluster.list",
0036                      const char *outname = "outputest.root")
0037 {
0038
0039
0040
0041   Fun4AllServer *se = Fun4AllServer::instance();
0042   int verbosity = 0;
0043
0044   se->Verbosity(verbosity);
0045   recoConsts *rc = recoConsts::instance();
0046
0047   PHG4CentralityReco *cent = new PHG4CentralityReco();
0048   cent->Verbosity(0);
0049   cent->GetCalibrationParameters().ReadFromFile("centrality", "xml", 0, 0, string(getenv("CALIBRATIONROOT")) + string("/Centrality/"));
0050   se->registerSubsystem( cent );
0051
0052   HIJetReco();
0053
0054   JetNconstituents *myJetNconst = new JetNconstituents("AntiKt_Tower_r04_Sub1", outname);
0055
0056   myJetNconst->setPtRange(5, 100);
0057   myJetNconst->setEtaRange(-1.1, 1.1);
0058   se->registerSubsystem(myJetNconst);
0059
0060   Fun4AllInputManager *in2 = new Fun4AllDstInputManager("DSTcalo");
0061   in2->AddListFile(filelistcalo,1);
0062   se->registerInputManager(in2);
0063
0064
0065   se->run(-1);
0066   se->End();
0067
0068   gSystem->Exit(0);
0069   return 0;
0070
0071 }