CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
Masks.cc
Go to the documentation of this file.
1 // $Id: Masks.cc,v 1.18 2012/04/27 13:46:04 yiiyama Exp $
2 
11 #include <sstream>
12 #include <iomanip>
13 
14 #include "FWCore/Framework/interface/ESHandle.h"
15 #include "CondFormats/DataRecord/interface/EcalDQMChannelStatusRcd.h"
16 #include "CondFormats/DataRecord/interface/EcalDQMTowerStatusRcd.h"
17 
18 #include "DataFormats/EcalDetId/interface/EBDetId.h"
19 #include "DataFormats/EcalDetId/interface/EEDetId.h"
20 
21 #include "DQM/EcalCommon/interface/Numbers.h"
22 
23 #include "DQM/EcalCommon/interface/Masks.h"
24 
25 //-------------------------------------------------------------------------
26 
27 const EcalDQMChannelStatus* Masks::channelStatus = 0;
28 const EcalDQMTowerStatus* Masks::towerStatus = 0;
29 
30 bool Masks::init = false;
31 
32 //-------------------------------------------------------------------------
33 
34 void Masks::initMasking( const edm::EventSetup& setup, bool verbose ) {
35 
36  if ( Masks::init ) return;
37 
38  if ( verbose ) std::cout << "Initializing EcalDQMChannelStatus and EcalDQMTowerStatus ..." << std::endl;
39 
40  Masks::init = true;
41 
42  if ( setup.find( edm::eventsetup::EventSetupRecordKey::makeKey<EcalDQMChannelStatusRcd>() ) ) {
43  edm::ESHandle<EcalDQMChannelStatus> handle;
44  setup.get<EcalDQMChannelStatusRcd>().get(handle);
45  if ( handle.isValid() ) Masks::channelStatus = handle.product();
46  }
47 
48  if ( setup.find( edm::eventsetup::EventSetupRecordKey::makeKey<EcalDQMTowerStatusRcd>() ) ) {
49  edm::ESHandle<EcalDQMTowerStatus> handle;
50  setup.get<EcalDQMTowerStatusRcd>().get(handle);
51  if ( handle.isValid() ) Masks::towerStatus = handle.product();
52  }
53 
54  if ( verbose ) std::cout << "done." << std::endl;
55 
56 }
57 
58 //-------------------------------------------------------------------------
59 
60 bool Masks::maskChannel( int ism, int i1, int i2, uint32_t bits, const EcalSubdetector subdet ) throw( cms::Exception ) {
61 
62  bool mask = false;
63 
64  if ( subdet == EcalBarrel ) {
65 
66  int jsm = Numbers::iSM(ism, EcalBarrel);
67  int ic = 20*(i1-1)+(i2-1)+1;
68 
69  EBDetId id(jsm, ic, EBDetId::SMCRYSTALMODE);
70  if ( Masks::channelStatus ) {
71  EcalDQMChannelStatus::const_iterator it = Masks::channelStatus->find( id.rawId() );
72  if ( it != Masks::channelStatus->end() ) mask |= it->getStatusCode() & bits;
73  }
74  if ( Masks::towerStatus ) {
75  EcalDQMTowerStatus::const_iterator it = Masks::towerStatus->find( id.tower().rawId() );
76  if ( it != Masks::towerStatus->end() ) mask |= it->getStatusCode() & bits;
77  }
78 
79  } else if ( subdet == EcalEndcap ) {
80 
81  int jx = i1 + Numbers::ix0EE(ism);
82  int jy = i2 + Numbers::iy0EE(ism);
83 
84  if ( ism >= 1 && ism <= 9 ) jx = 101 - jx;
85 
86  if ( Numbers::validEE(ism, jx, jy) ) {
87  EEDetId id(jx, jy, (ism>=1&&ism<=9)?-1:+1);
88  if ( Masks::channelStatus ) {
89  EcalDQMChannelStatus::const_iterator it = Masks::channelStatus->find( id.rawId() );
90  if ( it != Masks::channelStatus->end() ) mask |= it->getStatusCode() & bits;
91  }
92  if ( Masks::towerStatus ) {
93  EcalDQMTowerStatus::const_iterator it = Masks::towerStatus->find( id.sc().rawId() );
94  if ( it != Masks::towerStatus->end() ) mask |= it->getStatusCode() & bits;
95  }
96  }
97 
98  } else {
99 
100  std::ostringstream s;
101  s << "Invalid subdetector: subdet = " << subdet;
102  throw( cms::Exception( s.str() ) );
103 
104  }
105 
106  return ( mask );
107 
108 }
109 
110 //-------------------------------------------------------------------------
111 
112 bool Masks::maskPn( int ism, int i1, uint32_t bits, const EcalSubdetector subdet ) throw( cms::Exception ) {
113 
114  bool mask = false;
115 
116  if ( subdet == EcalBarrel ) {
117 
118  // EB-03
119  if ( ism == 3 && i1 == 1 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
120  if ( ism == 3 && i1 == 1 && (bits & (1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR)) ) mask = true;
121 
122  // EB-07
123  if ( ism == 7 && i1 == 4 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
124  if ( ism == 7 && i1 == 4 && (bits & (1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR)) ) mask = true;
125 
126  // EB-15
127  if ( ism == 15 && i1 == 6 && (bits & (1 << EcalDQMStatusHelper::TT_SIZE_ERROR)) ) mask = true;
128  if ( ism == 15 && i1 == 6 && (bits & (1 << EcalDQMStatusHelper::STATUS_FLAG_ERROR)) ) mask = true;
129  if ( ism == 15 && i1 == 6 && (bits & (1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR)) ) mask = true;
130  if ( ism == 15 && i1 == 6 && (bits & (1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR)) ) mask = true;
131  if ( ism == 15 && i1 == 6 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
132  if ( ism == 15 && i1 == 7 && (bits & (1 << EcalDQMStatusHelper::TT_SIZE_ERROR)) ) mask = true;
133  if ( ism == 15 && i1 == 7 && (bits & (1 << EcalDQMStatusHelper::STATUS_FLAG_ERROR)) ) mask = true;
134  if ( ism == 15 && i1 == 7 && (bits & (1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR)) ) mask = true;
135  if ( ism == 15 && i1 == 7 && (bits & (1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR)) ) mask = true;
136  if ( ism == 15 && i1 == 7 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
137  if ( ism == 15 && i1 == 8 && (bits & (1 << EcalDQMStatusHelper::TT_SIZE_ERROR)) ) mask = true;
138  if ( ism == 15 && i1 == 8 && (bits & (1 << EcalDQMStatusHelper::STATUS_FLAG_ERROR)) ) mask = true;
139  if ( ism == 15 && i1 == 8 && (bits & (1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR)) ) mask = true;
140  if ( ism == 15 && i1 == 8 && (bits & (1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR)) ) mask = true;
141  if ( ism == 15 && i1 == 8 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
142  if ( ism == 15 && i1 == 9 && (bits & (1 << EcalDQMStatusHelper::TT_SIZE_ERROR)) ) mask = true;
143  if ( ism == 15 && i1 == 9 && (bits & (1 << EcalDQMStatusHelper::STATUS_FLAG_ERROR)) ) mask = true;
144  if ( ism == 15 && i1 == 9 && (bits & (1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR)) ) mask = true;
145  if ( ism == 15 && i1 == 9 && (bits & (1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR)) ) mask = true;
146  if ( ism == 15 && i1 == 9 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
147  if ( ism == 15 && i1 == 10 && (bits & (1 << EcalDQMStatusHelper::TT_SIZE_ERROR)) ) mask = true;
148  if ( ism == 15 && i1 == 10 && (bits & (1 << EcalDQMStatusHelper::STATUS_FLAG_ERROR)) ) mask = true;
149  if ( ism == 15 && i1 == 10 && (bits & (1 << EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR)) ) mask = true;
150  if ( ism == 15 && i1 == 10 && (bits & (1 << EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR)) ) mask = true;
151  if ( ism == 15 && i1 == 10 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
152 
153  // EB+06
154  if ( ism == 24 && (bits & (1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR)) ) mask = true;
155 
156  // EB+07
157  if ( ism == 25 && i1 == 4 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
158  if ( ism == 25 && i1 == 4 && (bits & (1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR)) ) mask = true;
159 
160  // EB+12
161  if ( ism == 30 && i1 == 9 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
162  if ( ism == 30 && i1 == 9 && (bits & (1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR)) ) mask = true;
163 
164  // EB+15
165  if ( ism == 15 && i1 == 3 && (bits & (1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR)) ) mask = true;
166  if ( ism == 15 && i1 == 4 && (bits & (1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR)) ) mask = true;
167  if ( ism == 15 && i1 == 5 && (bits & (1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR)) ) mask = true;
168 
169  } else if ( subdet == EcalEndcap ) {
170 
171  // EE-02
172  if ( ism == 5 && i1 == 3 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
173  if ( ism == 5 && i1 == 3 && (bits & (1 << EcalDQMStatusHelper::LED_MEAN_ERROR)) ) mask = true;
174  if ( ism == 5 && i1 == 10 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
175  if ( ism == 5 && i1 == 10 && (bits & (1 << EcalDQMStatusHelper::LED_MEAN_ERROR)) ) mask = true;
176 
177  // EE-03
178  if ( ism == 6 && i1 == 5 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
179  if ( ism == 6 && i1 == 5 && (bits & (1 << EcalDQMStatusHelper::LED_MEAN_ERROR)) ) mask = true;
180 
181  // EE-07
182  if ( ism == 1 && i1 == 4 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
183  if ( ism == 1 && i1 == 4 && (bits & (1 << EcalDQMStatusHelper::LED_MEAN_ERROR)) ) mask = true;
184  if ( ism == 1 && i1 == 10 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
185  if ( ism == 1 && i1 == 10 && (bits & (1 << EcalDQMStatusHelper::LED_MEAN_ERROR)) ) mask = true;
186 
187  // EE-08
188  if ( ism == 2 && i1 == 1 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
189  if ( ism == 2 && i1 == 1 && (bits & (1 << EcalDQMStatusHelper::LED_MEAN_ERROR)) ) mask = true;
190  if ( ism == 2 && i1 == 4 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
191  if ( ism == 2 && i1 == 4 && (bits & (1 << EcalDQMStatusHelper::LED_MEAN_ERROR)) ) mask = true;
192  if ( ism == 2 && i1 == 4 && (bits & (1 << EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR)) ) mask = true;
193  if ( ism == 2 && i1 == 10 && (bits & (1 << EcalDQMStatusHelper::LASER_MEAN_ERROR)) ) mask = true;
194  if ( ism == 2 && i1 == 10 && (bits & (1 << EcalDQMStatusHelper::LED_MEAN_ERROR)) ) mask = true;
195 
196  } else {
197 
198  std::ostringstream s;
199  s << "Invalid subdetector: subdet = " << subdet;
200  throw( cms::Exception( s.str() ) );
201 
202  }
203 
204  return ( mask );
205 
206 }
207 
208 //-------------------------------------------------------------------------
209 
EcalDQMStatusHelper::LED_MEAN_ERROR
static const int LED_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:40
Numbers::validEE
static bool validEE(const unsigned ism, const unsigned ix, const unsigned iy)
Definition: Numbers.cc:832
EcalDQMStatusHelper::LASER_MEAN_ERROR
static const int LASER_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:35
Masks::maskChannel
static bool maskChannel(int ism, int i1, int i2, uint32_t bits, const EcalSubdetector subdet)
Definition: Masks.cc:60
Masks::init
static bool init
Definition: Masks.h:40
Numbers.h
Some &quot;id&quot; conversions.
Masks::towerStatus
static const EcalDQMTowerStatus * towerStatus
Definition: Masks.h:45
EcalDQMStatusHelper::PEDESTAL_HIGH_GAIN_MEAN_ERROR
static const int PEDESTAL_HIGH_GAIN_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:20
bits
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision bits
Definition: EventSelector-behavior.doc:35
EcalDQMStatusHelper::TT_SIZE_ERROR
static const int TT_SIZE_ERROR
Definition: EcalDQMStatusHelper.h:16
Numbers::ix0EE
static int ix0EE(const unsigned ism)
Definition: Numbers.cc:773
EcalCondTowerObjectContainer::find
const_iterator find(uint32_t rawId) const
Definition: EcalCondTowerObjectContainer.h:57
Masks::maskPn
static bool maskPn(int ism, int i1, uint32_t bits, const EcalSubdetector subdet)
Definition: Masks.cc:112
Numbers::iy0EE
static int iy0EE(const unsigned ism)
Definition: Numbers.cc:812
EcalDQMStatusHelper::STATUS_FLAG_ERROR
static const int STATUS_FLAG_ERROR
Definition: EcalDQMStatusHelper.h:45
edm::EventSetup::find
const eventsetup::EventSetupRecord * find(const eventsetup::EventSetupRecordKey &) const
Definition: EventSetup.cc:90
Masks.h
channel masking
patZpeak.handle
tuple handle
Definition: patZpeak.py:22
end
#define end
Definition: vmac.h:38
Masks::initMasking
static void initMasking(const edm::EventSetup &setup, bool verbose=false)
Definition: Masks.cc:34
EcalDQMStatusHelper::PEDESTAL_LOW_GAIN_MEAN_ERROR
static const int PEDESTAL_LOW_GAIN_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:18
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
EcalCondObjectContainer::const_iterator
std::vector< Item >::const_iterator const_iterator
Definition: EcalCondObjectContainer.h:15
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
Numbers::iSM
static unsigned iSM(const unsigned ism, const EcalSubdetector subdet)
Definition: Numbers.cc:246
Masks::channelStatus
static const EcalDQMChannelStatus * channelStatus
Definition: Masks.h:44
EcalCondObjectContainer::find
const_iterator find(uint32_t rawId) const
Definition: EcalCondObjectContainer.h:58
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
EcalCondTowerObjectContainer::const_iterator
std::vector< Item >::const_iterator const_iterator
Definition: EcalCondTowerObjectContainer.h:18
EcalDQMStatusHelper::TESTPULSE_HIGH_GAIN_MEAN_ERROR
static const int TESTPULSE_HIGH_GAIN_MEAN_ERROR
Definition: EcalDQMStatusHelper.h:30
edm::ESHandleBase::isValid
bool isValid() const
Definition: ESHandle.h:37
EBDetId::SMCRYSTALMODE
static const int SMCRYSTALMODE
Definition: EBDetId.h:168
EcalSubdetector
EcalSubdetector
Definition: EcalSubdetector.h:12
ecalpyutils::ism
int ism(int ieta, int iphi)
Definition: EcalPyUtils.cc:56
HcalObjRepresent::setup
void setup(std::vector< TH2F > &depth, std::string name, std::string units="")
Definition: HcalObjRepresent.h:628
Generated for CMSSW Reference Manual by doxygen 1.8.5