ColdQCDENC/SEnergyCorrelator/RunENCCalculation.cxx

0001 /// ---------------------------------------------------------------------------
0002 /*! \file    RunENCCalculation.cxx
0003  *  \authors Derek Anderson
0004  *  \date    02.02.2023
0005  *
0006  *  A ROOT macro to run the SEnergyCorrelator module.
0007  */
0008 /// ---------------------------------------------------------------------------
0009
0010 #ifndef RUNENCCALCULATION_CXX
0011 #define RUNENCCALCULATION_CXX
0012
0013 // c++ utilities
0014 #include <string>
0015 #include <vector>
0016 #include <cstdlib>
0017 #include <utility>
0018 // module definition
0019 #include <senergycorrelator/SEnergyCorrelator.h>
0020 // macro options
0021 #include "EnergyCorrelatorOptions.h"
0022
0023 using namespace std;
0024 using namespace SColdQcdCorrelatorAnalysis;
0025
0026 // load libraries
0027 R__LOAD_LIBRARY(libsenergycorrelator.so)
0028 R__LOAD_LIBRARY(libscorrelatorutilities.so)
0029
0030
0031
0032 // macro body -----------------------------------------------------------------
0033
0034 void RunENCCalculation(const bool doBatch = false, const int verbosity = 0) {
0035
0036   // input/output files
0037   const vector<string> inFiles = {
0038     "../SCorrelatorJetTreeMaker/output/intermediate_merge/pp200py8jet10run15_towerSubAntiKtR04_2024sep28/correlatorJetTree.pp200py8jet10run15_towerSubAntiKtR04_0000to0099.d29m9y2024.root"
0039   };
0040   const pair<string, string> outFiles = {
0041     "twoPoint.pp200py8jet10run15_antiKtTowerR04_packageCalcOnly.refactor_addGlobalCalc_reco.d6m11y2024.root",
0042     "twoPoint.pp200py8jet10run15_antiKtTowerR04_packageCalcOnly.refactor_addGlobalCalc_true.d6m11y2024.root"
0043   };
0044
0045   // get module configurations
0046   SEnergyCorrelatorConfig cfg_reco = EnergyCorrelatorOptions::GetRecoConfig(inFiles, outFiles.first, doBatch, verbosity);
0047   SEnergyCorrelatorConfig cfg_true = EnergyCorrelatorOptions::GetTruthConfig(inFiles, outFiles.second, doBatch, verbosity);
0048
0049   // do correlator calculation on reco jets
0050   SEnergyCorrelator* recoCorrelator = new SEnergyCorrelator(cfg_reco);
0051   recoCorrelator -> Init();
0052   recoCorrelator -> Analyze();
0053   recoCorrelator -> End();
0054
0055   // do correlator calculation on truth jets
0056   SEnergyCorrelator* trueCorrelator = new SEnergyCorrelator(cfg_true);
0057   trueCorrelator -> Init();
0058   trueCorrelator -> Analyze();
0059   trueCorrelator -> End();
0060   return;
0061
0062 }
0063
0064 #endif
0065
0066 // end ------------------------------------------------------------------------