CMS 3D CMS Logo
MESetBinningUtils.cc
Go to the documentation of this file.
1 #include "DQM/EcalCommon/interface/MESetBinningUtils.h"
2 
3 #include "DQM/EcalCommon/interface/EcalDQMCommonUtils.h"
4 
5 #include "FWCore/MessageLogger/interface/MessageLogger.h"
6 #include "FWCore/ParameterSet/interface/ParameterSet.h"
7 #include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
8 
9 #include "TPRegexp.h"
10 
11 namespace ecaldqm {
12  namespace binning {
13  AxisSpecs getBinning(ObjectType _otype, BinningType _btype, bool _isMap, int _axis, unsigned _iME) {
14  if (_otype >= nObjType || _btype >= unsigned(nPresetBinnings))
15  return AxisSpecs(); // you are on your own
16 
17  switch (_otype) {
18  case kEB:
19  return getBinningEB_(_btype, _isMap, _axis);
20  case kEE:
21  return getBinningEE_(_btype, _isMap, 0, _axis);
22  case kEEm:
23  return getBinningEE_(_btype, _isMap, -1, _axis);
24  case kEEp:
25  return getBinningEE_(_btype, _isMap, 1, _axis);
26  case kSM:
27  return getBinningSM_(_btype, _isMap, _iME, _axis);
28  case kEBSM:
29  return getBinningSM_(_btype, _isMap, _iME + 9, _axis);
30  case kEESM:
31  if (_iME <= kEEmHigh)
32  return getBinningSM_(_btype, _isMap, _iME, _axis);
33  else
34  return getBinningSM_(_btype, _isMap, _iME + nEBDCC, _axis);
35  case kSMMEM:
36  return getBinningSMMEM_(_btype, _isMap, _iME, _axis);
37  case kEBSMMEM:
38  return getBinningSMMEM_(_btype, _isMap, _iME + nEEDCCMEM / 2, _axis);
39  case kEESMMEM:
40  if (_iME <= kEEmHigh)
41  return getBinningSMMEM_(_btype, _isMap, _iME, _axis);
42  else
43  return getBinningSMMEM_(_btype, _isMap, _iME + nEBDCC, _axis);
44  case kEcal:
45  return getBinningEcal_(_btype, _isMap, _axis);
46  case kMEM:
47  return getBinningMEM_(_btype, _isMap, -1, _axis);
48  case kEBMEM:
49  return getBinningMEM_(_btype, _isMap, EcalBarrel, _axis);
50  case kEEMEM:
51  return getBinningMEM_(_btype, _isMap, EcalEndcap, _axis);
52  default:
53  return AxisSpecs();
54  }
55  }
56 
57  int findBin1D(ObjectType _otype, BinningType _btype, const DetId &_id) {
58  switch (_otype) {
59  case kSM:
60  case kEBSM:
61  case kEESM:
62  if (_btype == kSuperCrystal)
63  return towerId(_id);
64  else if (_btype == kTriggerTower) {
65  unsigned tccid(tccId(_id));
66  if (tccid <= 36 || tccid >= 73) { // EE
67  unsigned bin(ttId(_id));
68  bool outer((tccid >= 19 && tccid <= 36) || (tccid >= 73 && tccid <= 90));
69  // For the following, the constants nTTInner and nTTOuter are defined in
70  // EcalDQMCommonUtils.h.
71  if (outer)
72  bin += 2 * nTTInner; // For outer TCCs, sets bin number to increment
73  // by twice the number of TTs in inner TCCs, because numbering of bins
74  // is in the order (inner1, inner2, outer1, outer2).
75  // ("inner"" := closer to the beam)
76  bin += (tccid % 2) * (outer ? nTTOuter : nTTInner); // Yields x-axis bin number
77  // in the format above; TTs in even-numbered TCCs are filled in inner1
78  // or outer1, and those in odd-numbered TCC are filled in inner2 and
79  // outer2.
80  return bin;
81  } else
82  return ttId(_id);
83  } else
84  break;
85  case kEcal:
86  if (_btype == kDCC)
87  return dccId(_id);
88  else if (_btype == kTCC)
89  return tccId(_id);
90  else
91  break;
92  case kEB:
93  if (_btype == kDCC)
94  return dccId(_id) - 9;
95  else if (_btype == kTCC)
96  return tccId(_id) - 36;
97  else
98  break;
99  case kEEm:
100  if (_btype == kDCC)
101  return dccId(_id);
102  else if (_btype == kTCC)
103  return tccId(_id);
104  else
105  break;
106  case kEEp:
107  if (_btype == kDCC)
108  return dccId(_id) - 45;
109  else if (_btype == kTCC)
110  return tccId(_id) - 72;
111  else
112  break;
113  case kEE:
114  if (_btype == kDCC) {
115  int bin(dccId(_id));
116  if (bin >= 46)
117  bin -= 36;
118  return bin;
119  } else if (_btype == kTCC) {
120  int bin(tccId(_id));
121  if (bin >= 72)
122  bin -= 36;
123  return bin;
124  } else
125  break;
126  case kSMMEM:
127  case kEBSMMEM:
128  case kEESMMEM:
129  if (_btype == kCrystal)
130  return EcalPnDiodeDetId(_id).iPnId();
131  else
132  break;
133  default:
134  break;
135  }
136 
137  return 0;
138  }
139 
140  int findBin1D(ObjectType _otype, BinningType _btype, const EcalElectronicsId &_id) {
141  switch (_otype) {
142  case kSM:
143  case kEBSM:
144  case kEESM:
145  if (_btype == kSuperCrystal)
146  return towerId(_id);
147  else if (_btype == kTriggerTower) {
148  unsigned tccid(tccId(_id));
149  if (tccid <= 36 || tccid >= 73) { // EE
150  unsigned bin(ttId(_id));
151  bool outer((tccid >= 19 && tccid <= 36) || (tccid >= 73 && tccid <= 90));
152  // For the following, the constants nTTInner and nTTOuter are defined in
153  // EcalDQMCommonUtils.h.
154  if (outer)
155  bin += 2 * nTTInner; // For outer TCCs, sets bin number to increment
156  // by twice the number of TTs in inner TCCs, because numbering of bins
157  // is in the order (inner1, inner2, outer1, outer2).
158  // ("inner"" := closer to the beam)
159  bin += (tccid % 2) * (outer ? nTTOuter : nTTInner); // Yields x-axis bin number
160  // in the format above; TTs in even-numbered TCCs are filled in inner1
161  // or outer1, and those in odd-numbered TCC are filled in inner2 and
162  // outer2.
163  return bin;
164  } else
165  return ttId(_id);
166  } else
167  break;
168  case kEcal:
169  if (_btype == kDCC)
170  return dccId(_id);
171  else if (_btype == kTCC)
172  return tccId(_id);
173  else
174  break;
175  case kEB:
176  if (_btype == kDCC)
177  return dccId(_id) - 9;
178  else if (_btype == kTCC)
179  return tccId(_id) - 36;
180  else
181  break;
182  case kEEm:
183  if (_btype == kDCC)
184  return dccId(_id);
185  else if (_btype == kTCC)
186  return tccId(_id);
187  else
188  break;
189  case kEEp:
190  if (_btype == kDCC)
191  return dccId(_id) - 45;
192  else if (_btype == kTCC)
193  return tccId(_id) - 72;
194  else
195  break;
196  case kEE:
197  if (_btype == kDCC) {
198  int bin(dccId(_id));
199  if (bin >= 46)
200  bin -= 36;
201  return bin;
202  } else if (_btype == kTCC) {
203  int bin(tccId(_id));
204  if (bin >= 72)
205  bin -= 36;
206  return bin;
207  } else
208  break;
209  default:
210  break;
211  }
212 
213  return 0;
214  }
215 
216  int findBin1D(ObjectType _otype, BinningType _btype, int _dcctccid) {
217  if (_otype == kEcal && _btype == kDCC)
218  return _dcctccid;
219  else if (_otype == kEcal && _btype == kTCC)
220  return _dcctccid;
221  if (_otype == kEB && _btype == kDCC)
222  return _dcctccid - 9;
223  else if (_otype == kEB && _btype == kTCC)
224  return _dcctccid - 36;
225  else if (_otype == kEEm && _btype == kDCC)
226  return _dcctccid;
227  else if (_otype == kEEm && _btype == kTCC)
228  return _dcctccid;
229  else if (_otype == kEEp && _btype == kDCC)
230  return _dcctccid - 45;
231  else if (_otype == kEEp && _btype == kTCC)
232  return _dcctccid - 72;
233  else if (_otype == kEE && _btype == kDCC)
234  return _dcctccid <= 9 ? _dcctccid : _dcctccid - 36;
235  else if (_otype == kEE && _btype == kTCC)
236  return _dcctccid <= 36 ? _dcctccid : _dcctccid - 36;
237 
238  return 0;
239  }
240 
241  int findBin2D(ObjectType _otype, BinningType _btype, const DetId &_id) {
242  if (_otype >= nObjType || _btype >= unsigned(nPresetBinnings))
243  return 0;
244 
245  switch (_btype) {
246  case kCrystal:
247  return findBinCrystal_(_otype, _id);
248  break;
249  case kTriggerTower:
250  return findBinTriggerTower_(_otype, _id);
251  break;
252  case kSuperCrystal:
253  return findBinSuperCrystal_(_otype, _id);
254  break;
255  case kPseudoStrip:
256  return findBinPseudoStrip_(_otype, _id);
257  break;
258  case kRCT:
259  return findBinRCT_(_otype, _id);
260  break;
261  default:
262  return 0;
263  }
264  }
265 
266  int findBin2D(ObjectType _otype, BinningType _btype, const EcalElectronicsId &_id) {
267  if (_otype >= nObjType || _btype >= unsigned(nPresetBinnings))
268  return 0;
269 
270  switch (_btype) {
271  case kCrystal:
272  return findBinCrystal_(_otype, _id);
273  break;
274  case kSuperCrystal:
275  return findBinSuperCrystal_(_otype, _id);
276  break;
277  default:
278  return 0;
279  }
280  }
281 
282  int findBin2D(ObjectType _otype, BinningType _btype, int _dccid) {
283  if (_otype != kEcal || _btype != kDCC)
284  return 0;
285 
286  int nbinsX(9);
287  unsigned iDCC(_dccid - 1);
288  int xbin(0);
289  if (iDCC <= kEEmHigh || iDCC >= kEEpLow)
290  xbin = (iDCC + 6) % nbinsX + 1;
291  else
292  xbin = iDCC % nbinsX + 1;
293  int ybin(6 - iDCC / nbinsX);
294 
295  return (nbinsX + 2) * ybin + xbin;
296  }
297 
298  unsigned findPlotIndex(ObjectType _otype, const DetId &_id) {
299  if (getNObjects(_otype) == 1)
300  return 0;
301 
302  switch (_otype) {
303  case kEcal3P:
304  if (_id.subdetId() == EcalBarrel)
305  return 1;
306  else if (_id.subdetId() == EcalEndcap && zside(_id) > 0)
307  return 2;
308  else if (_id.subdetId() == EcalTriggerTower) {
309  if (!isEndcapTTId(_id))
310  return 1;
311  else {
312  if (zside(_id) > 0)
313  return 2;
314  else
315  return 0;
316  }
317  } else
318  return 0;
319 
320  case kEcal2P:
321  if (_id.subdetId() == EcalBarrel)
322  return 1;
323  else if (_id.subdetId() == EcalTriggerTower && !isEndcapTTId(_id))
324  return 1;
325  else
326  return 0;
327 
328  case kEE2P:
329  if (zside(_id) > 0)
330  return 1;
331  else
332  return 0;
333 
334  case kMEM2P:
335  if (_id.subdetId() == EcalLaserPnDiode) {
336  unsigned iDCC(dccId(_id) - 1);
337  if (iDCC >= kEBmLow && iDCC <= kEBpHigh)
338  return 1;
339  else
340  return 0;
341  } else
342  return -1;
343 
344  default:
345  return findPlotIndex(_otype, dccId(_id));
346  }
347  }
348 
349  unsigned findPlotIndex(ObjectType _otype, const EcalElectronicsId &_id) {
350  if (getNObjects(_otype) == 1)
351  return 0;
352 
353  return findPlotIndex(_otype, _id.dccId());
354  }
355 
356  unsigned findPlotIndex(ObjectType _otype, int _dcctccid, BinningType _btype /* = kDCC*/) {
357  if (getNObjects(_otype) == 1)
358  return 0;
359 
360  int iSM(_dcctccid - 1);
361 
362  switch (_otype) {
363  case kSM:
364  if (_btype == kPseudoStrip) {
365  iSM = iSM <= kEEmTCCHigh ? (iSM + 1) % 18 / 2
366  : iSM >= kEEpTCCLow ? (iSM + 1 - 72) % 18 / 2 + 45 : (iSM + 1) - kEEmTCCHigh;
367  return iSM;
368  } else
369  return iSM;
370 
371  case kEBSM:
372  return iSM - 9;
373 
374  case kEESM:
375  if (iSM <= kEEmHigh)
376  return iSM;
377  else
378  return iSM - nEBDCC;
379 
380  case kSMMEM:
381  return memDCCIndex(_dcctccid);
382 
383  case kEBSMMEM:
384  return memDCCIndex(_dcctccid) - nEEDCCMEM / 2;
385 
386  case kEESMMEM:
387  if (iSM <= kEEmHigh)
388  return memDCCIndex(_dcctccid);
389  else
390  return memDCCIndex(_dcctccid) - nEBDCC;
391 
392  case kEcal2P:
393  if (_btype == kDCC) {
394  if (iSM <= kEEmHigh || iSM >= kEEpLow)
395  return 0;
396  else
397  return 1;
398  } else if (_btype == kTCC) {
399  if (iSM <= kEEmTCCHigh || iSM >= kEEpTCCLow)
400  return 0;
401  else
402  return 1;
403  } else {
404  if (iSM == EcalBarrel - 1)
405  return 1;
406  else
407  return 0;
408  }
409 
410  case kEcal3P:
411  if (_btype == kDCC) {
412  if (iSM <= kEEmHigh)
413  return 0;
414  else if (iSM <= kEBpHigh)
415  return 1;
416  else
417  return 2;
418  } else if (_btype == kTCC) {
419  if (iSM <= kEEmTCCHigh)
420  return 0;
421  else if (iSM <= kEBTCCHigh)
422  return 1;
423  else
424  return 2;
425  } else {
426  if (iSM == -EcalEndcap - 1)
427  return 0;
428  else if (iSM == EcalBarrel - 1)
429  return 1;
430  else
431  return 2;
432  }
433 
434  case kEE2P:
435  if (_btype == kDCC) {
436  if (iSM >= kEEpLow)
437  return 1;
438  else
439  return 0;
440  } else {
441  if (iSM >= kEEpTCCLow)
442  return 1;
443  else
444  return 0;
445  }
446 
447  case kMEM2P:
448  if (_btype == kDCC) {
449  if (iSM <= kEEmHigh || iSM >= kEEpLow)
450  return 0;
451  else
452  return 1;
453  } else if (_btype == kTCC)
454  return -1;
455  else {
456  if (iSM == kEB)
457  return 1;
458  else
459  return 0;
460  }
461  default:
462  return -1;
463  }
464  }
465 
466  ObjectType getObject(ObjectType _otype, unsigned _iObj) {
467  switch (_otype) {
468  case kEcal3P:
469  switch (_iObj) {
470  case 0:
471  return kEEm;
472  case 1:
473  return kEB;
474  case 2:
475  return kEEp;
476  default:
477  return nObjType;
478  }
479  case kEcal2P:
480  switch (_iObj) {
481  case 0:
482  return kEE;
483  case 1:
484  return kEB;
485  default:
486  return nObjType;
487  }
488  case kEE2P:
489  switch (_iObj) {
490  case 0:
491  return kEEm;
492  case 1:
493  return kEEp;
494  default:
495  return nObjType;
496  }
497  case kMEM2P:
498  switch (_iObj) {
499  case 0:
500  return kEEMEM;
501  case 1:
502  return kEBMEM;
503  default:
504  return nObjType;
505  }
506  default:
507  return _otype;
508  }
509  }
510 
511  unsigned getNObjects(ObjectType _otype) {
512  switch (_otype) {
513  case kSM:
514  return nDCC;
515  case kEBSM:
516  return nEBDCC;
517  case kEESM:
518  return nEEDCC;
519  case kSMMEM:
520  return nDCCMEM;
521  case kEBSMMEM:
522  return nEBDCC;
523  case kEESMMEM:
524  return nEEDCCMEM;
525  case kEcal2P:
526  return 2;
527  case kEcal3P:
528  return 3;
529  case kEE2P:
530  return 2;
531  case kMEM2P:
532  return 2;
533  default:
534  return 1;
535  }
536  }
537 
538  bool isValidIdBin(ObjectType _otype, BinningType _btype, unsigned _iME, int _bin) {
539  if (_otype == kEEm || _otype == kEEp) {
540  if (_btype == kCrystal || _btype == kTriggerTower)
541  return EEDetId::validDetId(_bin % 102, _bin / 102, 1);
542  else if (_btype == kSuperCrystal)
543  return EcalScDetId::validDetId(_bin % 22, _bin / 22, 1);
544  } else if (_otype == kEE) {
545  if (_btype == kCrystal || _btype == kTriggerTower) {
546  int ix(_bin % 202);
547  if (ix > 100)
548  ix = (ix - 100) % 101;
549  return EEDetId::validDetId(ix, _bin / 202, 1);
550  } else if (_btype == kSuperCrystal) {
551  int ix(_bin % 42);
552  if (ix > 20)
553  ix = (ix - 20) % 21;
554  return EcalScDetId::validDetId(ix, _bin / 42, 1);
555  }
556  } else if (_otype == kSM || _otype == kEESM) {
557  unsigned iSM(_iME);
558  if (_otype == kEESM && iSM > kEEmHigh)
559  iSM += nEBDCC;
560 
561  if (iSM >= kEBmLow && iSM <= kEBpHigh)
562  return true;
563 
564  if (_btype == kCrystal || _btype == kTriggerTower) {
565  int nX(nEESMX);
566  if (iSM == kEEm02 || iSM == kEEm08 || iSM == kEEp02 || iSM == kEEp08)
567  nX = nEESMXExt;
568  if (iSM == kEEm01 || iSM == kEEm05 || iSM == kEEm09 || iSM == kEEp01 || iSM == kEEp05 || iSM == kEEp09)
569  nX = nEESMXRed;
570  int ix(_bin % (nX + 2) + xlow_(iSM));
571  int iy(_bin / (nX + 2) + ylow_(iSM));
572  int z(iSM <= kEEmHigh ? -1 : 1);
573  return EEDetId::validDetId(ix, iy, 1) && iSM == dccId(EEDetId(ix, iy, z)) - 1;
574  } else if (_btype == kSuperCrystal) {
575  int nX(nEESMX / 5);
576  if (iSM == kEEm02 || iSM == kEEm08 || iSM == kEEp02 || iSM == kEEp08)
577  nX = nEESMXExt / 5;
578  if (iSM == kEEm01 || iSM == kEEm05 || iSM == kEEm09 || iSM == kEEp01 || iSM == kEEp05 || iSM == kEEp09)
579  nX = nEESMXRed / 5;
580  int ix(_bin % (nX + 2) + xlow_(iSM) / 5);
581  int iy(_bin / (nX + 2) + ylow_(iSM) / 5);
582  int z(iSM <= kEEmHigh ? -1 : 1);
583  return EcalScDetId::validDetId(ix, iy, z) && iSM == dccId(EcalScDetId(ix, iy, z)) - 1;
584  }
585  }
586 
587  return true;
588  }
589 
590  std::string channelName(uint32_t _rawId, BinningType _btype /* = kDCC*/) {
591  // assume the following IDs for respective binning types:
592  // Crystal: EcalElectronicsId
593  // TriggerTower: EcalTriggerElectronicsId (pstrip and channel ignored)
594  // SuperCrystal: EcalElectronicsId (strip and crystal ignored)
595  // TCC: TCC ID
596  // DCC: DCC ID
597 
598  std::stringstream ss;
599 
600  switch (_btype) {
601  case kCrystal: {
602  // EB-03 DCC 12 CCU 12 strip 3 xtal 1 (EB ieta -13 iphi 60) (TCC 39 TT 12
603  // pstrip 3 chan 1)
604  EcalElectronicsId eid(_rawId);
605  if (eid.towerId() >= 69)
606  ss << smName(eid.dccId()) << " DCC " << eid.dccId() << " CCU " << eid.towerId() << " PN " << eid.xtalId();
607  else {
608  ss << smName(eid.dccId()) << " DCC " << eid.dccId() << " CCU " << eid.towerId() << " strip "
609  << eid.stripId() << " xtal " << eid.xtalId();
610 
611  if (eid.dccId() >= kEBmLow + 1 && eid.dccId() <= kEBpHigh + 1) {
612  EBDetId ebid(getElectronicsMap()->getDetId(eid));
613  ss << " (EB ieta " << std::showpos << ebid.ieta() << std::noshowpos << " iphi " << ebid.iphi() << ")";
614  } else {
615  EEDetId eeid(getElectronicsMap()->getDetId(eid));
616  ss << " (EE ix " << eeid.ix() << " iy " << eeid.iy() << ")";
617  }
618  EcalTriggerElectronicsId teid(getElectronicsMap()->getTriggerElectronicsId(eid));
619  ss << " (TCC " << teid.tccId() << " TT " << teid.ttId() << " pstrip " << teid.pseudoStripId() << " chan "
620  << teid.channelId() << ")";
621  break;
622  }
623  }
624  [[fallthrough]];
625  case kTriggerTower: {
626  // EB-03 DCC 12 TCC 18 TT 3
627  EcalTriggerElectronicsId teid(_rawId);
628  EcalElectronicsId eid(getElectronicsMap()->getElectronicsId(teid));
629  ss << smName(eid.dccId()) << " DCC " << eid.dccId() << " TCC " << teid.tccId() << " TT " << teid.ttId();
630  break;
631  }
632  case kSuperCrystal: {
633  // EB-03 DCC 12 CCU 18 (EBTT ieta -13 iphi 60)
634  EcalElectronicsId eid(_rawId);
635  ss << smName(eid.dccId()) << " DCC " << eid.dccId() << " CCU " << eid.towerId();
636  if (eid.dccId() >= kEBmLow + 1 && eid.dccId() <= kEBpHigh + 1) {
637  EcalTrigTowerDetId ttid(EBDetId(getElectronicsMap()->getDetId(eid)).tower());
638  ss << " (EBTT ieta " << std::showpos << ttid.ieta() << std::noshowpos << " iphi " << ttid.iphi() << ")";
639  } else {
640  EcalScDetId scid(EEDetId(getElectronicsMap()->getDetId(eid)).sc());
641  ss << " (EESC ix " << scid.ix() << " iy " << scid.iy() << ")";
642  }
643  break;
644  }
645  case kTCC: {
646  // EB-03 TCC 12
647  int tccid(_rawId - nDCC);
648  int dccid(getElectronicsMap()->DCCid(getElectronicsMap()->getTrigTowerDetId(tccid, 1)));
649  ss << smName(dccid) << " TCC " << (_rawId - nDCC);
650  break;
651  }
652  case kDCC:
653  ss << smName(_rawId);
654  break;
655  default:
656  break;
657  }
658 
659  return ss.str();
660  }
661 
662  uint32_t idFromName(std::string const &_name) {
663  TString name(_name);
664  TPRegexp re(
665  "(EB|EE)([+-][0-9][0-9])(?: TCC ([0-9]+)| DCC ([0-9]+) (CCU|TCC) "
666  "([0-9]+)(?: (TT|strip|PN) ([0-9]+)(?: xtal ([0-9]+)|)|)|)");
667  // 1 2 3 4 5 6 7 8 9
668  uint32_t rawId(0);
669 
670  TObjArray *matches(re.MatchS(name));
671  matches->SetOwner(true);
672  if (matches->GetEntries() == 0)
673  return 0;
674  else if (matches->GetEntries() == 3) {
675  TString subdet(static_cast<TObjString *>(matches->At(1))->GetString());
676  if (subdet == "EB") {
677  int dccid(static_cast<TObjString *>(matches->At(2))->GetString().Atoi());
678  unsigned offset(0);
679  if (dccid < 0) {
680  dccid *= -1;
681  offset = kEEmLow;
682  } else
683  offset = kEEpLow;
684  rawId = (dccid + 2) % 9 + 1 + offset;
685  } else {
686  int dccid(static_cast<TObjString *>(matches->At(2))->GetString().Atoi());
687  if (dccid < 0)
688  dccid *= -1;
689  else
690  dccid += 18;
691  rawId = kEBmLow + dccid;
692  }
693  } else if (matches->GetEntries() == 4)
694  rawId = static_cast<TObjString *>(matches->At(3))->GetString().Atoi() + nDCC;
695  else {
696  TString subtype(static_cast<TObjString *>(matches->At(5))->GetString());
697  if (subtype == "TCC") {
698  int tccid(static_cast<TObjString *>(matches->At(6))->GetString().Atoi());
699  int ttid(static_cast<TObjString *>(matches->At(8))->GetString().Atoi());
700  rawId = EcalTriggerElectronicsId(tccid, ttid, 1, 1).rawId();
701  } else {
702  int dccid(static_cast<TObjString *>(matches->At(4))->GetString().Atoi());
703  int towerid(static_cast<TObjString *>(matches->At(6))->GetString().Atoi());
704  if (matches->GetEntries() == 7)
705  rawId = EcalElectronicsId(dccid, towerid, 1, 1).rawId();
706  else {
707  TString chType(static_cast<TObjString *>(matches->At(7))->GetString());
708  int stripOrPNid(static_cast<TObjString *>(matches->At(8))->GetString().Atoi());
709  if (chType == "PN")
710  rawId = EcalElectronicsId(dccid, towerid, 1, stripOrPNid).rawId();
711  else if (chType == "strip") {
712  int xtalid(static_cast<TObjString *>(matches->At(9))->GetString().Atoi());
713  rawId = EcalElectronicsId(dccid, towerid, stripOrPNid, xtalid).rawId();
714  }
715  // case "TT" is already taken care of
716  }
717  }
718  }
719 
720  delete matches;
721 
722  return rawId;
723  }
724 
725  uint32_t idFromBin(ObjectType _otype, BinningType _btype, unsigned _iME, int _bin) {
726  if (_otype == kEB) {
727  if (_btype == kCrystal) {
728  int ieta(_bin / 362 - 86);
729  if (ieta >= 0)
730  ++ieta;
731  return EBDetId(ieta, _bin % 362);
732  } else if (_btype == kTriggerTower || _btype == kSuperCrystal) {
733  int ieta(_bin / 74 - 17);
734  int z(1);
735  if (ieta <= 0) {
736  z = -1;
737  ieta = -ieta + 1;
738  }
739  return EcalTrigTowerDetId(z, EcalBarrel, ieta, (_bin % 74 + 69) % 72 + 1);
740  }
741  } else if (_otype == kEEm || _otype == kEEp) {
742  int z(_otype == kEEm ? -1 : 1);
743  if (_btype == kCrystal || _btype == kTriggerTower)
744  return EEDetId(_bin % 102, _bin / 102, z).rawId();
745  else if (_btype == kSuperCrystal)
746  return EcalScDetId(_bin % 22, _bin / 22, z).rawId();
747  } else if (_otype == kEE) {
748  if (_btype == kCrystal || _btype == kTriggerTower) {
749  int ix(_bin % 202);
750  int z(ix > 100 ? 1 : -1);
751  if (z > 0)
752  ix = (ix - 100) % 101;
753  return EEDetId(ix, _bin / 202, z).rawId();
754  } else if (_btype == kSuperCrystal) {
755  int ix(_bin % 42);
756  int z(ix > 20 ? 1 : -1);
757  if (z > 0)
758  ix = (ix - 20) % 21;
759  return EcalScDetId(ix, _bin / 42, z).rawId();
760  }
761  } else if (_otype == kSM || _otype == kEBSM || _otype == kEESM) {
762  unsigned iSM(_iME);
763  if (_otype == kEBSM)
764  iSM += 9;
765  else if (_otype == kEESM && iSM > kEEmHigh)
766  iSM += nEBDCC;
767 
768  int z(iSM <= kEBmHigh ? -1 : 1);
769 
770  if ((iSM >= kEBmLow && iSM <= kEBmHigh) || (iSM >= kEBpLow && iSM <= kEBpHigh)) {
771  if (_btype == kCrystal) {
772  int iphi(((iSM - 9) % 18) * 20 + (z < 0 ? _bin / 87 : 21 - _bin / 87));
773  int ieta((_bin % 87) * z);
774  return EBDetId(ieta, iphi).rawId();
775  } else if (_btype == kTriggerTower || _btype == kSuperCrystal) {
776  int iphi((((iSM - 9) % 18) * 4 + (z < 0 ? _bin / 19 : 5 - _bin / 19) + 69) % 72 + 1);
777  int ieta(_bin % 19);
778  return EcalTrigTowerDetId(z, EcalBarrel, ieta, iphi).rawId();
779  }
780  } else {
781  if (_btype == kCrystal || _btype == kTriggerTower) {
782  int nX(nEESMX);
783  if (iSM == kEEm02 || iSM == kEEm08 || iSM == kEEp02 || iSM == kEEp08)
784  nX = nEESMXExt;
785  if (iSM == kEEm01 || iSM == kEEm05 || iSM == kEEm09 || iSM == kEEp01 || iSM == kEEp05 || iSM == kEEp09)
786  nX = nEESMXRed;
787  return EEDetId(_bin % (nX + 2) + xlow_(iSM), _bin / (nX + 2) + ylow_(iSM), z).rawId();
788  } else if (_btype == kSuperCrystal) {
789  int nX(nEESMX / 5);
790  if (iSM == kEEm02 || iSM == kEEm08 || iSM == kEEp02 || iSM == kEEp08)
791  nX = nEESMXExt / 5;
792  if (iSM == kEEm01 || iSM == kEEm05 || iSM == kEEm09 || iSM == kEEp01 || iSM == kEEp05 || iSM == kEEp09)
793  nX = nEESMXRed / 5;
794  return EcalScDetId(_bin % (nX + 2) + xlow_(iSM) / 5, _bin / (nX + 2) + ylow_(iSM) / 5, z).rawId();
795  }
796  }
797  }
798 
799  return 0;
800  }
801 
802  AxisSpecs formAxis(edm::ParameterSet const &_axisParams) {
803  AxisSpecs axis;
804 
805  if (_axisParams.existsAs<std::vector<double>>("edges", false)) {
806  std::vector<double> const &vEdges(_axisParams.getUntrackedParameter<std::vector<double>>("edges"));
807  axis.nbins = vEdges.size() - 1;
808  axis.edges = new float[vEdges.size()];
809  std::copy(vEdges.begin(), vEdges.end(), axis.edges);
810  } else {
811  axis.nbins = _axisParams.getUntrackedParameter<int>("nbins");
812  axis.low = _axisParams.getUntrackedParameter<double>("low");
813  bool highSet(_axisParams.existsAs<double>("high", false));
814  bool perBinSet(_axisParams.existsAs<double>("unitsPerBin", false));
815  if (highSet) {
816  if (perBinSet)
817  edm::LogWarning("EcalDQM") << "Maximum and bin width both set in an axis; using the former";
818  axis.high = _axisParams.getUntrackedParameter<double>("high");
819  } else if (perBinSet)
820  axis.high = axis.low + _axisParams.getUntrackedParameter<double>("unitsPerBin") * axis.nbins;
821  else
822  axis.high = 0.;
823  }
824 
825  if (_axisParams.existsAs<std::vector<std::string>>("labels", false)) {
826  std::vector<std::string> const &labels(_axisParams.getUntrackedParameter<std::vector<std::string>>("labels"));
827  if (int(labels.size()) == axis.nbins) {
828  axis.labels = new std::string[axis.nbins];
829  for (int iB(0); iB != axis.nbins; ++iB)
830  axis.labels[iB] = labels[iB];
831  }
832  }
833 
834  axis.title = _axisParams.getUntrackedParameter<std::string>("title");
835 
836  return axis;
837  }
838 
839  void fillAxisDescriptions(edm::ParameterSetDescription &_desc) {
840  _desc.addUntracked<std::string>("title", "");
841  _desc.addUntracked<int>("nbins", 0);
842  _desc.addUntracked<double>("low", 0.);
843  _desc.addOptionalNode(edm::ParameterDescription<double>("high", 0., false) ^
844  edm::ParameterDescription<double>("unitsPerBin", 0., false),
845  false);
846  _desc.addOptionalUntracked<std::vector<double>>("edges");
847  _desc.addOptionalUntracked<std::vector<std::string>>("labels");
848  }
849 
850  ObjectType translateObjectType(std::string const &_otypeName) {
851  if (_otypeName == "EB")
852  return kEB;
853  else if (_otypeName == "EE")
854  return kEE;
855  else if (_otypeName == "EEm")
856  return kEEm;
857  else if (_otypeName == "EEp")
858  return kEEp;
859  else if (_otypeName == "SM")
860  return kSM;
861  else if (_otypeName == "EBSM")
862  return kEBSM;
863  else if (_otypeName == "EESM")
864  return kEESM;
865  else if (_otypeName == "SMMEM")
866  return kSMMEM;
867  else if (_otypeName == "EBSMMEM")
868  return kEBSMMEM;
869  else if (_otypeName == "EESMMEM")
870  return kEESMMEM;
871  else if (_otypeName == "Ecal")
872  return kEcal;
873  else if (_otypeName == "MEM")
874  return kMEM;
875  else if (_otypeName == "EBMEM")
876  return kEBMEM;
877  else if (_otypeName == "EEMEM")
878  return kEEMEM;
879  else if (_otypeName == "Ecal2P")
880  return kEcal2P;
881  else if (_otypeName == "Ecal3P")
882  return kEcal3P;
883  else if (_otypeName == "EE2P")
884  return kEE2P;
885  else if (_otypeName == "MEM2P")
886  return kMEM2P;
887  else if (_otypeName == "Channel")
888  return kChannel;
889  else if (_otypeName == "None")
890  return nObjType;
891 
892  throw cms::Exception("InvalidConfiguration") << "No object type " << _otypeName << " defined";
893  }
894 
895  BinningType translateBinningType(std::string const &_btypeName) {
896  if (_btypeName == "Crystal")
897  return kCrystal;
898  else if (_btypeName == "TriggerTower")
899  return kTriggerTower;
900  else if (_btypeName == "SuperCrystal")
901  return kSuperCrystal;
902  else if (_btypeName == "PseudoStrip")
903  return kPseudoStrip;
904  else if (_btypeName == "TCC")
905  return kTCC;
906  else if (_btypeName == "DCC")
907  return kDCC;
908  else if (_btypeName == "ProjEta")
909  return kProjEta;
910  else if (_btypeName == "ProjPhi")
911  return kProjPhi;
912  else if (_btypeName == "RCT")
913  return kRCT;
914  else if (_btypeName == "User")
915  return kUser;
916  else if (_btypeName == "Report")
917  return kReport;
918  else if (_btypeName == "Trend")
919  return kTrend;
920 
921  throw cms::Exception("InvalidConfiguration") << "No binning type " << _btypeName << " defined";
922  }
923 
924  dqm::reco::MonitorElement::Kind translateKind(std::string const &_kindName) {
925  if (_kindName == "REAL")
926  return dqm::reco::MonitorElement::Kind::REAL;
927  else if (_kindName == "TH1F")
928  return dqm::reco::MonitorElement::Kind::TH1F;
929  else if (_kindName == "TProfile")
930  return dqm::reco::MonitorElement::Kind::TPROFILE;
931  else if (_kindName == "TH2F")
932  return dqm::reco::MonitorElement::Kind::TH2F;
933  else if (_kindName == "TProfile2D")
934  return dqm::reco::MonitorElement::Kind::TPROFILE2D;
935  else
936  return dqm::reco::MonitorElement::Kind::INVALID;
937  }
938  } // namespace binning
939 } // namespace ecaldqm
EcalScDetId::validDetId
static bool validDetId(int ix, int iy, int iz)
Definition: EcalScDetId.cc:59
ecaldqm::binning::findBinCrystal_
int findBinCrystal_(ObjectType, DetId const &, int=-1)
ecaldqm::memDCCIndex
unsigned memDCCIndex(unsigned)
Definition: EcalDQMCommonUtils.cc:48
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
ecaldqm::binning::findBinRCT_
int findBinRCT_(ObjectType, DetId const &)
Definition: MESetBinningUtils2.cc:523
ecaldqm::binning::findBinTriggerTower_
int findBinTriggerTower_(ObjectType, DetId const &)
Definition: MESetBinningUtils2.cc:539
EEDetId::ix
int ix() const
Definition: EEDetId.h:77
edm::ParameterSet::existsAs
bool existsAs(std::string const &parameterName, bool trackiness=true) const
checks if a parameter exists as a given type
Definition: ParameterSet.h:160
EcalTriggerElectronicsId::tccId
int tccId() const
get the DCC (Ecal Local DCC value not global one) id
Definition: EcalTriggerElectronicsId.h:31
EcalElectronicsId::xtalId
int xtalId() const
get the channel id
Definition: EcalElectronicsId.h:37
edm::ParameterSetDescription::addUntracked
ParameterDescriptionBase * addUntracked(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:100
EcalElectronicsId::stripId
int stripId() const
get the tower id
Definition: EcalElectronicsId.h:35
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
ecaldqm::binning::isValidIdBin
bool isValidIdBin(ObjectType, BinningType, unsigned, int)
Definition: MESetBinningUtils.cc:538
EcalElectronicsId
Ecal readout channel identification [32:20] Unused (so far) [19:13] DCC id [12:6] tower [5:3] strip [...
Definition: EcalElectronicsId.h:18
ecaldqm::binning::getNObjects
unsigned getNObjects(ObjectType)
Definition: MESetBinningUtils.cc:511
ecaldqm::getElectronicsMap
EcalElectronicsMapping const * getElectronicsMap()
Definition: EcalDQMCommonUtils.cc:438
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
ecaldqm::tccId
unsigned tccId(DetId const &)
EcalElectronicsId::towerId
int towerId() const
get the tower id
Definition: EcalElectronicsId.h:33
ecaldqm::zside
int zside(DetId const &)
ecaldqm::binning::getObject
ObjectType getObject(ObjectType, unsigned)
Definition: MESetBinningUtils.cc:466
ecaldqm::ttId
unsigned ttId(DetId const &)
EcalTriggerElectronicsId::channelId
int channelId() const
get the channel id
Definition: EcalTriggerElectronicsId.h:37
EcalTrigTowerDetId::ieta
int ieta() const
get the tower ieta
Definition: EcalTrigTowerDetId.h:44
ecaldqm::binning::translateObjectType
ObjectType translateObjectType(std::string const &)
Definition: MESetBinningUtils.cc:850
EBDetId::iphi
int iphi() const
get the crystal iphi
Definition: EBDetId.h:51
ecaldqm::binning::getBinning
AxisSpecs getBinning(ObjectType, BinningType, bool, int, unsigned)
Definition: MESetBinningUtils.cc:13
ecaldqm::isEndcapTTId
bool isEndcapTTId(DetId const &)
EcalPnDiodeDetId::iPnId
int iPnId() const
get the PnId
Definition: EcalPnDiodeDetId.h:46
ecaldqm::binning::findBinSuperCrystal_
int findBinSuperCrystal_(ObjectType, DetId const &, int=-1)
ecaldqm::binning::translateKind
dqm::legacy::MonitorElement::Kind translateKind(std::string const &)
Definition: MESetBinningUtils.cc:924
EcalScDetId::ix
int ix() const
Definition: EcalScDetId.h:70
EcalTriggerElectronicsId::ttId
int ttId() const
get the tower id
Definition: EcalTriggerElectronicsId.h:33
ecaldqm::binning::getBinningSMMEM_
AxisSpecs getBinningSMMEM_(BinningType, bool, unsigned, int)
Definition: MESetBinningUtils2.cc:306
DetId::subdetId
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:48
edm::ParameterSetDescription::addOptionalNode
ParameterDescriptionNode * addOptionalNode(ParameterDescriptionNode const &node, bool writeToCfi)
Definition: ParameterSetDescription.cc:50
ecaldqm::binning::getBinningEB_
AxisSpecs getBinningEB_(BinningType, bool, int)
Definition: MESetBinningUtils2.cc:89
hlt_dqm_clientPB-live_cfg.nbinsX
nbinsX
Definition: hlt_dqm_clientPB-live_cfg.py:51
photons_cff._id
_id
Definition: photons_cff.py:52
EEDetId::iy
int iy() const
Definition: EEDetId.h:83
EBDetId::ieta
int ieta() const
get the crystal ieta
Definition: EBDetId.h:49
ecaldqm::binning::idFromName
uint32_t idFromName(std::string const &)
Definition: MESetBinningUtils.cc:662
EcalTriggerElectronicsId::pseudoStripId
int pseudoStripId() const
get the tower id
Definition: EcalTriggerElectronicsId.h:35
ecaldqm::binning::getBinningSM_
AxisSpecs getBinningSM_(BinningType, bool, unsigned, int)
Definition: MESetBinningUtils2.cc:225
ecaldqm::binning::findBin1D
int findBin1D(ObjectType, BinningType, DetId const &)
EcalElectronicsId::dccId
int dccId() const
get the DCC (Ecal Local DCC value not global one) id
Definition: EcalElectronicsId.h:31
ecaldqm::towerId
unsigned towerId(DetId const &)
EcalScDetId::iy
int iy() const
Definition: EcalScDetId.h:76
ecaldqm::binning::channelName
std::string channelName(uint32_t, BinningType _btype=kDCC)
Definition: MESetBinningUtils.cc:590
EcalElectronicsId::rawId
uint32_t rawId() const
Definition: EcalElectronicsId.h:28
DetId
Definition: DetId.h:17
EcalTrigTowerDetId::iphi
int iphi() const
get the tower iphi
Definition: EcalTrigTowerDetId.h:52
ecaldqm::smName
std::string smName(unsigned)
Definition: EcalDQMCommonUtils.cc:289
ecaldqm::binning::findBin2D
int findBin2D(ObjectType, BinningType, DetId const &)
ecaldqm::binning::findPlotIndex
unsigned findPlotIndex(ObjectType, DetId const &)
EEDetId::validDetId
static bool validDetId(int crystal_ix, int crystal_iy, int iz)
Definition: EEDetId.h:248
ecaldqm::binning::getBinningEE_
AxisSpecs getBinningEE_(BinningType, bool, int, int)
Definition: MESetBinningUtils2.cc:154
ecaldqm::binning::fillAxisDescriptions
void fillAxisDescriptions(edm::ParameterSetDescription &)
Definition: MESetBinningUtils.cc:839
hgcalTowerProducer_cfi.tower
tower
Definition: hgcalTowerProducer_cfi.py:3
SummaryClient_cfi.labels
labels
Definition: SummaryClient_cfi.py:61
ecaldqm::binning::findBinPseudoStrip_
int findBinPseudoStrip_(ObjectType, DetId const &)
Definition: MESetBinningUtils2.cc:582
edm::ParameterSetDescription::addOptionalUntracked
ParameterDescriptionBase * addOptionalUntracked(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:110
ecaldqm::binning::translateBinningType
BinningType translateBinningType(std::string const &)
Definition: MESetBinningUtils.cc:895
ecaldqm::binning::formAxis
AxisSpecs formAxis(edm::ParameterSet const &)
Definition: MESetBinningUtils.cc:802
ecaldqm::binning::getBinningMEM_
AxisSpecs getBinningMEM_(BinningType, bool, int, int)
Definition: MESetBinningUtils2.cc:403
ecaldqm::dccId
unsigned dccId(DetId const &)
ecaldqm::binning::getBinningEcal_
AxisSpecs getBinningEcal_(BinningType, bool, int)
Definition: MESetBinningUtils2.cc:332
patCandidatesForDimuonsSequences_cff.matches
matches
Definition: patCandidatesForDimuonsSequences_cff.py:131
EcalTriggerElectronicsId
Ecal trigger electronics identification [32:20] Unused (so far) [19:13] TCC id [12:6] TT id [5:3] pse...
Definition: EcalTriggerElectronicsId.h:18
hltrates_dqm_sourceclient-live_cfg.offset
offset
Definition: hltrates_dqm_sourceclient-live_cfg.py:79
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
ecaldqm::binning::idFromBin
uint32_t idFromBin(ObjectType, BinningType, unsigned, int)
Definition: MESetBinningUtils.cc:725
Generated for CMSSW Reference Manual by doxygen 1.8.11