HcalLogicalMap.cc

Go to the documentation of this file.

#include "DataFormats/HcalDetId/interface/HcalGenericDetId.h"
#include "DataFormats/HcalDetId/interface/HcalOtherDetId.h"
#include "Geometry/CaloTopology/interface/HcalTopology.h"
#include "CondFormats/HcalObjects/interface/HcalLogicalMap.h"

#include <string>
#include <cstring>
#include <iostream>
#include <fstream>
#include <sstream>
#include <cstdio>
#include <ctime>
#include <memory>

using namespace std;

HcalLogicalMap::HcalLogicalMap(const HcalTopology* topo,
                               std::vector<HBHEHFLogicalMapEntry>& HBHEHFEntries,
                               std::vector<HOHXLogicalMapEntry>& HOHXEntries,
                               std::vector<CALIBLogicalMapEntry>& CALIBEntries,
                               std::vector<ZDCLogicalMapEntry>& ZDCEntries,
                               std::vector<HTLogicalMapEntry>& HTEntries,
                               std::vector<uint32_t>& LinearIndex2Entry,
                               std::vector<uint32_t>& HbHash2Entry,
                               std::vector<uint32_t>& HeHash2Entry,
                               std::vector<uint32_t>& HfHash2Entry,
                               std::vector<uint32_t>& HtHash2Entry,
                               std::vector<uint32_t>& HoHash2Entry,
                               std::vector<uint32_t>& HxCalibHash2Entry,
                               std::vector<uint32_t>& ZdcHash2Entry)
    : topo_(topo)

{
  HBHEHFEntries_.resize(HBHEHFEntries.size());
  HOHXEntries_.resize(HOHXEntries.size());
  CALIBEntries_.resize(CALIBEntries.size());
  ZDCEntries_.resize(ZDCEntries.size());
  HTEntries_.resize(HTEntries.size());

  LinearIndex2Entry_.resize(LinearIndex2Entry.size());
  HbHash2Entry_.resize(HbHash2Entry.size());
  HeHash2Entry_.resize(HeHash2Entry.size());
  HfHash2Entry_.resize(HfHash2Entry.size());
  HtHash2Entry_.resize(HtHash2Entry.size());
  HoHash2Entry_.resize(HoHash2Entry.size());
  HxCalibHash2Entry_.resize(HxCalibHash2Entry.size());
  ZdcHash2Entry_.resize(ZdcHash2Entry.size());

  copy(HBHEHFEntries.begin(), HBHEHFEntries.end(), HBHEHFEntries_.begin());
  copy(HOHXEntries.begin(), HOHXEntries.end(), HOHXEntries_.begin());
  copy(CALIBEntries.begin(), CALIBEntries.end(), CALIBEntries_.begin());
  copy(ZDCEntries.begin(), ZDCEntries.end(), ZDCEntries_.begin());
  copy(HTEntries.begin(), HTEntries.end(), HTEntries_.begin());
  copy(LinearIndex2Entry.begin(), LinearIndex2Entry.end(), LinearIndex2Entry_.begin());
  copy(HbHash2Entry.begin(), HbHash2Entry.end(), HbHash2Entry_.begin());
  copy(HeHash2Entry.begin(), HeHash2Entry.end(), HeHash2Entry_.begin());
  copy(HfHash2Entry.begin(), HfHash2Entry.end(), HfHash2Entry_.begin());
  copy(HtHash2Entry.begin(), HtHash2Entry.end(), HtHash2Entry_.begin());
  copy(HoHash2Entry.begin(), HoHash2Entry.end(), HoHash2Entry_.begin());
  copy(HxCalibHash2Entry.begin(), HxCalibHash2Entry.end(), HxCalibHash2Entry_.begin());
  copy(ZdcHash2Entry.begin(), ZdcHash2Entry.end(), ZdcHash2Entry_.begin());
}

HcalLogicalMap::~HcalLogicalMap() {}

uint32_t HcalLogicalMap::makeEntryNumber(bool isvalid, int vectorid, int entry) {
  uint32_t answer = 0;
  answer |= isvalid;
  answer |= (vectorid & 0x3) << 1;
  answer |= entry << 3;
  return answer;
}

void HcalLogicalMap::printMap(unsigned int mapIOV) {
  using namespace std;

  static FILE* HBEFmap;
  static FILE* HOXmap;
  static FILE* CALIBmap;
  static FILE* ZDCmap;
  static FILE* HTmap;

  char tempbuff[30];

  stringstream mystream;
  string HBEFmapstr, HOXmapstr, CALIBmapstr, ZDCmapstr, HTmapstr;
  string date;
  string IOVlabel;

  time_t myTime;
  time(&myTime);

  strftime(tempbuff, 128, "%d.%b.%Y", localtime(&myTime));
  mystream << tempbuff;
  date = mystream.str();

  mystream.str("");
  if (mapIOV == 1)
    IOVlabel = "A";
  else if (mapIOV == 2)
    IOVlabel = "B";
  else if (mapIOV == 3)
    IOVlabel = "C";
  else if (mapIOV == 4)
    IOVlabel = "D";
  else
    IOVlabel = "E";

  HBEFmapstr = "./HCALmapHBEF_" + IOVlabel + ".txt";
  HOXmapstr = "./HCALmapHO_" + IOVlabel + ".txt";
  CALIBmapstr = "./HCALmapCALIB_" + IOVlabel + ".txt";
  ZDCmapstr = "./ZDCmap_" + IOVlabel + ".txt";
  HTmapstr = "./HCALmapHT_" + IOVlabel + ".txt";

  //  HBEFmapstr  = "./HCALmapHBEF_"+date+".txt";
  //  HOXmapstr   = "./HCALmapHO_"+date+"_"+IOVlabel+".txt";
  //  CALIBmapstr = "./HCALmapCALIB_"+date+".txt";
  //  ZDCmapstr   = "./ZDCmap_"+date+".txt";
  //  HTmapstr    = "./HCALmapHT_"+date+".txt";

  HBEFmap = fopen(HBEFmapstr.c_str(), "w");
  HOXmap = fopen(HOXmapstr.c_str(), "w");
  CALIBmap = fopen(CALIBmapstr.c_str(), "w");
  ZDCmap = fopen(ZDCmapstr.c_str(), "w");
  HTmap = fopen(HTmapstr.c_str(), "w");
  /**********************/

  if (HBEFmap) {
    fprintf(HBEFmap, "## file created %s ##\n", date.c_str());
    printHBEFMap(HBEFmap);
  } else
    cout << HBEFmapstr << " not found!" << endl;

  if (HOXmap) {
    fprintf(HOXmap, "## file created %s ##\n", date.c_str());
    printHOXMap(HOXmap);
  } else
    cout << HOXmapstr << " not found!" << endl;

  if (CALIBmap) {
    fprintf(CALIBmap, "## file created %s ##\n", date.c_str());
    printCalibMap(CALIBmap);
  } else
    cout << CALIBmapstr << " not found!" << endl;

  if (ZDCmap) {
    fprintf(ZDCmap, "## file created %s ##\n", date.c_str());
    printZDCMap(ZDCmap);
  } else
    cout << ZDCmapstr << " not found!" << endl;

  if (HTmap) {
    fprintf(HTmap, "## file created %s ##\n", date.c_str());
    printHTMap(HTmap);
  } else
    cout << HTmapstr << " not found!" << endl;
}

/**************/
std::unique_ptr<HcalElectronicsMap> HcalLogicalMap::generateHcalElectronicsMap() {
  HcalElectronicsMapAddons::Helper theemapHelper;

  for (std::vector<HBHEHFLogicalMapEntry>::iterator it = HBHEHFEntries_.begin(); it != HBHEHFEntries_.end(); ++it) {
    theemapHelper.mapEId2chId(it->getHcalElectronicsId(), it->getDetId());
  }
  for (std::vector<HOHXLogicalMapEntry>::iterator it = HOHXEntries_.begin(); it != HOHXEntries_.end(); ++it) {
    theemapHelper.mapEId2chId(it->getHcalElectronicsId(), it->getDetId());
  }
  for (std::vector<CALIBLogicalMapEntry>::iterator it = CALIBEntries_.begin(); it != CALIBEntries_.end(); ++it) {
    theemapHelper.mapEId2chId(it->getHcalElectronicsId(), it->getDetId());
  }
  for (std::vector<ZDCLogicalMapEntry>::iterator it = ZDCEntries_.begin(); it != ZDCEntries_.end(); ++it) {
    theemapHelper.mapEId2chId(it->getHcalElectronicsId(), it->getDetId());
  }
  for (std::vector<HTLogicalMapEntry>::iterator it = HTEntries_.begin(); it != HTEntries_.end(); ++it) {
    theemapHelper.mapEId2tId(it->getHcalTrigElectronicsId(), it->getDetId());
  }
  return std::make_unique<HcalElectronicsMap>(theemapHelper);
}
/**************/

void HcalLogicalMap::printHBEFMap(FILE* hbefmapfile) {
  int titlecounter = 0;

  for (std::vector<HBHEHFLogicalMapEntry>::iterator it = HBHEHFEntries_.begin(); it != HBHEHFEntries_.end(); ++it) {
    titlecounter = titlecounter % 21;
    if (titlecounter == 0) {
      fprintf(hbefmapfile, "#   side    eta    phi   dphi  depth    det     rbx  wedge     rm  pixel   qie    adc");
      fprintf(hbefmapfile, "  rm_fi  fi_ch  crate    htr   fpga  htr_fi  dcc_sl  spigo    dcc    slb  slbin  slbin2");
      fprintf(hbefmapfile, "           slnam    rctcra rctcar rctcon               rctnam     fedid\n");
    }
    titlecounter++;
    fprintf(hbefmapfile, "%s", it->printLMapLine());
  }
}

void HcalLogicalMap::printHOXMap(FILE* hoxmapfile) {
  int titlecounter = 0;

  for (std::vector<HOHXLogicalMapEntry>::iterator it = HOHXEntries_.begin(); it != HOHXEntries_.end(); ++it) {
    titlecounter = titlecounter % 21;
    if (titlecounter == 0) {
      fprintf(hoxmapfile, "#   side    eta    phi   dphi  depth    det     rbx  sector    rm  pixel   qie    adc");
      fprintf(hoxmapfile, "  rm_fi  fi_ch let_code  crate    htr   fpga  htr_fi  dcc_sl  spigo    dcc  fedid\n");
    }
    titlecounter++;
    fprintf(hoxmapfile, "%s", it->printLMapLine());
  }
}

void HcalLogicalMap::printCalibMap(FILE* calibmapfile) {
  int titlecounter = 0;

  for (std::vector<CALIBLogicalMapEntry>::iterator it = CALIBEntries_.begin(); it != CALIBEntries_.end(); ++it) {
    titlecounter = titlecounter % 21;
    if (titlecounter == 0) {
      fprintf(calibmapfile, "#   side    eta    phi   dphi    det     rbx  sector  rm  rm_fi ");
      fprintf(calibmapfile, " fi_ch  crate  htr  fpga  htr_fi  dcc_sl  spigo  dcc  fedid  ch_type      name\n");
    }
    titlecounter++;
    fprintf(calibmapfile, "%s", it->printLMapLine());
  }
}

void HcalLogicalMap::printZDCMap(FILE* zdcmapfile) {
  int titlecounter = 0;

  for (std::vector<ZDCLogicalMapEntry>::iterator it = ZDCEntries_.begin(); it != ZDCEntries_.end(); ++it) {
    titlecounter = titlecounter % 21;
    if (titlecounter == 0) {
      fprintf(zdcmapfile, "#  side  x  y  dx  depth     det  det_ch  cable  rm  qie ");
      fprintf(zdcmapfile, " adc  rm_fi  fi_ch  crate  htr  fpga  htr_fi  dcc_sl  spigo  dcc  fedid\n");
    }
    titlecounter++;
    fprintf(zdcmapfile, "%s", it->printLMapLine());
  }
}

void HcalLogicalMap::printHTMap(FILE* htmapfile) {
  int titlecounter = 0;

  for (std::vector<HTLogicalMapEntry>::iterator it = HTEntries_.begin(); it != HTEntries_.end(); ++it) {
    titlecounter = titlecounter % 21;
    if (titlecounter == 0) {
      fprintf(htmapfile, "#  side  eta  phi  dphi  depth  det   wedge  crate");
      fprintf(htmapfile, "  htr  fpga  dcc_sl  spigo  dcc  slb  slbin  slbin2  nDat    ");
      fprintf(htmapfile, "     slnam  rctcra  rctcar  rctcon            rctnam  fedid\n");
    }
    titlecounter++;
    fprintf(htmapfile, "%s", it->printLMapLine());
  }
}

const DetId HcalLogicalMap::getDetId(const HcalElectronicsId& eid) {
  const uint32_t entry = LinearIndex2Entry_.at(eid.linearIndex());
  if (!(entry & 0x1))
    return DetId(0);
  if (((entry >> 1) & 0x3) == 0)
    return HBHEHFEntries_.at(entry >> 3).getDetId();
  if (((entry >> 1) & 0x3) == 1)
    return HOHXEntries_.at(entry >> 3).getDetId();
  if (((entry >> 1) & 0x3) == 2)
    return CALIBEntries_.at(entry >> 3).getDetId();
  if (((entry >> 1) & 0x3) == 3)
    return ZDCEntries_.at(entry >> 3).getDetId();
  return DetId(0);
}

const HcalFrontEndId HcalLogicalMap::getHcalFrontEndId(const DetId& did) {
  const HcalGenericDetId hgdi(did);

  const HcalGenericDetId::HcalGenericSubdetector hgsd = hgdi.genericSubdet();
  if (hgsd == HcalGenericDetId::HcalGenBarrel) {
    const int hashedId = topo_->detId2denseIdHB(did);
    const uint32_t entry = HbHash2Entry_.at(hashedId) - 1;
    return HBHEHFEntries_.at(entry).getHcalFrontEndId();
  } else if (hgsd == HcalGenericDetId::HcalGenEndcap) {
    const int hashedId = topo_->detId2denseIdHE(did);
    const uint32_t entry = HeHash2Entry_.at(hashedId) - 1;
    return HBHEHFEntries_.at(entry).getHcalFrontEndId();
  } else if (hgsd == HcalGenericDetId::HcalGenForward) {
    const int hashedId = topo_->detId2denseIdHF(did);
    const uint32_t entry = HfHash2Entry_.at(hashedId) - 1;
    return HBHEHFEntries_.at(entry).getHcalFrontEndId();
  } else if (hgsd == HcalGenericDetId::HcalGenOuter) {
    const int hashedId = topo_->detId2denseIdHO(did);
    const uint32_t entry = HoHash2Entry_.at(hashedId) - 1;
    return HOHXEntries_.at(entry).getHcalFrontEndId();
  } else if (hgsd == HcalGenericDetId::HcalGenCalibration) {
    HcalCalibDetId hcid(did);
    if (hcid.calibFlavor() == HcalCalibDetId::HOCrosstalk) {
      const int hashedId = topo_->detId2denseIdCALIB(did);
      const uint32_t entry = HxCalibHash2Entry_.at(hashedId) - 1;
      return HOHXEntries_.at(entry).getHcalFrontEndId();
    } else if (hcid.calibFlavor() == HcalCalibDetId::CalibrationBox) {
      const int hashedId = topo_->detId2denseIdCALIB(did);
      const uint32_t entry = HxCalibHash2Entry_.at(hashedId) - 1;
      return CALIBEntries_.at(entry).getHcalFrontEndId();
    }
  }
  return HcalFrontEndId(0);
}

void HcalLogicalMap::checkIdFunctions() {
  int HBHEHF_EID_pass = 0;
  int HBHEHF_EID_fail = 0;
  int HOHX_EID_pass = 0;
  int HOHX_EID_fail = 0;
  int CALIB_EID_pass = 0;
  int CALIB_EID_fail = 0;
  int ZDC_EID_pass = 0;
  int ZDC_EID_fail = 0;

  int HBHEHF_FEID_pass = 0;
  int HBHEHF_FEID_fail = 0;
  int HOHX_FEID_pass = 0;
  int HOHX_FEID_fail = 0;
  int CALIB_FEID_pass = 0;
  int CALIB_FEID_fail = 0;

  cout << "\nRunning the id function checker..." << endl;

  for (std::vector<HBHEHFLogicalMapEntry>::iterator it = HBHEHFEntries_.begin(); it != HBHEHFEntries_.end(); ++it) {
    const HcalElectronicsId heid = it->getHcalElectronicsId();
    const DetId did0 = it->getDetId();
    const DetId did1 = getDetId(heid);
    if (did0 == did1)
      HBHEHF_EID_pass++;
    else
      HBHEHF_EID_fail++;

    const HcalFrontEndId hfeid0 = it->getHcalFrontEndId();
    const HcalFrontEndId hfeid1 = getHcalFrontEndId(did0);
    if (hfeid0 == hfeid1)
      HBHEHF_FEID_pass++;
    else
      HBHEHF_FEID_fail++;
  }
  for (std::vector<HOHXLogicalMapEntry>::iterator it = HOHXEntries_.begin(); it != HOHXEntries_.end(); ++it) {
    const HcalElectronicsId heid = it->getHcalElectronicsId();
    const DetId did0 = it->getDetId();
    const DetId did1 = getDetId(heid);
    if (did0 == did1)
      HOHX_EID_pass++;
    else
      HOHX_EID_fail++;

    const HcalFrontEndId hfeid0 = it->getHcalFrontEndId();
    const HcalFrontEndId hfeid1 = getHcalFrontEndId(did0);
    if (hfeid0 == hfeid1)
      HOHX_FEID_pass++;
    else
      HOHX_FEID_fail++;
  }
  for (std::vector<CALIBLogicalMapEntry>::iterator it = CALIBEntries_.begin(); it != CALIBEntries_.end(); ++it) {
    const HcalElectronicsId heid = it->getHcalElectronicsId();
    const DetId did0 = it->getDetId();
    const DetId did1 = getDetId(heid);
    if (did0 == did1)
      CALIB_EID_pass++;
    else
      CALIB_EID_fail++;

    const HcalFrontEndId hfeid0 = it->getHcalFrontEndId();
    const HcalFrontEndId hfeid1 = getHcalFrontEndId(did0);
    if (hfeid0 == hfeid1)
      CALIB_FEID_pass++;
    else
      CALIB_FEID_fail++;
  }
  for (std::vector<ZDCLogicalMapEntry>::iterator it = ZDCEntries_.begin(); it != ZDCEntries_.end(); ++it) {
    const HcalElectronicsId heid = it->getHcalElectronicsId();
    const DetId did0 = it->getDetId();
    const DetId did1 = getDetId(heid);
    if (did0 == did1)
      ZDC_EID_pass++;
    else
      ZDC_EID_fail++;
  }

  cout << "Checking detIds from electronics ids..." << endl;
  cout << "HBHEHF EID (pass,fail) = (" << HBHEHF_EID_pass << "," << HBHEHF_EID_fail << ")" << endl;
  cout << "HOHX EID (pass,fail) = (" << HOHX_EID_pass << "," << HOHX_EID_fail << ")" << endl;
  cout << "CALIB EID (pass,fail) = (" << CALIB_EID_pass << "," << CALIB_EID_fail << ")" << endl;
  cout << "ZDC EID (pass,fail) = (" << ZDC_EID_pass << "," << ZDC_EID_fail << ")" << endl;
  cout << endl;
  cout << "Checking frontEndIds from electronics ids..." << endl;
  cout << "HBHEHF FEID (pass,fail) = (" << HBHEHF_FEID_pass << "," << HBHEHF_FEID_fail << ")" << endl;
  cout << "HOHX FEID (pass,fail) = (" << HOHX_FEID_pass << "," << HOHX_FEID_fail << ")" << endl;
  cout << "CALIB FEID (pass,fail) = (" << CALIB_FEID_pass << "," << CALIB_FEID_fail << ")" << endl;
}

void HcalLogicalMap::checkHashIds() {
  std::vector<int> HB_Hashes_;       // index 0
  std::vector<int> HE_Hashes_;       // index 1
  std::vector<int> HF_Hashes_;       // index 2
  std::vector<int> HO_Hashes_;       // index 3
  std::vector<int> CALIBHX_Hashes_;  // index 4
  std::vector<int> ZDC_Hashes_;      // index 5
  std::vector<int> HT_Hashes_;       // index 6

  int numfails[7] = {0, 0, 0, 0, 0, 0, 0};
  int numpass[7] = {0, 0, 0, 0, 0, 0, 0};
  int numnotdense[7] = {0, 0, 0, 0, 0, 0, 0};

  cout << "\nRunning the hash checker for detIds..." << endl;
  for (std::vector<HBHEHFLogicalMapEntry>::iterator it = HBHEHFEntries_.begin(); it != HBHEHFEntries_.end(); ++it) {
    if (it->getDetId().subdetId() == HcalBarrel) {
      HB_Hashes_.push_back(topo_->detId2denseIdHB(it->getDetId()));
    } else if (it->getDetId().subdetId() == HcalEndcap) {
      HE_Hashes_.push_back(topo_->detId2denseIdHE(it->getDetId()));
    } else if (it->getDetId().subdetId() == HcalForward) {
      HF_Hashes_.push_back(topo_->detId2denseIdHF(it->getDetId()));
    }
  }
  for (std::vector<HOHXLogicalMapEntry>::iterator it = HOHXEntries_.begin(); it != HOHXEntries_.end(); ++it) {
    if (HcalGenericDetId(it->getDetId().rawId()).isHcalCalibDetId()) {
      CALIBHX_Hashes_.push_back(topo_->detId2denseIdCALIB(it->getDetId()));
    } else {
      HO_Hashes_.push_back(topo_->detId2denseIdHO(it->getDetId()));
    }
  }
  for (std::vector<CALIBLogicalMapEntry>::iterator it = CALIBEntries_.begin(); it != CALIBEntries_.end(); ++it) {
    CALIBHX_Hashes_.push_back(topo_->detId2denseIdCALIB(it->getDetId()));
  }
  for (std::vector<ZDCLogicalMapEntry>::iterator it = ZDCEntries_.begin(); it != ZDCEntries_.end(); ++it) {
    ZDC_Hashes_.push_back(HcalZDCDetId(it->getDetId()).denseIndex());
  }
  for (std::vector<HTLogicalMapEntry>::iterator it = HTEntries_.begin(); it != HTEntries_.end(); ++it) {
    HT_Hashes_.push_back(topo_->detId2denseIdHT(it->getDetId()));
  }

  sort(HB_Hashes_.begin(), HB_Hashes_.end());
  sort(HE_Hashes_.begin(), HE_Hashes_.end());
  sort(HF_Hashes_.begin(), HF_Hashes_.end());
  sort(HO_Hashes_.begin(), HO_Hashes_.end());
  sort(CALIBHX_Hashes_.begin(), CALIBHX_Hashes_.end());
  sort(ZDC_Hashes_.begin(), ZDC_Hashes_.end());
  sort(HT_Hashes_.begin(), HT_Hashes_.end());

  for (unsigned int i = 0; i < HB_Hashes_.size() - 1; i++) {
    int diff = HB_Hashes_.at(i + 1) - HB_Hashes_.at(i);
    if (diff == 0)
      numfails[0]++;
    else if (diff > 1)
      numnotdense[0]++;
    else
      numpass[0]++;
  }
  for (unsigned int i = 0; i < HE_Hashes_.size() - 1; i++) {
    int diff = HE_Hashes_.at(i + 1) - HE_Hashes_.at(i);
    if (diff == 0)
      numfails[1]++;
    else if (diff > 1)
      numnotdense[1]++;
    else
      numpass[1]++;
  }
  for (unsigned int i = 0; i < HF_Hashes_.size() - 1; i++) {
    int diff = HF_Hashes_.at(i + 1) - HF_Hashes_.at(i);
    if (diff == 0)
      numfails[2]++;
    else if (diff > 1)
      numnotdense[2]++;
    else
      numpass[2]++;
  }
  for (unsigned int i = 0; i < HO_Hashes_.size() - 1; i++) {
    int diff = HO_Hashes_.at(i + 1) - HO_Hashes_.at(i);
    if (diff == 0)
      numfails[3]++;
    else if (diff > 1)
      numnotdense[3]++;
    else
      numpass[3]++;
  }
  for (unsigned int i = 0; i < CALIBHX_Hashes_.size() - 1; i++) {
    int diff = CALIBHX_Hashes_.at(i + 1) - CALIBHX_Hashes_.at(i);
    if (diff == 0)
      numfails[4]++;
    else if (diff > 1)
      numnotdense[4]++;
    else
      numpass[4]++;
  }
  for (unsigned int i = 0; i < ZDC_Hashes_.size() - 1; i++) {
    int diff = ZDC_Hashes_.at(i + 1) - ZDC_Hashes_.at(i);
    if (diff == 0)
      numfails[5]++;
    else if (diff > 1)
      numnotdense[5]++;
    else
      numpass[5]++;
  }
  for (unsigned int i = 0; i < HT_Hashes_.size() - 1; i++) {
    int diff = HT_Hashes_.at(i + 1) - HT_Hashes_.at(i);
    if (diff == 0)
      numfails[6]++;
    else if (diff > 1)
      numnotdense[6]++;
    else
      numpass[6]++;
  }
  cout << "HB HashIds (pass, collisions, non-dense) = (" << numpass[0] << "," << numfails[0] << "," << numnotdense[0]
       << ")" << endl;
  cout << "HE HashIds (pass, collisions, non-dense) = (" << numpass[1] << "," << numfails[1] << "," << numnotdense[1]
       << ")" << endl;
  cout << "HF HashIds (pass, collisions, non-dense) = (" << numpass[2] << "," << numfails[2] << "," << numnotdense[2]
       << ")" << endl;
  cout << "HO HashIds (pass, collisions, non-dense) = (" << numpass[3] << "," << numfails[3] << "," << numnotdense[3]
       << ")" << endl;
  cout << "CALIB/HX HashIds (pass, collisions, non-dense) = (" << numpass[4] << "," << numfails[4] << ","
       << numnotdense[4] << ")" << endl;
  cout << "ZDC HashIds (pass, collisions, non-dense) = (" << numpass[5] << "," << numfails[5] << "," << numnotdense[5]
       << ")" << endl;
  cout << "HT HashIds (pass, collisions, non-dense) = (" << numpass[6] << "," << numfails[6] << "," << numnotdense[6]
       << ")" << endl;
}

void HcalLogicalMap::checkElectronicsHashIds() {
  std::vector<int> Electronics_Hashes_;

  int numfails = 0;
  int numpass = 0;
  int numnotdense = 0;

  cout << "\nRunning the hash checker for electronics Ids..." << endl;
  for (std::vector<HBHEHFLogicalMapEntry>::iterator it = HBHEHFEntries_.begin(); it != HBHEHFEntries_.end(); ++it) {
    Electronics_Hashes_.push_back((it->getHcalElectronicsId()).linearIndex());
  }
  for (std::vector<ZDCLogicalMapEntry>::iterator it = ZDCEntries_.begin(); it != ZDCEntries_.end(); ++it) {
    Electronics_Hashes_.push_back((it->getHcalElectronicsId()).linearIndex());
  }
  for (std::vector<CALIBLogicalMapEntry>::iterator it = CALIBEntries_.begin(); it != CALIBEntries_.end(); ++it) {
    Electronics_Hashes_.push_back((it->getHcalElectronicsId()).linearIndex());
  }
  for (std::vector<HOHXLogicalMapEntry>::iterator it = HOHXEntries_.begin(); it != HOHXEntries_.end(); ++it) {
    Electronics_Hashes_.push_back((it->getHcalElectronicsId()).linearIndex());
  }
  for (std::vector<HTLogicalMapEntry>::iterator it = HTEntries_.begin(); it != HTEntries_.end(); ++it) {
    Electronics_Hashes_.push_back((it->getHcalTrigElectronicsId()).linearIndex());
  }

  sort(Electronics_Hashes_.begin(), Electronics_Hashes_.end());

  for (unsigned int i = 0; i < Electronics_Hashes_.size() - 1; i++) {
    int diff = Electronics_Hashes_.at(i + 1) - Electronics_Hashes_.at(i);
    if (diff == 0)
      numfails++;
    else if (diff > 1)
      numnotdense++;
    else
      numpass++;
  }
  cout << "Electronics Id linearIndex (pass, collisions, nondense) = (" << numpass << "," << numfails << ","
       << numnotdense << ")" << endl;
}