myjetanalysis/macro/Fun4All_JetAna_MDC.C

0001 #ifndef MACRO_FUN4ALLJETANA_C
0002 #define MACRO_FUN4ALLJETANA_C
0003
0004 #include <G4_Jets.C>
0005
0006 #include <fun4all/Fun4AllDstInputManager.h>
0007 #include <fun4all/Fun4AllInputManager.h>
0008 #include <fun4all/Fun4AllServer.h>
0009 #include <fun4all/SubsysReco.h>
0010
0011
0012 // here you need your package name (set in configure.ac)
0013 // this suppresses a cling warning when reading an eigen include
0014 #pragma GCC diagnostic push
0015 #pragma GCC diagnostic ignored "-Wundefined-internal"
0016 #include <myjetanalysis/MyJetAnalysis.h>
0017 #pragma GCC diagnostic pop
0018
0019 R__LOAD_LIBRARY(libfun4all.so)
0020 R__LOAD_LIBRARY(libmyjetanalysis.so)
0021
0022 void Fun4All_JetAna_MDC(const int nevnt = 0, const std::string &inlist1 = "dst_tracks.list",  const std::string &inlist2 = "dst_calo_cluster.list", const std::string &inlist3 = "dst_truth.list")
0023 {
0024   gSystem->Load("libmyjetanalysis");
0025   gSystem->Load("libg4dst");
0026
0027   Jet_Reco();
0028
0029   Fun4AllServer *se = Fun4AllServer::instance();
0030   // uncomment if you want the event counter printed
0031   //se->Verbosity(1);
0032   MyJetAnalysis *myJetAnalysis = new MyJetAnalysis("AntiKt_Tower_r04", "AntiKt_Truth_r04", "myjetanalysis.root");
0033   //myJetAnalysis->Verbosity(10);
0034   // change lower pt and eta cut to make them visible
0035   myJetAnalysis->setPtRange(1, 100);
0036   myJetAnalysis->setEtaRange(-1.1, 1.1);
0037   myJetAnalysis->use_initial_vertex(true);
0038   se->registerSubsystem(myJetAnalysis);
0039
0040   Fun4AllInputManager *in = new Fun4AllDstInputManager("DSTin1");
0041   in->AddListFile(inlist1);
0042   se->registerInputManager(in);
0043   in = new Fun4AllDstInputManager("DSTin2");
0044   in->AddListFile(inlist2);
0045   se->registerInputManager(in);
0046   in = new Fun4AllDstInputManager("DSTin3");
0047   in->AddListFile(inlist3);
0048   se->registerInputManager(in);
0049
0050   se->run(nevnt);
0051   se->End();
0052   delete se;
0053   gSystem->Exit(0);
0054 }
0055
0056 #endif