TowerRhoAna/macro/Fun4All_jetQA.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 <jetbackground/DetermineTowerRho.h>
0009
0010 #include <jetbase/TowerJetInput.h>
0011
0012 #include <towerrhoana/TowerRhoAna.h>
0013
0014
0015 #include <HIJetReco.C>
0016
0017 R__LOAD_LIBRARY(libfun4all.so)
0018 R__LOAD_LIBRARY(libg4jets.so)
0019 R__LOAD_LIBRARY(libjetbackground.so)
0020 R__LOAD_LIBRARY(libg4dst.so)
0021
0022 R__LOAD_LIBRARY(libjetbase.so)
0023 R__LOAD_LIBRARY(libTowerRhoAna.so)
0024
0025 #endif
0026
0027 void Fun4All_jetQA(const int nEvents = -1,
0028            const char *filelistcalo = "dst_calo_cluster.list",
0029            const char *filelistglobal = "dst_global.list",
0030            const char *outname = "test.root")
0031 {
0032
0033   //make the fun4all server
0034   Fun4AllServer *se = Fun4AllServer::instance();
0035
0036   //set the verbosity (increase to print more output for debugging)
0037   int verbosity = 0;
0038   se->Verbosity(verbosity);
0039   Enable::VERBOSITY = verbosity;
0040
0041   //run the jet reconstruction
0042   HIJetReco();
0043
0044   std::vector<TowerJetInput*> background_inputs =
0045   {
0046     new TowerJetInput(Jet::CEMC_TOWERINFO_RETOWER),
0047     new TowerJetInput(Jet::HCALIN_TOWERINFO),
0048     new TowerJetInput(Jet::HCALOUT_TOWERINFO)
0049   };
0050
0051
0052   DetermineTowerRho *dtr = new DetermineTowerRho();
0053   dtr->add_method(TowerRho::Method::AREA); // do area method
0054   dtr->add_method(TowerRho::Method::MULT); // do multiplicy method
0055   dtr->set_algo(Jet::ALGO::KT); // default = KT
0056   dtr->set_par(0.4); // default = 0.4
0057   dtr->set_tower_abs_eta(1.1); // default=1.1
0058   dtr->set_jet_abs_eta(0.7); //default is tower_abs_eta - R
0059   dtr->set_omit_nhardest(2); // default=2
0060   dtr->set_ghost_area(0.01); // default=0.01
0061   dtr->set_jet_min_pT(0.0); // default=0.0
0062   for(auto input : background_inputs) dtr->add_tower_input(input); // add inputs to estimate rho
0063   dtr->Verbosity(verbosity);
0064   se->registerSubsystem(dtr);
0065
0066
0067   // TowerRhoAna
0068   TowerRhoAna *ana = new TowerRhoAna(outname);
0069   ana->Verbosity(verbosity);
0070   ana->add_area_rho_node("TowerRho_AREA"); // add rho node to output tree
0071   ana->add_mult_rho_node("TowerRho_MULT"); // add rho node to output tree
0072   se->registerSubsystem(ana);
0073
0074   //add inputs
0075   //add any additional needed inputs here
0076   Fun4AllInputManager *incalo = new Fun4AllDstInputManager("DSTcalo");
0077   incalo->AddListFile(filelistcalo,1);
0078   se->registerInputManager(incalo);
0079
0080   Fun4AllInputManager *inglobal = new Fun4AllDstInputManager("DSTglobal");
0081   inglobal->AddListFile(filelistglobal,1);
0082   se->registerInputManager(inglobal);
0083
0084   //run
0085   se->run(nEvents);
0086   se->End();
0087
0088   gSystem->Exit(0);
0089   return 0;
0090 }