CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
SiStripBadAPVAlgorithmFromClusterOccupancy.cc
Go to the documentation of this file.
1 #include "CalibTracker/SiStripQuality/interface/SiStripBadAPVAlgorithmFromClusterOccupancy.h"
2 
3 
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 SiStripBadAPVAlgorithmFromClusterOccupancy::SiStripBadAPVAlgorithmFromClusterOccupancy(const edm::ParameterSet& iConfig, const TrackerTopology* theTopo):
14  lowoccupancy_(0),
15  highoccupancy_(100),
16  absolutelow_(0),
17  numberiterations_(2),
18  Nevents_(0),
19  occupancy_(0),
20  OutFileName_("Occupancy.root"),
21  UseInputDB_(iConfig.getUntrackedParameter<bool>("UseInputDB",false)),
22  tTopo(theTopo)
23  {
24  minNevents_=Nevents_*occupancy_;
25  }
26 
27 SiStripBadAPVAlgorithmFromClusterOccupancy::~SiStripBadAPVAlgorithmFromClusterOccupancy(){
28  LogTrace("SiStripBadAPVAlgorithmFromClusterOccupancy")<<"[SiStripBadAPVAlgorithmFromClusterOccupancy::~SiStripBadAPVAlgorithmFromClusterOccupancy] "<<std::endl;
29 }
30 
31 void SiStripBadAPVAlgorithmFromClusterOccupancy::extractBadAPVs(SiStripQuality* siStripQuality,HistoMap& DM, edm::ESHandle<SiStripQuality>& inSiStripQuality){
32 
33  LogTrace("SiStripBadAPVAlgorithmFromClusterOccupancy")<<"[SiStripBadAPVAlgorithmFromClusterOccupancy::extractBadAPVs] "<<std::endl;
34 
35  if (WriteOutputFile_==true){
36  f = new TFile(OutFileName_.c_str(),"RECREATE");
37  f->cd();
38 
39  apvtree = new TTree("moduleOccupancy","tree");
40 
41  apvtree->Branch("DetRawId", &detrawid, "DetRawId/I");
42  apvtree->Branch("SubDetId", &subdetid, "SubDetId/I");
43  apvtree->Branch("Layer_Ring", &layer_ring, "Layer_Ring/I");
44  apvtree->Branch("Disc", &disc, "Disc/I");
45  apvtree->Branch("IsBack", &isback, "IsBack/I");
46  apvtree->Branch("IsExternalString", &isexternalstring, "IsExternalString/I");
47  apvtree->Branch("IsZMinusSide", &iszminusside, "IsZMinusSide/I");
48  apvtree->Branch("RodStringPetal", &rodstringpetal, "RodStringPetal/I");
49  apvtree->Branch("IsStereo", &isstereo, "IsStereo/I");
50  apvtree->Branch("ModuleNumber", &module_number, "ModuleNumber/I");
51  apvtree->Branch("NumberOfStrips", &number_strips, "NumberOfStrips/I");
52  apvtree->Branch("APVGlobalPositionX", &global_position_x, "APVGlobalPositionX/F");
53  apvtree->Branch("APVGlobalPositionY", &global_position_y, "APVGlobalPositionY/F");
54  apvtree->Branch("APVGlobalPositionZ", &global_position_z, "APVGlobalPositionZ/F");
55  apvtree->Branch("APVNumber", &apv_number, "APVNumber/I");
56  apvtree->Branch("APVAbsoluteOccupancy", &apvAbsoluteOccupancy, "apvAbsoluteOccupancy/D");
57  apvtree->Branch("APVMedianOccupancy", &apvMedianOccupancy, "apvMedianOccupancy/D");
58  }
59 
60  HistoMap::iterator it=DM.begin();
61  HistoMap::iterator itEnd=DM.end();
62  std::vector<unsigned int> badStripList;
63  uint32_t detid;
64  for (;it!=itEnd;++it){
65 
66  Apv APV;
67 
68  for (int apv=0; apv<6; apv++)
69  {
70  APV.apvMedian[apv] = 0;
71  APV.apvabsoluteOccupancy[apv] = 0;
72 
73  for (int strip=0; strip<128; strip++)
74  {
75  stripOccupancy[apv][strip] = 0;
76  stripWeight[apv][strip] = 0;
77  }
78  }
79 
80  pHisto phisto;
81  phisto._th1f = it->second.get();
82  phisto._NEntries = (int)phisto._th1f->GetEntries();
83  phisto._NBins = phisto._th1f->GetNbinsX();
84 
85  number_strips = (int)phisto._NBins;
86  number_apvs = number_strips/128;
87  APV.numberApvs = number_apvs;
88 
89  for (int apv=0; apv<number_apvs; apv++)
90  {
91  for (int strip=0; strip<128; strip++)
92  {
93  stripOccupancy[apv][strip] = phisto._th1f->GetBinContent((apv*128)+strip+1); // Remember: Bin=0 is underflow bin!
94  stripWeight[apv][strip] = 1;
95  APV.apvabsoluteOccupancy[apv] += phisto._th1f->GetBinContent((apv*128)+strip+1); // Remember: Bin=0 is underflow bin!
96  }
97  }
98 
99  for (int apv=0; apv<number_apvs; apv++)
100  {
101  APV.apvMedian[apv] = TMath::Median(128,stripOccupancy[apv],stripWeight[apv]);
102  }
103 
104  detid=it->first;
105  DetId detectorId=DetId(detid);
106 
107  if (edm::isDebugEnabled())
108  LogTrace("SiStripBadAPV") << "Analyzing detid " << detid<< std::endl;
109 
110  detrawid = detid;
111  APV.detrawId = detrawid;
112  subdetid = detectorId.subdetId();
113  switch (detectorId.subdetId())
114  {
115  case StripSubdetector::TIB :
116  layer_ring = tTopo->tibLayer(detrawid);
117  disc = -1;
118  isstereo = tTopo->tibIsStereo(detrawid);
119  isback = -1;
120  if (tTopo->tibIsExternalString(detrawid)) isexternalstring = 1;
121  else isexternalstring = 0;
122  if (tTopo->tibIsZMinusSide(detrawid)) iszminusside = 1;
123  else iszminusside = 0;
124  rodstringpetal = tTopo->tibString(detrawid);
125  module_number = tTopo->tibModule(detrawid);
126  APV.modulePosition = module_number;
127 
128  if (layer_ring == 1) medianValues_TIB_Layer1.push_back(APV);
129  else if (layer_ring == 2) medianValues_TIB_Layer2.push_back(APV);
130  else if (layer_ring == 3) medianValues_TIB_Layer3.push_back(APV);
131  else if (layer_ring == 4) medianValues_TIB_Layer4.push_back(APV);
132  break;
133 
134  case StripSubdetector::TID :
135  layer_ring = tTopo->tidRing(detrawid);
136  disc = tTopo->tidWheel(detrawid);
137  isstereo = tTopo->tidIsStereo(detrawid);
138  if (tTopo->tidIsBackRing(detrawid)) isback = 1;
139  else isback = 0;
140  if (tTopo->tidIsZMinusSide(detrawid)) iszminusside = 1;
141  else iszminusside = 0;
142  isexternalstring = -1;
143  rodstringpetal = -1;
144  module_number = tTopo->tidModule(detrawid);
145  APV.modulePosition = layer_ring;
146 
147  if (iszminusside==0)
148  {
149  if (disc==1) medianValues_TIDPlus_Disc1.push_back(APV);
150  else if (disc==2) medianValues_TIDPlus_Disc2.push_back(APV);
151  else if (disc==3) medianValues_TIDPlus_Disc3.push_back(APV);
152  }
153  else if (iszminusside==1)
154  {
155  if (disc==1) medianValues_TIDMinus_Disc1.push_back(APV);
156  else if (disc==2) medianValues_TIDMinus_Disc2.push_back(APV);
157  else if (disc==3) medianValues_TIDMinus_Disc3.push_back(APV);
158  }
159  break;
160 
161  case StripSubdetector::TOB :
162  layer_ring = tTopo->tobLayer(detrawid);
163  disc = -1;
164  isstereo = tTopo->tobIsStereo(detrawid);
165  isback = -1;
166  if (tTopo->tobIsZMinusSide(detrawid)) iszminusside = 1;
167  else iszminusside = 0;
168  isexternalstring = -1;
169  rodstringpetal = tTopo->tobRod(detrawid);
170  module_number = tTopo->tobModule(detrawid);
171  APV.modulePosition = module_number;
172 
173  if (layer_ring == 1) medianValues_TOB_Layer1.push_back(APV);
174  else if (layer_ring == 2) medianValues_TOB_Layer2.push_back(APV);
175  else if (layer_ring == 3) medianValues_TOB_Layer3.push_back(APV);
176  else if (layer_ring == 4) medianValues_TOB_Layer4.push_back(APV);
177  else if (layer_ring == 5) medianValues_TOB_Layer5.push_back(APV);
178  else if (layer_ring == 6) medianValues_TOB_Layer6.push_back(APV);
179  break;
180 
181  case StripSubdetector::TEC :
182  layer_ring = tTopo->tecRing(detrawid);
183  disc = tTopo->tecWheel(detrawid);
184  isstereo = tTopo->tecIsStereo(detrawid);
185  if (tTopo->tecIsBackPetal(detrawid)) isback = 1;
186  else isback = 0;
187  if (tTopo->tecIsZMinusSide(detrawid)) iszminusside = 1;
188  else iszminusside = 0;
189  isexternalstring = -1;
190  rodstringpetal = tTopo->tecPetalNumber(detrawid);
191  module_number = tTopo->tecModule(detrawid);
192  APV.modulePosition = layer_ring;
193 
194  if (iszminusside==0)
195  {
196  if (disc==1) medianValues_TECPlus_Disc1.push_back(APV);
197  else if (disc==2) medianValues_TECPlus_Disc2.push_back(APV);
198  else if (disc==3) medianValues_TECPlus_Disc3.push_back(APV);
199  else if (disc==4) medianValues_TECPlus_Disc4.push_back(APV);
200  else if (disc==5) medianValues_TECPlus_Disc5.push_back(APV);
201  else if (disc==6) medianValues_TECPlus_Disc6.push_back(APV);
202  else if (disc==7) medianValues_TECPlus_Disc7.push_back(APV);
203  else if (disc==8) medianValues_TECPlus_Disc8.push_back(APV);
204  else if (disc==9) medianValues_TECPlus_Disc9.push_back(APV);
205  }
206  else if (iszminusside==1)
207  {
208  if (disc==1) medianValues_TECMinus_Disc1.push_back(APV);
209  else if (disc==2) medianValues_TECMinus_Disc2.push_back(APV);
210  else if (disc==3) medianValues_TECMinus_Disc3.push_back(APV);
211  else if (disc==4) medianValues_TECMinus_Disc4.push_back(APV);
212  else if (disc==5) medianValues_TECMinus_Disc5.push_back(APV);
213  else if (disc==6) medianValues_TECMinus_Disc6.push_back(APV);
214  else if (disc==7) medianValues_TECMinus_Disc7.push_back(APV);
215  else if (disc==8) medianValues_TECMinus_Disc8.push_back(APV);
216  else if (disc==9) medianValues_TECMinus_Disc9.push_back(APV);
217  }
218  break;
219 
220  default :
221  std::cout << "### Detector does not belong to TIB, TID, TOB or TEC !? ###" << std::endl;
222  std::cout << "### DetRawId: " << detrawid << " ###" << std::endl;
223  }
224 
225  const StripGeomDetUnit* theStripDet = dynamic_cast<const StripGeomDetUnit*>( (TkGeom->idToDet(detectorId)) );
226  const StripTopology* theStripTopol = dynamic_cast<const StripTopology*>( &(theStripDet->specificTopology()) );
227 
228  for (int apv=0; apv<number_apvs; apv++)
229  {
230  apv_number = apv+1;
231  apvMedianOccupancy = APV.apvMedian[apv];
232  apvAbsoluteOccupancy = APV.apvabsoluteOccupancy[apv];
233 
234  LocalPoint pos_strip_local = theStripTopol->localPosition((apv*128));
235  GlobalPoint pos_strip_global = (TkGeom->idToDet(detectorId))->surface().toGlobal(pos_strip_local);
236 
237  global_position_x = pos_strip_global.x();
238  global_position_y = pos_strip_global.y();
239  global_position_z = pos_strip_global.z();
240 
241  if (WriteOutputFile_==true) apvtree->Fill();
242  }
243 
244  } // end loop on modules
245 
246  // Calculate Mean and RMS for each Layer
247  CalculateMeanAndRMS(medianValues_TIB_Layer1,MeanAndRms_TIB_Layer1,numberiterations_);
248  CalculateMeanAndRMS(medianValues_TIB_Layer2,MeanAndRms_TIB_Layer2,numberiterations_);
249  CalculateMeanAndRMS(medianValues_TIB_Layer3,MeanAndRms_TIB_Layer3,numberiterations_);
250  CalculateMeanAndRMS(medianValues_TIB_Layer4,MeanAndRms_TIB_Layer4,numberiterations_);
251 
252  CalculateMeanAndRMS(medianValues_TOB_Layer1,MeanAndRms_TOB_Layer1,numberiterations_);
253  CalculateMeanAndRMS(medianValues_TOB_Layer2,MeanAndRms_TOB_Layer2,numberiterations_);
254  CalculateMeanAndRMS(medianValues_TOB_Layer3,MeanAndRms_TOB_Layer3,numberiterations_);
255  CalculateMeanAndRMS(medianValues_TOB_Layer4,MeanAndRms_TOB_Layer4,numberiterations_);
256  CalculateMeanAndRMS(medianValues_TOB_Layer5,MeanAndRms_TOB_Layer5,numberiterations_);
257  CalculateMeanAndRMS(medianValues_TOB_Layer6,MeanAndRms_TOB_Layer6,numberiterations_);
258 
259  CalculateMeanAndRMS(medianValues_TIDPlus_Disc1,MeanAndRms_TIDPlus_Disc1,numberiterations_);
260  CalculateMeanAndRMS(medianValues_TIDPlus_Disc2,MeanAndRms_TIDPlus_Disc2,numberiterations_);
261  CalculateMeanAndRMS(medianValues_TIDPlus_Disc3,MeanAndRms_TIDPlus_Disc3,numberiterations_);
262  CalculateMeanAndRMS(medianValues_TIDMinus_Disc1,MeanAndRms_TIDMinus_Disc1,numberiterations_);
263  CalculateMeanAndRMS(medianValues_TIDMinus_Disc2,MeanAndRms_TIDMinus_Disc2,numberiterations_);
264  CalculateMeanAndRMS(medianValues_TIDMinus_Disc3,MeanAndRms_TIDMinus_Disc3,numberiterations_);
265 
266  CalculateMeanAndRMS(medianValues_TECPlus_Disc1,MeanAndRms_TECPlus_Disc1,numberiterations_);
267  CalculateMeanAndRMS(medianValues_TECPlus_Disc2,MeanAndRms_TECPlus_Disc2,numberiterations_);
268  CalculateMeanAndRMS(medianValues_TECPlus_Disc3,MeanAndRms_TECPlus_Disc3,numberiterations_);
269  CalculateMeanAndRMS(medianValues_TECPlus_Disc4,MeanAndRms_TECPlus_Disc4,numberiterations_);
270  CalculateMeanAndRMS(medianValues_TECPlus_Disc5,MeanAndRms_TECPlus_Disc5,numberiterations_);
271  CalculateMeanAndRMS(medianValues_TECPlus_Disc6,MeanAndRms_TECPlus_Disc6,numberiterations_);
272  CalculateMeanAndRMS(medianValues_TECPlus_Disc7,MeanAndRms_TECPlus_Disc7,numberiterations_);
273  CalculateMeanAndRMS(medianValues_TECPlus_Disc8,MeanAndRms_TECPlus_Disc8,numberiterations_);
274  CalculateMeanAndRMS(medianValues_TECPlus_Disc9,MeanAndRms_TECPlus_Disc9,numberiterations_);
275 
276  CalculateMeanAndRMS(medianValues_TECMinus_Disc1,MeanAndRms_TECMinus_Disc1,numberiterations_);
277  CalculateMeanAndRMS(medianValues_TECMinus_Disc2,MeanAndRms_TECMinus_Disc2,numberiterations_);
278  CalculateMeanAndRMS(medianValues_TECMinus_Disc3,MeanAndRms_TECMinus_Disc3,numberiterations_);
279  CalculateMeanAndRMS(medianValues_TECMinus_Disc4,MeanAndRms_TECMinus_Disc4,numberiterations_);
280  CalculateMeanAndRMS(medianValues_TECMinus_Disc5,MeanAndRms_TECMinus_Disc5,numberiterations_);
281  CalculateMeanAndRMS(medianValues_TECMinus_Disc6,MeanAndRms_TECMinus_Disc6,numberiterations_);
282  CalculateMeanAndRMS(medianValues_TECMinus_Disc7,MeanAndRms_TECMinus_Disc7,numberiterations_);
283  CalculateMeanAndRMS(medianValues_TECMinus_Disc8,MeanAndRms_TECMinus_Disc8,numberiterations_);
284  CalculateMeanAndRMS(medianValues_TECMinus_Disc9,MeanAndRms_TECMinus_Disc9,numberiterations_);
285 
286  pQuality=siStripQuality;
287  badStripList.clear();
288 
289  // Analyze the APV Occupancy for hot APVs
290  AnalyzeOccupancy(siStripQuality,medianValues_TIB_Layer1,MeanAndRms_TIB_Layer1,badStripList,inSiStripQuality);
291  AnalyzeOccupancy(siStripQuality,medianValues_TIB_Layer2,MeanAndRms_TIB_Layer2,badStripList,inSiStripQuality);
292  AnalyzeOccupancy(siStripQuality,medianValues_TIB_Layer3,MeanAndRms_TIB_Layer3,badStripList,inSiStripQuality);
293  AnalyzeOccupancy(siStripQuality,medianValues_TIB_Layer4,MeanAndRms_TIB_Layer4,badStripList,inSiStripQuality);
294 
295  AnalyzeOccupancy(siStripQuality,medianValues_TOB_Layer1,MeanAndRms_TOB_Layer1,badStripList,inSiStripQuality);
296  AnalyzeOccupancy(siStripQuality,medianValues_TOB_Layer2,MeanAndRms_TOB_Layer2,badStripList,inSiStripQuality);
297  AnalyzeOccupancy(siStripQuality,medianValues_TOB_Layer3,MeanAndRms_TOB_Layer3,badStripList,inSiStripQuality);
298  AnalyzeOccupancy(siStripQuality,medianValues_TOB_Layer4,MeanAndRms_TOB_Layer4,badStripList,inSiStripQuality);
299  AnalyzeOccupancy(siStripQuality,medianValues_TOB_Layer5,MeanAndRms_TOB_Layer5,badStripList,inSiStripQuality);
300  AnalyzeOccupancy(siStripQuality,medianValues_TOB_Layer6,MeanAndRms_TOB_Layer6,badStripList,inSiStripQuality);
301 
302  AnalyzeOccupancy(siStripQuality,medianValues_TIDPlus_Disc1,MeanAndRms_TIDPlus_Disc1,badStripList,inSiStripQuality);
303  AnalyzeOccupancy(siStripQuality,medianValues_TIDPlus_Disc2,MeanAndRms_TIDPlus_Disc2,badStripList,inSiStripQuality);
304  AnalyzeOccupancy(siStripQuality,medianValues_TIDPlus_Disc3,MeanAndRms_TIDPlus_Disc3,badStripList,inSiStripQuality);
305  AnalyzeOccupancy(siStripQuality,medianValues_TIDMinus_Disc1,MeanAndRms_TIDMinus_Disc1,badStripList,inSiStripQuality);
306  AnalyzeOccupancy(siStripQuality,medianValues_TIDMinus_Disc2,MeanAndRms_TIDMinus_Disc2,badStripList,inSiStripQuality);
307  AnalyzeOccupancy(siStripQuality,medianValues_TIDMinus_Disc3,MeanAndRms_TIDMinus_Disc3,badStripList,inSiStripQuality);
308 
309  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc1,MeanAndRms_TECPlus_Disc1,badStripList,inSiStripQuality);
310  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc2,MeanAndRms_TECPlus_Disc2,badStripList,inSiStripQuality);
311  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc3,MeanAndRms_TECPlus_Disc3,badStripList,inSiStripQuality);
312  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc4,MeanAndRms_TECPlus_Disc4,badStripList,inSiStripQuality);
313  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc5,MeanAndRms_TECPlus_Disc5,badStripList,inSiStripQuality);
314  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc6,MeanAndRms_TECPlus_Disc6,badStripList,inSiStripQuality);
315  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc7,MeanAndRms_TECPlus_Disc7,badStripList,inSiStripQuality);
316  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc8,MeanAndRms_TECPlus_Disc8,badStripList,inSiStripQuality);
317  AnalyzeOccupancy(siStripQuality,medianValues_TECPlus_Disc9,MeanAndRms_TECPlus_Disc9,badStripList,inSiStripQuality);
318 
319  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc1,MeanAndRms_TECMinus_Disc1,badStripList,inSiStripQuality);
320  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc2,MeanAndRms_TECMinus_Disc2,badStripList,inSiStripQuality);
321  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc3,MeanAndRms_TECMinus_Disc3,badStripList,inSiStripQuality);
322  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc4,MeanAndRms_TECMinus_Disc4,badStripList,inSiStripQuality);
323  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc5,MeanAndRms_TECMinus_Disc5,badStripList,inSiStripQuality);
324  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc6,MeanAndRms_TECMinus_Disc6,badStripList,inSiStripQuality);
325  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc7,MeanAndRms_TECMinus_Disc7,badStripList,inSiStripQuality);
326  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc8,MeanAndRms_TECMinus_Disc8,badStripList,inSiStripQuality);
327  AnalyzeOccupancy(siStripQuality,medianValues_TECMinus_Disc9,MeanAndRms_TECMinus_Disc9,badStripList,inSiStripQuality);
328 
329  siStripQuality->fillBadComponents();
330 
331  if (WriteOutputFile_==true){
332  f->cd();
333  apvtree->Write();
334  f->Close();
335  }
336 
337  LogTrace("SiStripBadAPV") << ss.str() << std::endl;
338 }
339 
340 
341 void SiStripBadAPVAlgorithmFromClusterOccupancy::CalculateMeanAndRMS(const std::vector<Apv>& a, std::pair<double,double>* MeanRMS, int number_iterations)
342 {
343  Double_t tot[7], tot2[7];
344  Double_t n[7];
345 
346  Double_t Mean[7] = {0};
347  Double_t Rms[7] = {1000,1000,1000,1000,1000,1000,1000};
348 
349  int Moduleposition;
350 
351  for (int i=0; i<number_iterations; i++)
352  {
353  for (int j=0; j<7; j++)
354  {
355  n[j] = 0;
356  tot[j] = 0;
357  tot2[j] = 0;
358  }
359 
360  for (uint32_t it=0; it<a.size(); it++)
361  {
362  Moduleposition = (a[it].modulePosition)-1;
363 
364  for (int apv=0; apv<a[it].numberApvs; apv++)
365  {
366  if (i>0)
367  {
368  if (a[it].apvMedian[apv]<(Mean[Moduleposition]-3*Rms[Moduleposition]) || (a[it].apvMedian[apv]>(Mean[Moduleposition]+5*Rms[Moduleposition])))
369  {
370  continue;
371  }
372  }
373  tot[Moduleposition] += a[it].apvMedian[apv];
374  tot2[Moduleposition] += (a[it].apvMedian[apv])*(a[it].apvMedian[apv]);
375  n[Moduleposition]++;
376  }
377  }
378 
379  for (int j=0; j<7; j++)
380  {
381  if (n[j]!=0)
382  {
383  Mean[j] = tot[j]/n[j];
384  Rms[j] = TMath::Sqrt(TMath::Abs(tot2[j]/n[j] -Mean[j]*Mean[j]));
385  }
386  }
387  }
388 
389  for (int j=0; j<7; j++)
390  {
391  MeanRMS[j] = std::make_pair(Mean[j],Rms[j]);
392  }
393 
394 }
395 
396 void SiStripBadAPVAlgorithmFromClusterOccupancy::AnalyzeOccupancy(SiStripQuality* quality, std::vector<Apv>& medianValues, std::pair<double,double>* MeanAndRms, std::vector<unsigned int>& BadStripList, edm::ESHandle<SiStripQuality>& InSiStripQuality)
397 {
398  int Moduleposition;
399  uint32_t Detid;
400 
401  for (uint32_t it=0; it<medianValues.size(); it++)
402  {
403  Moduleposition = (medianValues[it].modulePosition)-1;
404  Detid = medianValues[it].detrawId;
405 
406  for (int apv=0; apv<medianValues[it].numberApvs; apv++)
407  {
408  if(UseInputDB_)
409  {
410  if(InSiStripQuality->IsApvBad(Detid,apv) )
411  {
412  continue;//if the apv is already flagged as bad, continue.
413  }
414  }
415  if (medianValues[it].apvMedian[apv] > minNevents_)
416  {
417  if ((medianValues[it].apvMedian[apv]>(MeanAndRms[Moduleposition].first+highoccupancy_*MeanAndRms[Moduleposition].second)) && (medianValues[it].apvMedian[apv]>absolutelow_))
418  BadStripList.push_back(pQuality->encode((apv*128),128,0));
419  }
420  else if (medianValues[it].apvMedian[apv]<(MeanAndRms[Moduleposition].first-lowoccupancy_*MeanAndRms[Moduleposition].second) && (MeanAndRms[Moduleposition].first>2 || medianValues[it].apvabsoluteOccupancy[apv]==0))
421  {
422  BadStripList.push_back(pQuality->encode((apv*128),128,0));
423  std::cout << "Dead APV! DetId: " << medianValues[it].detrawId << ", APV number: " << apv+1 << ", APVMedian: " << medianValues[it].apvMedian[apv] << ", Mean: " << MeanAndRms[Moduleposition].first << ", RMS: " << MeanAndRms[Moduleposition].second << ", LowThreshold: " << lowoccupancy_ << ", Mean-Low*RMS: " << (MeanAndRms[Moduleposition].first-lowoccupancy_*MeanAndRms[Moduleposition].second) << std::endl;
424  }
425  }
426  if (BadStripList.begin()!=BadStripList.end())
427  {
428  quality->compact(Detid,BadStripList);
429  SiStripQuality::Range range(BadStripList.begin(),BadStripList.end());
430  quality->put(Detid,range);
431  }
432  BadStripList.clear();
433  }
434 }
435 
436 void SiStripBadAPVAlgorithmFromClusterOccupancy::setMinNumOfEvents()
437 {
438  minNevents_=occupancy_*Nevents_;
439 }
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:122
TrackerTopology::tibString
unsigned int tibString(const DetId &id) const
Definition: TrackerTopology.h:324
TrackerTopology::tidRing
unsigned int tidRing(const DetId &id) const
Definition: TrackerTopology.h:203
TrackerTopology::tobIsStereo
bool tobIsStereo(const DetId &id) const
Definition: TrackerTopology.h:241
SiStripBadAPVAlgorithmFromClusterOccupancy::AnalyzeOccupancy
void AnalyzeOccupancy(SiStripQuality *, std::vector< Apv > &, std::pair< double, double > *, std::vector< unsigned int > &, edm::ESHandle< SiStripQuality > &)
Definition: SiStripBadAPVAlgorithmFromClusterOccupancy.cc:396
SiStripBadAPVAlgorithmFromClusterOccupancy::CalculateMeanAndRMS
void CalculateMeanAndRMS(const std::vector< Apv > &, std::pair< double, double > *, int)
Definition: SiStripBadAPVAlgorithmFromClusterOccupancy.cc:341
TrackerTopology::tecRing
unsigned int tecRing(const DetId &id) const
ring id
Definition: TrackerTopology.h:200
SiStripBadAPVAlgorithmFromClusterOccupancy::SiStripBadAPVAlgorithmFromClusterOccupancy
SiStripBadAPVAlgorithmFromClusterOccupancy(const edm::ParameterSet &, const TrackerTopology *)
Definition: SiStripBadAPVAlgorithmFromClusterOccupancy.cc:13
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
TrackerTopology::tidWheel
unsigned int tidWheel(const DetId &id) const
Definition: TrackerTopology.h:183
StripGeomDetUnit::specificTopology
virtual const StripTopology & specificTopology() const
Returns a reference to the strip proxy topology.
Definition: StripGeomDetUnit.cc:18
TrackerTopology::tecIsZMinusSide
bool tecIsZMinusSide(const DetId &id) const
Definition: TrackerTopology.h:238
TrackerTopology::tidIsStereo
bool tidIsStereo(const DetId &id) const
Definition: TrackerTopology.h:244
TrackerTopology::tidIsZMinusSide
bool tidIsZMinusSide(const DetId &id) const
Definition: TrackerTopology.h:235
sistrip::extrainfo::occupancy_
static const char occupancy_[]
Definition: ConstantsForDqm.h:60
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:248
TrackerTopology::tecIsStereo
bool tecIsStereo(const DetId &id) const
Definition: TrackerTopology.h:242
TrackerTopology::tibIsExternalString
bool tibIsExternalString(const DetId &id) const
Definition: TrackerTopology.h:340
TrackerTopology::tibIsZMinusSide
bool tibIsZMinusSide(const DetId &id) const
Definition: TrackerTopology.h:232
TrackerTopology::tidModule
unsigned int tidModule(const DetId &id) const
Definition: TrackerTopology.h:148
TrackerTopology::tidIsBackRing
bool tidIsBackRing(const DetId &id) const
Definition: TrackerTopology.h:352
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
SiStripQuality::compact
void compact(unsigned int &, std::vector< unsigned int > &)
Definition: SiStripQuality.cc:336
Abs
T Abs(T a)
Definition: MathUtil.h:49
TrackerTopology::tobIsZMinusSide
bool tobIsZMinusSide(const DetId &id) const
Definition: TrackerTopology.h:229
j
int j
Definition: DBlmapReader.cc:9
SiStripQuality::fillBadComponents
void fillBadComponents()
Definition: SiStripQuality.cc:530
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
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
TrackerTopology::tibModule
unsigned int tibModule(const DetId &id) const
Definition: TrackerTopology.h:145
TrackerTopology::tecModule
unsigned int tecModule(const DetId &id) const
Definition: TrackerTopology.h:142
TrackerTopology::tecIsBackPetal
bool tecIsBackPetal(const DetId &id) const
Definition: TrackerTopology.h:317
DetId
Definition: DetId.h:18
TrackerTopology::tibIsStereo
bool tibIsStereo(const DetId &id) const
Definition: TrackerTopology.h:243
TrackerTopology::tobModule
unsigned int tobModule(const DetId &id) const
Definition: TrackerTopology.h:139
a
double a
Definition: hdecay.h:121
SiStripBadStrip::Range
std::pair< ContainerIterator, ContainerIterator > Range
Definition: SiStripBadStrip.h:54
SiStripBadAPVAlgorithmFromClusterOccupancy::extractBadAPVs
void extractBadAPVs(SiStripQuality *, HistoMap &, edm::ESHandle< SiStripQuality > &)
Definition: SiStripBadAPVAlgorithmFromClusterOccupancy.cc:31
HLT_25ns14e33_v1_cff.quality
tuple quality
Definition: HLT_25ns14e33_v1_cff.py:33032
StripTopology::localPosition
virtual LocalPoint localPosition(float strip) const =0
TrackerTopology::tecPetalNumber
unsigned int tecPetalNumber(const DetId &id) const
Definition: TrackerTopology.h:209
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
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:175
SiStripBadStrip::put
bool put(const uint32_t &detID, const InputVector &vect)
Definition: SiStripBadStrip.h:66
citk::for
for(const auto &isodef:isoDefs)
Definition: CITKPFIsolationSumProducer.cc:77
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
TrackerTopology::tecWheel
unsigned int tecWheel(const DetId &id) const
Definition: TrackerTopology.h:180
SiStripBadStrip::encode
unsigned int encode(const unsigned short &first, const unsigned short &NconsecutiveBadStrips, const unsigned short &flag=0)
Definition: SiStripBadStrip.h:88
TrackerTopology::tobLayer
unsigned int tobLayer(const DetId &id) const
Definition: TrackerTopology.h:119
TrackerGeometry::idToDet
virtual const TrackerGeomDet * idToDet(DetId) const
Definition: TrackerGeometry.cc:178
Generated for CMSSW Reference Manual by doxygen 1.8.5