JS-Jet/JetUESize-v2Psi2Centrality/Fun4All_UE_v2_Psi2.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 <g4jets/FastJetAlgo.h>
0016 #include <g4jets/JetReco.h>
0017 #include <g4jets/TowerJetInput.h>
0018 #include <g4jets/TruthJetInput.h>
0019
0020 #include <g4centrality/PHG4CentralityReco.h>
0021
0022 #include <jetbackground/CopyAndSubtractJets.h>
0023 #include <jetbackground/DetermineTowerBackground.h>
0024 #include <jetbackground/FastJetAlgoSub.h>
0025 #include <jetbackground/RetowerCEMC.h>
0026 #include <jetbackground/SubtractTowers.h>
0027 #include <jetbackground/SubtractTowersCS.h>
0028
0029 #include <HIJetReco.C>
0030
0031 #include <uevsetacentrality/UEvsEtaCentrality.h>
0032
0033
0034 R__LOAD_LIBRARY(libfun4all.so)
0035 R__LOAD_LIBRARY(libg4jets.so)
0036 R__LOAD_LIBRARY(libjetbackground.so)
0037 R__LOAD_LIBRARY(libUEvsEtaCentrality.so)
0038 R__LOAD_LIBRARY(libg4centrality.so)
0039 R__LOAD_LIBRARY(libg4dst.so)
0040
0041
0042 #endif
0043
0044
0045 void  Fun4All_UE_v2_Psi2(const char *filelistbbc = "dst_bbc_g4hit.list",
0046                      const char *filelistcalo = "dst_calo_cluster.list",
0047                      const char *outname = "outputest.root")
0048 {
0049
0050
0051   Fun4AllServer *se = Fun4AllServer::instance();
0052   int verbosity = 0;
0053
0054   se->Verbosity(verbosity);
0055   recoConsts *rc = recoConsts::instance();
0056
0057   PHG4CentralityReco *cent = new PHG4CentralityReco();
0058   cent->Verbosity(0);
0059   cent->GetCalibrationParameters().ReadFromFile("centrality", "xml", 0, 0, string(getenv("CALIBRATIONROOT")) + string("/Centrality/"));
0060   se->registerSubsystem( cent );
0061
0062
0063   HIJetReco();
0064
0065   UEvsEtaCentrality *myJetVal = new UEvsEtaCentrality();
0066   se->registerSubsystem(myJetVal);
0067
0068   Fun4AllInputManager *intrue = new Fun4AllDstInputManager("DSTbbc");
0069   intrue->AddListFile(filelistbbc,1);
0070   se->registerInputManager(intrue);
0071
0072   Fun4AllInputManager *in2 = new Fun4AllDstInputManager("DSTcalo");
0073   in2->AddListFile(filelistcalo,1);
0074   se->registerInputManager(in2);
0075
0076
0077   se->run(-1);
0078   se->End();
0079   se->PrintTimer();
0080   gSystem->Exit(0);
0081   return 0;
0082
0083 }