CaloTPGTranscoderULUTs.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    CaloTPGTranscoderULUTs
// Class:      CaloTPGTranscoderULUTs
//
//
// Original Author:  Jeremiah Mans
//         Created:  Fri Sep 15 11:49:44 CDT 2006
//
//

// system include files
#include <memory>

// user include files

#include "FWCore/Framework/interface/ModuleFactory.h"
#include "FWCore/Framework/interface/ESProducer.h"
#include "FWCore/Framework/interface/EventSetupRecordIntervalFinder.h"

#include "CalibCalorimetry/CaloTPG/interface/CaloTPGTranscoderULUT.h"
#include "CalibFormats/CaloTPG/interface/CaloTPGRecord.h"
#include "FWCore/Framework/interface/ValidityInterval.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "FWCore/Framework/interface/ESHandle.h"
#include "CondFormats/HcalObjects/interface/HcalLutMetadata.h"
#include "Geometry/HcalTowerAlgo/interface/HcalTrigTowerGeometry.h"
#include "Geometry/Records/interface/CaloGeometryRecord.h"
#include "Geometry/Records/interface/HcalRecNumberingRecord.h"

//
// class decleration
//

class CaloTPGTranscoderULUTs : public edm::ESProducer {
public:
  CaloTPGTranscoderULUTs(const edm::ParameterSet&);
  ~CaloTPGTranscoderULUTs() override;

  typedef std::unique_ptr<CaloTPGTranscoder> ReturnType;

  ReturnType produce(const CaloTPGRecord&);

private:
  // ----------member data ---------------------------
  const edm::FileInPath hfilename1_;
  const edm::FileInPath hfilename2_;
  const bool read_Ascii_Compression;
  const bool read_Ascii_RCT;
  const std::vector<int> ietal;
  const std::vector<int> ietah;
  const std::vector<int> ZS;
  const std::vector<int> LUTfactor;
  const bool linearLUTs_;
  const double nominal_gain;
  const double RCTLSB;
  const int NCTScaleShift;
  const int RCTScaleShift;
  const double lsbQIE8;
  const double lsbQIE11;
  edm::ESGetToken<HcalLutMetadata, HcalLutMetadataRcd> lutMetadataToken;
  edm::ESGetToken<HcalTrigTowerGeometry, CaloGeometryRecord> theTrigTowerGeometryToken;
  edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> topoToken;
};

//
// constants, enums and typedefs
//

//
// static data member definitions
//

//
// constructors and destructor
//
CaloTPGTranscoderULUTs::CaloTPGTranscoderULUTs(const edm::ParameterSet& iConfig)
    : hfilename1_(iConfig.getParameter<edm::FileInPath>("hcalLUT1")),
      hfilename2_(iConfig.getParameter<edm::FileInPath>("hcalLUT2")),
      read_Ascii_Compression(iConfig.getParameter<bool>("read_Ascii_Compression_LUTs")),
      read_Ascii_RCT(iConfig.getParameter<bool>("read_Ascii_RCT_LUTs")),
      ietal(iConfig.getParameter<std::vector<int>>("ietaLowerBound")),
      ietah(iConfig.getParameter<std::vector<int>>("ietaUpperBound")),
      ZS(iConfig.getParameter<std::vector<int>>("ZS")),
      LUTfactor(iConfig.getParameter<std::vector<int>>("LUTfactor")),
      linearLUTs_(iConfig.getParameter<bool>("linearLUTs")),
      nominal_gain(iConfig.getParameter<double>("nominal_gain")),
      RCTLSB(iConfig.getParameter<double>("RCTLSB")),
      NCTScaleShift(iConfig.getParameter<edm::ParameterSet>("tpScales")
                        .getParameter<edm::ParameterSet>("HF")
                        .getParameter<int>("NCTShift")),
      RCTScaleShift(iConfig.getParameter<edm::ParameterSet>("tpScales")
                        .getParameter<edm::ParameterSet>("HF")
                        .getParameter<int>("RCTShift")),
      lsbQIE8(iConfig.getParameter<edm::ParameterSet>("tpScales")
                  .getParameter<edm::ParameterSet>("HBHE")
                  .getParameter<double>("LSBQIE8")),
      lsbQIE11(iConfig.getParameter<edm::ParameterSet>("tpScales")
                   .getParameter<edm::ParameterSet>("HBHE")
                   .getParameter<double>("LSBQIE11")) {
  auto cc = setWhatProduced(this);
  lutMetadataToken = cc.consumes();
  theTrigTowerGeometryToken = cc.consumes();
  topoToken = cc.consumes();
}

CaloTPGTranscoderULUTs::~CaloTPGTranscoderULUTs() {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
}

//
// member functions
//

// ------------ method called to produce the data  ------------
CaloTPGTranscoderULUTs::ReturnType CaloTPGTranscoderULUTs::produce(const CaloTPGRecord& iRecord) {
  using namespace edm::es;
  std::string file1 = "";
  std::string file2 = "";
  if (read_Ascii_RCT && read_Ascii_Compression) {
    edm::LogInfo("Level1") << "Using " << hfilename1_.fullPath() << " & " << hfilename2_.fullPath()
                           << " for CaloTPGTranscoderULUTs HCAL initialization";
    //std::unique_ptr<CaloTPGTranscoder> pTCoder(new CaloTPGTranscoderULUT(hfilename1_.fullPath(), hfilename2_.fullPath()));
    //return pTCoder;
    file1 = hfilename1_.fullPath();
    file2 = hfilename2_.fullPath();
  } else if (read_Ascii_RCT && !read_Ascii_Compression) {
    edm::LogInfo("Level1") << "Using analytical compression and " << hfilename2_.fullPath()
                           << " RCT decompression for CaloTPGTranscoderULUTs HCAL initialization";
    //std::unique_ptr<CaloTPGTranscoder> pTCoder(new CaloTPGTranscoderULUT("", hfilename2_.fullPath()));
    //return pTCoder;
    file2 = hfilename2_.fullPath();
  } else if (read_Ascii_Compression && !read_Ascii_RCT) {
    edm::LogInfo("Level1") << "Using ASCII compression tables " << hfilename1_.fullPath()
                           << " and automatic RCT decompression for CaloTPGTranscoderULUTs HCAL initialization";
    //std::unique_ptr<CaloTPGTranscoder> pTCoder(new CaloTPGTranscoderULUT(hfilename1_.fullPath(),""));
    //return pTCoder;
    file1 = hfilename1_.fullPath();
  } else {
    edm::LogInfo("Level1")
        << "Using analytical compression and RCT decompression for CaloTPGTranscoderULUTs HCAL initialization";
    //std::unique_ptr<CaloTPGTranscoder> pTCoder(new CaloTPGTranscoderULUT());
    //return pTCoder;
  }
  //std::unique_ptr<CaloTPGTranscoder> pTCoder(new CaloTPGTranscoderULUT(ietal, ietah, ZS, LUTfactor, RCTLSB, nominal_gain, file1, file2));

  const auto& lutMetadata = iRecord.get(lutMetadataToken);
  const auto& theTrigTowerGeometry = iRecord.get(theTrigTowerGeometryToken);
  const auto& topoRecord = iRecord.getRecord<HcalLutMetadataRcd>();
  const auto& topo = topoRecord.get(topoToken);

  HcalLutMetadata fullLut{lutMetadata};
  fullLut.setTopo(&topo);

  auto pTCoder = std::make_unique<CaloTPGTranscoderULUT>(file1, file2);
  pTCoder->setup(fullLut, theTrigTowerGeometry, NCTScaleShift, RCTScaleShift, lsbQIE8, lsbQIE11, linearLUTs_);
  return pTCoder;
}

//define this as a plug-in
DEFINE_FWK_EVENTSETUP_MODULE(CaloTPGTranscoderULUTs);