CondTools/JetMET/plugins/JetCorrectorDBWriter.cc - master - githubmirror/cmssw

CondTools/JetMET/plugins/JetCorrectorDBWriter.cc
93 строки3 KB
Lauren Hay
Move DB related codes to CondTools
02 дек 2021, 01:59
02 дек 2021, 01:592fee424
// Author: Benedikt Hegner
// Email:  benedikt.hegner@cern.ch

#include <memory>
#include <string>
#include <fstream>
#include <iostream>
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/EDMException.h"
#include "CondCore/DBOutputService/interface/PoolDBOutputService.h"
#include "CondFormats/JetMETObjects/interface/JetCorrectorParameters.h"

class JetCorrectorDBWriter : public edm::one::EDAnalyzer<> {
public:
  JetCorrectorDBWriter(const edm::ParameterSet&);
  void beginJob() override;
  void analyze(const edm::Event&, const edm::EventSetup&) override {}
  void endJob() override {}
  ~JetCorrectorDBWriter() override {}

private:
  std::string era;
  std::string algo;
  std::string path;
  std::string inputTxtFile;
  std::string payloadTag;
};

// Constructor
JetCorrectorDBWriter::JetCorrectorDBWriter(const edm::ParameterSet& pSet) {
  era = pSet.getUntrackedParameter<std::string>("era");
  algo = pSet.getUntrackedParameter<std::string>("algo");
  path = pSet.getUntrackedParameter<std::string>("path");
  //payloadTag = "JetCorrectorParametersCollection_"+era+"_"+algo;
  payloadTag = algo;
}

// Begin Job
void JetCorrectorDBWriter::beginJob() {
  //  std::string path("CondFormats/JetMETObjects/data/");

  JetCorrectorParametersCollection payload;
  std::cout << "Starting to import payload " << payloadTag << " from text files." << std::endl;
  for (int i = 0; i < JetCorrectorParametersCollection::N_LEVELS; ++i) {
    std::string append("_");
    std::string ilev =
        JetCorrectorParametersCollection::findLabel(static_cast<JetCorrectorParametersCollection::Level_t>(i));
    append += ilev;
    append += "_";
    append += algo;
    append += ".txt";
    inputTxtFile = path + era + append;
    try {
      edm::FileInPath fip(inputTxtFile);
      std::cout << "Opened file " << inputTxtFile << std::endl;
      // create the parameter object from file
      std::vector<std::string> sections;
      JetCorrectorParametersCollection::getSections(fip.fullPath(), sections);
      if (sections.empty()) {
        payload.push_back(i, JetCorrectorParameters(fip.fullPath(), ""));
      } else {
        for (std::vector<std::string>::const_iterator isectbegin = sections.begin(),
                                                      isectend = sections.end(),
                                                      isect = isectbegin;
             isect != isectend;
             ++isect) {
          payload.push_back(i, JetCorrectorParameters(fip.fullPath(), *isect), ilev + "_" + *isect);
          std::cout << "Added " << ilev + "_" + *isect << " to record " << i << std::endl;
        }
      }
      std::cout << "Added record " << i << std::endl;
    } catch (edm::Exception const&) {
      std::cout << "Did not find JEC file " << inputTxtFile << std::endl;
    }
  }

  std::cout << "Opening PoolDBOutputService" << std::endl;

  // now write it into the DB
  edm::Service<cond::service::PoolDBOutputService> s;
  if (s.isAvailable()) {
    std::cout << "Setting up payload tag " << payloadTag << std::endl;
    s->writeOneIOV(payload, s->beginOfTime(), payloadTag);
  }
  std::cout << "Wrote in CondDB payload label: " << payloadTag << std::endl;
}

DEFINE_FWK_MODULE(JetCorrectorDBWriter);