CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EcalDQMStatusReader.cc
Go to the documentation of this file.
1 /*
2  * \file EcalDQMStatusReader.cc
3  *
4  * \author G. Della Ricca
5  *
6 */
7 
8 #include <fstream>
9 #include <iostream>
10 #include <string>
11 #include <cstring>
12 #include <time.h>
13 #include <unistd.h>
14 
15 #include "FWCore/Framework/interface/ESHandle.h"
16 #include "CondFormats/DataRecord/interface/EcalDQMChannelStatusRcd.h"
17 #include "CondFormats/DataRecord/interface/EcalDQMTowerStatusRcd.h"
18 
19 #include "DataFormats/EcalDetId/interface/EBDetId.h"
20 #include "DataFormats/EcalDetId/interface/EEDetId.h"
21 #include "DataFormats/EcalDetId/interface/EcalTrigTowerDetId.h"
22 #include "DataFormats/EcalDetId/interface/EcalScDetId.h"
23 
24 #include "DQM/EcalCommon/interface/Numbers.h"
25 
26 #include "DQM/EcalCommon/interface/EcalDQMStatusDictionary.h"
27 
28 #include "DQM/EcalCommon/interface/EcalDQMStatusReader.h"
29 
30 EcalDQMStatusReader::EcalDQMStatusReader(const edm::ParameterSet& ps) {
31 
32  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
33 
34 }
35 
36 void EcalDQMStatusReader::beginRun(const edm::Run& r, const edm::EventSetup& c) {
37 
38  Numbers::initGeometry(c, verbose_);
39 
40  std::vector<EcalDQMStatusDictionary::codeDef> dictionary;
41  EcalDQMStatusDictionary::getDictionary( dictionary );
42 
43  const EcalDQMChannelStatus* channelStatus = 0;
44  if ( c.find( edm::eventsetup::EventSetupRecordKey::makeKey< EcalDQMChannelStatusRcd >() ) ) {
45  edm::ESHandle< EcalDQMChannelStatus > handle;
46  c.get< EcalDQMChannelStatusRcd >().get(handle);
47  if ( handle.isValid() ) channelStatus = handle.product();
48  }
49 
50  const EcalDQMTowerStatus* towerStatus = 0;
51  if ( c.find( edm::eventsetup::EventSetupRecordKey::makeKey< EcalDQMTowerStatusRcd >() ) ) {
52  edm::ESHandle< EcalDQMTowerStatus > handle;
53  c.get< EcalDQMTowerStatusRcd >().get(handle);
54  if ( handle.isValid() ) towerStatus = handle.product();
55  }
56 
57  // barrel
58  for ( int ism=1; ism<=36; ism++ ) {
59  for ( int ic=1; ic<=1700; ic++ ) {
60  int jsm = Numbers::iSM(ism, EcalBarrel);
61  EBDetId id(jsm, ic, EBDetId::SMCRYSTALMODE);
62  if ( channelStatus ) {
63  EcalDQMChannelStatus::const_iterator it = channelStatus->find( id.rawId() );
64  if ( it != channelStatus->end() ) {
65  if ( it->getStatusCode() != 0 ) {
66  if ( verbose_ ) std::cout << "# EB:channel, ic=" << ic << " hi=" << id.hashedIndex() << " status=" << it->getStatusCode() << std::endl;
67  std::vector<EcalDQMStatusDictionary::codeDef> codes;
68  EcalDQMStatusDictionary::getCodes( codes, it->getStatusCode() );
69  for ( unsigned int i=0; i<codes.size(); i++ ) {
70  std::cout << "Crystal " << Numbers::sEB(ism) << " " << ic << " " << codes[i].desc << std::endl;
71  }
72  }
73  }
74  }
75  }
76  }
77 
78  for ( int ix=1; ix<=17; ix++ ) {
79  for ( int iy=1; iy<=72; iy++ ) {
80  if ( EcalTrigTowerDetId::validDetId(+1, EcalBarrel, ix, iy) ) {
81  EcalTrigTowerDetId id(+1, EcalBarrel, ix, iy);
82  if ( towerStatus ) {
83  EcalDQMTowerStatus::const_iterator it = towerStatus->find( id.rawId() );
84  if ( it != towerStatus->end() ) {
85  if ( it->getStatusCode() != 0 ) {
86  if ( verbose_ ) std::cout << "# EB:tower, tt=" << Numbers::iTT(id) << " hi=" << id.hashedIndex() << " status=" << it->getStatusCode() << std::endl;
87  std::vector<EcalDQMStatusDictionary::codeDef> codes;
88  EcalDQMStatusDictionary::getCodes( codes, it->getStatusCode() );
89  for ( unsigned int i=0; i<codes.size(); i++ ) {
90  std::cout << "TT " << Numbers::sEB(Numbers::iSM(id)) << " " << Numbers::iTT(id) << " " << codes[i].desc << std::endl;
91  }
92  }
93  }
94  }
95  }
96  if ( EcalTrigTowerDetId::validDetId(-1, EcalBarrel, ix, iy) ) {
97  EcalTrigTowerDetId id(-1, EcalBarrel, ix, iy);
98  if ( towerStatus ) {
99  EcalDQMTowerStatus::const_iterator it = towerStatus->find( id.rawId() );
100  if ( it != towerStatus->end() ) {
101  if ( it->getStatusCode() != 0 ) {
102  if ( verbose_ ) std::cout << "# EB:tower, tt=" << Numbers::iTT(id) << " hi=" << id.hashedIndex() << " status=" << it->getStatusCode() << std::endl;
103  std::vector<EcalDQMStatusDictionary::codeDef> codes;
104  EcalDQMStatusDictionary::getCodes( codes, it->getStatusCode() );
105  for ( unsigned int i=0; i<codes.size(); i++ ) {
106  std::cout << "TT " << Numbers::sEB(Numbers::iSM(id)) << " " << Numbers::iTT(id) << " " << codes[i].desc << std::endl;
107  }
108  }
109  }
110  }
111  }
112  }
113  }
114 
115  // endcap
116  for ( int ix=1; ix<=100; ix++ ) {
117  for ( int iy=1; iy<=100; iy++ ) {
118  if ( EEDetId::validDetId(ix, iy, +1) ) {
119  EEDetId id(ix, iy, +1);
120  if ( channelStatus ) {
121  EcalDQMChannelStatus::const_iterator it = channelStatus->find( id.rawId() );
122  if ( it != channelStatus->end() ) {
123  if ( it->getStatusCode() != 0 ) {
124  if ( verbose_ ) std::cout << "# EE:channel, " << Numbers::indexEE(Numbers::iSM(id), ix, iy) << " hi=" << id.hashedIndex() << " " << it->getStatusCode() << std::endl;
125  std::vector<EcalDQMStatusDictionary::codeDef> codes;
126  EcalDQMStatusDictionary::getCodes( codes, it->getStatusCode() );
127  for ( unsigned int i=0; i<codes.size(); i++ ) {
128  std::cout << "Crystal " << Numbers::sEE(Numbers::iSM(id)) << " " << Numbers::indexEE(Numbers::iSM(id), ix, iy) << " " << codes[i].desc << std::endl;
129  }
130  }
131  }
132  }
133  }
134  if ( EEDetId::validDetId(ix, iy, -1) ) {
135  EEDetId id(ix, iy, -1);
136  if ( channelStatus ) {
137  EcalDQMChannelStatus::const_iterator it = channelStatus->find( id.rawId() );
138  if ( it != channelStatus->end() ) {
139  if ( it->getStatusCode() != 0 ) {
140  if ( verbose_ ) std::cout << "# EE:channel, " << Numbers::indexEE(Numbers::iSM(id), ix, iy) << " hi=" << id.hashedIndex() << " " << it->getStatusCode() << std::endl;
141  std::vector<EcalDQMStatusDictionary::codeDef> codes;
142  EcalDQMStatusDictionary::getCodes( codes, it->getStatusCode() );
143  for ( unsigned int i=0; i<codes.size(); i++ ) {
144  std::cout << "Crystal " << Numbers::sEE(Numbers::iSM(id)) << " " << Numbers::indexEE(Numbers::iSM(id), ix, iy) << " " << codes[i].desc << std::endl;
145  }
146  }
147  }
148  }
149  }
150  }
151  }
152 
153  for ( int ix=1; ix<=20; ix++ ) {
154  for ( int iy=1; iy<=20; iy++ ) {
155  if ( EcalScDetId::validDetId(ix, iy, +1) ) {
156  EcalScDetId id(ix, iy, +1);
157  if ( towerStatus ) {
158  EcalDQMTowerStatus::const_iterator it = towerStatus->find( id.rawId() );
159  if ( it != towerStatus->end() ) {
160  if ( it->getStatusCode() != 0 ) {
161  if ( verbose_ ) std::cout << "# EE:tower, " << Numbers::iSC(id) << " hi=" << id.hashedIndex() << " " << it->getStatusCode() << std::endl;
162  std::vector<EcalDQMStatusDictionary::codeDef> codes;
163  EcalDQMStatusDictionary::getCodes( codes, it->getStatusCode() );
164  for ( unsigned int i=0; i<codes.size(); i++ ) {
165  std::cout << "TT " << Numbers::sEE(Numbers::iSM(id)) << " " << Numbers::iSC(id) << " " << codes[i].desc << std::endl;
166  }
167  }
168  }
169  }
170  }
171  if ( EcalScDetId::validDetId(ix, iy, -1) ) {
172  EcalScDetId id(ix, iy, -1);
173  if ( towerStatus ) {
174  EcalDQMTowerStatus::const_iterator it = towerStatus->find( id.rawId() );
175  if ( it != towerStatus->end() ) {
176  if ( it->getStatusCode() != 0 ) {
177  if ( verbose_ ) std::cout << "# EE:tower, " << Numbers::iSC(id) << " hi=" << id.hashedIndex() << " " << it->getStatusCode() << std::endl;
178  std::vector<EcalDQMStatusDictionary::codeDef> codes;
179  EcalDQMStatusDictionary::getCodes( codes, it->getStatusCode() );
180  for ( unsigned int i=0; i<codes.size(); i++ ) {
181  std::cout << "TT " << Numbers::sEE(Numbers::iSM(id)) << " " << Numbers::iSC(id) << " " << codes[i].desc << std::endl;
182  }
183  }
184  }
185  }
186  }
187  }
188  }
189 
190 }
191 
EcalScDetId::validDetId
static bool validDetId(int ix, int iy, int iz)
Definition: EcalScDetId.cc:64
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
Numbers::sEE
static std::string sEE(const unsigned ism)
Definition: Numbers.cc:223
i
int i
Definition: DBlmapReader.cc:9
EcalDQMStatusDictionary::getCodes
static void getCodes(std::vector< codeDef > &codeVec, uint32_t code)
Definition: EcalDQMStatusDictionary.h:30
Numbers.h
Some &quot;id&quot; conversions.
Numbers::sEB
static std::string sEB(const unsigned ism)
Definition: Numbers.cc:91
EcalDQMStatusReader::beginRun
void beginRun(const edm::Run &r, const edm::EventSetup &c)
Definition: EcalDQMStatusReader.cc:36
EcalCondTowerObjectContainer::find
const_iterator find(uint32_t rawId) const
Definition: EcalCondTowerObjectContainer.h:57
edm::EventSetup::find
const eventsetup::EventSetupRecord * find(const eventsetup::EventSetupRecordKey &) const
Definition: EventSetup.cc:90
EcalDQMStatusDictionary::getDictionary
static void getDictionary(std::vector< codeDef > &dict)
Definition: EcalDQMStatusDictionary.h:25
Numbers::iSC
static unsigned iSC(const EcalScDetId &id)
Definition: Numbers.cc:434
patZpeak.handle
tuple handle
Definition: patZpeak.py:22
Numbers::initGeometry
static void initGeometry(const edm::EventSetup &setup, bool verbose=false)
Definition: Numbers.cc:47
EEDetId::validDetId
static bool validDetId(int crystal_ix, int crystal_iy, int iz)
Definition: EEDetId.h:248
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
EcalCondObjectContainer< EcalDQMStatusCode >::const_iterator
std::vector< Item >::const_iterator const_iterator
Definition: EcalCondObjectContainer.h:15
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
EcalTrigTowerDetId::validDetId
static bool validDetId(int iz, EcalSubdetector sd, int i, int j)
check if a valid index combination
Definition: EcalTrigTowerDetId.cc:97
Numbers::indexEE
static unsigned indexEE(const unsigned ism, const unsigned ix, const unsigned iy)
Definition: Numbers.cc:709
Numbers::iSM
static unsigned iSM(const unsigned ism, const EcalSubdetector subdet)
Definition: Numbers.cc:243
EcalCondObjectContainer::find
const_iterator find(uint32_t rawId) const
Definition: EcalCondObjectContainer.h:58
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
EcalCondObjectContainer::end
const_iterator end() const
Definition: EcalCondObjectContainer.h:88
EcalDQMStatusReader::EcalDQMStatusReader
EcalDQMStatusReader(const edm::ParameterSet &ps)
Definition: EcalDQMStatusReader.cc:30
edm::ESHandleBase::isValid
bool isValid() const
Definition: ESHandle.h:37
EBDetId::SMCRYSTALMODE
static const int SMCRYSTALMODE
Definition: EBDetId.h:167
ecalpyutils::ism
int ism(int ieta, int iphi)
Definition: EcalPyUtils.cc:56
Numbers::iTT
static unsigned iTT(const unsigned ism, const EcalSubdetector subdet, const unsigned i1, const unsigned i2)
Definition: Numbers.cc:482
edm::Run
Definition: Run.h:41
Generated for CMSSW Reference Manual by doxygen 1.8.5