Operations.
92 LogVerbatim(
"NoiseDBValidation") <<
"[DTnoiseDBValidation] Parameters initialization";
99 dbe_->book1D(
"noisyCellDiff",
"percentual (wrt the previous db) total number of noisy cells", 1, 0.5, 1.5);
101 wheelHisto_ =
dbe_->book1D(
"wheelOccupancy",
"percentual noisy cells occupancy per wheel", 5, -2.5, 2.5);
107 stationHisto_ =
dbe_->book1D(
"stationOccupancy",
"percentual noisy cells occupancy per station", 4, 0.5, 4.5);
112 sectorHisto_ =
dbe_->book1D(
"sectorOccupancy",
"percentual noisy cells occupancy per sector", 12, 0.5, 12.5);
125 layerHisto_ =
dbe_->book1D(
"layerOccupancy",
"percentual noisy cells occupancy per layer", 3, 0.5, 3.5);
135 map<int, int> sectMap;
137 map<int, int> layerMap;
142 DTWireId wireId((*noise).first.wheelId,
143 (*noise).first.stationId,
144 (*noise).first.sectorId,
146 (*noise).first.layerId,
147 (*noise).first.cellId);
148 LogVerbatim(
"NoiseDBValidation") <<
"Ref. noisy wire: " << wireId;
154 DTWireId wireId((*noise).first.wheelId,
155 (*noise).first.stationId,
156 (*noise).first.sectorId,
158 (*noise).first.layerId,
159 (*noise).first.cellId);
160 LogVerbatim(
"NoiseDBValidation") <<
"Valid. noisy wire: " << wireId;
163 whMap[(*noise).first.wheelId]++;
164 stMap[(*noise).first.stationId]++;
165 sectMap[(*noise).first.sectorId]++;
169 if ((*noise).first.cellId <= 10)
171 if ((*noise).first.cellId > 10 && (*noise).first.cellId < (lastWire - 10))
173 if ((*noise).first.cellId >= (lastWire - 10))
176 const DTChamberId chId = wireId.layerId().superlayerId().chamberId();
179 int binNumber = 4 * (wireId.superLayer() - 1) + wireId.layer();
188 for (map<int, int>::const_iterator
wheel = whMap.begin();
wheel != whMap.end();
wheel++) {
194 for (map<int, int>::const_iterator sector = sectMap.begin(); sector != sectMap.end(); sector++) {
197 for (map<int, int>::const_iterator layer = layerMap.begin(); layer != layerMap.end(); layer++) {
MonitorElement * wheelHisto_
const_iterator end() const
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
void bookHisto(const DTChamberId &)
const DTStatusFlag * noiseRefMap_
const DTStatusFlag * noiseMap_
int lastChannel() const
Returns the wire number of the last wire.
def binNumber(station, sl)
std::map< DTChamberId, MonitorElement * > noiseHistoMap_
Abs< T >::type abs(const T &t)
std::vector< std::pair< DTStatusFlagId, DTStatusFlagData > >::const_iterator const_iterator
Access methods to data.
MonitorElement * diffHisto_
edm::ESHandle< DTGeometry > dtGeom_
MonitorElement * stationHisto_
const DTLayer * layer(const DTLayerId &id) const
Return a layer given its id.
MonitorElement * layerHisto_
const_iterator begin() const
MonitorElement * sectorHisto_