diff --git a/TIMBER/Framework/Corrections/JetCorrections/AK4corrections.cc b/TIMBER/Framework/Corrections/JetCorrections/AK4corrections.cc deleted file mode 100644 index e69de29..0000000 diff --git a/TIMBER/Framework/Corrections/JetCorrections/AK8corrections.cc b/TIMBER/Framework/Corrections/JetCorrections/AK8corrections.cc deleted file mode 100644 index e69de29..0000000 diff --git a/TIMBER/Framework/Corrections/JetCorrections/JetRecalibrator.cc b/TIMBER/Framework/Corrections/JetCorrections/JetRecalibrator.cc deleted file mode 100644 index 68bcb80..0000000 --- a/TIMBER/Framework/Corrections/JetCorrections/JetRecalibrator.cc +++ /dev/null @@ -1,28 +0,0 @@ -using namespace std; - -class JetRecalibrator -{ -private: - string globalTag; - string jetFlavor; - bool doResidualJECs; - string jecPath; - int upToLevel; - bool calculateSeparateCorrections; - bool calculateType1METCorrection; - float type1MET_jetPtThreshold; - float type1MET_skipEMfractionThreshold; - bool type1MET_skipMuons; - -public: - JetRecalibrator(): - upToLevel(3),calculateSeparateCorrections(false), - calculateType1METCorrection(false), - type1MET_jetPtThreshold(15.0),type1MET_skipEMfractionThreshold(0.9),type1MET_skipMuons(true) - { - - - } - ~JetRecalibrator(); - -}; \ No newline at end of file diff --git a/TIMBER/Framework/Corrections/JetCorrections/JetSmearer.cc b/TIMBER/Framework/Corrections/JetCorrections/JetSmearer.cc deleted file mode 100644 index e69de29..0000000 diff --git a/TIMBER/Framework/Corrections/in_progress/pileup.cc b/TIMBER/Framework/Corrections/in_progress/pileup.cc deleted file mode 100644 index 2fd2a0a..0000000 --- a/TIMBER/Framework/Corrections/in_progress/pileup.cc +++ /dev/null @@ -1,29 +0,0 @@ -// NOT USABLE CURRENTLY WITHOUT PROPPER FILES -using namespace ROOT::VecOps; - -namespace analyzer { - puWeightNom - TFile *puWeightNom = new TFile("JHUanalyzer/data/pileup/...","READ"); - TFile *puWeightUp = new TFile("JHUanalyzer/data/pileup/...up","READ"); - TFile *puWeightDown = new TFile("JHUanalyzer/data/pileup/...down","READ"); - - std::vector PUWeightLookup(int nvtx) { - std::vector weight; - double weightNom = 1; - double weightUp = 1; - double weightDown = 1; - - weightNom *= puWeightNom->GetBinContent(puWeightNom->FindBin(nvtx)); - weightUp *= puWeightUp->GetBinContent(puWeightUp->FindBin(nvtx)); - weightDown *= puWeightDown->GetBinContent(puWeightDown->FindBin(nvtx)); - - weight.push_back(weightNom); - weight.push_back(weightUp); - weight.push_back(weightDown); - return weight; - } - - std::vector PUWeightCalc(...) { - ... - } -} \ No newline at end of file diff --git a/TIMBER/Framework/Corrections/in_progress/tagging.cc b/TIMBER/Framework/Corrections/in_progress/tagging.cc deleted file mode 100644 index db5c95d..0000000 --- a/TIMBER/Framework/Corrections/in_progress/tagging.cc +++ /dev/null @@ -1,84 +0,0 @@ -#include -#include -#include "../PhysicsObjects/GenParticleTree.h" -#include "TLorentzVector.h" -#include "CondFormats/BTauObjects/interface/BTagCalibration.h" -#include "CondTools/BTau/interface/BTagCalibrationReader.h" - -namespace analyzer { - int TopMergeChecker() { - - } - bool WMatchChecker() { - - } - - float TopTaggingSF(char taggername[], int year, char mergecategory[], TLorentzVector jet) { - // If tau32 only (yes, == 0 is correct because of how strcmp returns) - if (strcmp(taggername,"tau32") == 0){ - if - } else if (strcmp(taggername,"tau32+sjbtag")){ - - } - } - float WTaggingSF(char taggername[], int year, char purity[], float pt, float eta) { - - } - float bTaggingSF(char taggername[], int year, char category[], float pt) { - if (strcmp(taggername,"DeepCSV") == 0){ - std::string csv_name = ; - } - - - BTagCalibration calib("csvv1", "CSVV1.csv"); - BTagCalibrationReader reader(BTagEntry::OP_LOOSE, // operating point - "central", // central sys type - {"up", "down"}); // other sys types - - reader.load(calib, // calibration instance - BTagEntry::FLAV_B, // btag flavour - "incl") // measurement type - - double jet_scalefactor = reader.eval_auto_bounds( - "central", - BTagEntry::FLAV_B, - b_jet.eta(), // absolute value of eta - b_jet.pt() - ); - double jet_scalefactor_up = reader.eval_auto_bounds( - "up", BTagEntry::FLAV_B, b_jet.eta(), b_jet.pt()); - double jet_scalefactor_do = reader.eval_auto_bounds( - "down", BTagEntry::FLAV_B, b_jet.eta(), b_jet.pt()); - } -} - -int main(char particle[], char tagger[], char wp[] int year, TLorentzVector jet){ - if (strcmp(particle,"top") == 0) { - if (strcmp(tagger,"tau32") == 0) { - - } else if (strcmp(tagger,"tau32+sjbtag") == 0) { - - } else { - printf ("Top tagger not supported. Quitting..."); - exit(EXIT_FAILURE); - } - } else if (strcmp(particle,"W") == 0) { - if strcmp(tagger,"tau21") { - - } else { - printf ("W tagger not supported. Quitting..."); - exit(EXIT_FAILURE); - } - - } else if (strcmp(particle,"bottom") == 0) { - if strcmp(tagger,"DeepCSV") { - - } else { - printf ("Bottom tagger not supported. Quitting..."); - exit(EXIT_FAILURE); - } - - } - - return 0; -} \ No newline at end of file diff --git a/TIMBER/Framework/Corrections/TopPt_weight.cc b/TIMBER/Framework/in_development/TopPt_weight.cc similarity index 100% rename from TIMBER/Framework/Corrections/TopPt_weight.cc rename to TIMBER/Framework/in_development/TopPt_weight.cc diff --git a/TIMBER/Framework/Corrections/Trigger_weight.cc b/TIMBER/Framework/in_development/Trigger_weight.cc similarity index 100% rename from TIMBER/Framework/Corrections/Trigger_weight.cc rename to TIMBER/Framework/in_development/Trigger_weight.cc