JetUnderlyingEvent/macro/Fun4All_JetVal.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 #include <jetunderlyingevent/JetUnderlyingEvent.h>
0021
0022
0023 R__LOAD_LIBRARY(libfun4all.so)
0024 R__LOAD_LIBRARY(libg4jets.so)
0025 R__LOAD_LIBRARY(libjetbase.so)
0026 R__LOAD_LIBRARY(libjetbackground.so)
0027 R__LOAD_LIBRARY(JetUnderlyingEvent.so)
0028 R__LOAD_LIBRARY(libg4centrality.so)
0029 R__LOAD_LIBRARY(libg4dst.so)
0030
0031
0032 #endif
0033
0034 void Fun4All_JetVal(const char *filelisttruth = "dst_03_27/dst_truth_jet.list",
0035             const char *filelistcalo = "dst_03_27/dst_calo_cluster.list",
0036             const char *filelistglobal = "dst_03_27/dst_global.list",
0037             const char *filelistg4hit = "dst_03_27/dst_bbc_g4hit.list",
0038             const char *outname = "outputest.root")
0039 {
0040
0041
0042   Fun4AllServer *se = Fun4AllServer::instance();
0043   int verbosity = 0;
0044
0045   se->Verbosity(verbosity);
0046   recoConsts *rc = recoConsts::instance();
0047   rc->set_StringFlag("CDB_GLOBALTAG", "ProdA_2023");
0048
0049   PHG4CentralityReco *cent = new PHG4CentralityReco();
0050   cent->Verbosity(verbosity);
0051   cent->GetCalibrationParameters().ReadFromFile("centrality", "xml", 0, 0, string(getenv("CALIBRATIONROOT")) + string("/Centrality/"));
0052   se->registerSubsystem( cent );
0053 <<<<<<< HEAD:JS-Jet/JetValidation/macro/Fun4All_JetVal.C
0054
0055   Enable::VERBOSITY = verbosity;
0056   HIJetReco();
0057
0058   JetValidation *myJetVal = new JetValidation("AntiKt_Tower_r04_Sub1", "AntiKt_Truth_r04", outname);
0059 =======
0060 >>>>>>> 07b79f5351fbab0f10c9a719c1cb49f741419ec1:JS-Jet/JetUnderlyingEvent/macro/Fun4All_JetVal.C
0061
0062   HIJetReco();
0063
0064
0065   JetUnderlyingEvent *myJetVal = new JetUnderlyingEvent("name",outname);
0066   se->registerSubsystem(myJetVal);
0067
0068   Fun4AllInputManager *intrue = new Fun4AllDstInputManager("DSTtruth");
0069   intrue->AddListFile(filelisttruth,1);
0070   se->registerInputManager(intrue);
0071
0072   Fun4AllInputManager *in2 = new Fun4AllDstInputManager("DSTcalo");
0073   in2->AddListFile(filelistcalo,1);
0074   se->registerInputManager(in2);
0075
0076   Fun4AllInputManager *in3 = new Fun4AllDstInputManager("DSTglobal");
0077   in3->AddListFile(filelistglobal,1);
0078   se->registerInputManager(in3);
0079
0080   se->run(-1);
0081
0082   se->End();
0083
0084   gSystem->Exit(0);
0085   return 0;
0086
0087 }