DijetThirdJetCorrection/macros/RunThirdJetSpectra.C

0001
0002 #pragma once
0003 #if ROOT_VERSION_CODE >= ROOT_VERSION(6,00,0)
0004 #include <string>
0005 #include <sstream>
0006 #include <vector>
0007 #include <utility>
0008 #include <map>
0009 #include <algorithm>
0010 #include <fun4all/Fun4AllBase.h>
0011 #include <fun4all/Fun4AllUtils.h>
0012 #include <fun4all/Fun4AllServer.h>
0013 #include <fun4all/Fun4AllInputManager.h>
0014 #include <phhepmc/Fun4AllHepMCInputManager.h>
0015 #include <fun4all/Fun4AllDstInputManager.h>
0016 #include <fun4all/SubsysReco.h>
0017 #include <ffamodules/CDBInterface.h>
0018 #include <hepmcjettrigger/HepMCJetTrigger.h>
0019 #include <thirdjetspectra/ThirdJetSpectra.h>
0020
0021
0022 R__LOAD_LIBRARY(libfun4all.so);
0023 R__LOAD_LIBRARY(libffamodules.so);
0024 R__LOAD_LIBRARY(libHepMCJetTrigger.so);
0025 R__LOAD_LIBRARY(libThirdJetSpectra.so);
0026 R__LOAD_LIBRARY(libphhepmc.so);
0027
0028 int RunThirdJetSpectra(std::string filename, std::string segment)
0029 {
0030     //run jet spectra using some parameters
0031     Fun4AllServer* se = Fun4AllServer::instance();
0032     Fun4AllHepMCInputManager* in = new Fun4AllHepMCInputManager("in");
0033     se->registerInputManager(in);
0034     se->fileopen(in->Name().c_str(), filename);
0035     ThirdJetSpectra* spec=new ThirdJetSpectra(segment);
0036     HepMCJetTrigger* trig=new HepMCJetTrigger(30.);
0037     se->registerSubsystem(trig);
0038     se->registerSubsystem(spec);
0039     se->run(); //for test purposes
0040     spec->Print();
0041     return 0;
0042 }
0043 #endif