CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EcalDQMStatusWriter.cc
Go to the documentation of this file.
1 /*
2  * \file EcalDQMStatusWriter.cc
3  *
4  * $Date: 2010/08/30 13:14:08 $
5  * $Revision: 1.19 $
6  * \author G. Della Ricca
7  *
8 */
9 
10 #include <fstream>
11 #include <iostream>
12 #include <string>
13 #include <cstring>
14 #include <time.h>
15 #include <unistd.h>
16 
17 #include "FWCore/ServiceRegistry/interface/Service.h"
18 #include "CondCore/DBOutputService/interface/PoolDBOutputService.h"
19 
20 #include "DQM/EcalCommon/interface/Numbers.h"
21 
22 #include "DataFormats/EcalDetId/interface/EBDetId.h"
23 #include "DataFormats/EcalDetId/interface/EEDetId.h"
24 #include "DataFormats/EcalDetId/interface/EcalTrigTowerDetId.h"
25 #include "DataFormats/EcalDetId/interface/EcalScDetId.h"
26 
27 #include "DQM/EcalCommon/interface/EcalDQMStatusDictionary.h"
28 
29 #include "DQM/EcalCommon/interface/EcalDQMStatusWriter.h"
30 
31 EcalDQMStatusWriter::EcalDQMStatusWriter(const edm::ParameterSet& ps) {
32 
33  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
34 
35  typedef std::vector<edm::ParameterSet> Parameters;
36  Parameters toPut = ps.getParameter<Parameters>("toPut");
37 
38  for ( Parameters::iterator itToPut=toPut.begin(); itToPut!=toPut.end(); itToPut++ ) {
39  inpFileName_.push_back(itToPut->getUntrackedParameter<std::string>("inputFile"));
40  objectName_.push_back(itToPut->getUntrackedParameter<std::string>("conditionType"));
41  since_.push_back(itToPut->getUntrackedParameter<unsigned int>("since"));
42  }
43 
44 }
45 
46 void EcalDQMStatusWriter::beginRun(const edm::Run& r, const edm::EventSetup& c) {
47 
48  Numbers::initGeometry(c, verbose_);
49 
50  edm::Service<cond::service::PoolDBOutputService> dbservice;
51  if ( !dbservice.isAvailable() ){
52  std::cout << "PoolDBOutputService is unavailable" << std::endl;
53  return;
54  }
55 
56  for ( unsigned int i=0; i<objectName_.size(); i++ ) {
57 
58  bool toAppend;
59  cond::Time_t newTime;
60 
61  if ( dbservice->isNewTagRequest( objectName_[i]+std::string("Rcd") ) ) {
62  // This is the first object for this tag.
63  // Append mode should be off.
64  // newTime is the end of this new objects IOV.
65  toAppend = false;
66  newTime = dbservice->beginOfTime();
67  } else {
68  // There should already be an object in the DB for this tag.
69  // Append IOV mode should be on.
70  // newTime is the beginning of this new objects IOV.
71  toAppend = true;
72  newTime = (cond::Time_t) since_[i];
73  }
74 
75  std::cout << "Reading " << objectName_[i] << " from file and writing to DB with newTime " << newTime << std::endl;
76 
77  if (objectName_[i] == "EcalDQMChannelStatus") {
78 
79  EcalDQMChannelStatus* status = readEcalDQMChannelStatusFromFile(inpFileName_[i].c_str());
80 
81  if ( !toAppend ) {
82  dbservice->createNewIOV<EcalDQMChannelStatus>(status, newTime, dbservice->endOfTime(), "EcalDQMChannelStatusRcd");
83  } else {
84  dbservice->appendSinceTime<EcalDQMChannelStatus>(status, newTime, "EcalDQMChannelStatusRcd");
85  }
86 
87  } else if (objectName_[i] == "EcalDQMTowerStatus") {
88 
89  EcalDQMTowerStatus* status = readEcalDQMTowerStatusFromFile(inpFileName_[i].c_str());
90 
91  if ( !toAppend ) {
92  dbservice->createNewIOV<EcalDQMTowerStatus>(status, newTime, dbservice->endOfTime(), "EcalDQMTowerStatusRcd");
93  } else {
94  dbservice->appendSinceTime<EcalDQMTowerStatus>(status, newTime, "EcalDQMTowerStatusRcd");
95  }
96 
97  } else {
98 
99  std::cout << "Object " << objectName_[i] << " is not supported by this program." << std::endl;
100 
101  }
102 
103  }
104 
105 }
106 
107 EcalDQMChannelStatus* EcalDQMStatusWriter::readEcalDQMChannelStatusFromFile(const char* inputFile) {
108 
109  EcalDQMChannelStatus* status = new EcalDQMChannelStatus();
110 
111  std::vector<EcalDQMStatusDictionary::codeDef> dictionary;
112  EcalDQMStatusDictionary::getDictionary( dictionary );
113 
114  // barrel
115  for ( int ism=1; ism<=36; ism++ ) {
116  int jsm = Numbers::iSM(ism, EcalBarrel);
117  for ( int ic=1; ic<=1700; ic++ ) {
118  EBDetId id(jsm, ic, EBDetId::SMCRYSTALMODE);
119  status->setValue(id, 0);
120  }
121  }
122 
123  // endcap
124  for ( int ix=1; ix<=100; ix++ ) {
125  for ( int iy=1; iy<=100; iy++ ) {
126  if ( EEDetId::validDetId(ix, iy, +1) ) {
127  EEDetId id(ix, iy, +1);
128  status->setValue(id, 0);
129  }
130  if ( EEDetId::validDetId(ix, iy, -1) ) {
131  EEDetId id(ix, iy, -1);
132  status->setValue(id, 0);
133  }
134  }
135  }
136 
137  std::cout << "Reading channel status from file " << inputFile << std::endl;
138  FILE *ifile = fopen( inputFile ,"r" );
139 
140  if ( !ifile ) throw cms::Exception ("Cannot open file") ;
141 
142  char line[256];
143 
144  int ii = 0;
145  while ( fgets(line, 255, ifile) ) {
146 
147  std::stringstream stream;
148 
149  ii++;
150  std::string key;
151  stream << line;
152  if ( verbose_ ) std::cout << line;
153  stream >> key;
154 
155  if ( key.size() == 0 || strcmp(key.c_str(), " ") == 0 || strncmp(key.c_str(), "#", 1) == 0 ) {
156 
157  // skip comments
158 
159  } else if ( strcmp(key.c_str(), "EB") == 0 ) {
160 
161  int index;
162  uint32_t code;
163  stream >> index >> code;
164 
165  EBDetId id = EBDetId::unhashIndex(index);
166  code = convert(code);
167 
168  EcalDQMChannelStatus::const_iterator it = status->find(id);
169  if ( it != status->end() ) code |= it->getStatusCode();
170  status->setValue(id, code);
171 
172  } else if ( strcmp(key.c_str(), "EE") == 0 ) {
173 
174  int ix, iy, iz;
175  uint32_t code;
176  stream >> ix >> iy >> iz >> code;
177 
178  EEDetId id(ix, iy, iz);
179  code = convert(code);
180 
181  EcalDQMChannelStatus::const_iterator it = status->find(id);
182  if ( it != status->end() ) code |= it->getStatusCode();
183  status->setValue(id, code);
184 
185  } else if ( strcmp(key.c_str(), "Crystal") == 0 ) {
186 
187  std::string module;
188  stream >> module;
189 
190  if ( strncmp(module.c_str(), "EB+", 3) == 0 || strncmp(module.c_str(), "EB-", 3) == 0 ) {
191 
192 #if 0
193  int ie, ip;
194  std::string token;
195  stream >> ie >> ip >> token;
196 
197  int sm = atoi( module.substr(2, module.size()-2).c_str() );
198 
199  int ism = (sm>=1&&sm<=18) ? sm+18 : -sm;
200  int jsm = Numbers::iSM(ism, EcalBarrel);
201  int ic = 20*(ie-1)+(ip-1)+1;
202 #else
203  int ic;
204  std::string token;
205  stream >> ic >> token;
206 
207  int sm = atoi( module.substr(2, module.size()-2).c_str() );
208 
209  int ism = (sm>=1&&sm<=18) ? sm+18 : -sm;
210  int jsm = Numbers::iSM(ism, EcalBarrel);
211 #endif
212 
213  EBDetId id(jsm, ic, EBDetId::SMCRYSTALMODE);
214  uint32_t code = 0;
215  for ( unsigned int i=0; i<dictionary.size(); i++ ) {
216  if ( strcmp(token.c_str(), dictionary[i].desc) == 0 ) {
217  code = dictionary[i].code;
218  }
219  }
220  if ( code == 0 ) {
221  std::cout << " --> not found in the dictionary: " << token << std::endl;
222  continue;
223  }
224 
225  EcalDQMChannelStatus::const_iterator it = status->find(id);
226  if ( it != status->end() ) code |= it->getStatusCode();
227  status->setValue(id, code);
228 
229  } else if ( strncmp(module.c_str(), "EE+", 3) == 0 || strncmp(module.c_str(), "EE-", 3) == 0 ) {
230 
231 #if 0
232  int jx, jy;
233  std::string token;
234  stream >> jx >> jy >> token;
235 #else
236  int index;
237  std::string token;
238  stream >> index >> token;
239  int jx = index/1000;
240  int jy = index%1000;
241 #endif
242 
243  int sm = atoi( module.substr(2, module.size()-2).c_str() );
244 
245  int ism = 0;
246  if( sm == -99 ) sm = -1;
247  if( sm == +99 ) sm = +1;
248  switch ( sm ) {
249  case -7: ism = 1; break;
250  case -8: ism = 2; break;
251  case -9: ism = 3; break;
252  case -1: ism = 4; break;
253  case -2: ism = 5; break;
254  case -3: ism = 6; break;
255  case -4: ism = 7; break;
256  case -5: ism = 8; break;
257  case -6: ism = 9; break;
258  case +7: ism = 10; break;
259  case +8: ism = 11; break;
260  case +9: ism = 12; break;
261  case +1: ism = 13; break;
262  case +2: ism = 14; break;
263  case +3: ism = 15; break;
264  case +4: ism = 16; break;
265  case +5: ism = 17; break;
266  case +6: ism = 18; break;
267  }
268 
269  EEDetId id(jx, jy, (ism>=1&&ism<=9)?-1:+1);
270  uint32_t code = 0;
271  for ( unsigned int i=0; i<dictionary.size(); i++ ) {
272  if ( strcmp(token.c_str(), dictionary[i].desc) == 0 ) {
273  code = dictionary[i].code;
274  }
275  }
276  if ( code == 0 ) {
277  std::cout << " --> not found in the dictionary: " << token << std::endl;
278  continue;
279  }
280 
281  EcalDQMChannelStatus::const_iterator it = status->find(id);
282  if ( it != status->end() ) code |= it->getStatusCode();
283  status->setValue(id, code);
284 
285  } else {
286 
287  std:: cout << "--> unknown token at line #" << ii << " : " << line;
288 
289  }
290 
291  } else if ( strcmp(key.c_str(), "TT") == 0 ) {
292 
293  // see readEcalDQMTowerStatusFromFile
294 
295  } else if ( strcmp(key.c_str(), "PN") == 0 || strcmp(key.c_str(), "MemCh") == 0 || strcmp(key.c_str(), "MemTT") == 0 ) {
296 
297  std::cout << "--> unsupported key at line #" << ii << " : " << line;
298 
299  } else {
300 
301  std:: cout << "--> skipped line #" << ii << " : " << line;
302 
303  }
304 
305  }
306 
307  fclose(ifile);
308 
309  return status;
310 
311 }
312 
313 EcalDQMTowerStatus* EcalDQMStatusWriter::readEcalDQMTowerStatusFromFile(const char* inputFile) {
314 
315  EcalDQMTowerStatus* status = new EcalDQMTowerStatus();
316 
317  std::vector<EcalDQMStatusDictionary::codeDef> dictionary;
318  EcalDQMStatusDictionary::getDictionary( dictionary );
319 
320  // barrel
321  for ( int ix=1; ix<=17; ix++ ) {
322  for ( int iy=1; iy<=72; iy++ ) {
323  if ( EcalTrigTowerDetId::validDetId(+1, EcalBarrel, ix, iy) ) {
324  EcalTrigTowerDetId id(+1, EcalBarrel, ix, iy);
325  status->setValue(id, 0);
326  }
327  if ( EcalTrigTowerDetId::validDetId(-1, EcalBarrel, ix, iy) ) {
328  EcalTrigTowerDetId id(-1, EcalBarrel, ix, iy);
329  status->setValue(id, 0);
330  }
331  }
332  }
333 
334  // endcap
335  for ( int ix=1; ix<=20; ix++ ) {
336  for ( int iy=1; iy<=20; iy++ ) {
337  if ( EcalScDetId::validDetId(ix, iy, +1) ) {
338  EcalScDetId id(ix, iy, +1);
339  status->setValue(id, 0);
340  }
341  if ( EcalScDetId::validDetId(ix, iy, -1) ) {
342  EcalScDetId id(ix, iy, -1);
343  status->setValue(id, 0);
344  }
345  }
346  }
347 
348  std::cout << "Reading tower status from file " << inputFile << std::endl;
349  FILE *ifile = fopen( inputFile ,"r" );
350 
351  if ( !ifile ) throw cms::Exception ("Cannot open file") ;
352 
353  char line[256];
354 
355  int ii = 0;
356  while ( fgets(line, 255, ifile) ) {
357 
358  std::stringstream stream;
359 
360  ii++;
361  std::string key;
362  stream << line;
363  if ( verbose_ ) std::cout << line;
364  stream >> key;
365 
366  if ( key.size() == 0 || strcmp(key.c_str(), " ") == 0 || strncmp(key.c_str(), "#", 1) == 0 ) {
367 
368  // skip comments
369 
370  } else if ( strcmp(key.c_str(), "EB") == 0 ) {
371 
372  // see readEcalDQMChannelStatusFromFile
373 
374  } else if ( strcmp(key.c_str(), "EE") == 0 ) {
375 
376  // see readEcalDQMChannelStatusFromFile
377 
378  } else if ( strcmp(key.c_str(), "Crystal") == 0 ) {
379 
380  // see readEcalDQMChannelStatusFromFile
381 
382  } else if ( strcmp(key.c_str(), "TT") == 0 ) {
383 
384  std::string module;
385  stream >> module;
386 
387  if ( strncmp(module.c_str(), "EB+", 3) == 0 || strncmp(module.c_str(), "EB-", 3) == 0 ) {
388 
389  int itt;
390  std::string token;
391  stream >> itt >> token;
392 
393  if ( itt >= 1 && itt <= 68 ) {
394 
395  int sm = atoi( module.substr(2, module.size()-2).c_str() );
396 
397  int iet = (itt-1)/4+1;
398  int ipt = (itt-1)%4+1;
399 
400  if ( sm<0 ) {
401  ipt = ipt+(std::abs(sm)-1)*4-2;
402  if ( ipt < 1 ) ipt = ipt+72;
403  if ( ipt > 72 ) ipt = ipt-72;
404  } else {
405  ipt = (5-ipt)+(std::abs(sm)-1)*4-2;
406  if ( ipt < 1 ) ipt = ipt+72;
407  if ( ipt > 72 ) ipt = ipt-72;
408  }
409 
410  EcalTrigTowerDetId id((sm<0)?-1:+1, EcalBarrel, iet, ipt);
411  uint32_t code = 0;
412  for ( unsigned int i=0; i<dictionary.size(); i++ ) {
413  if ( strcmp(token.c_str(), dictionary[i].desc) == 0 ) {
414  code = dictionary[i].code;
415  }
416  }
417  if ( code == 0 ) {
418  std::cout << " --> not found in the dictionary: " << token << std::endl;
419  continue;
420  }
421 
422  EcalDQMTowerStatus::const_iterator it = status->find(id);
423  if ( it != status->end() ) code |= it->getStatusCode();
424  status->setValue(id, code);
425 
426  } else {
427 
428  std::cout << "--> unsupported configuration at line #" << ii << " : " << line;
429 
430  }
431 
432  } else if ( strncmp(module.c_str(), "EE+", 3) == 0 || strncmp(module.c_str(), "EE-", 3) == 0 ) {
433 
434  int isc;
435  std::string token;
436  stream >> isc >> token;
437 
438  if ( isc >= 1 && isc <= 68 ) {
439 
440  int sm = atoi( module.substr(2, module.size()-2).c_str() );
441 
442  int ism = 0;
443  switch ( sm ) {
444  case -7: ism = 1; break;
445  case -8: ism = 2; break;
446  case -9: ism = 3; break;
447  case -1: ism = 4; break;
448  case -2: ism = 5; break;
449  case -3: ism = 6; break;
450  case -4: ism = 7; break;
451  case -5: ism = 8; break;
452  case -6: ism = 9; break;
453  case +7: ism = 10; break;
454  case +8: ism = 11; break;
455  case +9: ism = 12; break;
456  case +1: ism = 13; break;
457  case +2: ism = 14; break;
458  case +3: ism = 15; break;
459  case +4: ism = 16; break;
460  case +5: ism = 17; break;
461  case +6: ism = 18; break;
462  }
463 
464  int idcc = (ism>=1&&ism<=9) ? ism : ism-9+45;
465 
466  std::vector<DetId>* crystals = Numbers::crystals(idcc, isc);
467 
468  for ( unsigned int i=0; i<crystals->size(); i++ ) {
469 
470  EcalScDetId id = ((EEDetId) (*crystals)[i]).sc();
471  uint32_t code = 0;
472  for ( unsigned int i=0; i<dictionary.size(); i++ ) {
473  if ( strcmp(token.c_str(), dictionary[i].desc) == 0 ) {
474  code = dictionary[i].code;
475  }
476  }
477  if ( code == 0 ) {
478  std::cout << " --> not found in the dictionary: " << token << std::endl;
479  continue;
480  }
481 
482  EcalDQMTowerStatus::const_iterator it = status->find(id);
483  if ( it != status->end() ) code |= it->getStatusCode();
484  status->setValue(id, code);
485 
486  }
487 
488  } else {
489 
490  std::cout << "--> unsupported configuration at line #" << ii << " : " << line;
491 
492  }
493 
494  } else {
495 
496  std:: cout << "--> unknown token at line #" << ii << " : " << line;
497 
498  }
499 
500  } else if ( strcmp(key.c_str(), "PN") == 0 || strcmp(key.c_str(), "MemCh") == 0 || strcmp(key.c_str(), "MemTT") == 0 ) {
501 
502  std::cout << "--> unsupported key at line #" << ii << " : " << line;
503 
504  } else {
505 
506  std:: cout << "--> skipped line #" << ii << " : " << line;
507 
508  }
509 
510  }
511 
512  fclose(ifile);
513 
514  return status;
515 
516 }
517 
518 uint32_t EcalDQMStatusWriter::convert(uint32_t chStatus) {
519 
520  if ( chStatus == 1 || (chStatus >= 8 && chStatus <= 12 )) {
521 
522  chStatus = 0;
523  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR;
524  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_MEAN_ERROR;
525  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR;
526  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_RMS_ERROR;
527  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_RMS_ERROR;
528  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_RMS_ERROR;
529  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR;
530  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR;
531  chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_MEAN_ERROR;
532  chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_MEAN_ERROR;
533  chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR;
534  chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_RMS_ERROR;
535  chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_RMS_ERROR;
536  chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_RMS_ERROR;
537  chStatus |= 1 << EcalDQMStatusHelper::LASER_MEAN_ERROR;
538  chStatus |= 1 << EcalDQMStatusHelper::LASER_RMS_ERROR;
539  chStatus |= 1 << EcalDQMStatusHelper::LASER_TIMING_MEAN_ERROR;
540  chStatus |= 1 << EcalDQMStatusHelper::LASER_TIMING_RMS_ERROR;
541  chStatus |= 1 << EcalDQMStatusHelper::LED_MEAN_ERROR;
542  chStatus |= 1 << EcalDQMStatusHelper::LED_RMS_ERROR;
543  chStatus |= 1 << EcalDQMStatusHelper::LED_TIMING_MEAN_ERROR;
544  chStatus |= 1 << EcalDQMStatusHelper::LED_TIMING_RMS_ERROR;
545 
546  } else if ( chStatus == 2 ) {
547 
548  chStatus = 0;
549  chStatus |= 1 << EcalDQMStatusHelper::LASER_MEAN_ERROR;
550  chStatus |= 1 << EcalDQMStatusHelper::LASER_RMS_ERROR;
551  chStatus |= 1 << EcalDQMStatusHelper::LASER_TIMING_MEAN_ERROR;
552  chStatus |= 1 << EcalDQMStatusHelper::LASER_TIMING_RMS_ERROR;
553  chStatus |= 1 << EcalDQMStatusHelper::LED_MEAN_ERROR;
554  chStatus |= 1 << EcalDQMStatusHelper::LED_RMS_ERROR;
555  chStatus |= 1 << EcalDQMStatusHelper::LED_TIMING_MEAN_ERROR;
556  chStatus |= 1 << EcalDQMStatusHelper::LED_TIMING_RMS_ERROR;
557 
558  } else if ( chStatus == 3 ) {
559 
560  chStatus = 0;
561  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR;
562  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_MEAN_ERROR;
563  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR;
564  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_RMS_ERROR;
565  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_RMS_ERROR;
566  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_RMS_ERROR;
567  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR;
568  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR;
569  chStatus |= 1 << EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_WARNING;
570 
571  } else if ( chStatus == 4 ) {
572 
573  chStatus = 0;
574  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR;
575  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_MEAN_ERROR;
576  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR;
577  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_RMS_ERROR;
578  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_RMS_ERROR;
579  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_RMS_ERROR;
580  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR;
581  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR;
582 
583  } else if ( chStatus == 13 || chStatus == 14 ) {
584 
585  chStatus = 0;
586  chStatus |= 1 << EcalDQMStatusHelper::STATUS_FLAG_ERROR;
587 
588  } else if ( (chStatus&0x3f) == 13 || (chStatus&0x7f) == 13 ||
589  (chStatus&0x3f) == 14 || (chStatus&0x7f) == 14 ) {
590 
591  chStatus = 0;
592  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR;
593  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_MEAN_ERROR;
594  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR;
595  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_RMS_ERROR;
596  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_RMS_ERROR;
597  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_RMS_ERROR;
598  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR;
599  chStatus |= 1 << EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR;
600  chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_MEAN_ERROR;
601  chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_MEAN_ERROR;
602  chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR;
603  chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_RMS_ERROR;
604  chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_RMS_ERROR;
605  chStatus |= 1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_RMS_ERROR;
606  chStatus |= 1 << EcalDQMStatusHelper::LASER_MEAN_ERROR;
607  chStatus |= 1 << EcalDQMStatusHelper::LASER_RMS_ERROR;
608  chStatus |= 1 << EcalDQMStatusHelper::LASER_TIMING_MEAN_ERROR;
609  chStatus |= 1 << EcalDQMStatusHelper::LASER_TIMING_RMS_ERROR;
610  chStatus |= 1 << EcalDQMStatusHelper::LED_MEAN_ERROR;
611  chStatus |= 1 << EcalDQMStatusHelper::LED_RMS_ERROR;
612  chStatus |= 1 << EcalDQMStatusHelper::LED_TIMING_MEAN_ERROR;
613  chStatus |= 1 << EcalDQMStatusHelper::LED_TIMING_RMS_ERROR;
614  chStatus |= 1 << EcalDQMStatusHelper::STATUS_FLAG_ERROR;
615 
616  }
617 
618  return( chStatus );
619 
620 }
621 
EcalScDetId::validDetId
static bool validDetId(int ix, int iy, int iz)
Definition: EcalScDetId.cc:63
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
EcalDQMStatusHelper::LED_MEAN_ERROR
static const int LED_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:40
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_RMS_ERROR
static const int PEDESTAL_MIDDLE_GAIN_RMS_ERROR
Definition: EcalDQMStatusHelper.h:22
EcalDQMStatusHelper::LASER_MEAN_ERROR
static const int LASER_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:35
EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_RMS_ERROR
static const int TESTPULSE_LOW_GAIN_RMS_ERROR
Definition: EcalDQMStatusHelper.h:31
EcalDQMStatusHelper::PHYSICS_BAD_CHANNEL_WARNING
static const int PHYSICS_BAD_CHANNEL_WARNING
Definition: EcalDQMStatusHelper.h:47
EcalDQMStatusWriter::since_
std::vector< unsigned long long > since_
Definition: EcalDQMStatusWriter.h:50
Numbers.h
Some &quot;id&quot; conversions.
EcalDQMStatusHelper::LASER_TIMING_RMS_ERROR
static const int LASER_TIMING_RMS_ERROR
Definition: EcalDQMStatusHelper.h:38
EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR
static const int PEDESTAL_HIGH_GAIN_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:20
EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR
static const int PEDESTAL_ONLINE_HIGH_GAIN_RMS_ERROR
Definition: EcalDQMStatusHelper.h:26
abs
#define abs(x)
Definition: mlp_lapack.h:159
EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_MEAN_ERROR
static const int TESTPULSE_MIDDLE_GAIN_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:29
EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_RMS_ERROR
static const int PEDESTAL_LOW_GAIN_RMS_ERROR
Definition: EcalDQMStatusHelper.h:21
EcalDQMStatusWriter::readEcalDQMTowerStatusFromFile
EcalDQMTowerStatus * readEcalDQMTowerStatusFromFile(const char *)
Definition: EcalDQMStatusWriter.cc:313
cuy.ii
int ii
Definition: cuy.py:588
EcalCondTowerObjectContainer::setValue
void setValue(const uint32_t id, const Item &item)
Definition: EcalCondTowerObjectContainer.h:84
EcalCondTowerObjectContainer::find
const_iterator find(uint32_t rawId) const
Definition: EcalCondTowerObjectContainer.h:57
EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_RMS_ERROR
static const int PEDESTAL_HIGH_GAIN_RMS_ERROR
Definition: EcalDQMStatusHelper.h:23
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
EcalDQMStatusHelper::STATUS_FLAG_ERROR
static const int STATUS_FLAG_ERROR
Definition: EcalDQMStatusHelper.h:45
EcalCondObjectContainer::setValue
void setValue(const uint32_t id, const Item &item)
Definition: EcalCondObjectContainer.h:93
EcalDQMStatusHelper::TESTPULSE_MIDDLE_GAIN_RMS_ERROR
static const int TESTPULSE_MIDDLE_GAIN_RMS_ERROR
Definition: EcalDQMStatusHelper.h:32
Numbers::crystals
static std::vector< DetId > * crystals(const EcalTrigTowerDetId &id)
Definition: Numbers.cc:585
cond::service::PoolDBOutputService::appendSinceTime
void appendSinceTime(T *payloadObj, cond::Time_t sinceTime, const std::string &recordName, bool withlogging=false)
Definition: PoolDBOutputService.h:155
EcalDQMChannelStatus
EcalDQMChannelStatusMap EcalDQMChannelStatus
Definition: EcalDQMChannelStatus.h:10
cond::Time_t
unsigned long long Time_t
Definition: Time.h:16
Parameters
vector< ParameterSet > Parameters
Definition: HLTMuonPlotter.cc:34
EcalDQMStatusDictionary::getDictionary
static void getDictionary(std::vector< codeDef > &dict)
Definition: EcalDQMStatusDictionary.h:27
cond::service::PoolDBOutputService::isNewTagRequest
bool isNewTagRequest(const std::string &recordName)
Definition: PoolDBOutputService.cc:106
EcalDQMStatusWriter::convert
uint32_t convert(uint32_t c)
Definition: EcalDQMStatusWriter.cc:518
EcalDQMTowerStatus
EcalDQMTowerStatusMap EcalDQMTowerStatus
Definition: EcalDQMTowerStatus.h:13
EcalDQMStatusWriter::readEcalDQMChannelStatusFromFile
EcalDQMChannelStatus * readEcalDQMChannelStatusFromFile(const char *)
Definition: EcalDQMStatusWriter.cc:107
edm::Service::isAvailable
bool isAvailable() const
Definition: Service.h:47
compare_using_db.ifile
tuple ifile
Definition: compare_using_db.py:253
EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_RMS_ERROR
static const int TESTPULSE_HIGH_GAIN_RMS_ERROR
Definition: EcalDQMStatusHelper.h:33
EcalDQMStatusHelper::LED_TIMING_MEAN_ERROR
static const int LED_TIMING_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:42
EcalDQMStatusWriter::beginRun
void beginRun(const edm::Run &r, const edm::EventSetup &c)
Definition: EcalDQMStatusWriter.cc:46
cond::service::PoolDBOutputService::createNewIOV
void createNewIOV(T *firstPayloadObj, cond::Time_t firstSinceTime, cond::Time_t firstTillTime, const std::string &recordName, bool withlogging=false)
Definition: PoolDBOutputService.h:129
analyzePatCOC_cfg.inputFile
tuple inputFile
Definition: analyzePatCOC_cfg.py:6
EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR
static const int PEDESTAL_LOW_GAIN_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:18
EcalDQMStatusWriter::inpFileName_
std::vector< std::string > inpFileName_
Definition: EcalDQMStatusWriter.h:49
EcalDQMStatusHelper::TESTPULSE_LOW_GAIN_MEAN_ERROR
static const int TESTPULSE_LOW_GAIN_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:28
Numbers::initGeometry
static void initGeometry(const edm::EventSetup &setup, bool verbose=false)
Definition: Numbers.cc:50
EcalDQMStatusWriter::objectName_
std::vector< std::string > objectName_
Definition: EcalDQMStatusWriter.h:48
convertXMLtoSQLite_cfg.toPut
tuple toPut
Definition: convertXMLtoSQLite_cfg.py:21
EEDetId::validDetId
static bool validDetId(int crystal_ix, int crystal_iy, int iz)
Definition: EEDetId.h:249
EcalDQMStatusHelper::LED_TIMING_RMS_ERROR
static const int LED_TIMING_RMS_ERROR
Definition: EcalDQMStatusHelper.h:43
EcalCondObjectContainer::const_iterator
std::vector< Item >::const_iterator const_iterator
Definition: EcalCondObjectContainer.h:15
EcalDQMStatusHelper::PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR
static const int PEDESTAL_ONLINE_HIGH_GAIN_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:25
EcalDQMStatusHelper::LASER_TIMING_MEAN_ERROR
static const int LASER_TIMING_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:37
EBDetId::unhashIndex
static EBDetId unhashIndex(int hi)
get a DetId from a compact index for arrays
Definition: EBDetId.h:115
EcalTrigTowerDetId::validDetId
static bool validDetId(int iz, EcalSubdetector sd, int i, int j)
check if a valid index combination
Definition: EcalTrigTowerDetId.cc:97
combine.key
list key
Definition: combine.py:13
Numbers::iSM
static unsigned iSM(const unsigned ism, const EcalSubdetector subdet)
Definition: Numbers.cc:246
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
EcalDQMStatusHelper::LED_RMS_ERROR
static const int LED_RMS_ERROR
Definition: EcalDQMStatusHelper.h:41
EcalDQMStatusWriter::EcalDQMStatusWriter
EcalDQMStatusWriter(const edm::ParameterSet &ps)
Definition: EcalDQMStatusWriter.cc:31
EcalCondTowerObjectContainer::const_iterator
std::vector< Item >::const_iterator const_iterator
Definition: EcalCondTowerObjectContainer.h:18
EcalDQMStatusHelper::PEDESTAL_MIDDLE_GAIN_MEAN_ERROR
static const int PEDESTAL_MIDDLE_GAIN_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:19
EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR
static const int TESTPULSE_HIGH_GAIN_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:30
ntuplemaker.status
tuple status
Definition: ntuplemaker.py:245
EcalDQMStatusHelper::LASER_RMS_ERROR
static const int LASER_RMS_ERROR
Definition: EcalDQMStatusHelper.h:36
EBDetId::SMCRYSTALMODE
static const int SMCRYSTALMODE
Definition: EBDetId.h:168
ecalpyutils::ism
int ism(int ieta, int iphi)
Definition: EcalPyUtils.cc:56
edm::Run
Definition: Run.h:36
Generated for CMSSW Reference Manual by doxygen 1.8.5