CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
SiStripQualityStatistics.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // Package: SiStripQualityStatistics
4 // Class: SiStripQualityStatistics
5 //
13 //
14 // Original Author: Domenico GIORDANO
15 // Created: Wed Oct 3 12:11:10 CEST 2007
16 //
17 //
18 #include "CalibTracker/Records/interface/SiStripQualityRcd.h"
19 
20 #include "CalibTracker/SiStripQuality/plugins/SiStripQualityStatistics.h"
21 
22 #include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
23 #include "Geometry/Records/interface/IdealGeometryRecord.h"
24 
25 #include <iostream>
26 #include <iomanip>
27 #include <stdio.h>
28 #include <sys/time.h>
29 
30 
31 SiStripQualityStatistics::SiStripQualityStatistics( const edm::ParameterSet& iConfig ):
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(0),tkMapFullIOVs(0)
38 {
39  reader = new SiStripDetInfoFileReader(fp_.fullPath());
40 
41  tkMapFullIOVs=new TrackerMap( "BadComponents" );
42  tkhisto=0;
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 }
47 
48 void SiStripQualityStatistics::endJob(){
49 
50  std::string filename=TkMapFileName_;
51  if (filename!=""){
52  tkMapFullIOVs->save(false,0,0,filename.c_str());
53  filename.erase(filename.begin()+filename.find("."),filename.end());
54  tkMapFullIOVs->print(false,0,0,filename.c_str());
55 
56  if(saveTkHistoMap_){
57  tkhisto->save(filename+".root");
58  tkhisto->saveAsCanvas(filename+"_Canvas.root","E");
59  }
60  }
61 }
62 
63 void SiStripQualityStatistics::analyze( const edm::Event& e, const edm::EventSetup& iSetup){
64  //Retrieve tracker topology from geometry
65  edm::ESHandle<TrackerTopology> tTopoHandle;
66  iSetup.get<IdealGeometryRecord>().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!=NULL)
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 
266  std::string filetype=TkMapFileName_,filename=TkMapFileName_;
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.c_str());
273  filename.erase(filename.begin()+filename.find("."),filename.end());
274  tkMap->print(true,0,0,filename.c_str());
275  }
276 }
277 
278 
279 void SiStripQualityStatistics::SetBadComponents(int i, int component,SiStripQuality::BadComponent& BC){
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 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:39
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:41
i
int i
Definition: DBlmapReader.cc:9
TrackerTopology::tibLayer
unsigned int tibLayer(const DetId &id) const
Definition: TrackerTopology.h:122
SiStripQualityStatistics::SetBadComponents
void SetBadComponents(int, int, SiStripQuality::BadComponent &)
Definition: SiStripQualityStatistics.cc:279
SiStripDetInfoFileReader::getNumberOfApvsAndStripLength
const std::pair< unsigned short, double > getNumberOfApvsAndStripLength(uint32_t detId) const
Definition: SiStripDetInfoFileReader.cc:120
edm::EventBase::luminosityBlock
edm::LuminosityBlockNumber_t luminosityBlock() const
Definition: EventBase.h:63
NULL
#define NULL
Definition: scimark2.h:8
TrackerTopology::tidWheel
unsigned int tidWheel(const DetId &id) const
Definition: TrackerTopology.h:183
SiStripBadStrip::RegistryIterator
Registry::const_iterator RegistryIterator
Definition: SiStripBadStrip.h:56
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
TkHistoMap::saveAsCanvas
void saveAsCanvas(std::string filename, std::string options="", std::string mode="RECREATE")
Definition: TkHistoMap.cc:252
SiStripQualityStatistics::analyze
void analyze(const edm::Event &, const edm::EventSetup &)
Definition: SiStripQualityStatistics.cc:63
TkHistoMap::save
void save(std::string filename)
Definition: TkHistoMap.cc:57
TrackerTopology::tidSide
unsigned int tidSide(const DetId &id) const
Definition: TrackerTopology.h:166
SiStripDetInfoFileReader::getAllDetIds
const std::vector< uint32_t > & getAllDetIds() const
Definition: SiStripDetInfoFileReader.h:56
TkHistoMap::fill
void fill(uint32_t &detid, float value)
Definition: TkHistoMap.cc:180
TrackerMap::save
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:698
j
int j
Definition: DBlmapReader.cc:9
DetId::subdetId
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:37
TrackerMap::fillc
void fillc(int idmod, int RGBcode)
Definition: TrackerMap.h:103
SiStripQualityStatistics::ssV
std::stringstream ssV[4][19]
Definition: SiStripQualityStatistics.h:48
OfflineSiStripConfigDb_cff.SiStripDetInfoFileReader
tuple SiStripDetInfoFileReader
Definition: OfflineSiStripConfigDb_cff.py:20
DetId
Definition: DetId.h:18
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:60
SiStripBadStrip::Range
std::pair< ContainerIterator, ContainerIterator > Range
Definition: SiStripBadStrip.h:54
TrackerMap::print
void print(bool print_total=true, float minval=0., float maxval=0., std::string s="svgmap")
Definition: TrackerMap.cc:2432
lut2db_cfg.filename
tuple filename
Definition: lut2db_cfg.py:20
SiStripQualityStatistics::reader
SiStripDetInfoFileReader * reader
Definition: SiStripQualityStatistics.h:51
SiStripQualityStatistics::SiStripQualityStatistics
SiStripQualityStatistics(const edm::ParameterSet &)
Definition: SiStripQualityStatistics.cc:31
edm::FileInPath::fullPath
std::string fullPath() const
Definition: FileInPath.cc:165
relval_steps.k
list k
Definition: relval_steps.py:1469
TrackerTopology::tecWheel
unsigned int tecWheel(const DetId &id) const
Definition: TrackerTopology.h:180
edm::Timestamp::value
TimeValue_t value() const
Definition: Timestamp.h:56
edm::EventBase::time
edm::Timestamp time() const
Definition: EventBase.h:61
TrackerMap::fill
void fill(int layer, int ring, int nmod, float x)
Definition: TrackerMap.cc:2776
TrackerTopology::tobLayer
unsigned int tobLayer(const DetId &id) const
Definition: TrackerTopology.h:119
TrackerTopology::tecSide
unsigned int tecSide(const DetId &id) const
Definition: TrackerTopology.h:158
Generated for CMSSW Reference Manual by doxygen 1.8.5