CondTools/Hcal/plugins/HcalQIEDataCheck.cc - master - githubmirror/cmssw

CondTools/Hcal/plugins/HcalQIEDataCheck.cc
177 строк7 KB
Sunanda
Use ESGetToken for accessing objects from EventSetup in CondTools/Hcal
16 янв 2021, 00:26
16 янв 2021, 00:264f3b1f2
#include "CondTools/Hcal/interface/HcalQIEDataCheck.h"

HcalQIEDataCheck::HcalQIEDataCheck(edm::ParameterSet const& ps) {
  outfile = ps.getUntrackedParameter<std::string>("outFile", "null");
  dumprefs = ps.getUntrackedParameter<std::string>("dumpRefQIEsTo", "null");
  dumpupdate = ps.getUntrackedParameter<std::string>("dumpUpdateQIEsTo", "null");
  checkemapflag = ps.getUntrackedParameter<bool>("checkEmap", false);
  validateflag = ps.getUntrackedParameter<bool>("validateQIEs", false);
  //  epsilon = ps.getUntrackedParameter<double>("deltaQIE",0);
  m_tok1 = esConsumes<HcalQIEData, HcalQIEDataRcd>(edm::ESInputTag("", "update"));
  m_tok2 = esConsumes<HcalQIEData, HcalQIEDataRcd>(edm::ESInputTag("", "reference"));
  m_tokmap = esConsumes<HcalElectronicsMap, HcalElectronicsMapRcd>(edm::ESInputTag("", "reference"));
}

HcalQIEDataCheck::~HcalQIEDataCheck() {}

void HcalQIEDataCheck::analyze(const edm::Event& ev, const edm::EventSetup& es) {
  using namespace edm::eventsetup;

  const HcalQIEData* myNewQIEs = &es.getData(m_tok1);
  const HcalQIEData* myRefQIEs = &es.getData(m_tok2);
  const HcalElectronicsMap* myRefEMap = &es.getData(m_tokmap);

  if (dumpupdate != "null") {
    std::ofstream outStream(dumpupdate.c_str());
    std::cout << "--- Dumping QIEs - update ---" << std::endl;
    HcalDbASCIIIO::dumpObject(outStream, (*myNewQIEs));
  }
  if (dumprefs != "null") {
    std::ofstream outStream2(dumprefs.c_str());
    std::cout << "--- Dumping QIEs - reference ---" << std::endl;
    HcalDbASCIIIO::dumpObject(outStream2, (*myRefQIEs));
  }

  // first get the list of all channels from the update
  std::vector<DetId> listNewChan = myNewQIEs->getAllChannels();

  HcalQIEData* resultQIEs = new HcalQIEData(myRefQIEs->topo());
  std::vector<DetId> listRefChan = myRefQIEs->getAllChannels();
  std::vector<DetId>::iterator cell;

  if (validateflag) {
    for (std::vector<DetId>::iterator it = listRefChan.begin(); it != listRefChan.end(); it++) {
      DetId mydetid = *it;
      cell = std::find(listNewChan.begin(), listNewChan.end(), mydetid);
      if (cell == listNewChan.end())  // not present in new list
      {
        throw cms::Exception("DataDoesNotMatch") << "Value not found in reference" << std::endl;
      } else  // present in new list
      {
        const HcalQIECoder* first = myNewQIEs->getCoder(mydetid);
        const HcalQIECoder* second = myRefQIEs->getCoder(mydetid);
        {
          bool failflag = false;
          if (first->offset(0, 0) != second->offset(0, 0))
            failflag = true;
          if (first->offset(0, 1) != second->offset(0, 1))
            failflag = true;
          if (first->offset(0, 2) != second->offset(0, 2))
            failflag = true;
          if (first->offset(0, 3) != second->offset(0, 3))
            failflag = true;
          if (first->offset(1, 0) != second->offset(1, 0))
            failflag = true;
          if (first->offset(1, 1) != second->offset(1, 1))
            failflag = true;
          if (first->offset(1, 2) != second->offset(1, 2))
            failflag = true;
          if (first->offset(1, 3) != second->offset(1, 3))
            failflag = true;
          if (first->offset(2, 0) != second->offset(2, 0))
            failflag = true;
          if (first->offset(2, 1) != second->offset(2, 1))
            failflag = true;
          if (first->offset(2, 2) != second->offset(2, 2))
            failflag = true;
          if (first->offset(2, 3) != second->offset(2, 3))
            failflag = true;
          if (first->offset(3, 0) != second->offset(3, 0))
            failflag = true;
          if (first->offset(3, 1) != second->offset(3, 1))
            failflag = true;
          if (first->offset(3, 2) != second->offset(3, 2))
            failflag = true;
          if (first->offset(3, 3) != second->offset(3, 3))
            failflag = true;
          if (first->slope(0, 0) != second->slope(0, 0))
            failflag = true;
          if (first->slope(0, 1) != second->slope(0, 1))
            failflag = true;
          if (first->slope(0, 2) != second->slope(0, 2))
            failflag = true;
          if (first->slope(0, 3) != second->slope(0, 3))
            failflag = true;
          if (first->slope(1, 0) != second->slope(1, 0))
            failflag = true;
          if (first->slope(1, 1) != second->slope(1, 1))
            failflag = true;
          if (first->slope(1, 2) != second->slope(1, 2))
            failflag = true;
          if (first->slope(1, 3) != second->slope(1, 3))
            failflag = true;
          if (first->slope(2, 0) != second->slope(2, 0))
            failflag = true;
          if (first->slope(2, 1) != second->slope(2, 1))
            failflag = true;
          if (first->slope(2, 2) != second->slope(2, 2))
            failflag = true;
          if (first->slope(2, 3) != second->slope(2, 3))
            failflag = true;
          if (first->slope(3, 0) != second->slope(3, 0))
            failflag = true;
          if (first->slope(3, 1) != second->slope(3, 1))
            failflag = true;
          if (first->slope(3, 2) != second->slope(3, 2))
            failflag = true;
          if (first->slope(3, 3) != second->slope(3, 3))
            failflag = true;
          if (failflag)
            throw cms::Exception("DataDoesNotMatch") << "Values are do not match";
        }
        listNewChan.erase(cell);  // fix 25.02.08
      }
    }
    std::cout << "These are identical" << std::endl;
  }

  //  if(epsilon!=0){
  //implement compare qies -- different epsilon for slope and offset?
  //  }

  if (outfile != "null") {
    for (std::vector<DetId>::iterator it = listRefChan.begin(); it != listRefChan.end(); it++) {
      DetId mydetid = *it;
      cell = std::find(listNewChan.begin(), listNewChan.end(), mydetid);
      if (cell == listNewChan.end())  // not present in new list
      {
        const HcalQIECoder* myCoder = myRefQIEs->getCoder(mydetid);
        std::cout << "o";
        resultQIEs->addCoder(*myCoder);
      } else  // present in new list
      {
        const HcalQIECoder* myCoder = myNewQIEs->getCoder(mydetid);
        std::cout << "n";
        resultQIEs->addCoder(*myCoder);
        listNewChan.erase(cell);  // fix 25.02.08
      }
    }
    for (std::vector<DetId>::iterator it = listNewChan.begin(); it != listNewChan.end(); it++)  // fix 25.02.08
    {
      DetId mydetid = *it;
      const HcalQIECoder* myCoder = myNewQIEs->getCoder(mydetid);
      std::cout << "N";
      resultQIEs->addCoder(*myCoder);
    }

    std::ofstream outStream3(outfile.c_str());
    std::cout << "--- Dumping QIEs - the combined ones ---" << std::endl;
    resultQIEs->sort();
    HcalDbASCIIIO::dumpObject(outStream3, (*resultQIEs));
  }

  std::cout << std::endl;
  if (checkemapflag) {
    std::vector<DetId> listResult = resultQIEs->getAllChannels();
    // get the e-map list of channels
    std::vector<HcalGenericDetId> listEMap = myRefEMap->allPrecisionId();
    for (std::vector<HcalGenericDetId>::const_iterator it = listEMap.begin(); it != listEMap.end(); it++) {
      DetId mydetid = DetId(it->rawId());
      if (std::find(listResult.begin(), listResult.end(), mydetid) == listResult.end()) {
        std::cout << "Conditions not found for DetId = " << HcalGenericDetId(it->rawId()) << std::endl;
      }
    }
  }
}

DEFINE_FWK_MODULE(HcalQIEDataCheck);