JetValidation/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 #include <HIJetReco.C>
0018 #include <JetValidation.h>
0019
0020 R__LOAD_LIBRARY(libfun4all.so)
0021 R__LOAD_LIBRARY(libg4jets.so)
0022 R__LOAD_LIBRARY(libjetbackground.so)
0023 R__LOAD_LIBRARY(libJetValidation.so)
0024 R__LOAD_LIBRARY(libg4centrality.so)
0025 R__LOAD_LIBRARY(libg4dst.so)
0026
0027
0028 #endif
0029 void Fun4All_JetVal(const char *filelisttruth = "dst_truth_jet.list",
0030                     const char *filelistcalo = "dst_calo_cluster.list",
0031             const char *filelistglobal = "dst_global.list",
0032             const char *outname = "outputest.root")
0033 {
0034
0035
0036   Fun4AllServer *se = Fun4AllServer::instance();
0037   int verbosity = 0;
0038
0039   se->Verbosity(verbosity);
0040   recoConsts *rc = recoConsts::instance();
0041
0042   PHG4CentralityReco *cent = new PHG4CentralityReco();
0043   cent->Verbosity(verbosity);
0044   cent->GetCalibrationParameters().ReadFromFile("centrality", "xml", 0, 0, string(getenv("CALIBRATIONROOT")) + string("/Centrality/"));
0045   se->registerSubsystem( cent );
0046
0047   Enable::VERBOSITY = verbosity;
0048   HIJetReco();
0049
0050   JetValidation *myJetVal = new JetValidation("AntiKt_Tower_r04_Sub1", "AntiKt_Truth_r04", outname);
0051
0052   myJetVal->setPtRange(5, 100);
0053   myJetVal->setEtaRange(-1.1, 1.1);
0054   myJetVal->doUnsub(1);
0055   myJetVal->doTruth(1);
0056   myJetVal->doSeeds(1);
0057   se->registerSubsystem(myJetVal);
0058
0059   Fun4AllInputManager *intrue = new Fun4AllDstInputManager("DSTtruth");
0060   intrue->AddListFile(filelisttruth,1);
0061   se->registerInputManager(intrue);
0062
0063   Fun4AllInputManager *in2 = new Fun4AllDstInputManager("DSTcalo");
0064   in2->AddListFile(filelistcalo,1);
0065   se->registerInputManager(in2);
0066
0067   Fun4AllInputManager *in3 = new Fun4AllDstInputManager("DSTglobal");
0068   in3->AddListFile(filelistglobal,1);
0069   se->registerInputManager(in3);
0070
0071   se->run(-1);
0072   se->End();
0073
0074   gSystem->Exit(0);
0075   return 0;
0076
0077 }