CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy.cc
Go to the documentation of this file.
1 #include "CalibTracker/SiStripQuality/interface/SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy.h"
2 
3 #include "FWCore/ServiceRegistry/interface/Service.h"
4 #include "FWCore/MessageLogger/interface/MessageLogger.h"
5 #include "DataFormats/DetId/interface/DetId.h"
6 #include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
7 #include "Geometry/Records/interface/IdealGeometryRecord.h"
8 #include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"
9 #include "Geometry/CommonTopologies/interface/StripTopology.h"
10 #include "CalibFormats/SiStripObjects/interface/SiStripQuality.h"
11 
12 
13 SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy(const edm::ParameterSet& iConfig, const TrackerTopology* theTopo):
14  lowoccupancy_(0),
15  highoccupancy_(100),
16  absolutelow_(0),
17  numberiterations_(2),
18  Nevents_(0),
19  absolute_occupancy_(0),
20  OutFileName_("Occupancy.root"),
21  DQMOutfileName_("DQMOutput"),
22  UseInputDB_(iConfig.getUntrackedParameter<bool>("UseInputDB",false)),
23  tTopo(theTopo)
24  {
25  minNevents_=Nevents_*absolute_occupancy_;
26  }
27 
28 SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::~SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy(){
29  LogTrace("SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy")<<"[SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::~SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy] "<<std::endl;
30 }
31 
32 void SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::extractBadAPVSandStrips(SiStripQuality* siStripQuality,HistoMap& DM, edm::ESHandle<SiStripQuality>& inSiStripQuality){
33 
34  LogTrace("SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy")<<"[SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::extractBadAPVs] "<<std::endl;
35 
36  if (WriteOutputFile_== true)
37  {
38  f = new TFile(OutFileName_.c_str(),"RECREATE");
39  f->cd();
40 
41  apvtree = new TTree("moduleOccupancy","tree");
42 
43  apvtree->Branch("DetRawId", &detrawid, "DetRawId/I");
44  apvtree->Branch("SubDetId", &subdetid, "SubDetId/I");
45  apvtree->Branch("Layer_Ring", &layer_ring, "Layer_Ring/I");
46  apvtree->Branch("Disc", &disc, "Disc/I");
47  apvtree->Branch("IsBack", &isback, "IsBack/I");
48  apvtree->Branch("IsExternalString", &isexternalstring, "IsExternalString/I");
49  apvtree->Branch("IsZMinusSide", &iszminusside, "IsZMinusSide/I");
50  apvtree->Branch("RodStringPetal", &rodstringpetal, "RodStringPetal/I");
51  apvtree->Branch("IsStereo", &isstereo, "IsStereo/I");
52  apvtree->Branch("ModuleNumber", &module_number, "ModuleNumber/I");
53  apvtree->Branch("NumberOfStrips", &number_strips, "NumberOfStrips/I");
54  apvtree->Branch("APVGlobalPositionX", &global_position_x, "APVGlobalPositionX/F");
55  apvtree->Branch("APVGlobalPositionY", &global_position_y, "APVGlobalPositionY/F");
56  apvtree->Branch("APVGlobalPositionZ", &global_position_z, "APVGlobalPositionZ/F");
57  apvtree->Branch("APVNumber", &apv_number, "APVNumber/I");
58  apvtree->Branch("APVAbsoluteOccupancy", &apvAbsoluteOccupancy, "apvAbsoluteOccupancy/I");
59  apvtree->Branch("APVMedianOccupancy", &apvMedianOccupancy, "apvMedianOccupancy/D");
60  apvtree->Branch("IsBad", &isBad, "IsBad/I");
61 
62  striptree = new TTree("stripOccupancy","tree");
63 
64  striptree->Branch("DetRawId", &detrawid, "DetRawId/I");
65  striptree->Branch("SubDetId", &subdetid, "SubDetId/I");
66  striptree->Branch("Layer_Ring", &layer_ring, "Layer_Ring/I");
67  striptree->Branch("Disc", &disc, "Disc/I");
68  striptree->Branch("IsBack", &isback, "IsBack/I");
69  striptree->Branch("IsExternalString", &isexternalstring, "IsExternalString/I");
70  striptree->Branch("IsZMinusSide", &iszminusside, "IsZMinusSide/I");
71  striptree->Branch("RodStringPetal", &rodstringpetal, "RodStringPetal/I");
72  striptree->Branch("IsStereo", &isstereo, "IsStereo/I");
73  striptree->Branch("ModulePosition", &module_number, "ModulePosition/I");
74  striptree->Branch("NumberOfStrips", &number_strips, "NumberOfStrips/I");
75  striptree->Branch("StripNumber", &strip_number, "StripNumber/I");
76  striptree->Branch("APVChannel", &apv_channel, "APVChannel/I");
77  striptree->Branch("StripGlobalPositionX", &strip_global_position_x, "StripGlobalPositionX/F");
78  striptree->Branch("StripGlobalPositionY", &strip_global_position_y, "StripGlobalPositionY/F");
79  striptree->Branch("StripGlobalPositionZ", &strip_global_position_z, "StripGlobalPositionZ/F");
80  striptree->Branch("IsHot", &isHot, "IsHot/I");
81  striptree->Branch("HotStripsPerAPV", &hotStripsPerAPV, "HotStripsPerAPV/I");
82  striptree->Branch("HotStripsPerModule", &hotStripsPerModule,"HotStripsPerModule/I");
83  striptree->Branch("StripOccupancy", &singleStripOccupancy, "StripOccupancy/D");
84  striptree->Branch("StripHits", &stripHits, "StripHits/I");
85  striptree->Branch("PoissonProb", &poissonProb, "PoissonProb/D");
86  }
87 
88  HistoMap::iterator it=DM.begin();
89  HistoMap::iterator itEnd=DM.end();
90  std::vector<unsigned int> badStripList;
91  uint32_t detid;
92  for (;it!=itEnd;++it){
93 
94  Apv APV;
95 
96  for (int apv=0; apv<6; apv++)
97  {
98  APV.apvMedian[apv] = 0;
99  APV.apvabsoluteOccupancy[apv] = 0;
100  APV.NEntries[apv] = 0;
101  APV.NEmptyBins[apv] = 0;
102 
103  for (int strip=0; strip<128; strip++)
104  {
105  stripOccupancy[apv][strip] = 0;
106  stripWeight[apv][strip] = 0;
107  }
108  }
109 
110  number_strips = (int)((it->second.get())->GetNbinsX());
111  number_apvs = number_strips/128;
112  APV.numberApvs = number_apvs;
113 
114  for (int apv=0; apv<number_apvs; apv++)
115  {
116  APV.th1f[apv] = new TH1F("tmp","tmp",128,0.5,128.5);
117  int NumberEntriesPerAPV=0;
118 
119  for (int strip=0; strip<128; strip++)
120  {
121  stripOccupancy[apv][strip] = (it->second.get())->GetBinContent((apv*128)+strip+1); // Remember: Bin=0 is underflow bin!
122  stripWeight[apv][strip] = 1;
123  APV.apvabsoluteOccupancy[apv] += (it->second.get())->GetBinContent((apv*128)+strip+1); // Remember: Bin=0 is underflow bin!
124  APV.th1f[apv]->SetBinContent(strip+1,(it->second.get())->GetBinContent((apv*128)+strip+1));
125  NumberEntriesPerAPV += (int)(it->second.get())->GetBinContent((apv*128)+strip+1);
126  }
127 
128  APV.th1f[apv]->SetEntries(NumberEntriesPerAPV);
129  APV.NEntries[apv]=(int)APV.th1f[apv]->GetEntries();
130  }
131 
132  for (int apv=0; apv<number_apvs; apv++)
133  {
134  APV.apvMedian[apv] = TMath::Median(128,stripOccupancy[apv],stripWeight[apv]);
135  }
136 
137  detid=it->first;
138  DetId detectorId=DetId(detid);
139 
140  if (edm::isDebugEnabled())
141  LogTrace("SiStripBadAPV") << "Analyzing detid " << detid<< std::endl;
142 
143  detrawid = detid;
144  APV.detrawId = detrawid;
145  subdetid = detectorId.subdetId();
146 
147  switch (detectorId.subdetId())
148  {
149  case StripSubdetector::TIB :
150  layer_ring = tTopo->tibLayer(detrawid);
151  module_number = tTopo->tibModule(detrawid);
152  APV.modulePosition = module_number;
153 
154  if (layer_ring == 1) medianValues_TIB_Layer1.push_back(APV);
155  else if (layer_ring == 2) medianValues_TIB_Layer2.push_back(APV);
156  else if (layer_ring == 3) medianValues_TIB_Layer3.push_back(APV);
157  else if (layer_ring == 4) medianValues_TIB_Layer4.push_back(APV);
158  break;
159 
160  case StripSubdetector::TID :
161  layer_ring = tTopo->tidRing(detrawid);
162  disc = tTopo->tidWheel(detrawid);
163  APV.modulePosition = layer_ring;
164 
165  if (tTopo->tidIsZMinusSide(detrawid)) iszminusside = 1;
166  else iszminusside = 0;
167 
168  if (iszminusside==0)
169  {
170  if (disc==1) medianValues_TIDPlus_Disc1.push_back(APV);
171  else if (disc==2) medianValues_TIDPlus_Disc2.push_back(APV);
172  else if (disc==3) medianValues_TIDPlus_Disc3.push_back(APV);
173  }
174  else if (iszminusside==1)
175  {
176  if (disc==1) medianValues_TIDMinus_Disc1.push_back(APV);
177  else if (disc==2) medianValues_TIDMinus_Disc2.push_back(APV);
178  else if (disc==3) medianValues_TIDMinus_Disc3.push_back(APV);
179  }
180  break;
181 
182  case StripSubdetector::TOB :
183  layer_ring = tTopo->tobLayer(detrawid);
184  module_number = tTopo->tobModule(detrawid);
185  APV.modulePosition = module_number;
186 
187  if (layer_ring == 1) medianValues_TOB_Layer1.push_back(APV);
188  else if (layer_ring == 2) medianValues_TOB_Layer2.push_back(APV);
189  else if (layer_ring == 3) medianValues_TOB_Layer3.push_back(APV);
190  else if (layer_ring == 4) medianValues_TOB_Layer4.push_back(APV);
191  else if (layer_ring == 5) medianValues_TOB_Layer5.push_back(APV);
192  else if (layer_ring == 6) medianValues_TOB_Layer6.push_back(APV);
193  break;
194 
195  case StripSubdetector::TEC :
196  layer_ring = tTopo->tecRing(detrawid);
197  disc = tTopo->tecWheel(detrawid);
198  APV.modulePosition = layer_ring;
199 
200  if (tTopo->tecIsZMinusSide(detrawid)) iszminusside = 1;
201  else iszminusside = 0;
202 
203  if (iszminusside==0)
204  {
205  if (disc==1) medianValues_TECPlus_Disc1.push_back(APV);
206  else if (disc==2) medianValues_TECPlus_Disc2.push_back(APV);
207  else if (disc==3) medianValues_TECPlus_Disc3.push_back(APV);
208  else if (disc==4) medianValues_TECPlus_Disc4.push_back(APV);
209  else if (disc==5) medianValues_TECPlus_Disc5.push_back(APV);
210  else if (disc==6) medianValues_TECPlus_Disc6.push_back(APV);
211  else if (disc==7) medianValues_TECPlus_Disc7.push_back(APV);
212  else if (disc==8) medianValues_TECPlus_Disc8.push_back(APV);
213  else if (disc==9) medianValues_TECPlus_Disc9.push_back(APV);
214  }
215  else if (iszminusside==1)
216  {
217  if (disc==1) medianValues_TECMinus_Disc1.push_back(APV);
218  else if (disc==2) medianValues_TECMinus_Disc2.push_back(APV);
219  else if (disc==3) medianValues_TECMinus_Disc3.push_back(APV);
220  else if (disc==4) medianValues_TECMinus_Disc4.push_back(APV);
221  else if (disc==5) medianValues_TECMinus_Disc5.push_back(APV);
222  else if (disc==6) medianValues_TECMinus_Disc6.push_back(APV);
223  else if (disc==7) medianValues_TECMinus_Disc7.push_back(APV);
224  else if (disc==8) medianValues_TECMinus_Disc8.push_back(APV);
225  else if (disc==9) medianValues_TECMinus_Disc9.push_back(APV);
226  }
227  break;
228 
229  default :
230  std::cout << "### Detector does not belong to TIB, TID, TOB or TEC !? ###" << std::endl;
231  std::cout << "### DetRawId: " << detrawid << " ###" << std::endl;
232  }
233 
234  } // end loop on modules
235 
236  // Calculate Mean and RMS for each Layer
237  CalculateMeanAndRMS(medianValues_TIB_Layer1,MeanAndRms_TIB_Layer1,numberiterations_);
238  CalculateMeanAndRMS(medianValues_TIB_Layer2,MeanAndRms_TIB_Layer2,numberiterations_);
239  CalculateMeanAndRMS(medianValues_TIB_Layer3,MeanAndRms_TIB_Layer3,numberiterations_);
240  CalculateMeanAndRMS(medianValues_TIB_Layer4,MeanAndRms_TIB_Layer4,numberiterations_);
241 
242  CalculateMeanAndRMS(medianValues_TOB_Layer1,MeanAndRms_TOB_Layer1,numberiterations_);
243  CalculateMeanAndRMS(medianValues_TOB_Layer2,MeanAndRms_TOB_Layer2,numberiterations_);
244  CalculateMeanAndRMS(medianValues_TOB_Layer3,MeanAndRms_TOB_Layer3,numberiterations_);
245  CalculateMeanAndRMS(medianValues_TOB_Layer4,MeanAndRms_TOB_Layer4,numberiterations_);
246  CalculateMeanAndRMS(medianValues_TOB_Layer5,MeanAndRms_TOB_Layer5,numberiterations_);
247  CalculateMeanAndRMS(medianValues_TOB_Layer6,MeanAndRms_TOB_Layer6,numberiterations_);
248 
249  CalculateMeanAndRMS(medianValues_TIDPlus_Disc1,MeanAndRms_TIDPlus_Disc1,numberiterations_);
250  CalculateMeanAndRMS(medianValues_TIDPlus_Disc2,MeanAndRms_TIDPlus_Disc2,numberiterations_);
251  CalculateMeanAndRMS(medianValues_TIDPlus_Disc3,MeanAndRms_TIDPlus_Disc3,numberiterations_);
252  CalculateMeanAndRMS(medianValues_TIDMinus_Disc1,MeanAndRms_TIDMinus_Disc1,numberiterations_);
253  CalculateMeanAndRMS(medianValues_TIDMinus_Disc2,MeanAndRms_TIDMinus_Disc2,numberiterations_);
254  CalculateMeanAndRMS(medianValues_TIDMinus_Disc3,MeanAndRms_TIDMinus_Disc3,numberiterations_);
255 
256  CalculateMeanAndRMS(medianValues_TECPlus_Disc1,MeanAndRms_TECPlus_Disc1,numberiterations_);
257  CalculateMeanAndRMS(medianValues_TECPlus_Disc2,MeanAndRms_TECPlus_Disc2,numberiterations_);
258  CalculateMeanAndRMS(medianValues_TECPlus_Disc3,MeanAndRms_TECPlus_Disc3,numberiterations_);
259  CalculateMeanAndRMS(medianValues_TECPlus_Disc4,MeanAndRms_TECPlus_Disc4,numberiterations_);
260  CalculateMeanAndRMS(medianValues_TECPlus_Disc5,MeanAndRms_TECPlus_Disc5,numberiterations_);
261  CalculateMeanAndRMS(medianValues_TECPlus_Disc6,MeanAndRms_TECPlus_Disc6,numberiterations_);
262  CalculateMeanAndRMS(medianValues_TECPlus_Disc7,MeanAndRms_TECPlus_Disc7,numberiterations_);
263  CalculateMeanAndRMS(medianValues_TECPlus_Disc8,MeanAndRms_TECPlus_Disc8,numberiterations_);
264  CalculateMeanAndRMS(medianValues_TECPlus_Disc9,MeanAndRms_TECPlus_Disc9,numberiterations_);
265 
266  CalculateMeanAndRMS(medianValues_TECMinus_Disc1,MeanAndRms_TECMinus_Disc1,numberiterations_);
267  CalculateMeanAndRMS(medianValues_TECMinus_Disc2,MeanAndRms_TECMinus_Disc2,numberiterations_);
268  CalculateMeanAndRMS(medianValues_TECMinus_Disc3,MeanAndRms_TECMinus_Disc3,numberiterations_);
269  CalculateMeanAndRMS(medianValues_TECMinus_Disc4,MeanAndRms_TECMinus_Disc4,numberiterations_);
270  CalculateMeanAndRMS(medianValues_TECMinus_Disc5,MeanAndRms_TECMinus_Disc5,numberiterations_);
271  CalculateMeanAndRMS(medianValues_TECMinus_Disc6,MeanAndRms_TECMinus_Disc6,numberiterations_);
272  CalculateMeanAndRMS(medianValues_TECMinus_Disc7,MeanAndRms_TECMinus_Disc7,numberiterations_);
273  CalculateMeanAndRMS(medianValues_TECMinus_Disc8,MeanAndRms_TECMinus_Disc8,numberiterations_);
274  CalculateMeanAndRMS(medianValues_TECMinus_Disc9,MeanAndRms_TECMinus_Disc9,numberiterations_);
275 
276  pQuality=siStripQuality;
277  badStripList.clear();
278 
279  // Initialize the DQM output histograms
280  initializeDQMHistograms();
281 
282  // Analyze the Occupancy for both APVs and Strips
283  AnalyzeOccupancy(siStripQuality,medianValues_TIB_Layer1,MeanAndRms_TIB_Layer1,badStripList,inSiStripQuality);
284  AnalyzeOccupancy(siStripQuality,medianValues_TIB_Layer2,MeanAndRms_TIB_Layer2,badStripList,inSiStripQuality);
285  AnalyzeOccupancy(siStripQuality,medianValues_TIB_Layer3,MeanAndRms_TIB_Layer3,badStripList,inSiStripQuality);
286  AnalyzeOccupancy(siStripQuality,medianValues_TIB_Layer4,MeanAndRms_TIB_Layer4,badStripList,inSiStripQuality);
287 
288  AnalyzeOccupancy(siStripQuality,medianValues_TOB_Layer1,MeanAndRms_TOB_Layer1,badStripList,inSiStripQuality);
289  AnalyzeOccupancy(siStripQuality,medianValues_TOB_Layer2,MeanAndRms_TOB_Layer2,badStripList,inSiStripQuality);
290  AnalyzeOccupancy(siStripQuality,medianValues_TOB_Layer3,MeanAndRms_TOB_Layer3,badStripList,inSiStripQuality);
291  AnalyzeOccupancy(siStripQuality,medianValues_TOB_Layer4,MeanAndRms_TOB_Layer4,badStripList,inSiStripQuality);
292  AnalyzeOccupancy(siStripQuality,medianValues_TOB_Layer5,MeanAndRms_TOB_Layer5,badStripList,inSiStripQuality);
293  AnalyzeOccupancy(siStripQuality,medianValues_TOB_Layer6,MeanAndRms_TOB_Layer6,badStripList,inSiStripQuality);
294 
295  AnalyzeOccupancy(siStripQuality,medianValues_TIDPlus_Disc1,MeanAndRms_TIDPlus_Disc1,badStripList,inSiStripQuality);
296  AnalyzeOccupancy(siStripQuality,medianValues_TIDPlus_Disc2,MeanAndRms_TIDPlus_Disc2,badStripList,inSiStripQuality);
297  AnalyzeOccupancy(siStripQuality,medianValues_TIDPlus_Disc3,MeanAndRms_TIDPlus_Disc3,badStripList,inSiStripQuality);
298  AnalyzeOccupancy(siStripQuality,medianValues_TIDMinus_Disc1,MeanAndRms_TIDMinus_Disc1,badStripList,inSiStripQuality);
299  AnalyzeOccupancy(siStripQuality,medianValues_TIDMinus_Disc2,MeanAndRms_TIDMinus_Disc2,badStripList,inSiStripQuality);
300  AnalyzeOccupancy(siStripQuality,medianValues_TIDMinus_Disc3,MeanAndRms_TIDMinus_Disc3,badStripList,inSiStripQuality);
301 
302  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc1,MeanAndRms_TECPlus_Disc1,badStripList,inSiStripQuality);
303  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc2,MeanAndRms_TECPlus_Disc2,badStripList,inSiStripQuality);
304  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc3,MeanAndRms_TECPlus_Disc3,badStripList,inSiStripQuality);
305  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc4,MeanAndRms_TECPlus_Disc4,badStripList,inSiStripQuality);
306  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc5,MeanAndRms_TECPlus_Disc5,badStripList,inSiStripQuality);
307  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc6,MeanAndRms_TECPlus_Disc6,badStripList,inSiStripQuality);
308  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc7,MeanAndRms_TECPlus_Disc7,badStripList,inSiStripQuality);
309  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc8,MeanAndRms_TECPlus_Disc8,badStripList,inSiStripQuality);
310  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc9,MeanAndRms_TECPlus_Disc9,badStripList,inSiStripQuality);
311 
312  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc1,MeanAndRms_TECMinus_Disc1,badStripList,inSiStripQuality);
313  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc2,MeanAndRms_TECMinus_Disc2,badStripList,inSiStripQuality);
314  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc3,MeanAndRms_TECMinus_Disc3,badStripList,inSiStripQuality);
315  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc4,MeanAndRms_TECMinus_Disc4,badStripList,inSiStripQuality);
316  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc5,MeanAndRms_TECMinus_Disc5,badStripList,inSiStripQuality);
317  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc6,MeanAndRms_TECMinus_Disc6,badStripList,inSiStripQuality);
318  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc7,MeanAndRms_TECMinus_Disc7,badStripList,inSiStripQuality);
319  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc8,MeanAndRms_TECMinus_Disc8,badStripList,inSiStripQuality);
320  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc9,MeanAndRms_TECMinus_Disc9,badStripList,inSiStripQuality);
321 
322  siStripQuality->fillBadComponents();
323 
324  // Fill DQM histograms
325  for(unsigned int i = 0; i < subDetName.size(); i++)
326  {
327  projYDistanceVsStripNumber[i]->Add((TH1F*)distanceVsStripNumber[i]->ProjectionY());
328  pfxDistanceVsStripNumber[i]->Add(distanceVsStripNumber[i]->ProfileX(pfxDistanceVsStripNumber[i]->GetName(),1,998));
329  projYNHitsVsStripNumber[i]->Add(nHitsVsStripNumber[i]->ProjectionY());
330  projYNHitsGoodStripsVsStripNumber[i]->Add(nHitsGoodStripsVsStripNumber[i]->ProjectionY());
331  projYNHitsHotStripsVsStripNumber[i]->Add(nHitsHotStripsVsStripNumber[i]->ProjectionY());
332  projYOccupancyVsStripNumber[i]->Add(occupancyVsStripNumber[i]->ProjectionY());
333  projYOccupancyGoodStripsVsStripNumber[i]->Add(occupancyGoodStripsVsStripNumber[i]->ProjectionY());
334  projYOccupancyHotStripsVsStripNumber[i]->Add(occupancyHotStripsVsStripNumber[i]->ProjectionY());
335  pfxOccupancyVsStripNumber[i]->Add(occupancyVsStripNumber[i]->ProfileX(pfxOccupancyVsStripNumber[i]->GetName(),-8.,0.));
336  pfxOccupancyGoodStripsVsStripNumber[i]->Add(occupancyGoodStripsVsStripNumber[i]->ProfileX(pfxOccupancyGoodStripsVsStripNumber[i]->GetName(),-8.,0.));
337  pfxOccupancyHotStripsVsStripNumber[i]->Add(occupancyHotStripsVsStripNumber[i]->ProfileX(pfxOccupancyHotStripsVsStripNumber[i]->GetName(),-8.,0.));
338  projYPoissonProbVsStripNumber[i]->Add(poissonProbVsStripNumber[i]->ProjectionY());
339  projYPoissonProbGoodStripsVsStripNumber[i]->Add(poissonProbGoodStripsVsStripNumber[i]->ProjectionY());
340  projYPoissonProbHotStripsVsStripNumber[i]->Add(poissonProbHotStripsVsStripNumber[i]->ProjectionY());
341  pfxPoissonProbVsStripNumber[i]->Add(poissonProbVsStripNumber[i]->ProfileX(pfxPoissonProbVsStripNumber[i]->GetName(),-18., 0.));
342  pfxPoissonProbGoodStripsVsStripNumber[i]->Add(poissonProbGoodStripsVsStripNumber[i]->ProfileX(pfxPoissonProbGoodStripsVsStripNumber[i]->GetName(),-18., 0.));
343  pfxPoissonProbHotStripsVsStripNumber[i]->Add(poissonProbHotStripsVsStripNumber[i]->ProfileX(pfxPoissonProbHotStripsVsStripNumber[i]->GetName(),-18., 0.));
344  projXDistanceVsStripNumber[i]->Add(distanceVsStripNumber[i]->ProjectionX(projXDistanceVsStripNumber[i]->GetName(),1,998));
345 
346  }
347 
348  // Save output files
349 
350  if (WriteOutputFile_==true){
351  f->cd();
352  apvtree->Write();
353  striptree->Write();
354  f->Close();
355  }
356 
357  if (WriteDQMHistograms_==true){
358  dqmStore->cd();
359  dqmStore->save(DQMOutfileName_.c_str());
360  }
361 
362  LogTrace("SiStripBadAPV") << ss.str() << std::endl;
363 }
364 
365 
366 void SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::CalculateMeanAndRMS(const std::vector<Apv>& a, std::pair<double,double>* MeanRMS, int number_iterations)
367 {
368  Double_t tot[7], tot2[7];
369  Double_t n[7];
370 
371  Double_t Mean[7] = {0};
372  Double_t Rms[7] = {1000,1000,1000,1000,1000,1000,1000};
373 
374  int Moduleposition;
375 
376  for (int i=0; i<number_iterations; i++)
377  {
378  for (int j=0; j<7; j++)
379  {
380  n[j] = 0;
381  tot[j] = 0;
382  tot2[j] = 0;
383  }
384 
385  for (uint32_t it=0; it<a.size(); it++)
386  {
387  Moduleposition = (a[it].modulePosition)-1;
388 
389  for (int apv=0; apv<a[it].numberApvs; apv++)
390  {
391  if (i>0)
392  {
393  if (a[it].apvMedian[apv]<(Mean[Moduleposition]-3*Rms[Moduleposition]) || (a[it].apvMedian[apv]>(Mean[Moduleposition]+5*Rms[Moduleposition])))
394  {
395  continue;
396  }
397  }
398  tot[Moduleposition] += a[it].apvMedian[apv];
399  tot2[Moduleposition] += (a[it].apvMedian[apv])*(a[it].apvMedian[apv]);
400  n[Moduleposition]++;
401  }
402  }
403 
404  for (int j=0; j<7; j++)
405  {
406  if (n[j]!=0)
407  {
408  Mean[j] = tot[j]/n[j];
409  Rms[j] = TMath::Sqrt(TMath::Abs(tot2[j]/n[j] -Mean[j]*Mean[j]));
410  }
411  }
412  }
413 
414  for (int j=0; j<7; j++)
415  {
416  MeanRMS[j] = std::make_pair(Mean[j],Rms[j]);
417  }
418 
419 }
420 
421 void SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::AnalyzeOccupancy(SiStripQuality* quality, std::vector<Apv>& medianValues, std::pair<double,double>* MeanAndRms, std::vector<unsigned int>& BadStripList, edm::ESHandle<SiStripQuality>& InSiStripQuality)
422 {
423  int Moduleposition;
424  uint32_t Detid;
425 
426  for (uint32_t it=0; it<medianValues.size(); it++)
427  {
428  Moduleposition = (medianValues[it].modulePosition)-1;
429  Detid = medianValues[it].detrawId;
430 
431  setBasicTreeParameters(Detid);
432 
433  DetId DetectorId=DetId(Detid);
434  const StripGeomDetUnit* TheStripDet = dynamic_cast<const StripGeomDetUnit*>( (TkGeom->idToDet(DetectorId)) );
435  const StripTopology* TheStripTopol = dynamic_cast<const StripTopology*>( &(TheStripDet->specificTopology()) );
436 
437  //Analyze the occupancies
438  hotstripspermodule = 0;
439  vHotStripsInModule.clear();
440 
441  for (int apv=0; apv<medianValues[it].numberApvs; apv++)
442  {
443  double logMedianOccupancy = -1;
444  double logAbsoluteOccupancy = -1;
445 
446  for (int i=0; i<128; i++)
447  {
448  ishot[i] = 0;
449  stripoccupancy[i] = 0;
450  striphits[i] = 0;
451  poissonprob[i] = 0;
452  }
453 
454  number_strips = (medianValues[it].numberApvs)*128;
455  apv_number = apv+1;
456  apvMedianOccupancy = medianValues[it].apvMedian[apv];
457  apvAbsoluteOccupancy = medianValues[it].apvabsoluteOccupancy[apv];
458  isBad = 0;
459  hotstripsperapv[apv] = 0;
460 
461  LocalPoint pos_apv_local = TheStripTopol->localPosition((apv*128));
462  GlobalPoint pos_apv_global = (TkGeom->idToDet(DetectorId))->surface().toGlobal(pos_apv_local);
463 
464  global_position_x = pos_apv_global.x();
465  global_position_y = pos_apv_global.y();
466  global_position_z = pos_apv_global.z();
467 
468  if (apvMedianOccupancy>0) logMedianOccupancy = log10(apvMedianOccupancy);
469  if (apvAbsoluteOccupancy>0) logAbsoluteOccupancy = log10(apvAbsoluteOccupancy);
470 
471  //Fill the DQM histograms
472  unsigned int layer = 0;
473  if(subdetid==3 || subdetid==5)
474  layer=layer_ring;
475  else
476  layer=disc;
477 
478  // Fill histograms for all the tracker
479  medianVsAbsoluteOccupancy[0][0]->Fill(logAbsoluteOccupancy,logMedianOccupancy);
480  medianOccupancy[0][0]->Fill(logMedianOccupancy);
481  absoluteOccupancy[0][0]->Fill(logAbsoluteOccupancy);
482  // Fill summary histograms for each subdetector
483  medianVsAbsoluteOccupancy[subdetid-2][0]->Fill(logAbsoluteOccupancy,logMedianOccupancy);
484  medianOccupancy[subdetid-2][0]->Fill(logMedianOccupancy);
485  absoluteOccupancy[subdetid-2][0]->Fill(logAbsoluteOccupancy);
486  // Fill histograms for each layer/disk
487  medianVsAbsoluteOccupancy[subdetid-2][layer]->Fill(logAbsoluteOccupancy,logMedianOccupancy);
488  medianOccupancy[subdetid-2][layer]->Fill(logMedianOccupancy);
489  absoluteOccupancy[subdetid-2][layer]->Fill(logAbsoluteOccupancy);
490 
491  if(UseInputDB_)
492  {
493  if(InSiStripQuality->IsApvBad(Detid,apv) )
494  {
495  if (WriteOutputFile_==true)
496  {
497  apvtree->Fill();
498  for (int strip=0; strip<128; strip++)
499  {
500  strip_number = (apv*128)+strip+1;
501  apv_channel = apv+1;
502  isHot = ishot[strip];
503  singleStripOccupancy = stripoccupancy[strip];
504  stripHits = striphits[strip];
505  poissonProb = poissonprob[strip];
506 
507  hotStripsPerModule = hotstripspermodule;
508  hotStripsPerAPV = hotstripsperapv[apv];
509 
510  LocalPoint pos_strip_local = TheStripTopol->localPosition(strip);
511  GlobalPoint pos_strip_global = (TkGeom->idToDet(DetectorId))->surface().toGlobal(pos_strip_local);
512 
513  strip_global_position_x = pos_strip_global.x();
514  strip_global_position_y = pos_strip_global.y();
515  strip_global_position_z = pos_strip_global.z();
516  striptree->Fill();
517 
518  // Fill the strip DQM Plots
519  fillStripDQMHistograms();
520  }
521 
522  if(vHotStripsInModule.size()==1)
523  {
524  distance = 999;
525  distanceVsStripNumber[0]->Fill(vHotStripsInModule[0], distance);
526  distanceVsStripNumber[subdetid-2]->Fill(vHotStripsInModule[0], distance);
527  }
528  else if(vHotStripsInModule.size()>1)
529  {
530  for(unsigned int iVec = 0; iVec != vHotStripsInModule.size(); iVec++)
531  {
532  if(iVec==0)
533  distance = vHotStripsInModule[1] - vHotStripsInModule[0];
534  else if(iVec==vHotStripsInModule.size()-1)
535  {
536  distance = vHotStripsInModule[vHotStripsInModule.size()-1] - vHotStripsInModule[vHotStripsInModule.size() -2];
537  }
538  else if(vHotStripsInModule.size()>2)
539  {
540  distanceR = vHotStripsInModule[iVec + 1] - vHotStripsInModule[iVec];
541  distanceL = vHotStripsInModule[iVec] - vHotStripsInModule[iVec - 1];
542  distance = distanceL>distanceR?distanceR:distanceL;
543  }
544  else
545  {
546  std::cout << "ERROR! distance is never computed!!!\n";
547  }
548  distanceVsStripNumber[0]->Fill(vHotStripsInModule[iVec], distance);
549  distanceVsStripNumber[subdetid-2]->Fill(vHotStripsInModule[iVec], distance);
550  }
551  }
552 
553  }
554  continue;//if the apv is already flagged as bad, continue.
555  }
556  }
557 
558  if (medianValues[it].apvMedian[apv] > minNevents_)
559  {
560  if ((medianValues[it].apvMedian[apv]>(MeanAndRms[Moduleposition].first+highoccupancy_*MeanAndRms[Moduleposition].second)) && (medianValues[it].apvMedian[apv]>absolutelow_))
561  {
562  BadStripList.push_back(pQuality->encode((apv*128),128,0));
563  isBad = 1;
564  }
565  }
566  else if (medianValues[it].apvMedian[apv]<(MeanAndRms[Moduleposition].first-lowoccupancy_*MeanAndRms[Moduleposition].second) && (MeanAndRms[Moduleposition].first>2 || medianValues[it].apvabsoluteOccupancy[apv]==0))
567  {
568  BadStripList.push_back(pQuality->encode((apv*128),128,0));
569  isBad = 1;
570  }
571 
572  if (isBad!=1)
573  {
574  iterativeSearch(medianValues[it],BadStripList,apv);
575  }
576 
577  if (WriteOutputFile_==true)
578  {
579  apvtree->Fill();
580  for (int strip=0; strip<128; strip++)
581  {
582  strip_number = (apv*128)+strip+1;
583  apv_channel = apv+1;
584  isHot = ishot[strip];
585  singleStripOccupancy = stripoccupancy[strip];
586  stripHits = striphits[strip];
587  poissonProb = poissonprob[strip];
588 
589  hotStripsPerModule = hotstripspermodule;
590  hotStripsPerAPV = hotstripsperapv[apv];
591 
592  LocalPoint pos_strip_local = TheStripTopol->localPosition(strip);
593  GlobalPoint pos_strip_global = (TkGeom->idToDet(DetectorId))->surface().toGlobal(pos_strip_local);
594 
595  strip_global_position_x = pos_strip_global.x();
596  strip_global_position_y = pos_strip_global.y();
597  strip_global_position_z = pos_strip_global.z();
598  striptree->Fill();
599 
600  // Fill the strip DQM Plots
601  fillStripDQMHistograms();
602  }
603  if(vHotStripsInModule.size()==1)
604  {
605  distance = 999;
606  distanceVsStripNumber[0]->Fill(vHotStripsInModule[0], distance);
607  distanceVsStripNumber[subdetid-2]->Fill(vHotStripsInModule[0], distance);
608  }
609  else if(vHotStripsInModule.size()>1)
610  {
611  for(unsigned int iVec = 0; iVec != vHotStripsInModule.size(); iVec++)
612  {
613  if(iVec==0)
614  distance = vHotStripsInModule[1] - vHotStripsInModule[0];
615  else if(iVec==vHotStripsInModule.size()-1)
616  {
617  distance = vHotStripsInModule[vHotStripsInModule.size()-1] - vHotStripsInModule[vHotStripsInModule.size() -2];
618  }
619  else if(vHotStripsInModule.size()>2)
620  {
621  distanceR = vHotStripsInModule[iVec + 1] - vHotStripsInModule[iVec];
622  distanceL = vHotStripsInModule[iVec] - vHotStripsInModule[iVec - 1];
623  distance = distanceL>distanceR?distanceR:distanceL;
624  }
625  else
626  {
627  std::cout << "ERROR! distance is never computed!!!\n";
628  }
629  distanceVsStripNumber[0]->Fill(vHotStripsInModule[iVec], distance);
630  distanceVsStripNumber[subdetid-2]->Fill(vHotStripsInModule[iVec], distance);
631  }
632  }
633  }
634  }
635 
636  if (BadStripList.begin()!=BadStripList.end())
637  {
638  quality->compact(Detid,BadStripList);
639  SiStripQuality::Range range(BadStripList.begin(),BadStripList.end());
640  quality->put(Detid,range);
641  }
642  BadStripList.clear();
643  }
644 }
645 
646 void SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::iterativeSearch(Apv& histo,std::vector<unsigned int>& vect, int apv){
647  if (!histo.NEntries[apv] || histo.NEntries[apv] <=MinNumEntries_ || histo.NEntries[apv] <= minNevents_)
648  return;
649 
650  size_t startingSize=vect.size();
651  long double diff=1.-prob_;
652 
653  size_t Nbins = histo.th1f[apv]->GetNbinsX();
654  size_t ibinStart = 1;
655  size_t ibinStop = Nbins+1;
656  int MaxEntry = (int)histo.th1f[apv]->GetMaximum();
657 
658  std::vector<long double> vPoissonProbs(MaxEntry+1,0);
659  long double meanVal=1.*histo.NEntries[apv]/(1.*Nbins-histo.NEmptyBins[apv]);
660  evaluatePoissonian(vPoissonProbs,meanVal);
661 
662  for (size_t i=ibinStart; i<ibinStop; ++i){
663  unsigned int entries= (unsigned int)histo.th1f[apv]->GetBinContent(i);
664 
665  if (ishot[i-1]==0){
666  stripoccupancy[i-1] = entries/(double) Nevents_;
667  striphits[i-1] = entries;
668  poissonprob[i-1] = 1-vPoissonProbs[entries];
669  }
670 
671  if (entries<=MinNumEntriesPerStrip_ || entries <= minNevents_)
672  continue;
673 
674  if(diff<vPoissonProbs[entries]){
675  ishot[i-1] = 1;
676  hotstripspermodule++;
677  hotstripsperapv[apv]++;
678  histo.th1f[apv]->SetBinContent(i,0.);
679  histo.NEntries[apv]-=entries;
680  histo.NEmptyBins[apv]++;
681  if (edm::isDebugEnabled())
682  LogTrace("SiStripHotStrip")<< " rejecting strip " << (apv*128)+i-1 << " value " << entries << " diff " << diff << " prob " << vPoissonProbs[entries]<< std::endl;
683  vect.push_back(pQuality->encode((apv*128)+i-1,1,0));
684  }
685 
686  }
687  if (edm::isDebugEnabled())
688  LogTrace("SiStripHotStrip") << " [SiStripHotStripAlgorithmFromClusterOccupancy::iterativeSearch] Nbins="<< Nbins << " MaxEntry="<<MaxEntry << " meanVal=" << meanVal << " NEmptyBins="<<histo.NEmptyBins[apv]<< " NEntries=" << histo.NEntries[apv] << " thEntries " << histo.th1f[apv]->GetEntries()<< " startingSize " << startingSize << " vector.size " << vect.size() << std::endl;
689 
690  if (vect.size()!=startingSize)
691  iterativeSearch(histo,vect,apv);
692 }
693 
694 void SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::evaluatePoissonian(std::vector<long double>& vPoissonProbs, long double& meanVal){
695  for(size_t i=0;i<vPoissonProbs.size();++i){
696  vPoissonProbs[i]= (i==0)?TMath::Poisson(i,meanVal):vPoissonProbs[i-1]+TMath::Poisson(i,meanVal);
697  }
698 }
699 
700 void SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setBasicTreeParameters(int detid){
701  DetId DetectorID=DetId(detid);
702 
703  detrawid = detid;
704  subdetid = DetectorID.subdetId();
705 
706  if (SiStripDetId(detid).stereo() !=0 ) isstereo = 1; // It's a stereo module
707  else isstereo = 0; // It's an rphi module
708  switch (DetectorID.subdetId())
709  {
710  case StripSubdetector::TIB :
711  layer_ring = tTopo->tibLayer(detid);
712  disc = -1;
713  isback = -1;
714  if (tTopo->tibIsExternalString(detid)) isexternalstring = 1;
715  else isexternalstring = 0;
716  if (tTopo->tibIsZMinusSide(detid)) iszminusside = 1;
717  else iszminusside = 0;
718  rodstringpetal = tTopo->tibString(detid);
719  module_number = tTopo->tibModule(detid);
720 
721  break;
722 
723  case StripSubdetector::TID :
724  layer_ring = tTopo->tidRing(detid);
725  disc = tTopo->tidWheel(detid);
726  if (tTopo->tidIsBackRing(detid)) isback = 1;
727  else isback = 0;
728  if (tTopo->tidIsZMinusSide(detid)) iszminusside = 1;
729  else iszminusside = 0;
730  isexternalstring = -1;
731  rodstringpetal = -1;
732  module_number = tTopo->tidModule(detid);
733 
734  break;
735 
736  case StripSubdetector::TOB :
737  layer_ring = tTopo->tobLayer(detid);
738  disc = -1;
739  isback = -1;
740  if (tTopo->tobIsZMinusSide(detid)) iszminusside = 1;
741  else iszminusside = 0;
742  isexternalstring = -1;
743  rodstringpetal = tTopo->tobRod(detid);
744  module_number = tTopo->tobModule(detid);
745 
746  break;
747 
748  case StripSubdetector::TEC :
749  layer_ring = tTopo->tecRing(detid);
750  disc = tTopo->tecWheel(detid);
751  if (tTopo->tecIsBackPetal(detid)) isback = 1;
752  else isback = 0;
753  if (tTopo->tecIsZMinusSide(detid)) iszminusside = 1;
754  else iszminusside = 0;
755  isexternalstring = -1;
756  rodstringpetal = tTopo->tecPetalNumber(detid);
757  module_number = tTopo->tecModule(detid);
758 
759  break;
760 
761  default :
762  std::cout << "### Detector does not belong to TIB, TID, TOB or TEC !? ###" << std::endl;
763  std::cout << "### DetRawId: " << detid << " ###" << std::endl;
764  }
765 }
766 
767 void SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::setMinNumOfEvents()
768 {
769  minNevents_=absolute_occupancy_*Nevents_;
770 }
771 
772 void SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::initializeDQMHistograms()
773 {
774  oss.str("");
775  oss << 1; //runNumber
776 
777  dqmStore = edm::Service<DQMStore>().operator->();
778  dqmStore->setCurrentFolder("ChannelStatusPlots");
779 
780  // Initialize histograms
781  subDetName.push_back(""); subDetName.push_back("TIB"); subDetName.push_back("TID"); subDetName.push_back("TOB"); subDetName.push_back("TEC");
782  nLayers.push_back(0); nLayers.push_back(4); nLayers.push_back(3); nLayers.push_back(6); nLayers.push_back(9);
783  layerName.push_back(""); layerName.push_back("Layer"); layerName.push_back("Disk"); layerName.push_back("Layer"); layerName.push_back("Disk");
784 
785  std::string histoName;
786  std::string histoTitle;
787 
788  for(unsigned int i = 0; i < subDetName.size(); i++)
789  {
790  histoName = "distanceVsStripNumber" + subDetName[i];
791  histoTitle = "Distance between hot strips vs. strip number";
792  if(i!=0)
793  histoTitle += " in " + subDetName[i];
794  tmp = dqmStore->book2D(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5, 999, 0.5, 999.5);
795  distanceVsStripNumber.push_back(tmp->getTH2F());
796 
797  histoName = "pfxDistanceVsStripNumber" + subDetName[i];
798  tmp_prof = new TProfile(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5);
799  tmp = dqmStore->bookProfile(histoName.c_str(), tmp_prof);
800  pfxDistanceVsStripNumber.push_back(tmp->getTProfile());
801  pfxDistanceVsStripNumber[i]->GetXaxis()->SetTitle("Strip");
802  pfxDistanceVsStripNumber[i]->GetYaxis()->SetTitle("Distance");
803 
804  histoName = "projXDistanceVsStripNumber" + subDetName[i];
805  histoTitle = "Number of hot strips vs. strip number";
806  if(i!=0)
807  histoTitle += " in " + subDetName[i];
808  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5);
809  projXDistanceVsStripNumber.push_back(tmp->getTH1F());
810  projXDistanceVsStripNumber[i]->GetXaxis()->SetTitle("Strip");
811  projXDistanceVsStripNumber[i]->GetYaxis()->SetTitle("N_{hot}");
812 
813  histoName = "projYDistanceVsStripNumber" + subDetName[i];
814  histoTitle = "Distribution of distance between hot strips";
815  if(i!=0)
816  histoTitle += " in " + subDetName[i];
817  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 999, 0.5, 999.5);
818  projYDistanceVsStripNumber.push_back(tmp->getTH1F());
819  projYDistanceVsStripNumber[i]->GetXaxis()->SetTitle("Distance");
820  projYDistanceVsStripNumber[i]->GetYaxis()->SetTitle("N_{strips}");
821 
822  //
823  histoName = "occupancyVsStripNumber" + subDetName[i];
824  histoTitle = "Occupancy of strips vs. strip number";
825  if(i!=0)
826  histoTitle += " in " + subDetName[i];
827  tmp = dqmStore->book2D(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5, 1000, -8.,0.);
828  occupancyVsStripNumber.push_back(tmp->getTH2F());
829 
830  histoName = "pfxOccupancyVsStripNumber" + subDetName[i];
831  tmp_prof = new TProfile(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5);
832  tmp = dqmStore->bookProfile(histoName.c_str(), tmp_prof);
833  pfxOccupancyVsStripNumber.push_back(tmp->getTProfile());
834  pfxOccupancyVsStripNumber[i]->GetXaxis()->SetTitle("Strip");
835  pfxOccupancyVsStripNumber[i]->GetYaxis()->SetTitle("log_{10}(Occupancy)");
836 
837  histoName = "projYOccupancyVsStripNumber" + subDetName[i];
838  histoTitle = "Distribution of strip occupancy";
839  if(i!=0)
840  histoTitle += " in " + subDetName[i];
841  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 1000, -8., 0.);
842  projYOccupancyVsStripNumber.push_back(tmp->getTH1F());
843  projYOccupancyVsStripNumber[i]->GetXaxis()->SetTitle("log_{10}(Occupancy)");
844  projYOccupancyVsStripNumber[i]->GetYaxis()->SetTitle("N_{strips}");
845 
846  //
847  histoName = "occupancyHotStripsVsStripNumber" + subDetName[i];
848  histoTitle = "Occupancy of hot strips vs. strip number";
849  if(i!=0)
850  histoTitle += " in " + subDetName[i];
851  tmp = dqmStore->book2D(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5, 1000, -8., 0.);
852  occupancyHotStripsVsStripNumber.push_back(tmp->getTH2F());
853 
854  histoName = "pfxOccupancyHotStripsVsStripNumber" + subDetName[i];
855  tmp_prof = new TProfile(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5);
856  tmp = dqmStore->bookProfile(histoName.c_str(), tmp_prof);
857  pfxOccupancyHotStripsVsStripNumber.push_back(tmp->getTProfile());
858  pfxOccupancyHotStripsVsStripNumber[i]->GetXaxis()->SetTitle("Strip");
859  pfxOccupancyHotStripsVsStripNumber[i]->GetYaxis()->SetTitle("log_{10}(Occupancy)");
860 
861  histoName = "projYOccupancyHotStripsVsStripNumber" + subDetName[i];
862  histoTitle = "Distribution of hot strip occupancy";
863  if(i!=0)
864  histoTitle += " in " + subDetName[i];
865  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 1000, -8., 0.);
866  projYOccupancyHotStripsVsStripNumber.push_back(tmp->getTH1F());
867  projYOccupancyHotStripsVsStripNumber[i]->GetXaxis()->SetTitle("log_{10}(Occupancy)");
868  projYOccupancyHotStripsVsStripNumber[i]->GetYaxis()->SetTitle("N_{strips}");
869 
870  //
871  histoName = "occupancyGoodStripsVsStripNumber" + subDetName[i];
872  histoTitle = "Occupancy of good strips vs. strip number";
873  if(i!=0)
874  histoTitle += " in " + subDetName[i];
875  tmp = dqmStore->book2D(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5, 1000, -8., 0.);
876  occupancyGoodStripsVsStripNumber.push_back(tmp->getTH2F());
877 
878  histoName = "pfxOccupancyGoodStripsVsStripNumber" + subDetName[i];
879  tmp_prof = new TProfile(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5);
880  tmp = dqmStore->bookProfile(histoName.c_str(), tmp_prof);
881  pfxOccupancyGoodStripsVsStripNumber.push_back(tmp->getTProfile());
882  pfxOccupancyGoodStripsVsStripNumber[i]->GetXaxis()->SetTitle("Strip");
883  pfxOccupancyGoodStripsVsStripNumber[i]->GetYaxis()->SetTitle("log_{10}(Occupancy)");
884 
885  histoName = "projYOccupancyGoodStripsVsStripNumber" + subDetName[i];
886  histoTitle = "Distribution of good strip occupancy";
887  if(i!=0)
888  histoTitle += " in " + subDetName[i];
889  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 1000, -8., 0.);
890  projYOccupancyGoodStripsVsStripNumber.push_back(tmp->getTH1F());
891  projYOccupancyGoodStripsVsStripNumber[i]->GetXaxis()->SetTitle("log_{10}(Occupancy)");
892  projYOccupancyGoodStripsVsStripNumber[i]->GetYaxis()->SetTitle("N_{strips}");
893 
894  //
895  histoName = "poissonProbVsStripNumber" + subDetName[i];
896  histoTitle = "Poisson probability of strips vs. strip number";
897  if(i!=0)
898  histoTitle += " in " + subDetName[i];
899  tmp = dqmStore->book2D(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5, 1000, -18., 0.);
900  poissonProbVsStripNumber.push_back(tmp->getTH2F());
901 
902  histoName = "pfxPoissonProbVsStripNumber" + subDetName[i];
903  tmp_prof = new TProfile(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5);
904  tmp = dqmStore->bookProfile(histoName.c_str(), tmp_prof);
905  pfxPoissonProbVsStripNumber.push_back(tmp->getTProfile());
906  pfxPoissonProbVsStripNumber[i]->GetXaxis()->SetTitle("Strip");
907  pfxPoissonProbVsStripNumber[i]->GetYaxis()->SetTitle("log_{10}(Probability)");
908 
909  histoName = "projYPoissonProbVsStripNumber" + subDetName[i];
910  histoTitle = "Distribution of strip Poisson probability";
911  if(i!=0)
912  histoTitle += " in " + subDetName[i];
913  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 1000, -18., 0.);
914  projYPoissonProbVsStripNumber.push_back(tmp->getTH1F());
915  projYPoissonProbVsStripNumber[i]->GetXaxis()->SetTitle("log_{10}(Probability)");
916  projYPoissonProbVsStripNumber[i]->GetYaxis()->SetTitle("N_{strips}");
917 
918  //
919  histoName = "poissonProbHotStripsVsStripNumber" + subDetName[i];
920  histoTitle = "Poisson probability of hot strips vs. strip number";
921  if(i!=0)
922  histoTitle += " in " + subDetName[i];
923  tmp = dqmStore->book2D(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5, 1000, -18., 0.);
924  poissonProbHotStripsVsStripNumber.push_back(tmp->getTH2F());
925 
926  histoName = "pfxPoissonProbHotStripsVsStripNumber" + subDetName[i];
927  tmp_prof = new TProfile(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5);
928  tmp = dqmStore->bookProfile(histoName.c_str(), tmp_prof);
929  pfxPoissonProbHotStripsVsStripNumber.push_back(tmp->getTProfile());
930  pfxPoissonProbHotStripsVsStripNumber[i]->GetXaxis()->SetTitle("Strip");
931  pfxPoissonProbHotStripsVsStripNumber[i]->GetYaxis()->SetTitle("log_{10}(Probability)");
932 
933  histoName = "projYPoissonProbHotStripsVsStripNumber" + subDetName[i];
934  histoTitle = "Distribution of hot strip Poisson probability";
935  if(i!=0)
936  histoTitle += " in " + subDetName[i];
937  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 1000, -18., 0.);
938  projYPoissonProbHotStripsVsStripNumber.push_back(tmp->getTH1F());
939  projYPoissonProbHotStripsVsStripNumber[i]->GetXaxis()->SetTitle("log_{10}(Probability)");
940  projYPoissonProbHotStripsVsStripNumber[i]->GetYaxis()->SetTitle("N_{strips}");
941 
942  //
943  histoName = "poissonProbGoodStripsVsStripNumber" + subDetName[i];
944  histoTitle = "Poisson probability of good strips vs. strip number";
945  if(i!=0)
946  histoTitle += " in " + subDetName[i];
947  tmp = dqmStore->book2D(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5, 1000, -18., 0.);
948  poissonProbGoodStripsVsStripNumber.push_back(tmp->getTH2F());
949 
950  histoName = "pfxPoissonProbGoodStripsVsStripNumber" + subDetName[i];
951  tmp_prof = new TProfile(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5);
952  tmp = dqmStore->bookProfile(histoName.c_str(), tmp_prof);
953  pfxPoissonProbGoodStripsVsStripNumber.push_back(tmp->getTProfile());
954  pfxPoissonProbGoodStripsVsStripNumber[i]->GetXaxis()->SetTitle("Strip");
955  pfxPoissonProbGoodStripsVsStripNumber[i]->GetYaxis()->SetTitle("log_{10}(Probability)");
956 
957  histoName = "projYPoissonProbGoodStripsVsStripNumber" + subDetName[i];
958  histoTitle = "Distribution of good strip Poisson probability";
959  if(i!=0)
960  histoTitle += " in " + subDetName[i];
961  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 1000, -18., 0.);
962  projYPoissonProbGoodStripsVsStripNumber.push_back(tmp->getTH1F());
963  projYPoissonProbGoodStripsVsStripNumber[i]->GetXaxis()->SetTitle("log_{10}(Probability)");
964  projYPoissonProbGoodStripsVsStripNumber[i]->GetYaxis()->SetTitle("N_{strips}");
965 
966  //
967  histoName = "nHitsVsStripNumber" + subDetName[i];
968  histoTitle = "NHits in strips vs. strip number";
969  if(i!=0)
970  histoTitle += " in " + subDetName[i];
971  tmp = dqmStore->book2D(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5, 10000, -0.5, 9999.5);
972  nHitsVsStripNumber.push_back(tmp->getTH2F());
973 
974  histoName = "pfxNHitsVsStripNumber" + subDetName[i];
975  tmp_prof = new TProfile(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5);
976  tmp = dqmStore->bookProfile(histoName.c_str(), tmp_prof);
977  pfxNHitsVsStripNumber.push_back(tmp->getTProfile());
978 
979  histoName = "projXNHitsVsStripNumber" + subDetName[i];
980  histoTitle = "Cumulative nHits in strips vs. strip number";
981  if(i!=0)
982  histoTitle += " in " + subDetName[i];
983  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5);
984  projXNHitsVsStripNumber.push_back(tmp->getTH1F());
985 
986  histoName = "projYNHitsVsStripNumber" + subDetName[i];
987  histoTitle = "Distribution of nHits for all strips";
988  if(i!=0)
989  histoTitle += " in " + subDetName[i];
990  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 10000, -0.5, 9999.5);
991  projYNHitsVsStripNumber.push_back(tmp->getTH1F());
992  projYNHitsVsStripNumber[i]->GetXaxis()->SetTitle("N_{hits}");
993  projYNHitsVsStripNumber[i]->GetYaxis()->SetTitle("N_{strips}");
994 
995  //
996  histoName = "nHitsHotStripsVsStripNumber" + subDetName[i];
997  histoTitle = "NHits in hot strips vs. strip number";
998  if(i!=0)
999  histoTitle += " in " + subDetName[i];
1000  tmp = dqmStore->book2D(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5, 10000, -0.5, 9999.5);
1001  nHitsHotStripsVsStripNumber.push_back(tmp->getTH2F());
1002 
1003  histoName = "pfxNHitsHotStripsVsStripNumber" + subDetName[i];
1004  tmp_prof = new TProfile(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5);
1005  tmp = dqmStore->bookProfile(histoName.c_str(), tmp_prof);
1006  pfxNHitsHotStripsVsStripNumber.push_back(tmp->getTProfile());
1007 
1008  histoName = "projXNHitsHotStripsVsStripNumber" + subDetName[i];
1009  histoTitle = "Cumulative nHits in hot strips vs. strip number";
1010  if(i!=0)
1011  histoTitle += " in " + subDetName[i];
1012  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5);
1013  projXNHitsHotStripsVsStripNumber.push_back(tmp->getTH1F());
1014 
1015  histoName = "projYNHitsHotStripsVsStripNumber" + subDetName[i];
1016  histoTitle = "Distribution of nHits for hot strips";
1017  if(i!=0)
1018  histoTitle += " in " + subDetName[i];
1019  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 10000, -0.5, 9999.5);
1020  projYNHitsHotStripsVsStripNumber.push_back(tmp->getTH1F());
1021  projYNHitsHotStripsVsStripNumber[i]->GetXaxis()->SetTitle("N_{hits}");
1022  projYNHitsHotStripsVsStripNumber[i]->GetYaxis()->SetTitle("N_{strips}");
1023 
1024  //
1025  histoName = "nHitsGoodStripsVsStripNumber" + subDetName[i];
1026  histoTitle = "NHits in good strips vs. strip number";
1027  if(i!=0)
1028  histoTitle += " in " + subDetName[i];
1029  tmp = dqmStore->book2D(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5, 10000, -0.5, 9999.5);
1030  nHitsGoodStripsVsStripNumber.push_back(tmp->getTH2F());
1031 
1032  histoName = "pfxNHitsGoodStripsVsStripNumber" + subDetName[i];
1033  tmp_prof = new TProfile(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5);
1034  tmp = dqmStore->bookProfile(histoName.c_str(), tmp_prof);
1035  pfxNHitsGoodStripsVsStripNumber.push_back(tmp->getTProfile());
1036 
1037  histoName = "projXNHitsGoodStripsVsStripNumber" + subDetName[i];
1038  histoTitle = "Cumulative nHits in good strips vs. strip number";
1039  if(i!=0)
1040  histoTitle += " in " + subDetName[i];
1041  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 768, 0.5, 768.5);
1042  projXNHitsGoodStripsVsStripNumber.push_back(tmp->getTH1F());
1043 
1044  histoName = "projYNHitsGoodStripsVsStripNumber" + subDetName[i];
1045  histoTitle = "Distribution of nHits for good strips";
1046  if(i!=0)
1047  histoTitle += " in " + subDetName[i];
1048  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 10000, -0.5, 9999.5);
1049  projYNHitsGoodStripsVsStripNumber.push_back(tmp->getTH1F());
1050  projYNHitsGoodStripsVsStripNumber[i]->GetXaxis()->SetTitle("N_{hits}");
1051  projYNHitsGoodStripsVsStripNumber[i]->GetYaxis()->SetTitle("N_{strips}");
1052 
1053  for(unsigned int j = 0; j <= nLayers[i]; j++)
1054  {
1055  histoName = "medianVsAbsoluteOccupancy" + subDetName[i];
1056  if(j!=0)
1057  {
1058  oss.str("");
1059  oss << j;
1060  histoName += layerName[i] + oss.str();
1061  }
1062  histoTitle = "Median APV occupancy vs. absolute APV occupancy";
1063  if(i!=0)
1064  histoTitle += " in " + subDetName[i];
1065  if(j!=0)
1066  {
1067  histoTitle += " " + layerName[i] + " " + oss.str();
1068  }
1069  tmp = dqmStore->book2D(histoName.c_str(), histoTitle.c_str(), 1000, 0., 6., 1000, -1., 3.);
1070  medianVsAbsoluteOccupancy[i][j] = tmp->getTH2F();
1071  medianVsAbsoluteOccupancy[i][j]->Rebin2D(10,10);
1072  medianVsAbsoluteOccupancy[i][j]->GetXaxis()->SetTitle("log_{10}(Abs. Occupancy)");
1073  medianVsAbsoluteOccupancy[i][j]->GetYaxis()->SetTitle("log_{10}(Median Occupancy)");
1074  //
1075  histoName = "medianOccupancy" + subDetName[i];
1076  if(j!=0)
1077  {
1078  oss.str("");
1079  oss << j;
1080  histoName += layerName[i] + oss.str();
1081  }
1082  histoTitle = "Median APV occupancy";
1083  if(i!=0)
1084  histoTitle += " in " + subDetName[i];
1085  if(j!=0)
1086  {
1087  histoTitle += " " + layerName[i] + " " + oss.str();
1088  }
1089  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 1000, -1., 3.);
1090  medianOccupancy[i][j] = tmp->getTH1F();
1091  medianOccupancy[i][j]->GetXaxis()->SetTitle("log_{10}(Occupancy)");
1092  medianOccupancy[i][j]->GetYaxis()->SetTitle("APVs");
1093  //
1094  histoName = "absoluteOccupancy" + subDetName[i];
1095  if(j!=0)
1096  {
1097  oss.str("");
1098  oss << j;
1099  histoName += layerName[i] + oss.str();
1100  }
1101  histoTitle = "Absolute APV occupancy";
1102  if(i!=0)
1103  histoTitle += " in " + subDetName[i];
1104  if(j!=0)
1105  {
1106  histoTitle += " " + layerName[i] + " " + oss.str();
1107  }
1108  tmp = dqmStore->book1D(histoName.c_str(), histoTitle.c_str(), 1000, 0., 6.);
1109  absoluteOccupancy[i][j] = tmp->getTH1F();
1110  absoluteOccupancy[i][j]->GetXaxis()->SetTitle("log_{10}(Occupancy)");
1111  absoluteOccupancy[i][j]->GetYaxis()->SetTitle("APVs");
1112  }
1113  }
1114 }
1115 
1116 void SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::fillStripDQMHistograms()
1117 {
1118  double logStripOccupancy = -1;
1119  double logPoissonProb = -1;
1120 
1121  if (singleStripOccupancy>0) logStripOccupancy = log10(singleStripOccupancy);
1122  if (poissonProb>0) logPoissonProb = log10(fabs(poissonProb));
1123 
1124  occupancyVsStripNumber[0]->Fill(strip_number,logStripOccupancy);
1125  occupancyVsStripNumber[subdetid-2]->Fill(strip_number,logStripOccupancy);
1126  poissonProbVsStripNumber[0]->Fill(strip_number,logPoissonProb);
1127  poissonProbVsStripNumber[subdetid-2]->Fill(strip_number,logPoissonProb);
1128  nHitsVsStripNumber[0]->Fill(strip_number,stripHits);
1129  nHitsVsStripNumber[subdetid-2]->Fill(strip_number,stripHits);
1130 
1131  if(isHot)
1132  {
1133  vHotStripsInModule.push_back(strip_number);
1134  occupancyHotStripsVsStripNumber[0]->Fill(strip_number,logStripOccupancy);
1135  occupancyHotStripsVsStripNumber[subdetid-2]->Fill(strip_number,logStripOccupancy);
1136  poissonProbHotStripsVsStripNumber[0]->Fill(strip_number,logPoissonProb);
1137  poissonProbHotStripsVsStripNumber[subdetid-2]->Fill(strip_number,logPoissonProb);
1138  nHitsHotStripsVsStripNumber[0]->Fill(strip_number,stripHits);
1139  nHitsHotStripsVsStripNumber[subdetid-2]->Fill(strip_number,stripHits);
1140  }
1141  else
1142  {
1143  occupancyGoodStripsVsStripNumber[0]->Fill(strip_number,logStripOccupancy);
1144  occupancyGoodStripsVsStripNumber[subdetid-2]->Fill(strip_number,logStripOccupancy);
1145  poissonProbGoodStripsVsStripNumber[0]->Fill(strip_number,logPoissonProb);
1146  poissonProbGoodStripsVsStripNumber[subdetid-2]->Fill(strip_number,logPoissonProb);
1147  nHitsGoodStripsVsStripNumber[0]->Fill(strip_number,stripHits);
1148  nHitsGoodStripsVsStripNumber[subdetid-2]->Fill(strip_number,stripHits);
1149  }
1150 }
edm::isDebugEnabled
bool isDebugEnabled()
Definition: MessageLogger.cc:71
i
int i
Definition: DBlmapReader.cc:9
TrackerTopology::tibLayer
unsigned int tibLayer(const DetId &id) const
Definition: TrackerTopology.h:121
TrackerTopology::tibString
unsigned int tibString(const DetId &id) const
Definition: TrackerTopology.h:323
TrackerTopology::tidRing
unsigned int tidRing(const DetId &id) const
Definition: TrackerTopology.h:202
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:872
DQMStore::cd
void cd(void)
go to top directory (ie. root)
Definition: DQMStore.cc:561
SiStripMonitorClusterAlca_cfi.Nbins
tuple Nbins
Definition: SiStripMonitorClusterAlca_cfi.py:33
TrackerTopology::tecRing
unsigned int tecRing(const DetId &id) const
ring id
Definition: TrackerTopology.h:199
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
TrackerTopology::tidWheel
unsigned int tidWheel(const DetId &id) const
Definition: TrackerTopology.h:182
SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::CalculateMeanAndRMS
void CalculateMeanAndRMS(const std::vector< Apv > &, std::pair< double, double > *, int)
Definition: SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy.cc:366
StripGeomDetUnit::specificTopology
virtual const StripTopology & specificTopology() const
Returns a reference to the strip proxy topology.
Definition: StripGeomDetUnit.cc:19
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
TrackerTopology::tecIsZMinusSide
bool tecIsZMinusSide(const DetId &id) const
Definition: TrackerTopology.h:237
TrackerTopology::tidIsZMinusSide
bool tidIsZMinusSide(const DetId &id) const
Definition: TrackerTopology.h:234
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:250
TrackerTopology::tibIsExternalString
bool tibIsExternalString(const DetId &id) const
Definition: TrackerTopology.h:339
TrackerTopology::tibIsZMinusSide
bool tibIsZMinusSide(const DetId &id) const
Definition: TrackerTopology.h:231
TrackerTopology::tidModule
unsigned int tidModule(const DetId &id) const
Definition: TrackerTopology.h:147
TrackerTopology::tidIsBackRing
bool tidIsBackRing(const DetId &id) const
Definition: TrackerTopology.h:351
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
TrackerTopology::tobIsZMinusSide
bool tobIsZMinusSide(const DetId &id) const
Definition: TrackerTopology.h:228
j
int j
Definition: DBlmapReader.cc:9
DQMStore::save
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", const uint32_t run=0, SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
Definition: DQMStore.cc:2296
DQMStore::bookProfile
MonitorElement * bookProfile(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const char *option="s")
Definition: DQMStore.cc:1186
SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy
SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy(const edm::ParameterSet &, const TrackerTopology *)
Definition: SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy.cc:13
first
bool first
Definition: L1TdeRCT.cc:79
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
TrackerGeometry::idToDet
virtual const GeomDet * idToDet(DetId) const
Definition: TrackerGeometry.cc:137
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
TrackerTopology::tibModule
unsigned int tibModule(const DetId &id) const
Definition: TrackerTopology.h:144
TrackerTopology::tecModule
unsigned int tecModule(const DetId &id) const
Definition: TrackerTopology.h:141
SiStripDetId
Detector identifier class for the strip tracker.
Definition: SiStripDetId.h:17
TrackerTopology::tecIsBackPetal
bool tecIsBackPetal(const DetId &id) const
Definition: TrackerTopology.h:316
DetId
Definition: DetId.h:18
Range
PixelRecoRange< float > Range
Definition: PixelTripletHLTGenerator.cc:27
MonitorElement::getTH1F
TH1F * getTH1F(void) const
Definition: MonitorElement.cc:1498
SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::AnalyzeOccupancy
void AnalyzeOccupancy(SiStripQuality *, std::vector< Apv > &, std::pair< double, double > *, std::vector< unsigned int > &, edm::ESHandle< SiStripQuality > &)
Definition: SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy.cc:421
TrackerTopology::tobModule
unsigned int tobModule(const DetId &id) const
Definition: TrackerTopology.h:138
a
double a
Definition: hdecay.h:121
MonitorElement::getTProfile
TProfile * getTProfile(void) const
Definition: MonitorElement.cc:1554
StripTopology::localPosition
virtual LocalPoint localPosition(float strip) const =0
TrackerTopology::tecPetalNumber
unsigned int tecPetalNumber(const DetId &id) const
Definition: TrackerTopology.h:208
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::extractBadAPVSandStrips
void extractBadAPVSandStrips(SiStripQuality *, HistoMap &, edm::ESHandle< SiStripQuality > &)
Definition: SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy.cc:32
edm::false
volatile std::atomic< bool > shutdown_flag false
Definition: UnixSignalHandlers.cc:22
TrackerTopology::tobRod
unsigned int tobRod(const DetId &id) const
Definition: TrackerTopology.h:174
MonitorElement::getTH2F
TH2F * getTH2F(void) const
Definition: MonitorElement.cc:1522
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:1000
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
TrackerTopology::tecWheel
unsigned int tecWheel(const DetId &id) const
Definition: TrackerTopology.h:179
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:584
TrackerTopology::tobLayer
unsigned int tobLayer(const DetId &id) const
Definition: TrackerTopology.h:118
Generated for CMSSW Reference Manual by doxygen 1.8.5