TopoclusterJets/macro/Fun4All_EMJetVal.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/EMJetVal.h"
0019 //#include </sphenix/user/jamesj3j3/analysis/JS-Jet/JetValidation/src/EMJetVal.h>
0020 //#include </sphenix/user/jamesj3j3/analysis/JS-Jet/JetValidation/src/JetValidation.h>
0021
0022 R__LOAD_LIBRARY(libfun4all.so)
0023 R__LOAD_LIBRARY(libg4jets.so)
0024 R__LOAD_LIBRARY(libjetbackground.so)
0025 R__LOAD_LIBRARY(libJetValidation.so)
0026 R__LOAD_LIBRARY(libg4centrality.so)
0027 R__LOAD_LIBRARY(libg4dst.so)
0028
0029
0030 #endif
0031
0032
0033 void Fun4All_EMJetVal(const char *filelisttruth = "dst_truth_jet.list",
0034                     const char *filelistcalo = "dst_calo_cluster.list",
0035             const char *filelistglobal = "dst_global.list",
0036             const char *outname = "output.root")
0037 {
0038
0039
0040   Fun4AllServer *se = Fun4AllServer::instance();
0041   int verbosity = 0;
0042
0043   se->Verbosity(verbosity);
0044   recoConsts *rc = recoConsts::instance();
0045
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
0055   // JetValidation *myJetVal = new JetValidation("AntiKt_Tower_r04_Sub1", "AntiKt_Truth_r04", outname);
0056
0057   EMJetVal *myEMJetVal = new EMJetVal("AntiKt_Tower_r04_Sub1","AntiKt_Truth_r04", outname);
0058
0059   myEMJetVal->setPtRange(5, 100);
0060   myEMJetVal->setEtaRange(-1.1, 1.1);
0061   myEMJetVal->doUnsub(1);
0062   myEMJetVal->doTruth(0);
0063   myEMJetVal->doSeeds(0);
0064   //  se->registerSubsystem(myJetVal);
0065   se->registerSubsystem(myEMJetVal);
0066
0067
0068   Fun4AllInputManager *intrue = new Fun4AllDstInputManager("DSTtruth");
0069   //intrue->AddListFile(filelisttruth,1);
0070   // turn off for running data
0071   // se->registerInputManager(intrue);
0072
0073   Fun4AllInputManager *in2 = new Fun4AllDstInputManager("DSTcalo");
0074   in2->AddListFile(filelistcalo,1);
0075   // se->registerInputManager(in2);
0076   se->registerInputManager(in2);
0077
0078   Fun4AllInputManager *in3 = new Fun4AllDstInputManager("DSTglobal");
0079   // in3->AddListFile(filelistglobal,1);
0080   // turn off for running data
0081   // se->registerInputManager(in3);
0082
0083   se->run(10000);
0084   se->End();
0085
0086   gSystem->Exit(0);
0087   return 0;
0088
0089 }