CMS 3D CMS Logo

SiStripBadComponentsDQMServiceReader.cc
Go to the documentation of this file.
3 
8 
9 #include <iostream>
10 #include <cstdio>
11 #include <sys/time.h>
12 
13 #include <boost/lexical_cast.hpp>
14 
15 using namespace std;
16 
18  printdebug_(iConfig.getUntrackedParameter<bool>("printDebug",true)){}
19 
21 
23 {
24  //Retrieve tracker topology from geometry
26  iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
27  const TrackerTopology* const tTopo = tTopoHandle.product();
28 
29  uint32_t FedErrorMask = 1; // bit 0
30  uint32_t DigiErrorMask = 2; // bit 1
31  uint32_t ClusterErrorMask = 4; // bit 2
32 
33  edm::ESHandle<SiStripBadStrip> SiStripBadStrip_;
34  iSetup.get<SiStripBadStripRcd>().get(SiStripBadStrip_);
35  edm::LogInfo("SiStripBadComponentsDQMServiceReader") << "[SiStripBadComponentsDQMServiceReader::analyze] End Reading SiStripBadStrip" << std::endl;
36 
37  std::vector<uint32_t> detid;
38  SiStripBadStrip_->getDetIds(detid);
39 
40  std::stringstream ss;
41 
42  // ss << " detid" << " \t\t\t" << "FED error" << " \t" << "Digi test failed" << " \t" << "Cluster test failed" << std::endl;
43 
44  ss << "subdet layer stereo side \t detId \t\t Errors" << std::endl;
45 
46  for (size_t id=0;id<detid.size();id++) {
47  SiStripBadStrip::Range range=SiStripBadStrip_->getRange(detid[id]);
48 
49  for(int it=0;it<range.second-range.first;it++){
50  unsigned int value=(*(range.first+it));
51  ss << detIdToString(detid[id], tTopo) << "\t" << detid[id] << "\t";
52 
53  uint32_t flag = boost::lexical_cast<uint32_t>(SiStripBadStrip_->decode(value).flag);
54 
55  printError( ss, ((flag & FedErrorMask) == FedErrorMask), "Fed error, " );
56  printError( ss, ((flag & DigiErrorMask) == DigiErrorMask), "Digi error, " );
57  printError( ss, ((flag & ClusterErrorMask) == ClusterErrorMask), "Cluster error" );
58  ss << std::endl;
59 
60  if (printdebug_) {
61  ss << " firstBadStrip " << SiStripBadStrip_->decode(value).firstStrip << "\t "
62  << " NconsecutiveBadStrips " << SiStripBadStrip_->decode(value).range << "\t " // << std::endl;
63  << " flag " << SiStripBadStrip_->decode(value).flag << "\t "
64  << " packed integer " << std::hex << value << std::dec << "\t " << std::endl;
65  }
66  }
67  ss << std::endl;
68  }
69  edm::LogInfo("SiStripBadComponentsDQMServiceReader") << ss.str();
70 }
71 
72 void SiStripBadComponentsDQMServiceReader::printError( std::stringstream & ss, const bool error, const std::string & errorText)
73 {
74  if( error ) {
75  ss << errorText << "\t ";
76  }
77  else {
78  ss << "\t\t ";
79  }
80 }
81 
83 {
85  int layer = 0;
86  int stereo = 0;
87  int side = -1;
88 
89  // Using the operator[] if the element does not exist it is created with the default value. That is 0 for integral types.
90  switch (detid.subdetId()) {
92  {
93 
94  detector = "TIB";
95  layer = tTopo->tibLayer(detid.rawId());
96  stereo = tTopo->tibStereo(detid.rawId());
97  break;
98  }
100  {
101 
102  detector = "TOB";
103  layer = tTopo->tobLayer(detid.rawId());
104  stereo = tTopo->tobStereo(detid.rawId());
105  break;
106  }
108  {
109 
110  // is this module in TEC+ or TEC-?
111  side = tTopo->tecSide(detid.rawId());
112  detector = "TEC";
113  layer = tTopo->tecWheel(detid.rawId());
114  stereo = tTopo->tecStereo(detid.rawId());
115  break;
116  }
118  {
119 
120  // is this module in TID+ or TID-?
121  side = tTopo->tidSide(detid.rawId());
122  detector = "TID";
123  layer = tTopo->tidWheel(detid.rawId());
124  stereo = tTopo->tidStereo(detid.rawId());
125  break;
126  }
127  }
128  std::string name( detector + "\t" + std::to_string(layer) + "\t" + std::to_string(stereo) + "\t" );
129  if( side == 1 ) {
130  name += "-";
131  }
132  else if ( side == 2 ) {
133  name += "+";
134  }
135 // if( side != -1 ) {
136 // name += boost::lexical_cast<string>(side);
137 // }
138 
139  return name;
140 }
unsigned short range
void getDetIds(std::vector< uint32_t > &DetIds_) const
unsigned int tibLayer(const DetId &id) const
void printError(std::stringstream &ss, const bool error, const std::string &errorText)
uint32_t tobStereo(const DetId &id) const
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:47
unsigned int tidWheel(const DetId &id) const
unsigned int tidSide(const DetId &id) const
uint32_t tidStereo(const DetId &id) const
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:41
Definition: value.py:1
Definition: DetId.h:18
unsigned short firstStrip
std::string detIdToString(const DetId &detid, const TrackerTopology *tTopo)
const Range getRange(const uint32_t detID) const
std::pair< ContainerIterator, ContainerIterator > Range
T get() const
Definition: EventSetup.h:63
uint32_t tecStereo(const DetId &id) const
void analyze(const edm::Event &, const edm::EventSetup &) override
uint32_t tibStereo(const DetId &id) const
unsigned int tecWheel(const DetId &id) const
T const * product() const
Definition: ESHandle.h:86
data decode(const unsigned int &value) const
unsigned int tobLayer(const DetId &id) const
unsigned int tecSide(const DetId &id) const