SiStripBadStrip.cc

Go to the documentation of this file.

#include "CondFormats/SiStripObjects/interface/SiStripBadStrip.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "CondFormats/SiStripObjects/interface/SiStripDetSummary.h"

#include <algorithm>

bool SiStripBadStrip::put(const uint32_t& DetId, Range input) {
  // put in SiStripBadStrip::v_badstrips of DetId
  Registry::iterator p = std::lower_bound(indexes.begin(),indexes.end(),DetId,SiStripBadStrip::StrictWeakOrdering());
  if (p!=indexes.end() && p->detid==DetId){
    edm::LogError("SiStripBadStrip") << "[" << __PRETTY_FUNCTION__ << "] SiStripBadStrip for DetID " << DetId << " is already stored. Skippig this put" << std::endl;
    return false;
  }
  
  size_t sd= input.second-input.first;
  DetRegistry detregistry;
  detregistry.detid=DetId;
  detregistry.ibegin=v_badstrips.size();
  detregistry.iend=v_badstrips.size()+sd;
  indexes.insert(p,detregistry);

  v_badstrips.insert(v_badstrips.end(),input.first,input.second);
  return true;
}

const SiStripBadStrip::Range SiStripBadStrip::getRange(const uint32_t DetId) const {
  // get SiStripBadStrip Range of DetId
  
  RegistryIterator p = std::lower_bound(indexes.begin(),indexes.end(),DetId,SiStripBadStrip::StrictWeakOrdering());
  if (p==indexes.end()|| p->detid!=DetId) 
    return SiStripBadStrip::Range(v_badstrips.end(),v_badstrips.end()); 
  else {
    __builtin_prefetch((&v_badstrips.front())+p->ibegin);
    __builtin_prefetch((&v_badstrips.front())+p->ibegin+24);
    __builtin_prefetch((&v_badstrips.front())+p->iend-24);
    return SiStripBadStrip::Range(v_badstrips.begin()+p->ibegin,v_badstrips.begin()+p->iend);
  }
}

SiStripBadStrip::Range SiStripBadStrip::getRangeByPos(unsigned short pos) const {
  if (pos>indexes.size()) return Range(v_badstrips.end(),v_badstrips.end()); 
  auto p = indexes.begin()+pos;
  __builtin_prefetch((&v_badstrips.front())+p->ibegin);
  __builtin_prefetch((&v_badstrips.front())+p->ibegin+24);
  __builtin_prefetch((&v_badstrips.front())+p->iend-24);
  return Range(v_badstrips.begin()+p->ibegin,v_badstrips.begin()+p->iend);
}


void SiStripBadStrip::getDetIds(std::vector<uint32_t>& DetIds_) const {
  // returns vector of DetIds in map
  SiStripBadStrip::RegistryIterator begin = indexes.begin();
  SiStripBadStrip::RegistryIterator end   = indexes.end();
  for (SiStripBadStrip::RegistryIterator p=begin; p != end; ++p) {
    DetIds_.push_back(p->detid);
  }
}

void SiStripBadStrip::printSummary(std::stringstream & ss, const TrackerTopology* trackerTopo) const {
  SiStripDetSummary summaryBadModules{trackerTopo};
  SiStripDetSummary summaryBadStrips{trackerTopo};

  // Loop on the vector<DetRegistry> and take the bad modules and bad strips
  Registry::const_iterator it = indexes.begin();
  for( ; it!=indexes.end(); ++it ) {
    summaryBadModules.add(it->detid);
    summaryBadStrips.add(it->iend - it->ibegin);
  }
  ss << "Summary of bad modules in detector:" << std::endl;
  summaryBadModules.print(ss, false);
  ss << "Summary of bad strip in detectors:" << std::endl;
  summaryBadStrips.print(ss, false);
}

void SiStripBadStrip::printDebug(std::stringstream & ss, const TrackerTopology* /*trackerTopo*/) const {
  ss << "Printing all bad strips for all DetIds" << std::endl;
  // Loop on the vector<DetRegistry> and take the bad modules and bad strips
  Registry::const_iterator it = indexes.begin();
  for( ; it!=indexes.end(); ++it ) {
//    ss << "For DetId = " << it->detid << std::endl;
    SiStripBadStrip::Range range(getRange(it->detid));
    for( std::vector<unsigned int>::const_iterator badStrip = range.first;
         badStrip != range.second; ++badStrip ) {
      ss << "DetId="<<it->detid << " Strip=" << decode(*badStrip).firstStrip <<":"<<decode(*badStrip).range << " flag="<< decode(*badStrip).flag << std::endl;
    }
  }
}