SiPhase2BadStripChannelReader

Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94
// user include files
#include "CondFormats/DataRecord/interface/SiPhase2OuterTrackerCondDataRecords.h"
#include "CondFormats/SiStripObjects/interface/SiStripBadStrip.h"
#include "DataFormats/Phase2TrackerDigi/interface/Phase2TrackerDigi.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/global/EDAnalyzer.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"

//
//
// class decleration
//
class SiPhase2BadStripChannelReader : public edm::global::EDAnalyzer<> {
public:
  explicit SiPhase2BadStripChannelReader(const edm::ParameterSet&);
  ~SiPhase2BadStripChannelReader() override = default;
  void analyze(edm::StreamID, edm::Event const&, edm::EventSetup const&) const override;

  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

private:
  const bool printverbose_;
  const int printdebug_;
  const std::string label_;
  const edm::ESGetToken<SiStripBadStrip, SiPhase2OuterTrackerBadStripRcd> badStripToken_;
};

SiPhase2BadStripChannelReader::SiPhase2BadStripChannelReader(const edm::ParameterSet& iConfig)
    : printverbose_(iConfig.getUntrackedParameter<bool>("printVerbose", false)),
      printdebug_(iConfig.getUntrackedParameter<int>("printDebug", 1)),
      label_(iConfig.getUntrackedParameter<std::string>("label", "")),
      badStripToken_(esConsumes(edm::ESInputTag{"", label_})) {}

void SiPhase2BadStripChannelReader::analyze(edm::StreamID,
                                            edm::Event const& iEvent,
                                            edm::EventSetup const& iSetup) const {
  const auto& payload = &iSetup.getData(badStripToken_);

  edm::LogInfo("SiPhase2BadStripChannelReader")
      << "[SiPhase2BadStripChannelReader::analyze] End Reading SiStriBadStrip with label " << label_ << std::endl;

  const TrackerTopology* ttopo = nullptr;

  std::stringstream ss;
  if (printverbose_) {
    payload->printDebug(ss, ttopo);
  }

  std::vector<uint32_t> detIds;
  payload->getDetIds(detIds);

  int countMessages{0};
  for (const auto& d : detIds) {
    SiStripBadStrip::Range range = payload->getRange(d);
    for (std::vector<unsigned int>::const_iterator badChannel = range.first; badChannel != range.second; ++badChannel) {
      const auto& firstStrip = payload->decodePhase2(*badChannel).firstStrip;
      const auto& range = payload->decodePhase2(*badChannel).range;

      if (printverbose_) {
        ss << "DetId= " << d << " Channel=" << payload->decodePhase2(*badChannel).firstStrip << ":"
           << payload->decodePhase2(*badChannel).range << std::endl;
      }

      if (countMessages < printdebug_) {
        for (unsigned int index = 0; index < range; index++) {
          std::pair<int, int> badStrip = Phase2TrackerDigi::channelToPixel(firstStrip + index);
          ss << "DetId= " << d << " Channel= " << firstStrip + index << " -> strip (row,col)=(" << badStrip.first << ","
             << badStrip.second << ")" << std::endl;
        }
      }  // if doesn't exceed the maximum printout level
    }  // loop over the range
    countMessages++;
  }  // loop over the detids

  edm::LogInfo("SiPhase2BadStripChannelReader") << ss.str() << std::endl;
}

// ------------ method fills 'descriptions' with the allowed parameters for the module  ------------
void SiPhase2BadStripChannelReader::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  edm::ParameterSetDescription desc;
  desc.setComment("Module to read SiStripBadStrip Payloads through SiPhase2OuterTrackerBadStripRcd");
  desc.addUntracked<bool>("printVerbose", false)->setComment("if active, very verbose output");
  desc.addUntracked<int>("printDebug", 1)->setComment("tunes the amount of debug level print-outs");
  desc.addUntracked<std::string>("label", "")->setComment("label from which to read the payload");
  descriptions.addWithDefaultLabel(desc);
}

//define this as a plug-in
DEFINE_FWK_MODULE(SiPhase2BadStripChannelReader);