HcalQIEDataCheck.cc

Go to the documentation of this file.

#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);
}

HcalQIEDataCheck::~HcalQIEDataCheck() {}

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

  edm::ESHandle<HcalQIEData> newQIEs;
  es.get<HcalQIEDataRcd>().get("update",newQIEs);
  const HcalQIEData* myNewQIEs = newQIEs.product();

  edm::ESHandle<HcalQIEData> refQIEs;
  es.get<HcalQIEDataRcd>().get("reference",refQIEs);
  const HcalQIEData* myRefQIEs = refQIEs.product();

  edm::ESHandle<HcalElectronicsMap> refEMap;
  es.get<HcalElectronicsMapRcd>().get("reference",refEMap);
  const HcalElectronicsMap* myRefEMap = refEMap.product();

  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);