78 tkhisto = std::make_unique<TkHistoMap>(tkDetMapHandle.
product(),
"BadComp",
"BadComp",-1.);
95 for(
int j=0;j<19;++j){
109 std::vector<uint32_t>::const_iterator idet=detids.begin();
110 for(;idet!=detids.end();++idet){
111 ss <<
"detid " << (*idet) <<
" IsModuleUsable " << SiStripQuality_->
IsModuleUsable((*idet)) <<
"\n";
115 LogDebug(
"SiStripQualityStatistics") << ss.str() << std::endl;
120 for (
size_t i=0;
i<BC.size();++
i){
132 ( (BC[
i].BadApvs>>2)&0
x1 )+ ( (BC[
i].BadApvs>>1)&0
x1 ) + ( (BC[
i].BadApvs)&0
x1 );
184 uint32_t detid=rp->detid;
186 int subdet=-999;
int component=-999;
206 for(
int it=0;it<sqrange.second-sqrange.first;it++){
207 unsigned int range=SiStripQuality_->
decode( *(sqrange.first+it) ).
range;
216 edm::LogError(
"SiStripQualityStatistics") <<
"PROBLEM detid " << detid <<
" value " << percentage<< std::endl;
221 tkhisto->fill(detid,percentage);
230 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";
231 ss <<
"\n-----------------\nGlobal Info\n-----------------";
232 ss <<
"\nBadComponent \t Modules \tFibers \tApvs\tStrips\n----------------------------------------------------------------";
236 ss <<
"\nTID:\t\t\t"<<NBadComponent[1][0][0]<<
"\t"<<NBadComponent[1][0][1]<<
"\t"<<NBadComponent[1][0][2]<<
"\t"<<NBadComponent[1][0][3];
237 ss <<
"\nTOB:\t\t\t"<<NBadComponent[2][0][0]<<
"\t"<<NBadComponent[2][0][1]<<
"\t"<<NBadComponent[2][0][2]<<
"\t"<<NBadComponent[2][0][3];
238 ss <<
"\nTEC:\t\t\t"<<NBadComponent[3][0][0]<<
"\t"<<NBadComponent[3][0][1]<<
"\t"<<NBadComponent[3][0][2]<<
"\t"<<NBadComponent[3][0][3];
241 for (
int i=1;
i<5;++
i)
242 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];
244 for (
int i=1;
i<4;++
i)
245 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];
246 for (
int i=4;
i<7;++
i)
247 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];
249 for (
int i=1;
i<7;++
i)
250 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];
252 for (
int i=1;
i<10;++
i)
253 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];
254 for (
int i=10;
i<19;++
i)
255 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];
258 ss <<
"\n----------------------------------------------------------------\n\t\t Detid \tModules Fibers Apvs\n----------------------------------------------------------------";
259 for (
int i=1;
i<5;++
i)
260 ss <<
"\nTIB Layer " <<
i <<
" :" <<
ssV[0][
i].
str();
262 for (
int i=1;
i<4;++
i)
263 ss <<
"\nTID+ Disk " <<
i <<
" :" <<
ssV[1][
i].
str();
264 for (
int i=4;
i<7;++
i)
265 ss <<
"\nTID- Disk " <<
i-3 <<
" :" <<
ssV[1][
i].
str();
267 for (
int i=1;
i<7;++
i)
268 ss <<
"\nTOB Layer " <<
i <<
" :" <<
ssV[2][
i].
str();
270 for (
int i=1;
i<10;++
i)
271 ss <<
"\nTEC+ Disk " <<
i <<
" :" <<
ssV[3][
i].
str();
272 for (
int i=10;
i<19;++
i)
273 ss <<
"\nTEC- Disk " <<
i-9 <<
" :" <<
ssV[3][
i].
str();
276 edm::LogInfo(
"SiStripQualityStatistics") << ss.str() << std::endl;
279 std::stringstream sRun; sRun.str(
"");
280 sRun <<
"_Run_" << std::setw(6) << std::setfill(
'0')<< e.
id().
run() << std::setw(0) ;
283 filename.insert(filename.find(
"."),sRun.str());
285 filename.erase(filename.begin()+filename.find(
"."),filename.end());
std::string TkMapFileName_
EventNumber_t event() const
const std::vector< BadComponent > & getBadComponentList() const
unsigned int tibLayer(const DetId &id) const
void SetBadComponents(int, int, SiStripQuality::BadComponent &)
TrackerMap * tkMapFullIOVs
const std::pair< unsigned short, double > getNumberOfApvsAndStripLength(uint32_t detId) const
edm::LuminosityBlockNumber_t luminosityBlock() const
int NBadComponent[4][19][4]
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 save(bool print_total=true, float minval=0., float maxval=0., std::string s="svgmap.svg", int width=1500, int height=800)
int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum) ...
void fillc(int idmod, int RGBcode)
std::stringstream ssV[4][19]
ContainerIterator getDataVectorBegin() const
std::unique_ptr< TkHistoMap > tkhisto
unsigned long long m_cacheID_
RegistryIterator getRegistryVectorBegin() const
std::pair< ContainerIterator, ContainerIterator > Range
void print(bool print_total=true, float minval=0., float maxval=0., std::string s="svgmap")
SiStripDetInfoFileReader * reader
unsigned int tecWheel(const DetId &id) const
T const * product() const
TimeValue_t value() const
edm::Timestamp time() const
void fill(int layer, int ring, int nmod, float x)
data decode(const unsigned int &value) const
unsigned int tobLayer(const DetId &id) const
unsigned int tecSide(const DetId &id) const