Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-07-05 04:09:42

0001 #include "DQM/SiStripMonitorSummary/interface/SiStripQualityDQM.h"
0002 #include "DQMServices/Core/interface/DQMStore.h"
0003 #include "TCanvas.h"
0004 
0005 SiStripQualityDQM::SiStripQualityDQM(edm::ESGetToken<SiStripQuality, SiStripQualityRcd> token,
0006                                      edm::RunNumber_t iRun,
0007                                      edm::ParameterSet const &hPSet,
0008                                      edm::ParameterSet const &fPSet,
0009                                      const TrackerTopology *tTopo,
0010                                      const TkDetMap *tkDetMap)
0011     : SiStripBaseCondObjDQMGet<SiStripQuality, SiStripQualityRcd>{token, iRun, hPSet, fPSet, tTopo} {
0012   if (HistoMaps_On_) {
0013     Tk_HM_ = std::make_unique<TkHistoMap>(tkDetMap, "SiStrip/Histo_Map", "Quality_TkMap", 0.);
0014   }
0015 }
0016 
0017 SiStripQualityDQM::~SiStripQualityDQM() {}
0018 
0019 void SiStripQualityDQM::getActiveDetIds(const edm::EventSetup &eSetup) {
0020   getConditionObject(eSetup);
0021   condObj_->getDetIds(activeDetIds);
0022 }
0023 
0024 void SiStripQualityDQM::fillModMEs(const std::vector<uint32_t> &selectedDetIds) {
0025   ModMEs CondObj_ME;
0026   for (const auto det : selectedDetIds) {
0027     fillMEsForDet(CondObj_ME, det);
0028   }
0029 }
0030 
0031 void SiStripQualityDQM::fillMEsForDet(const ModMEs &_selModME_, uint32_t selDetId_) {
0032   ModMEs selModME_ = _selModME_;
0033   getModMEs(selModME_, selDetId_);
0034 
0035   const auto qualityRange = condObj_->getRange(selDetId_);
0036   int nStrip = detInfo_.getNumberOfApvsAndStripLength(selDetId_).first * 128;
0037 
0038   for (int istrip = 0; istrip < nStrip; ++istrip) {
0039     selModME_.ProfileDistr->Fill(istrip + 1, condObj_->IsStripBad(qualityRange, istrip) ? 0. : 1.);
0040 
0041   }  // istrip
0042 }
0043 
0044 void SiStripQualityDQM::fillSummaryMEs(const std::vector<uint32_t> &selectedDetIds) {
0045   for (const auto det : selectedDetIds) {
0046     fillMEsForLayer(/*SummaryMEsMap_,*/ det);
0047   }
0048 
0049   for (const auto &itm : SummaryMEsMap_) {
0050     ModMEs selME;
0051     selME = itm.second;
0052 
0053     if (hPSet_.getParameter<bool>("FillSummaryAtLayerLevel") &&
0054         fPSet_.getParameter<bool>("OutputSummaryAtLayerLevelAsImage")) {
0055       TCanvas c1("c1");
0056       selME.SummaryDistr->getTH1()->Draw();
0057       std::string name(selME.SummaryDistr->getTitle());
0058       name += ".png";
0059       c1.Print(name.c_str());
0060     }
0061   }
0062 }
0063 
0064 void SiStripQualityDQM::fillMEsForLayer(
0065     /* std::map<uint32_t, ModMEs> selMEsMap_,*/ uint32_t selDetId_) {
0066   float numberOfBadStrips = 0;
0067 
0068   SiStripHistoId hidmanager;
0069 
0070   if (hPSet_.getParameter<bool>("FillSummaryAtLayerLevel")) {
0071     std::string hSummary_description;
0072     hSummary_description = hPSet_.getParameter<std::string>("Summary_description");
0073 
0074     std::string hSummary_name;
0075 
0076     // ----
0077     int subDetId_ = ((selDetId_ >> 25) & 0x7);
0078 
0079     if (subDetId_ < 3 || subDetId_ > 6) {
0080       edm::LogError("SiStripQualityDQM") << "[SiStripQualityDQM::fillMEsForLayer] WRONG INPUT : no such "
0081                                             "subdetector type : "
0082                                          << subDetId_ << " no folder set!" << std::endl;
0083       return;
0084     }
0085     // ----
0086 
0087     hSummary_name = hidmanager.createHistoLayer(hSummary_description, "layer", getLayerNameAndId(selDetId_).first, "");
0088 
0089     const auto selMEsMapIter_ = SummaryMEsMap_.find(getLayerNameAndId(selDetId_).second);
0090 
0091     ModMEs selME_;
0092     if (selMEsMapIter_ != SummaryMEsMap_.end())
0093       selME_ = selMEsMapIter_->second;
0094 
0095     getSummaryMEs(selME_, selDetId_);
0096 
0097     std::vector<uint32_t> sameLayerDetIds_;
0098     sameLayerDetIds_.clear();
0099     sameLayerDetIds_ = GetSameLayerDetId(activeDetIds, selDetId_);
0100     // -----
0101     //  unsigned int iBin=0;
0102 
0103     //%%%%%%%%%%%%%%%%%%%
0104     //%%% FIXME: patch to fix the bug causing double counting on each layer
0105     //%%%%%%%%%%%%%%%%%%%%%%%
0106     if (std::binary_search(alreadyFilledLayers.begin(), alreadyFilledLayers.end(), sameLayerDetIds_[0]))
0107       return;
0108     alreadyFilledLayers.push_back(sameLayerDetIds_[0]);
0109     sort(alreadyFilledLayers.begin(), alreadyFilledLayers.end());
0110     //%%%%%%%%%%%%%%%%%%%%%%% END %%%%%%%%%%
0111 
0112     for (unsigned int i = 0; i < sameLayerDetIds_.size(); i++) {
0113       const auto qualityRange = condObj_->getRange(sameLayerDetIds_[i]);
0114       int nStrip = detInfo_.getNumberOfApvsAndStripLength(sameLayerDetIds_[i]).first * 128;
0115 
0116       numberOfBadStrips = 0;
0117 
0118       for (int istrip = 0; istrip < nStrip; ++istrip) {
0119         if (condObj_->IsStripBad(qualityRange, istrip)) {
0120           numberOfBadStrips++;
0121         }
0122       }
0123 
0124       float fr = 100 * float(numberOfBadStrips) / nStrip;
0125       selME_.SummaryDistr->Fill(i + 1, fr);
0126       if (fr > 20) {
0127         char c[9];
0128         sprintf(c, "%d", sameLayerDetIds_[i]);
0129         selME_.SummaryDistr->setBinLabel(i + 1, c);
0130       }
0131 
0132       // Fill the TkHistoMap with Quality output :
0133       if (HistoMaps_On_)
0134         Tk_HM_->setBinContent(sameLayerDetIds_[i], fr);
0135 
0136       // Fill the TkMap
0137       if (fPSet_.getParameter<bool>("TkMap_On") || hPSet_.getParameter<bool>("TkMap_On")) {
0138         fillTkMap(sameLayerDetIds_[i], fr);
0139       }
0140     }
0141   }  // if Fill ...
0142 }
0143 
0144 void SiStripQualityDQM::fillGrandSummaryMEs() {
0145   std::string hSummary_BadObjects_xTitle = hPSet_.getParameter<std::string>("Summary_BadObjects_histo_xTitle");
0146 
0147   std::string hSummary_BadModules_name = hPSet_.getParameter<std::string>("Summary_BadModules_histo_name");
0148   std::string hSummary_BadModules_yTitle = hPSet_.getParameter<std::string>("Summary_BadModules_histo_yTitle");
0149 
0150   std::string hSummary_BadFibers_name = hPSet_.getParameter<std::string>("Summary_BadFibers_histo_name");
0151   std::string hSummary_BadFibers_yTitle = hPSet_.getParameter<std::string>("Summary_BadFibers_histo_yTitle");
0152 
0153   std::string hSummary_BadApvs_name = hPSet_.getParameter<std::string>("Summary_BadApvs_histo_name");
0154   std::string hSummary_BadApvs_yTitle = hPSet_.getParameter<std::string>("Summary_BadApvs_histo_yTitle");
0155 
0156   std::string hSummary_BadStrips_name = hPSet_.getParameter<std::string>("Summary_BadStrips_histo_name");
0157   std::string hSummary_BadStrips_yTitle = hPSet_.getParameter<std::string>("Summary_BadStrips_histo_yTitle");
0158 
0159   int NchX = 34;
0160   double LowX = 0.5;
0161   double HighX = 34.5;
0162 
0163   MonitorElement *ME[4];
0164 
0165   DQMStore *dqmStore_ = edm::Service<DQMStore>().operator->();
0166 
0167   std::string FolderName = fPSet_.getParameter<std::string>("FolderName_For_QualityAndCabling_SummaryHistos");
0168 
0169   dqmStore_->setCurrentFolder(FolderName);
0170 
0171   ME[0] = dqmStore_->book1D(hSummary_BadModules_name, hSummary_BadModules_name, NchX, LowX, HighX);
0172   ME[0]->setAxisTitle(hSummary_BadObjects_xTitle, 1);
0173   ME[0]->setAxisTitle(hSummary_BadModules_yTitle, 2);
0174 
0175   ME[1] = dqmStore_->book1D(hSummary_BadFibers_name, hSummary_BadFibers_name, NchX, LowX, HighX);
0176   ME[1]->setAxisTitle(hSummary_BadObjects_xTitle, 1);
0177   ME[1]->setAxisTitle(hSummary_BadFibers_yTitle, 2);
0178 
0179   ME[2] = dqmStore_->book1D(hSummary_BadApvs_name, hSummary_BadApvs_name, NchX, LowX, HighX);
0180   ME[2]->setAxisTitle(hSummary_BadObjects_xTitle, 1);
0181   ME[2]->setAxisTitle(hSummary_BadApvs_yTitle, 2);
0182 
0183   ME[3] = dqmStore_->book1D(hSummary_BadStrips_name, hSummary_BadStrips_name, NchX, LowX, HighX);
0184   ME[3]->setAxisTitle(hSummary_BadObjects_xTitle, 1);
0185   ME[3]->setAxisTitle(hSummary_BadStrips_yTitle, 2);
0186 
0187   //==============================
0188 
0189   for (int i = 0; i < 4; ++i) {
0190     NTkBadComponent[i] = 0;
0191     for (int j = 0; j < 19; ++j) {
0192       ssV[i][j].str("");
0193       for (int k = 0; k < 4; ++k)
0194         NBadComponent[i][j][k] = 0;
0195     }
0196   }
0197 
0198   std::stringstream ss;
0199   ss.str("");
0200   for (const auto det : detInfo_.getAllDetIds()) {
0201     ss << "detid " << det << " IsModuleUsable " << condObj_->IsModuleUsable(det) << "\n";
0202   }
0203   LogDebug("SiStripQualityDQM") << ss.str() << std::endl;
0204 
0205   std::vector<SiStripQuality::BadComponent> BC = condObj_->getBadComponentList();
0206 
0207   for (size_t i = 0; i < BC.size(); ++i) {
0208     //&&&&&&&&&&&&&
0209     // Full Tk
0210     //&&&&&&&&&&&&&
0211 
0212     if (BC[i].BadModule)
0213       NTkBadComponent[0]++;
0214     if (BC[i].BadFibers)
0215       NTkBadComponent[1] += ((BC[i].BadFibers >> 2) & 0x1) + ((BC[i].BadFibers >> 1) & 0x1) + ((BC[i].BadFibers) & 0x1);
0216     if (BC[i].BadApvs)
0217       NTkBadComponent[2] += ((BC[i].BadApvs >> 5) & 0x1) + ((BC[i].BadApvs >> 4) & 0x1) + ((BC[i].BadApvs >> 3) & 0x1) +
0218                             ((BC[i].BadApvs >> 2) & 0x1) + ((BC[i].BadApvs >> 1) & 0x1) + ((BC[i].BadApvs) & 0x1);
0219 
0220     //&&&&&&&&&&&&&&&&&
0221     // Single SubSyste
0222     //&&&&&&&&&&&&&&&&&
0223 
0224     int component;
0225     SiStripDetId a(BC[i].detid);
0226     if (a.subdetId() == SiStripDetId::TIB) {
0227       //&&&&&&&&&&&&&&&&&
0228       // TIB
0229       //&&&&&&&&&&&&&&&&&
0230 
0231       component = tTopo_->tibLayer(BC[i].detid);
0232       SetBadComponents(0, component, BC[i]);
0233 
0234     } else if (a.subdetId() == SiStripDetId::TID) {
0235       //&&&&&&&&&&&&&&&&&
0236       // TID
0237       //&&&&&&&&&&&&&&&&&
0238 
0239       component = tTopo_->tidSide(BC[i].detid) == 2 ? tTopo_->tidWheel(BC[i].detid) : tTopo_->tidWheel(BC[i].detid) + 3;
0240       SetBadComponents(1, component, BC[i]);
0241 
0242     } else if (a.subdetId() == SiStripDetId::TOB) {
0243       //&&&&&&&&&&&&&&&&&
0244       // TOB
0245       //&&&&&&&&&&&&&&&&&
0246 
0247       component = tTopo_->tobLayer(BC[i].detid);
0248       SetBadComponents(2, component, BC[i]);
0249 
0250     } else if (a.subdetId() == SiStripDetId::TEC) {
0251       //&&&&&&&&&&&&&&&&&
0252       // TEC
0253       //&&&&&&&&&&&&&&&&&
0254 
0255       component = tTopo_->tecSide(BC[i].detid) == 2 ? tTopo_->tecWheel(BC[i].detid) : tTopo_->tecWheel(BC[i].detid) + 9;
0256       SetBadComponents(3, component, BC[i]);
0257     }
0258   }
0259 
0260   //&&&&&&&&&&&&&&&&&&
0261   // Single Strip Info
0262   //&&&&&&&&&&&&&&&&&&
0263 
0264   SiStripQuality::RegistryIterator rbegin = condObj_->getRegistryVectorBegin();
0265   SiStripQuality::RegistryIterator rend = condObj_->getRegistryVectorEnd();
0266 
0267   for (SiStripBadStrip::RegistryIterator rp = rbegin; rp != rend; ++rp) {
0268     uint32_t detid = rp->detid;
0269 
0270     int subdet = 0;
0271     int component = 0;
0272     SiStripDetId a(detid);
0273     if (a.subdetId() == 3) {
0274       subdet = 0;
0275       component = tTopo_->tibLayer(detid);
0276     } else if (a.subdetId() == 4) {
0277       subdet = 1;
0278       component = tTopo_->tidSide(detid) == 2 ? tTopo_->tidWheel(detid) : tTopo_->tidWheel(detid) + 3;
0279     } else if (a.subdetId() == 5) {
0280       subdet = 2;
0281       component = tTopo_->tobLayer(detid);
0282     } else if (a.subdetId() == 6) {
0283       subdet = 3;
0284       component = tTopo_->tecSide(detid) == 2 ? tTopo_->tecWheel(detid) : tTopo_->tecWheel(detid) + 9;
0285     }
0286 
0287     SiStripQuality::Range sqrange =
0288         SiStripQuality::Range(condObj_->getDataVectorBegin() + rp->ibegin, condObj_->getDataVectorBegin() + rp->iend);
0289 
0290     for (int it = 0; it < sqrange.second - sqrange.first; it++) {
0291       unsigned int range = condObj_->decode(*(sqrange.first + it)).range;
0292       NTkBadComponent[3] += range;
0293       NBadComponent[subdet][0][3] += range;
0294       NBadComponent[subdet][component][3] += range;
0295     }
0296   }
0297 
0298   //&&&&&&&&&&&&&&&&&&
0299   // printout
0300   //&&&&&&&&&&&&&&&&&&
0301 
0302   ss.str("");
0303   ss << "\n-----------------\nGlobal Info\n-----------------";
0304   ss << "\nBadComponent \t   Modules \tFibers "
0305         "\tApvs\tStrips\n------------------------------------------------------"
0306         "----------";
0307   ss << "\nTracker:\t\t" << NTkBadComponent[0] << "\t" << NTkBadComponent[1] << "\t" << NTkBadComponent[2] << "\t"
0308      << NTkBadComponent[3];
0309   ss << "\n";
0310   ss << "\nTIB:\t\t\t" << NBadComponent[0][0][0] << "\t" << NBadComponent[0][0][1] << "\t" << NBadComponent[0][0][2]
0311      << "\t" << NBadComponent[0][0][3];
0312   ss << "\nTID:\t\t\t" << NBadComponent[1][0][0] << "\t" << NBadComponent[1][0][1] << "\t" << NBadComponent[1][0][2]
0313      << "\t" << NBadComponent[1][0][3];
0314   ss << "\nTOB:\t\t\t" << NBadComponent[2][0][0] << "\t" << NBadComponent[2][0][1] << "\t" << NBadComponent[2][0][2]
0315      << "\t" << NBadComponent[2][0][3];
0316   ss << "\nTEC:\t\t\t" << NBadComponent[3][0][0] << "\t" << NBadComponent[3][0][1] << "\t" << NBadComponent[3][0][2]
0317      << "\t" << NBadComponent[3][0][3];
0318   ss << "\n";
0319 
0320   for (int i = 1; i < 5; ++i) {
0321     ss << "\nTIB Layer " << i << " :\t\t" << NBadComponent[0][i][0] << "\t" << NBadComponent[0][i][1] << "\t"
0322        << NBadComponent[0][i][2] << "\t" << NBadComponent[0][i][3];
0323     std::stringstream binlabel;
0324     binlabel << "TIB L " << i;
0325 
0326     for (int j = 0; j < 4; j++) {
0327       ME[j]->Fill(i, NBadComponent[0][i][j]);
0328       ME[j]->setBinLabel(i, binlabel.str());
0329     }
0330   }
0331   ss << "\n";
0332   for (int i = 1; i < 4; ++i) {
0333     ss << "\nTID+ Disk " << i << " :\t\t" << NBadComponent[1][i][0] << "\t" << NBadComponent[1][i][1] << "\t"
0334        << NBadComponent[1][i][2] << "\t" << NBadComponent[1][i][3];
0335     std::stringstream binlabel;
0336     binlabel << "TID+ D " << i;
0337 
0338     for (int j = 0; j < 4; j++) {
0339       ME[j]->Fill(i + 4, NBadComponent[1][i][j]);
0340       ME[j]->setBinLabel(i + 4, binlabel.str());
0341     }
0342   }
0343   for (int i = 4; i < 7; ++i) {
0344     ss << "\nTID- Disk " << i - 3 << " :\t\t" << NBadComponent[1][i][0] << "\t" << NBadComponent[1][i][1] << "\t"
0345        << NBadComponent[1][i][2] << "\t" << NBadComponent[1][i][3];
0346     std::stringstream binlabel;
0347     binlabel << "TID- D " << i - 3;
0348 
0349     for (int j = 0; j < 4; j++) {
0350       ME[j]->Fill(i + 4, NBadComponent[1][i][j]);
0351       ME[j]->setBinLabel(i + 4, binlabel.str());
0352     }
0353   }
0354   ss << "\n";
0355   for (int i = 1; i < 7; ++i) {
0356     ss << "\nTOB Layer " << i << " :\t\t" << NBadComponent[2][i][0] << "\t" << NBadComponent[2][i][1] << "\t"
0357        << NBadComponent[2][i][2] << "\t" << NBadComponent[2][i][3];
0358     std::stringstream binlabel;
0359     binlabel << "TOB L " << i;
0360 
0361     for (int j = 0; j < 4; j++) {
0362       ME[j]->Fill(i + 10, NBadComponent[2][i][j]);
0363       ME[j]->setBinLabel(i + 10, binlabel.str());
0364     }
0365   }
0366   ss << "\n";
0367   for (int i = 1; i < 10; ++i) {
0368     ss << "\nTEC+ Disk " << i << " :\t\t" << NBadComponent[3][i][0] << "\t" << NBadComponent[3][i][1] << "\t"
0369        << NBadComponent[3][i][2] << "\t" << NBadComponent[3][i][3];
0370     std::stringstream binlabel;
0371     binlabel << "TEC+ D " << i;
0372 
0373     for (int j = 0; j < 4; j++) {
0374       ME[j]->Fill(i + 16, NBadComponent[3][i][j]);
0375       ME[j]->setBinLabel(i + 16, binlabel.str());
0376     }
0377   }
0378   for (int i = 10; i < 19; ++i) {
0379     ss << "\nTEC- Disk " << i - 9 << " :\t\t" << NBadComponent[3][i][0] << "\t" << NBadComponent[3][i][1] << "\t"
0380        << NBadComponent[3][i][2] << "\t" << NBadComponent[3][i][3];
0381     std::stringstream binlabel;
0382     binlabel << "TEC- D " << i - 9;
0383 
0384     for (int j = 0; j < 4; j++) {
0385       ME[j]->Fill(i + 16, NBadComponent[3][i][j]);
0386       ME[j]->setBinLabel(i + 16, binlabel.str());
0387     }
0388   }
0389   ss << "\n";
0390 
0391   ss << "\n----------------------------------------------------------------"
0392         "\n\t\t   Detid  \tModules Fibers "
0393         "Apvs\n---------------------------------------------------------------"
0394         "-";
0395   for (int i = 1; i < 5; ++i)
0396     ss << "\nTIB Layer " << i << " :" << ssV[0][i].str();
0397   ss << "\n";
0398   for (int i = 1; i < 4; ++i)
0399     ss << "\nTID+ Disk " << i << " :" << ssV[1][i].str();
0400   for (int i = 4; i < 7; ++i)
0401     ss << "\nTID- Disk " << i - 3 << " :" << ssV[1][i].str();
0402   ss << "\n";
0403   for (int i = 1; i < 7; ++i)
0404     ss << "\nTOB Layer " << i << " :" << ssV[2][i].str();
0405   ss << "\n";
0406   for (int i = 1; i < 10; ++i)
0407     ss << "\nTEC+ Disk " << i << " :" << ssV[3][i].str();
0408   for (int i = 10; i < 19; ++i)
0409     ss << "\nTEC- Disk " << i - 9 << " :" << ssV[3][i].str();
0410 
0411   edm::LogInfo("SiStripQualityDQM") << ss.str() << std::endl;
0412 
0413   for (int i = 0; i < 4; i++) {
0414     TCanvas c1("c1");
0415     ME[i]->getTH1()->Draw();
0416     std::string name(ME[i]->getTitle());
0417     name += ".png";
0418     c1.Print(name.c_str());
0419   }
0420 }
0421 
0422 void SiStripQualityDQM::SetBadComponents(int i, int component, SiStripQuality::BadComponent &BC) {
0423   int napv = detInfo_.getNumberOfApvsAndStripLength(BC.detid).first;
0424 
0425   ssV[i][component] << "\n\t\t " << BC.detid << " \t " << BC.BadModule << " \t " << ((BC.BadFibers) & 0x1) << " ";
0426   if (napv == 4)
0427     ssV[i][component] << "x " << ((BC.BadFibers >> 1) & 0x1);
0428 
0429   if (napv == 6)
0430     ssV[i][component] << ((BC.BadFibers >> 1) & 0x1) << " " << ((BC.BadFibers >> 2) & 0x1);
0431   ssV[i][component] << " \t " << ((BC.BadApvs) & 0x1) << " " << ((BC.BadApvs >> 1) & 0x1) << " ";
0432   if (napv == 4)
0433     ssV[i][component] << "x x " << ((BC.BadApvs >> 2) & 0x1) << " " << ((BC.BadApvs >> 3) & 0x1);
0434   if (napv == 6)
0435     ssV[i][component] << ((BC.BadApvs >> 2) & 0x1) << " " << ((BC.BadApvs >> 3) & 0x1) << " "
0436                       << ((BC.BadApvs >> 4) & 0x1) << " " << ((BC.BadApvs >> 5) & 0x1) << " ";
0437 
0438   if (BC.BadApvs) {
0439     NBadComponent[i][0][2] += ((BC.BadApvs >> 5) & 0x1) + ((BC.BadApvs >> 4) & 0x1) + ((BC.BadApvs >> 3) & 0x1) +
0440                               ((BC.BadApvs >> 2) & 0x1) + ((BC.BadApvs >> 1) & 0x1) + ((BC.BadApvs) & 0x1);
0441     NBadComponent[i][component][2] += ((BC.BadApvs >> 5) & 0x1) + ((BC.BadApvs >> 4) & 0x1) +
0442                                       ((BC.BadApvs >> 3) & 0x1) + ((BC.BadApvs >> 2) & 0x1) +
0443                                       ((BC.BadApvs >> 1) & 0x1) + ((BC.BadApvs) & 0x1);
0444     //    tkMap->fillc(BC.detid,0xff0000);
0445   }
0446   if (BC.BadFibers) {
0447     NBadComponent[i][0][1] += ((BC.BadFibers >> 2) & 0x1) + ((BC.BadFibers >> 1) & 0x1) + ((BC.BadFibers) & 0x1);
0448     NBadComponent[i][component][1] +=
0449         ((BC.BadFibers >> 2) & 0x1) + ((BC.BadFibers >> 1) & 0x1) + ((BC.BadFibers) & 0x1);
0450     //    tkMap->fillc(BC.detid,0x0000ff);
0451   }
0452   if (BC.BadModule) {
0453     NBadComponent[i][0][0]++;
0454     NBadComponent[i][component][0]++;
0455     //    tkMap->fillc(BC.detid,0x0);
0456   }
0457 }