CMS 3D CMS Logo

List of all members | Public Member Functions | Private Member Functions | Private Attributes
SiStripQualityStatistics Class Reference

#include <CalibTracker/SiStripQuality/plugins/SiStripQualityStatistics.cc>

Inheritance diagram for SiStripQualityStatistics:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void endJob () override
 
 SiStripQualityStatistics (const edm::ParameterSet &)
 
 ~SiStripQualityStatistics () override
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
 ~EDAnalyzer () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfoconsumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBaseoperator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

void SetBadComponents (int, int, SiStripQuality::BadComponent &)
 

Private Attributes

std::string dataLabel_
 
edm::FileInPath fp_
 
unsigned long long m_cacheID_
 
int NBadComponent [4][19][4]
 
int NTkBadComponent [4]
 
SiStripDetInfoFileReaderreader
 
bool saveTkHistoMap_
 
std::stringstream ssV [4][19]
 
TkHistoMaptkhisto
 
TrackerMaptkMap
 
std::string TkMapFileName_
 
TrackerMaptkMapFullIOVs
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 24 of file SiStripQualityStatistics.h.

Constructor & Destructor Documentation

SiStripQualityStatistics::SiStripQualityStatistics ( const edm::ParameterSet iConfig)
explicit

Definition at line 31 of file SiStripQualityStatistics.cc.

References fp_, edm::FileInPath::fullPath(), reader, tkhisto, TkMapFileName_, and tkMapFullIOVs.

31  :
32  m_cacheID_(0),
33  dataLabel_(iConfig.getUntrackedParameter<std::string>("dataLabel","")),
34  TkMapFileName_(iConfig.getUntrackedParameter<std::string>("TkMapFileName","")),
35  fp_(iConfig.getUntrackedParameter<edm::FileInPath>("file",edm::FileInPath("CalibTracker/SiStripCommon/data/SiStripDetInfo.dat"))),
36  saveTkHistoMap_(iConfig.getUntrackedParameter<bool>("SaveTkHistoMap",true)),
37  tkMap(nullptr),tkMapFullIOVs(nullptr)
38 {
40 
41  tkMapFullIOVs=new TrackerMap( "BadComponents" );
42  tkhisto=nullptr;
43  if (TkMapFileName_!=""){
44  tkhisto =new TkHistoMap("BadComp","BadComp",-1.); //here the baseline (the value of the empty,not assigned bins) is put to -1 (default is zero)
45  }
46 }
T getUntrackedParameter(std::string const &, T const &) const
std::string fullPath() const
Definition: FileInPath.cc:184
SiStripDetInfoFileReader * reader
SiStripQualityStatistics::~SiStripQualityStatistics ( )
inlineoverride

Definition at line 28 of file SiStripQualityStatistics.h.

References analyze(), endJob(), and SetBadComponents().

28 {};

Member Function Documentation

void SiStripQualityStatistics::analyze ( const edm::Event e,
const edm::EventSetup iSetup 
)
override

Definition at line 63 of file SiStripQualityStatistics.cc.

References dataLabel_, SiStripBadStrip::decode(), edm::EventID::event(), corrVsCorr::filename, symbols::filetype, TkHistoMap::fill(), TrackerMap::fill(), TrackerMap::fillc(), edm::EventSetup::get(), SiStripDetInfoFileReader::getAllDetIds(), SiStripQuality::getBadComponentList(), SiStripBadStrip::getDataVectorBegin(), SiStripDetInfoFileReader::getNumberOfApvsAndStripLength(), SiStripBadStrip::getRegistryVectorBegin(), SiStripBadStrip::getRegistryVectorEnd(), mps_fire::i, edm::EventBase::id(), SiStripQuality::IsModuleUsable(), gen::k, LogDebug, edm::EventBase::luminosityBlock(), m_cacheID_, NBadComponent, NTkBadComponent, TrackerMap::print(), edm::ESHandle< T >::product(), SiStripBadStrip::data::range, reader, edm::EventID::run(), TrackerMap::save(), SetBadComponents(), ssV, harvestTrackValidationPlots::str, AlCaHLTBitMon_QueryRunRegistry::string, DetId::subdetId(), StripSubdetector::TEC, TrackerTopology::tecSide(), TrackerTopology::tecWheel(), StripSubdetector::TIB, TrackerTopology::tibLayer(), StripSubdetector::TID, TrackerTopology::tidSide(), TrackerTopology::tidWheel(), edm::EventBase::time(), tkhisto, tkMap, TkMapFileName_, tkMapFullIOVs, StripSubdetector::TOB, TrackerTopology::tobLayer(), edm::Timestamp::value(), and globals_cff::x1.

Referenced by ~SiStripQualityStatistics().

63  {
64  //Retrieve tracker topology from geometry
66  iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
67  const TrackerTopology* const tTopo = tTopoHandle.product();
68 
69  unsigned long long cacheID = iSetup.get<SiStripQualityRcd>().cacheIdentifier();
70 
71  std::stringstream ss;
72 
73  if (m_cacheID_ == cacheID)
74  return;
75 
76  m_cacheID_ = cacheID;
77 
78  edm::ESHandle<SiStripQuality> SiStripQuality_;
79  iSetup.get<SiStripQualityRcd>().get(dataLabel_,SiStripQuality_);
80 
81  for(int i=0;i<4;++i){
82  NTkBadComponent[i]=0;
83  for(int j=0;j<19;++j){
84  ssV[i][j].str("");
85  for(int k=0;k<4;++k)
86  NBadComponent[i][j][k]=0;
87  }
88  }
89 
90  if (tkMap)
91  delete tkMap;
92  tkMap=new TrackerMap( "BadComponents" );
93 
94 
95  ss.str("");
96  std::vector<uint32_t> detids=reader->getAllDetIds();
97  std::vector<uint32_t>::const_iterator idet=detids.begin();
98  for(;idet!=detids.end();++idet){
99  ss << "detid " << (*idet) << " IsModuleUsable " << SiStripQuality_->IsModuleUsable((*idet)) << "\n";
100  if (SiStripQuality_->IsModuleUsable((*idet)))
101  tkMap->fillc(*idet,0x00ff00);
102  }
103  LogDebug("SiStripQualityStatistics") << ss.str() << std::endl;
104 
105 
106  std::vector<SiStripQuality::BadComponent> BC = SiStripQuality_->getBadComponentList();
107 
108  for (size_t i=0;i<BC.size();++i){
109 
110  //&&&&&&&&&&&&&
111  //Full Tk
112  //&&&&&&&&&&&&&
113 
114  if (BC[i].BadModule)
115  NTkBadComponent[0]++;
116  if (BC[i].BadFibers)
117  NTkBadComponent[1]+= ( (BC[i].BadFibers>>2)&0x1 )+ ( (BC[i].BadFibers>>1)&0x1 ) + ( (BC[i].BadFibers)&0x1 );
118  if (BC[i].BadApvs)
119  NTkBadComponent[2]+= ( (BC[i].BadApvs>>5)&0x1 )+ ( (BC[i].BadApvs>>4)&0x1 ) + ( (BC[i].BadApvs>>3)&0x1 ) +
120  ( (BC[i].BadApvs>>2)&0x1 )+ ( (BC[i].BadApvs>>1)&0x1 ) + ( (BC[i].BadApvs)&0x1 );
121 
122  //&&&&&&&&&&&&&&&&&
123  //Single SubSyste
124  //&&&&&&&&&&&&&&&&&
125  int component;
126  DetId detectorId=DetId(BC[i].detid);
127  int subDet = detectorId.subdetId();
128  if ( subDet == StripSubdetector::TIB ){
129  //&&&&&&&&&&&&&&&&&
130  //TIB
131  //&&&&&&&&&&&&&&&&&
132 
133  component=tTopo->tibLayer(BC[i].detid);
134  SetBadComponents(0, component, BC[i]);
135 
136  } else if ( subDet == StripSubdetector::TID ) {
137  //&&&&&&&&&&&&&&&&&
138  //TID
139  //&&&&&&&&&&&&&&&&&
140 
141  component=tTopo->tidSide(BC[i].detid)==2?tTopo->tidWheel(BC[i].detid):tTopo->tidWheel(BC[i].detid)+3;
142  SetBadComponents(1, component, BC[i]);
143 
144  } else if ( subDet == StripSubdetector::TOB ) {
145  //&&&&&&&&&&&&&&&&&
146  //TOB
147  //&&&&&&&&&&&&&&&&&
148 
149  component=tTopo->tobLayer(BC[i].detid);
150  SetBadComponents(2, component, BC[i]);
151 
152  } else if ( subDet == StripSubdetector::TEC ) {
153  //&&&&&&&&&&&&&&&&&
154  //TEC
155  //&&&&&&&&&&&&&&&&&
156 
157  component=tTopo->tecSide(BC[i].detid)==2?tTopo->tecWheel(BC[i].detid):tTopo->tecWheel(BC[i].detid)+9;
158  SetBadComponents(3, component, BC[i]);
159 
160  }
161  }
162 
163  //&&&&&&&&&&&&&&&&&&
164  // Single Strip Info
165  //&&&&&&&&&&&&&&&&&&
166  float percentage=0;
167 
168  SiStripQuality::RegistryIterator rbegin = SiStripQuality_->getRegistryVectorBegin();
169  SiStripQuality::RegistryIterator rend = SiStripQuality_->getRegistryVectorEnd();
170 
171  for (SiStripBadStrip::RegistryIterator rp=rbegin; rp != rend; ++rp) {
172  uint32_t detid=rp->detid;
173 
174  int subdet=-999; int component=-999;
175  DetId detectorId=DetId(detid);
176  int subDet = detectorId.subdetId();
177  if ( subDet == StripSubdetector::TIB ){
178  subdet=0;
179  component=tTopo->tibLayer(detid);
180  } else if ( subDet == StripSubdetector::TID ) {
181  subdet=1;
182  component=tTopo->tidSide(detid)==2?tTopo->tidWheel(detid):tTopo->tidWheel(detid)+3;
183  } else if ( subDet == StripSubdetector::TOB ) {
184  subdet=2;
185  component=tTopo->tobLayer(detid);
186  } else if ( subDet == StripSubdetector::TEC ) {
187  subdet=3;
188  component=tTopo->tecSide(detid)==2?tTopo->tecWheel(detid):tTopo->tecWheel(detid)+9;
189  }
190 
191  SiStripQuality::Range sqrange = SiStripQuality::Range( SiStripQuality_->getDataVectorBegin()+rp->ibegin , SiStripQuality_->getDataVectorBegin()+rp->iend );
192 
193  percentage=0;
194  for(int it=0;it<sqrange.second-sqrange.first;it++){
195  unsigned int range=SiStripQuality_->decode( *(sqrange.first+it) ).range;
196  NTkBadComponent[3]+=range;
197  NBadComponent[subdet][0][3]+=range;
198  NBadComponent[subdet][component][3]+=range;
199  percentage+=range;
200  }
201  if(percentage!=0)
202  percentage/=128.*reader->getNumberOfApvsAndStripLength(detid).first;
203  if(percentage>1)
204  edm::LogError("SiStripQualityStatistics") << "PROBLEM detid " << detid << " value " << percentage<< std::endl;
205 
206  //------- Global Statistics on percentage of bad components along the IOVs ------//
207  tkMapFullIOVs->fill(detid,percentage);
208  if(tkhisto!=nullptr)
209  tkhisto->fill(detid,percentage);
210  }
211 
212 
213  //&&&&&&&&&&&&&&&&&&
214  // printout
215  //&&&&&&&&&&&&&&&&&&
216 
217  ss.str("");
218  ss << "\n-----------------\nNew IOV starting from run " << e.id().run() << " event " << e.id().event() << " lumiBlock " << e.luminosityBlock() << " time " << e.time().value() << " chacheID " << m_cacheID_ << "\n-----------------\n";
219  ss << "\n-----------------\nGlobal Info\n-----------------";
220  ss << "\nBadComponent \t Modules \tFibers \tApvs\tStrips\n----------------------------------------------------------------";
221  ss << "\nTracker:\t\t"<<NTkBadComponent[0]<<"\t"<<NTkBadComponent[1]<<"\t"<<NTkBadComponent[2]<<"\t"<<NTkBadComponent[3];
222  ss<< "\n";
223  ss << "\nTIB:\t\t\t"<<NBadComponent[0][0][0]<<"\t"<<NBadComponent[0][0][1]<<"\t"<<NBadComponent[0][0][2]<<"\t"<<NBadComponent[0][0][3];
224  ss << "\nTID:\t\t\t"<<NBadComponent[1][0][0]<<"\t"<<NBadComponent[1][0][1]<<"\t"<<NBadComponent[1][0][2]<<"\t"<<NBadComponent[1][0][3];
225  ss << "\nTOB:\t\t\t"<<NBadComponent[2][0][0]<<"\t"<<NBadComponent[2][0][1]<<"\t"<<NBadComponent[2][0][2]<<"\t"<<NBadComponent[2][0][3];
226  ss << "\nTEC:\t\t\t"<<NBadComponent[3][0][0]<<"\t"<<NBadComponent[3][0][1]<<"\t"<<NBadComponent[3][0][2]<<"\t"<<NBadComponent[3][0][3];
227  ss << "\n";
228 
229  for (int i=1;i<5;++i)
230  ss << "\nTIB Layer " << i << " :\t\t"<<NBadComponent[0][i][0]<<"\t"<<NBadComponent[0][i][1]<<"\t"<<NBadComponent[0][i][2]<<"\t"<<NBadComponent[0][i][3];
231  ss << "\n";
232  for (int i=1;i<4;++i)
233  ss << "\nTID+ Disk " << i << " :\t\t"<<NBadComponent[1][i][0]<<"\t"<<NBadComponent[1][i][1]<<"\t"<<NBadComponent[1][i][2]<<"\t"<<NBadComponent[1][i][3];
234  for (int i=4;i<7;++i)
235  ss << "\nTID- Disk " << i-3 << " :\t\t"<<NBadComponent[1][i][0]<<"\t"<<NBadComponent[1][i][1]<<"\t"<<NBadComponent[1][i][2]<<"\t"<<NBadComponent[1][i][3];
236  ss << "\n";
237  for (int i=1;i<7;++i)
238  ss << "\nTOB Layer " << i << " :\t\t"<<NBadComponent[2][i][0]<<"\t"<<NBadComponent[2][i][1]<<"\t"<<NBadComponent[2][i][2]<<"\t"<<NBadComponent[2][i][3];
239  ss << "\n";
240  for (int i=1;i<10;++i)
241  ss << "\nTEC+ Disk " << i << " :\t\t"<<NBadComponent[3][i][0]<<"\t"<<NBadComponent[3][i][1]<<"\t"<<NBadComponent[3][i][2]<<"\t"<<NBadComponent[3][i][3];
242  for (int i=10;i<19;++i)
243  ss << "\nTEC- Disk " << i-9 << " :\t\t"<<NBadComponent[3][i][0]<<"\t"<<NBadComponent[3][i][1]<<"\t"<<NBadComponent[3][i][2]<<"\t"<<NBadComponent[3][i][3];
244  ss<< "\n";
245 
246  ss << "\n----------------------------------------------------------------\n\t\t Detid \tModules Fibers Apvs\n----------------------------------------------------------------";
247  for (int i=1;i<5;++i)
248  ss << "\nTIB Layer " << i << " :" << ssV[0][i].str();
249  ss << "\n";
250  for (int i=1;i<4;++i)
251  ss << "\nTID+ Disk " << i << " :" << ssV[1][i].str();
252  for (int i=4;i<7;++i)
253  ss << "\nTID- Disk " << i-3 << " :" << ssV[1][i].str();
254  ss << "\n";
255  for (int i=1;i<7;++i)
256  ss << "\nTOB Layer " << i << " :" << ssV[2][i].str();
257  ss << "\n";
258  for (int i=1;i<10;++i)
259  ss << "\nTEC+ Disk " << i << " :" << ssV[3][i].str();
260  for (int i=10;i<19;++i)
261  ss << "\nTEC- Disk " << i-9 << " :" << ssV[3][i].str();
262 
263 
264  edm::LogInfo("SiStripQualityStatistics") << ss.str() << std::endl;
265 
267  std::stringstream sRun; sRun.str("");
268  sRun << "_Run_" << std::setw(6) << std::setfill('0')<< e.id().run() << std::setw(0) ;
269 
270  if (filename!=""){
271  filename.insert(filename.find("."),sRun.str());
272  tkMap->save(true,0,0,filename);
273  filename.erase(filename.begin()+filename.find("."),filename.end());
274  tkMap->print(true,0,0,filename);
275  }
276 }
#define LogDebug(id)
unsigned short range
RunNumber_t run() const
Definition: EventID.h:39
EventNumber_t event() const
Definition: EventID.h:41
const std::vector< BadComponent > & getBadComponentList() const
unsigned int tibLayer(const DetId &id) const
void SetBadComponents(int, int, SiStripQuality::BadComponent &)
const std::pair< unsigned short, double > getNumberOfApvsAndStripLength(uint32_t detId) const
edm::LuminosityBlockNumber_t luminosityBlock() const
Definition: EventBase.h:63
unsigned int tidWheel(const DetId &id) const
Registry::const_iterator RegistryIterator
RegistryIterator getRegistryVectorEnd() const
unsigned int tidSide(const DetId &id) const
const std::vector< uint32_t > & getAllDetIds() const
bool IsModuleUsable(const uint32_t &detid) const
void fill(uint32_t &detid, float value)
Definition: TkHistoMap.cc:180
void save(bool print_total=true, float minval=0., float maxval=0., std::string s="svgmap.svg", int width=1500, int height=800)
Definition: TrackerMap.cc:699
filetype
Definition: symbols.py:48
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:37
void fillc(int idmod, int RGBcode)
Definition: TrackerMap.h:109
std::stringstream ssV[4][19]
int k[5][pyjets_maxn]
ContainerIterator getDataVectorBegin() const
Definition: DetId.h:18
const T & get() const
Definition: EventSetup.h:55
RegistryIterator getRegistryVectorBegin() const
edm::EventID id() const
Definition: EventBase.h:60
std::pair< ContainerIterator, ContainerIterator > Range
void print(bool print_total=true, float minval=0., float maxval=0., std::string s="svgmap")
Definition: TrackerMap.cc:2442
SiStripDetInfoFileReader * reader
unsigned int tecWheel(const DetId &id) const
T const * product() const
Definition: ESHandle.h:86
TimeValue_t value() const
Definition: Timestamp.h:56
edm::Timestamp time() const
Definition: EventBase.h:61
void fill(int layer, int ring, int nmod, float x)
Definition: TrackerMap.cc:2786
data decode(const unsigned int &value) const
unsigned int tobLayer(const DetId &id) const
unsigned int tecSide(const DetId &id) const
void SiStripQualityStatistics::endJob ( void  )
overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 48 of file SiStripQualityStatistics.cc.

References corrVsCorr::filename, TrackerMap::print(), TkHistoMap::save(), TrackerMap::save(), TkHistoMap::saveAsCanvas(), saveTkHistoMap_, AlCaHLTBitMon_QueryRunRegistry::string, tkhisto, TkMapFileName_, and tkMapFullIOVs.

Referenced by o2olib.O2ORunMgr::executeJob(), and ~SiStripQualityStatistics().

48  {
49 
51  if (filename!=""){
52  tkMapFullIOVs->save(false,0,0,filename);
53  filename.erase(filename.begin()+filename.find("."),filename.end());
54  tkMapFullIOVs->print(false,0,0,filename);
55 
56  if(saveTkHistoMap_){
57  tkhisto->save(filename+".root");
58  tkhisto->saveAsCanvas(filename+"_Canvas.root","E");
59  }
60  }
61 }
void saveAsCanvas(std::string filename, std::string options="", std::string mode="RECREATE")
Definition: TkHistoMap.cc:252
void save(std::string filename)
Definition: TkHistoMap.cc:57
void save(bool print_total=true, float minval=0., float maxval=0., std::string s="svgmap.svg", int width=1500, int height=800)
Definition: TrackerMap.cc:699
void print(bool print_total=true, float minval=0., float maxval=0., std::string s="svgmap")
Definition: TrackerMap.cc:2442
void SiStripQualityStatistics::SetBadComponents ( int  i,
int  component,
SiStripQuality::BadComponent BC 
)
private

Definition at line 279 of file SiStripQualityStatistics.cc.

References SiStripQuality::BadComponent::BadApvs, SiStripQuality::BadComponent::BadFibers, SiStripQuality::BadComponent::BadModule, SiStripQuality::BadComponent::detid, TrackerMap::fillc(), SiStripDetInfoFileReader::getNumberOfApvsAndStripLength(), mps_fire::i, NBadComponent, reader, ssV, tkMap, and globals_cff::x1.

Referenced by analyze(), and ~SiStripQualityStatistics().

279  {
280 
281  int napv=reader->getNumberOfApvsAndStripLength(BC.detid).first;
282 
283  ssV[i][component] << "\n\t\t "
284  << BC.detid
285  << " \t " << BC.BadModule << " \t "
286  << ( (BC.BadFibers)&0x1 ) << " ";
287  if (napv==4)
288  ssV[i][component] << "x " <<( (BC.BadFibers>>1)&0x1 );
289 
290  if (napv==6)
291  ssV[i][component] << ( (BC.BadFibers>>1)&0x1 ) << " "
292  << ( (BC.BadFibers>>2)&0x1 );
293  ssV[i][component] << " \t "
294  << ( (BC.BadApvs)&0x1 ) << " "
295  << ( (BC.BadApvs>>1)&0x1 ) << " ";
296  if (napv==4)
297  ssV[i][component] << "x x " << ( (BC.BadApvs>>2)&0x1 ) << " "
298  << ( (BC.BadApvs>>3)&0x1 );
299  if (napv==6)
300  ssV[i][component] << ( (BC.BadApvs>>2)&0x1 ) << " "
301  << ( (BC.BadApvs>>3)&0x1 ) << " "
302  << ( (BC.BadApvs>>4)&0x1 ) << " "
303  << ( (BC.BadApvs>>5)&0x1 ) << " ";
304 
305  if (BC.BadApvs){
306  NBadComponent[i][0][2]+= ( (BC.BadApvs>>5)&0x1 )+ ( (BC.BadApvs>>4)&0x1 ) + ( (BC.BadApvs>>3)&0x1 ) +
307  ( (BC.BadApvs>>2)&0x1 )+ ( (BC.BadApvs>>1)&0x1 ) + ( (BC.BadApvs)&0x1 );
308  NBadComponent[i][component][2]+= ( (BC.BadApvs>>5)&0x1 )+ ( (BC.BadApvs>>4)&0x1 ) + ( (BC.BadApvs>>3)&0x1 ) +
309  ( (BC.BadApvs>>2)&0x1 )+ ( (BC.BadApvs>>1)&0x1 ) + ( (BC.BadApvs)&0x1 );
310  tkMap->fillc(BC.detid,0xff0000);
311  }
312  if (BC.BadFibers){
313  NBadComponent[i][0][1]+= ( (BC.BadFibers>>2)&0x1 )+ ( (BC.BadFibers>>1)&0x1 ) + ( (BC.BadFibers)&0x1 );
314  NBadComponent[i][component][1]+= ( (BC.BadFibers>>2)&0x1 )+ ( (BC.BadFibers>>1)&0x1 ) + ( (BC.BadFibers)&0x1 );
315  tkMap->fillc(BC.detid,0x0000ff);
316  }
317  if (BC.BadModule){
318  NBadComponent[i][0][0]++;
319  NBadComponent[i][component][0]++;
320  tkMap->fillc(BC.detid,0x0);
321  }
322 }
const std::pair< unsigned short, double > getNumberOfApvsAndStripLength(uint32_t detId) const
void fillc(int idmod, int RGBcode)
Definition: TrackerMap.h:109
std::stringstream ssV[4][19]
SiStripDetInfoFileReader * reader

Member Data Documentation

std::string SiStripQualityStatistics::dataLabel_
private

Definition at line 38 of file SiStripQualityStatistics.h.

Referenced by analyze().

edm::FileInPath SiStripQualityStatistics::fp_
private

Definition at line 40 of file SiStripQualityStatistics.h.

Referenced by SiStripQualityStatistics().

unsigned long long SiStripQualityStatistics::m_cacheID_
private

Definition at line 37 of file SiStripQualityStatistics.h.

Referenced by analyze().

int SiStripQualityStatistics::NBadComponent[4][19][4]
private

Definition at line 44 of file SiStripQualityStatistics.h.

Referenced by analyze(), and SetBadComponents().

int SiStripQualityStatistics::NTkBadComponent[4]
private

Definition at line 43 of file SiStripQualityStatistics.h.

Referenced by analyze().

SiStripDetInfoFileReader* SiStripQualityStatistics::reader
private

Definition at line 51 of file SiStripQualityStatistics.h.

Referenced by analyze(), SetBadComponents(), and SiStripQualityStatistics().

bool SiStripQualityStatistics::saveTkHistoMap_
private

Definition at line 41 of file SiStripQualityStatistics.h.

Referenced by endJob().

std::stringstream SiStripQualityStatistics::ssV[4][19]
private

Definition at line 48 of file SiStripQualityStatistics.h.

Referenced by analyze(), and SetBadComponents().

TkHistoMap* SiStripQualityStatistics::tkhisto
private

Definition at line 52 of file SiStripQualityStatistics.h.

Referenced by analyze(), endJob(), and SiStripQualityStatistics().

TrackerMap* SiStripQualityStatistics::tkMap
private

Definition at line 50 of file SiStripQualityStatistics.h.

Referenced by analyze(), and SetBadComponents().

std::string SiStripQualityStatistics::TkMapFileName_
private

Definition at line 39 of file SiStripQualityStatistics.h.

Referenced by analyze(), endJob(), and SiStripQualityStatistics().

TrackerMap * SiStripQualityStatistics::tkMapFullIOVs
private

Definition at line 50 of file SiStripQualityStatistics.h.

Referenced by analyze(), endJob(), and SiStripQualityStatistics().