CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HcalDetDiagLEDClient.cc
Go to the documentation of this file.
4 
8 
10 
11 #include <iostream>
12 
13 /*
14  * \file HcalDetDiagLEDClient.cc
15  *
16  * $Date: 2012/11/12 23:26:19 $
17  * $Revision: 1.11 $
18  * \author J. Temple
19  * \brief Hcal DetDiagLED Client class
20  */
21 
23  name_=myname; status=0;
24  needLogicalMap_=true;
25 }
26 
28  name_=myname;
29  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup",false);
30  debug_ = ps.getUntrackedParameter<int>("debug",0);
31  prefixME_ = ps.getUntrackedParameter<std::string>("subSystemFolder","Hcal/");
32  if (prefixME_.substr(prefixME_.size()-1,prefixME_.size())!="/")
33  prefixME_.append("/");
34  subdir_ = ps.getUntrackedParameter<std::string>("DetDiagLEDFolder","DetDiagLEDMonitor_Hcal/");
35  if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
36  subdir_.append("/");
37  subdir_=prefixME_+subdir_;
38 
39  validHtmlOutput_ = ps.getUntrackedParameter<bool>("DetDiagLED_validHtmlOutput",true);
40  cloneME_ = ps.getUntrackedParameter<bool>("cloneME", true);
41  badChannelStatusMask_ = ps.getUntrackedParameter<int>("DetDiagLED_BadChannelStatusMask",
42  ps.getUntrackedParameter<int>("BadChannelStatusMask",(1<<HcalChannelStatus::HcalCellDead)));
43  needLogicalMap_=true;
44 }
45 
47  if (debug_>2) std::cout <<"\tHcalDetDiagLEDClient::analyze()"<<std::endl;
49 }
50 
52 void HcalDetDiagLEDClient::updateChannelStatus(std::map<HcalDetId, unsigned int>& myqual){
53  // This gets called by HcalMonitorClient
54  // trigger primitives don't yet contribute to channel status (though they could...)
55  // see dead or hot cell code for an example
56 
57 }
58 
61  if (debug_>0){
62  std::cout <<"<HcalDetDiagLEDClient::beginJob()> Displaying dqmStore directory structure:"<<std::endl;
64  }
65 }
67 
69  if (!dqmStore_)
70  {
71  if (debug_>0) std::cout <<"<HcalDetDiagLEDClient::beginRun> dqmStore does not exist!"<<std::endl;
72  return;
73  }
75  problemnames_.clear();
76 
77  // Put the appropriate name of your problem summary here
78  ProblemCells=dqmStore_->book2D(" ProblemDetDiagLED",
79  " Problem DetDiagLED Rate for all HCAL;ieta;iphi",
80  85,-42.5,42.5,
81  72,0.5,72.5);
82  problemnames_.push_back(ProblemCells->getName());
83  if (debug_>1)
84  std::cout << "Tried to create ProblemCells Monitor Element in directory "<<subdir_<<" \t Failed? "<<(ProblemCells==0)<<std::endl;
85  dqmStore_->setCurrentFolder(subdir_+"problem_DetDiagLED");
87  ProblemCellsByDepth->setup(dqmStore_," Problem DetDiagLED Rate");
88  for (unsigned int i=0; i<ProblemCellsByDepth->depth.size();++i)
89  problemnames_.push_back(ProblemCellsByDepth->depth[i]->getName());
90  nevts_=0;
91 }
95 
97  if(status&2) return true;
98  return false;
99 }
101  if(status&1) return true;
102  return false;
103 }
104 bool HcalDetDiagLEDClient::hasOther_Temp(void){return false;}
105 bool HcalDetDiagLEDClient::test_enabled(void){return true;}
106 
108  std::string s=subdir_+"HcalDetDiagLEDMonitor Event Number";
109  MonitorElement *me = dqmStore_->get(s.c_str());
110  int n=0;
111  if ( me ) {
112  s = me->valueString();
113  sscanf((s.substr(2,s.length()-2)).c_str(), "%d", &n);
114  }
115  if(n<100) return false;
116  return true;
117 }
118 static void printTableHeader(ofstream& file,std::string header){
119  file << "</html><html xmlns=\"http://www.w3.org/1999/xhtml\">"<< std::endl;
120  file << "<head>"<< std::endl;
121  file << "<meta http-equiv=\"Content-Type\" content=\"text/html\"/>"<< std::endl;
122  file << "<title>"<< header <<"</title>"<< std::endl;
123  file << "<style type=\"text/css\">"<< std::endl;
124  file << " body,td{ background-color: #FFFFCC; font-family: arial, arial ce, helvetica; font-size: 12px; }"<< std::endl;
125  file << " td.s0 { font-family: arial, arial ce, helvetica; }"<< std::endl;
126  file << " td.s1 { font-family: arial, arial ce, helvetica; font-weight: bold; background-color: #FFC169; text-align: center;}"<< std::endl;
127  file << " td.s2 { font-family: arial, arial ce, helvetica; background-color: #eeeeee; }"<< std::endl;
128  file << " td.s3 { font-family: arial, arial ce, helvetica; background-color: #d0d0d0; }"<< std::endl;
129  file << " td.s4 { font-family: arial, arial ce, helvetica; background-color: #FFC169; }"<< std::endl;
130  file << "</style>"<< std::endl;
131  file << "<body>"<< std::endl;
132  file << "<table>"<< std::endl;
133 }
134 static void printTableLine(ofstream& file,int ind,HcalDetId& detid,HcalFrontEndId& lmap_entry,HcalElectronicsId &emap_entry, std::string comment=""){
135  if(ind==0){
136  file << "<tr>";
137  file << "<td class=\"s4\" align=\"center\">#</td>" << std::endl;
138  file << "<td class=\"s1\" align=\"center\">ETA</td>" << std::endl;
139  file << "<td class=\"s1\" align=\"center\">PHI</td>" << std::endl;
140  file << "<td class=\"s1\" align=\"center\">DEPTH</td>"<< std::endl;
141  file << "<td class=\"s1\" align=\"center\">RBX</td>" << std::endl;
142  file << "<td class=\"s1\" align=\"center\">RM</td>" << std::endl;
143  file << "<td class=\"s1\" align=\"center\">PIXEL</td>" << std::endl;
144  file << "<td class=\"s1\" align=\"center\">RM_FIBER</td>" << std::endl;
145  file << "<td class=\"s1\" align=\"center\">FIBER_CH</td>" << std::endl;
146  file << "<td class=\"s1\" align=\"center\">QIE</td>" << std::endl;
147  file << "<td class=\"s1\" align=\"center\">ADC</td>" << std::endl;
148  file << "<td class=\"s1\" align=\"center\">CRATE</td>" << std::endl;
149  file << "<td class=\"s1\" align=\"center\">DCC</td>" << std::endl;
150  file << "<td class=\"s1\" align=\"center\">SPIGOT</td>" << std::endl;
151  file << "<td class=\"s1\" align=\"center\">HTR_FIBER</td>" << std::endl;
152  file << "<td class=\"s1\" align=\"center\">HTR_SLOT</td>" << std::endl;
153  file << "<td class=\"s1\" align=\"center\">HTR_FPGA</td>" << std::endl;
154  if(comment[0]!=0) file << "<td class=\"s1\" align=\"center\">Comment</td>" << std::endl;
155  file << "</tr>" << std::endl;
156  }
157  std::string raw_class;
158  file << "<tr>"<< std::endl;
159  if((ind%2)==1){
160  raw_class="<td class=\"s2\" align=\"center\">";
161  }else{
162  raw_class="<td class=\"s3\" align=\"center\">";
163  }
164  file << "<td class=\"s4\" align=\"center\">" << ind+1 <<"</td>"<< std::endl;
165  file << raw_class<< detid.ieta()<<"</td>"<< std::endl;
166  file << raw_class<< detid.iphi()<<"</td>"<< std::endl;
167  file << raw_class<< detid.depth() <<"</td>"<< std::endl;
168  file << raw_class<< lmap_entry.rbx()<<"</td>"<< std::endl;
169  file << raw_class<< lmap_entry.rm() <<"</td>"<< std::endl;
170  file << raw_class<< lmap_entry.pixel()<<"</td>"<< std::endl;
171  file << raw_class<< lmap_entry.rmFiber() <<"</td>"<< std::endl;
172  file << raw_class<< lmap_entry.fiberChannel()<<"</td>"<< std::endl;
173  file << raw_class<< lmap_entry.qieCard() <<"</td>"<< std::endl;
174  file << raw_class<< lmap_entry.adc()<<"</td>"<< std::endl;
175  file << raw_class<< emap_entry.readoutVMECrateId()<<"</td>"<< std::endl;
176  file << raw_class<< emap_entry.dccid()<<"</td>"<< std::endl;
177  file << raw_class<< emap_entry.spigot()<<"</td>"<< std::endl;
178  file << raw_class<< emap_entry.fiberIndex()<<"</td>"<< std::endl;
179  file << raw_class<< emap_entry.htrSlot()<<"</td>"<< std::endl;
180  file << raw_class<< emap_entry.htrTopBottom()<<"</td>"<< std::endl;
181  if(comment[0]!=0) file << raw_class<< comment<<"</td>"<< std::endl;
182 }
183 static void printTableTail(ofstream& file){
184  file << "</table>"<< std::endl;
185  file << "</body>"<< std::endl;
186  file << "</html>"<< std::endl;
187 }
188 double HcalDetDiagLEDClient::get_channel_status(std::string subdet,int eta,int phi,int depth,int type){
189  int subdetval=-1;
190  if (subdet.compare("HB")==0) subdetval=(int)HcalBarrel;
191  else if (subdet.compare("HE")==0) subdetval=(int)HcalEndcap;
192  else if (subdet.compare("HO")==0) subdetval=(int)HcalOuter;
193  else if (subdet.compare("HF")==0) subdetval=(int)HcalForward;
194  else return -1.0;
195  int ietabin=CalcEtaBin(subdetval, eta, depth)+1;
196  if(type==1) return ChannelStatusMissingChannels[depth-1]->GetBinContent(ietabin,phi);
197  if(type==2) return ChannelStatusUnstableChannels[depth-1]->GetBinContent(ietabin,phi);
198  if(type==3) return ChannelStatusUnstableLEDsignal[depth-1]->GetBinContent(ietabin,phi);
199  if(type==4) return ChannelStatusLEDMean[depth-1]->GetBinContent(ietabin,phi);
200  if(type==5) return ChannelStatusLEDRMS[depth-1]->GetBinContent(ietabin,phi);
201  if(type==6) return ChannelStatusTimeMean[depth-1]->GetBinContent(ietabin,phi);
202  if(type==7) return ChannelStatusTimeRMS[depth-1]->GetBinContent(ietabin,phi);
203  return -1.0;
204 }
205 double HcalDetDiagLEDClient::get_energy(std::string subdet,int eta,int phi,int depth,int type){
206  int subdetval=-1;
207  if (subdet.compare("HB")==0) subdetval=(int)HcalBarrel;
208  else if (subdet.compare("HE")==0) subdetval=(int)HcalEndcap;
209  else if (subdet.compare("HO")==0) subdetval=(int)HcalOuter;
210  else if (subdet.compare("HF")==0) subdetval=(int)HcalForward;
211  else return -1.0;
212  int ietabin=CalcEtaBin(subdetval, eta, depth)+1;
213  if(type==1) return ChannelsLEDEnergy[depth-1]->GetBinContent(ietabin,phi);
214  if(type==2) return ChannelsLEDEnergyRef[depth-1]->GetBinContent(ietabin,phi);
215  return -1.0;
216 }
218 MonitorElement* me;
219 int MissingCnt=0;
220 int UnstableCnt=0;
221 int UnstableLEDCnt=0;
222 int BadTimingCnt=0;
223 int HBP[7]={0,0,0,0,0,0,0},newHBP[7]={0,0,0,0,0,0,0};
224 int HBM[7]={0,0,0,0,0,0,0},newHBM[7]={0,0,0,0,0,0,0};
225 int HEP[7]={0,0,0,0,0,0,0},newHEP[7]={0,0,0,0,0,0,0};
226 int HEM[7]={0,0,0,0,0,0,0},newHEM[7]={0,0,0,0,0,0,0};
227 int HFP[7]={0,0,0,0,0,0,0},newHFP[7]={0,0,0,0,0,0,0};
228 int HFM[7]={0,0,0,0,0,0,0},newHFM[7]={0,0,0,0,0,0,0};
229 int HO[7] ={0,0,0,0,0,0,0},newHO[7] ={0,0,0,0,0,0,0};
230 std::string subdet[4]={"HB","HE","HO","HF"};
231 
232  if (debug_>0) std::cout << "<HcalDetDiagLEDClient::htmlOutput> Preparing html output ..." << std::endl;
233  if(!dqmStore_) return;
235  std::vector<std::string> name = HcalEtaPhiHistNames();
236 
237  for(int i=0;i<4;++i){
241  std::string s;
242  s=subdir_+"channel status/"+name[i]+" Missing Channels";
243  me=dqmStore_->get(s.c_str());
244  if (me!=0) ChannelStatusMissingChannels[i]=HcalUtilsClient::getHisto<TH2F*>(me,cloneME_,ChannelStatusMissingChannels[i],debug_); else return;
245  s=subdir_+"channel status/"+name[i]+" Unstable Channels";
246  me=dqmStore_->get(s.c_str());
247  if (me!=0) ChannelStatusUnstableChannels[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, ChannelStatusUnstableChannels[i], debug_); else return;
248  s=subdir_+"channel status/"+name[i]+" Unstable LED";
249  me=dqmStore_->get(s.c_str());
250  if (me!=0) ChannelStatusUnstableLEDsignal[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, ChannelStatusUnstableLEDsignal[i], debug_); else return;
251  s=subdir_+"channel status/"+name[i]+" LED Mean";
252  me=dqmStore_->get(s.c_str());
253  if (me!=0) ChannelStatusLEDMean[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, ChannelStatusLEDMean[i], debug_); else return;
254  s=subdir_+"channel status/"+name[i]+" LED RMS";
255  me=dqmStore_->get(s.c_str());
256  if (me!=0) ChannelStatusLEDRMS[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, ChannelStatusLEDRMS[i], debug_); else return;
257  s=subdir_+"channel status/"+name[i]+" Time Mean";
258  me=dqmStore_->get(s.c_str());
259  if (me!=0) ChannelStatusTimeMean[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, ChannelStatusTimeMean[i], debug_); else return;
260  s=subdir_+"channel status/"+name[i]+" Time RMS";
261  me=dqmStore_->get(s.c_str());
262  if (me!=0) ChannelStatusTimeRMS[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, ChannelStatusTimeRMS[i], debug_); else return;
263  s=subdir_+"Summary Plots/"+name[i]+" Channel LED Energy";
264  me=dqmStore_->get(s.c_str());
265  if (me!=0) ChannelsLEDEnergy[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, ChannelsLEDEnergy[i], debug_); else return;
266  s=subdir_+"Summary Plots/"+name[i]+" Channel LED Energy Reference";
267  me=dqmStore_->get(s.c_str());
268  if (me!=0) ChannelsLEDEnergyRef[i]=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, ChannelsLEDEnergyRef[i], debug_); else return;
269  }
270 
271  // check how many problems we have:
272  for(int sd=0;sd<4;sd++){
273  int feta=0,teta=0,fdepth=0,tdepth=0;
274  if(sd==0){ feta=-16; teta=16 ;fdepth=1; tdepth=2; }
275  if(sd==1){ feta=-29; teta=29 ;fdepth=1; tdepth=3; }
276  if(sd==2){ feta=-15; teta=15 ;fdepth=4; tdepth=4; }
277  if(sd==3){ feta=-42; teta=42 ;fdepth=1; tdepth=2; }
278  for(int phi=1;phi<=72;phi++) for(int depth=fdepth;depth<=tdepth;depth++) for(int eta=feta;eta<=teta;eta++){
279  if(sd==3 && eta>-29 && eta<29) continue;
280  double problem[7]={0,0,0,0,0,0,0};
281  for(int i=0;i<6;i++){
282  problem[i] =get_channel_status(subdet[sd],eta,phi,depth,i+1);
283  if(problem[i]!=0){
285  if(sd==0)SD=HcalBarrel;
286  else if(sd==1) SD=HcalEndcap;
287  else if(sd==2) SD=HcalOuter;
288  else if(sd==3) SD=HcalForward;
289  HcalDetId hcalid(SD, eta, phi, depth);
290  if(sd==0){ if(eta>0){
291  HBP[i]++;
292  if(badstatusmap.find(hcalid)==badstatusmap.end())newHBP[i]++;
293  }else{
294  HBM[i]++;
295  if(badstatusmap.find(hcalid)==badstatusmap.end())newHBM[i]++;
296  }}
297  if(sd==1){ if(eta>0){
298  HEP[i]++;
299  if(badstatusmap.find(hcalid)==badstatusmap.end())newHEP[i]++;
300  }else{
301  HEM[i]++;
302  if(badstatusmap.find(hcalid)==badstatusmap.end())newHEM[i]++;
303  }}
304  if(sd==2){
305  HO[i]++;
306  if(badstatusmap.find(hcalid)==badstatusmap.end())newHO[i]++;
307  }
308  if(sd==3){ if(eta>0){
309  HFP[i]++;
310  if(badstatusmap.find(hcalid)==badstatusmap.end())newHFP[i]++;
311  }else{
312  HFM[i]++;
313  if(badstatusmap.find(hcalid)==badstatusmap.end())newHFM[i]++;
314  }}
315  }
316  }
317  }
318  }
319  // missing channels list
320  ofstream Missing;
321  Missing.open((htmlDir + "Missing.html").c_str());
322  printTableHeader(Missing,"Missing Channels list");
323  // Bad timing channels list
324  ofstream BadTiming;
325  BadTiming.open((htmlDir + "BadTiming.html").c_str());
326  printTableHeader(BadTiming,"Bad Timing Channels list");
327  // unstable channels list
328  ofstream Unstable;
329  Unstable.open((htmlDir + "Unstable.html").c_str());
330  printTableHeader(Unstable,"Low LED signal Channels list");
331  // unstable LED signal list
332  ofstream BadLED;
333  BadLED.open((htmlDir + "UnstableLED.html").c_str());
334  printTableHeader(BadLED,"Unstable LED signal channels list");
335 
336  for(int sd=0;sd<4;sd++){
337  int cnt=0;
338  if(sd==0 && ((HBM[0]+HBP[0])==0 || (HBM[0]+HBP[0])==(1296*2))) continue;
339  if(sd==1 && ((HEM[0]+HEP[0])==0 || (HEM[0]+HEP[0])==(1296*2))) continue;
340  if(sd==2 && ((HO[0])==0 || HO[0]==2160)) continue;
341  if(sd==3 && ((HFM[0]+HFP[0])==0 || (HFM[0]+HFP[0])==(864*2))) continue;
342  Missing << "<tr><td align=\"center\"><h3>"<< subdet[sd] <<"</h3></td></tr>" << std::endl;
343  int feta=0,teta=0,fdepth=0,tdepth=0;
344  if(sd==0){ feta=-16; teta=16 ;fdepth=1; tdepth=2; if(HBM[0]==1296) feta=0; if(HBP[0]==1296) teta=0;}
345  if(sd==1){ feta=-29; teta=29 ;fdepth=1; tdepth=3; if(HEM[0]==1296) feta=0; if(HEP[0]==1296) teta=0;}
346  if(sd==2){ feta=-15; teta=15 ;fdepth=4; tdepth=4; if(HO[0] ==2160) {feta=0; teta=0; }}
347  if(sd==3){ feta=-42; teta=42 ;fdepth=1; tdepth=2; if(HFM[0]==864) feta=0; if(HFP[0]==864) teta=0; }
348  for(int phi=1;phi<=72;phi++) for(int depth=fdepth;depth<=tdepth;depth++) for(int eta=feta;eta<=teta;eta++){
349  if(sd==3 && eta>-29 && eta<29) continue;
350  double missing =get_channel_status(subdet[sd],eta,phi,depth,1);
351  if(missing>0){
352  HcalDetId *detid=0;
353  if(sd==0) detid=new HcalDetId(HcalBarrel,eta,phi,depth);
354  if(sd==1) detid=new HcalDetId(HcalEndcap,eta,phi,depth);
355  if(sd==2) detid=new HcalDetId(HcalOuter,eta,phi,depth);
356  if(sd==3) detid=new HcalDetId(HcalForward,eta,phi,depth);
357  HcalFrontEndId lmap_entry=logicalMap_->getHcalFrontEndId(*detid);
358  HcalElectronicsId emap_entry=emap.lookup(*detid);
359  std::string s=" ";
360  if(badstatusmap.find(*detid)!=badstatusmap.end()){ s="Known problem"; }
361  printTableLine(Missing,cnt++,*detid,lmap_entry,emap_entry,s); MissingCnt++;
362  delete detid;
363  }
364  }
365  }
366 
367  for(int sd=0;sd<4;sd++){
368  int cnt=0;
369  if(sd==0 && (HBM[5]+HBP[5])==0) continue;
370  if(sd==1 && (HEM[5]+HEP[5])==0) continue;
371  if(sd==2 && (HO[5])==0) continue;
372  if(sd==3 && (HFM[5]+HFP[5])==0) continue;
373  BadTiming << "<tr><td align=\"center\"><h3>"<< subdet[sd] <<"</h3></td></tr>" << std::endl;
374  int feta=0,teta=0,fdepth=0,tdepth=0;
375  if(sd==0){ feta=-16; teta=16 ;fdepth=1; tdepth=2; if(HBM[0]==1296) feta=0; if(HBP[0]==1296) teta=0;}
376  if(sd==1){ feta=-29; teta=29 ;fdepth=1; tdepth=3; if(HEM[0]==1296) feta=0; if(HEP[0]==1296) teta=0;}
377  if(sd==2){ feta=-15; teta=15 ;fdepth=4; tdepth=4; if(HO[0] ==2160) {feta=0; teta=0; }}
378  if(sd==3){ feta=-42; teta=42 ;fdepth=1; tdepth=2; if(HFM[0]==864) feta=0; if(HFP[0]==864) teta=0; }
379  for(int phi=1;phi<=72;phi++) for(int depth=fdepth;depth<=tdepth;depth++) for(int eta=feta;eta<=teta;eta++){
380  if(sd==3 && eta>-29 && eta<29) continue;
381  double badtiming =get_channel_status(subdet[sd],eta,phi,depth,6);
382  if(badtiming!=0){
383  try{
384  char comment[100]; sprintf(comment,"Time-mean=%.1f\n",badtiming);
385  HcalDetId *detid=0;
386  if(sd==0) detid=new HcalDetId(HcalBarrel,eta,phi,depth);
387  if(sd==1) detid=new HcalDetId(HcalEndcap,eta,phi,depth);
388  if(sd==2) detid=new HcalDetId(HcalOuter,eta,phi,depth);
389  if(sd==3) detid=new HcalDetId(HcalForward,eta,phi,depth);
390  HcalFrontEndId lmap_entry=logicalMap_->getHcalFrontEndId(*detid);
391  HcalElectronicsId emap_entry=emap.lookup(*detid);
392  printTableLine(BadTiming,cnt++,*detid,lmap_entry,emap_entry,comment); BadTimingCnt++;
393  delete detid;
394  }catch(cms::Exception &){ continue;}
395  }
396  }
397  }
398 
399  for(int sd=0;sd<4;sd++){
400  int cnt=0;
401  if(sd==0 && (HBM[1]+HBP[1])==0) continue;
402  if(sd==1 && (HEM[1]+HEP[1])==0) continue;
403  if(sd==2 && (HO[1])==0) continue;
404  if(sd==3 && (HFM[1]+HFP[1])==0) continue;
405  Unstable << "<tr><td align=\"center\"><h3>"<< subdet[sd] <<"</h3></td></tr>" << std::endl;
406  int feta=0,teta=0,fdepth=0,tdepth=0;
407  if(sd==0){ feta=-16; teta=16 ;fdepth=1; tdepth=2;}
408  if(sd==1){ feta=-29; teta=29 ;fdepth=1; tdepth=3;}
409  if(sd==2){ feta=-15; teta=15 ;fdepth=4; tdepth=4;}
410  if(sd==3){ feta=-42; teta=42 ;fdepth=1; tdepth=2;}
411  for(int phi=1;phi<=72;phi++) for(int depth=fdepth;depth<=tdepth;depth++) for(int eta=feta;eta<=teta;eta++){
412  if(sd==3 && eta>-29 && eta<29) continue;
413  double unstable =get_channel_status(subdet[sd],eta,phi,depth,2);
414  if(unstable>0){
415  try{
416  char comment[100]; sprintf(comment,"%.3f%%\n",(1.0-unstable)*100.0);
417  HcalDetId *detid=0;
418  if(sd==0) detid=new HcalDetId(HcalBarrel,eta,phi,depth);
419  if(sd==1) detid=new HcalDetId(HcalEndcap,eta,phi,depth);
420  if(sd==2) detid=new HcalDetId(HcalOuter,eta,phi,depth);
421  if(sd==3) detid=new HcalDetId(HcalForward,eta,phi,depth);
422  HcalFrontEndId lmap_entry=logicalMap_->getHcalFrontEndId(*detid);
423  HcalElectronicsId emap_entry=emap.lookup(*detid);
424  printTableLine(Unstable,cnt++,*detid,lmap_entry,emap_entry,comment); UnstableCnt++;
425  delete detid;
426  }catch(cms::Exception &){ continue;}
427  }
428  }
429  }
430 
431  for(int sd=0;sd<4;sd++){
432  int cnt=0;
433  if(sd==0 && (HBM[2]+HBP[2])==0) continue;
434  if(sd==1 && (HEM[2]+HEP[2])==0) continue;
435  if(sd==2 && (HO[2])==0) continue;
436  if(sd==3 && (HFM[2]+HFP[2])==0) continue;
437  BadLED << "<tr><td align=\"center\"><h3>"<< subdet[sd] <<"</h3></td></tr>" << std::endl;
438  int feta=0,teta=0,fdepth=0,tdepth=0;
439  if(sd==0){ feta=-16; teta=16 ;fdepth=1; tdepth=2;}
440  if(sd==1){ feta=-29; teta=29 ;fdepth=1; tdepth=3;}
441  if(sd==2){ feta=-15; teta=15 ;fdepth=4; tdepth=4;}
442  if(sd==3){ feta=-42; teta=42 ;fdepth=1; tdepth=2;}
443  for(int phi=1;phi<=72;phi++) for(int depth=fdepth;depth<=tdepth;depth++) for(int eta=feta;eta<=teta;eta++){
444  if(sd==3 && eta>-29 && eta<29) continue;
445  double badled =get_channel_status(subdet[sd],eta,phi,depth,3);
446  if(badled>0){
447  try{
448  char comment[100]; sprintf(comment,"%.3f%%\n",(badled)*100.0);
449  HcalDetId *detid=0;
450  if(sd==0) detid=new HcalDetId(HcalBarrel,eta,phi,depth);
451  if(sd==1) detid=new HcalDetId(HcalEndcap,eta,phi,depth);
452  if(sd==2) detid=new HcalDetId(HcalOuter,eta,phi,depth);
453  if(sd==3) detid=new HcalDetId(HcalForward,eta,phi,depth);
454  HcalFrontEndId lmap_entry=logicalMap_->getHcalFrontEndId(*detid);
455  HcalElectronicsId emap_entry=emap.lookup(*detid);
456  std::string s=" ";
457  if(badstatusmap.find(*detid)!=badstatusmap.end()){ s="Known problem"; }
458  printTableLine(BadLED,cnt++,*detid,lmap_entry,emap_entry,s); UnstableLEDCnt++;
459  delete detid;
460  }catch(cms::Exception &){ continue;}
461  }
462  }
463  }
464  printTableTail(Missing);
465  Missing.close();
466  printTableTail(BadTiming);
467  BadTiming.close();
468  printTableTail(Unstable);
469  Unstable.close();
470  printTableTail(BadLED);
471  BadLED.close();
473  int ievt_ = -1,runNo=-1;
474  std::string ref_run;
475  std::string s=subdir_+"HcalDetDiagLEDMonitor Event Number";
476  me = dqmStore_->get(s.c_str());
477  if ( me ) {
478  s = me->valueString();
479  sscanf((s.substr(2,s.length()-2)).c_str(), "%d", &ievt_);
480  }
481  s=subdir_+"HcalDetDiagLEDMonitor Run Number";
482  me = dqmStore_->get(s.c_str());
483  if ( me ) {
484  s = me->valueString();
485  sscanf((s.substr(2,s.length()-2)).c_str(), "%d", &runNo);
486  }
487  s=subdir_+"HcalDetDiagLEDMonitor Reference Run";
488  me = dqmStore_->get(s.c_str());
489  if(me) {
490  std::string s=me->valueString();
491  char str[200];
492  sscanf((s.substr(2,s.length()-2)).c_str(), "%s", str);
493  ref_run=str;
494  }
495  TH1F *Energy=0,*Timing=0,*EnergyHF=0,*TimingHF=0,*EnergyRMS=0,*TimingRMS=0,*EnergyRMSHF=0,*TimingRMSHF=0;
496  TH2F *Time2Dhbhehf=0,*Time2Dho=0,*Energy2Dhbhehf=0,*Energy2Dho=0;
497  TH2F *HBPphi=0,*HBMphi=0,*HEPphi=0,*HEMphi=0,*HFPphi=0,*HFMphi=0,*HO0phi=0,*HO1Pphi=0,*HO2Pphi=0,*HO1Mphi=0,*HO2Mphi=0;
498 
499  s=subdir_+"Summary Plots/HBHEHO LED Energy Distribution"; me=dqmStore_->get(s.c_str());
500  if(me!=0) Energy=HcalUtilsClient::getHisto<TH1F*>(me, cloneME_, Energy, debug_); else return;
501  s=subdir_+"Summary Plots/HBHEHO LED Timing Distribution"; me=dqmStore_->get(s.c_str());
502  if(me!=0) Timing=HcalUtilsClient::getHisto<TH1F*>(me, cloneME_, Timing, debug_); else return;
503  s=subdir_+"Summary Plots/HBHEHO LED Energy RMS_div_Energy Distribution"; me=dqmStore_->get(s.c_str());
504  if(me!=0) EnergyRMS=HcalUtilsClient::getHisto<TH1F*>(me, cloneME_, EnergyRMS, debug_); else return;
505  s=subdir_+"Summary Plots/HBHEHO LED Timing RMS Distribution"; me=dqmStore_->get(s.c_str());
506  if(me!=0) TimingRMS=HcalUtilsClient::getHisto<TH1F*>(me, cloneME_, TimingRMS, debug_); else return;
507  s=subdir_+"Summary Plots/HF LED Energy Distribution"; me=dqmStore_->get(s.c_str());
508  if(me!=0) EnergyHF=HcalUtilsClient::getHisto<TH1F*>(me, cloneME_, EnergyHF, debug_); else return;
509  s=subdir_+"Summary Plots/HF LED Timing Distribution"; me=dqmStore_->get(s.c_str());
510  if(me!=0) TimingHF=HcalUtilsClient::getHisto<TH1F*>(me, cloneME_, TimingHF, debug_); else return;
511  s=subdir_+"Summary Plots/HF LED Energy RMS_div_Energy Distribution"; me=dqmStore_->get(s.c_str());
512  if(me!=0) EnergyRMSHF=HcalUtilsClient::getHisto<TH1F*>(me, cloneME_, EnergyRMSHF, debug_); else return;
513  s=subdir_+"Summary Plots/HF LED Timing RMS Distribution"; me=dqmStore_->get(s.c_str());
514  if(me!=0) TimingRMSHF=HcalUtilsClient::getHisto<TH1F*>(me, cloneME_, TimingRMSHF, debug_); else return;
515 
516  s=subdir_+"Summary Plots/LED Timing HBHEHF"; me=dqmStore_->get(s.c_str());
517  if(me!=0) Time2Dhbhehf=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, Time2Dhbhehf, debug_); else return;
518  s=subdir_+"Summary Plots/LED Timing HO"; me=dqmStore_->get(s.c_str());
519  if(me!=0) Time2Dho=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, Time2Dho, debug_); else return;
520  s=subdir_+"Summary Plots/LED Energy HBHEHF"; me=dqmStore_->get(s.c_str());
521  if(me!=0) Energy2Dhbhehf=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, Energy2Dhbhehf, debug_); else return;
522  s=subdir_+"Summary Plots/LED Energy HO"; me=dqmStore_->get(s.c_str());
523  if(me!=0) Energy2Dho=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, Energy2Dho, debug_); else return;
524 
525  s=subdir_+"Summary Plots/HBP Average over HPD LED Ref"; me=dqmStore_->get(s.c_str());
526  if(me!=0) HBPphi=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, HBPphi, debug_); else return;
527  s=subdir_+"Summary Plots/HBM Average over HPD LED Ref"; me=dqmStore_->get(s.c_str());
528  if(me!=0) HBMphi=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, HBMphi, debug_); else return;
529  s=subdir_+"Summary Plots/HEP Average over HPD LED Ref"; me=dqmStore_->get(s.c_str());
530  if(me!=0) HEPphi=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, HEPphi, debug_); else return;
531  s=subdir_+"Summary Plots/HEM Average over HPD LED Ref"; me=dqmStore_->get(s.c_str());
532  if(me!=0) HEMphi=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, HEMphi, debug_); else return;
533  s=subdir_+"Summary Plots/HFP Average over RM LED Ref"; me=dqmStore_->get(s.c_str());
534  if(me!=0) HFPphi=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, HFPphi, debug_); else return;
535  s=subdir_+"Summary Plots/HFM Average over RM LED Ref"; me=dqmStore_->get(s.c_str());
536  if(me!=0) HFMphi=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, HFMphi, debug_); else return;
537 
538  s=subdir_+"Summary Plots/HO0 Average over HPD LED Ref"; me=dqmStore_->get(s.c_str());
539  if(me!=0) HO0phi=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, HO0phi, debug_); else return;
540  s=subdir_+"Summary Plots/HO1P Average over HPD LED Ref"; me=dqmStore_->get(s.c_str());
541  if(me!=0) HO1Pphi=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, HO1Pphi, debug_); else return;
542  s=subdir_+"Summary Plots/HO2P Average over HPD LED Ref"; me=dqmStore_->get(s.c_str());
543  if(me!=0) HO2Pphi=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, HO2Pphi, debug_); else return;
544  s=subdir_+"Summary Plots/HO1M Average over HPD LED Ref"; me=dqmStore_->get(s.c_str());
545  if(me!=0) HO1Mphi=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, HO1Mphi, debug_); else return;
546  s=subdir_+"Summary Plots/HO2M Average over HPD LED Ref"; me=dqmStore_->get(s.c_str());
547  if(me!=0) HO2Mphi=HcalUtilsClient::getHisto<TH2F*>(me, cloneME_, HO2Mphi, debug_); else return;
548 
549  gROOT->SetBatch(true);
550  gStyle->SetCanvasColor(0);
551  gStyle->SetPadColor(0);
552  gStyle->SetOptStat(111110);
553  gStyle->SetPalette(1);
554 
555  TCanvas *can=new TCanvas("HcalDetDiagLEDClient","HcalDetDiagLEDClient",0,0,500,350);
556  can->SetGridy();
557  can->SetGridx();
558  can->cd();
559 
560  ofstream htmlFile;
561  std::string outfile=htmlDir+name_+".html";
562  htmlFile.open(outfile.c_str());
563  // html page header
564  htmlFile << "<!DOCTYPE html PUBLIC \"-//W3C//DTD HTML 4.01 Transitional//EN\"> " << std::endl;
565  htmlFile << "<html> " << std::endl;
566  htmlFile << "<head> " << std::endl;
567  htmlFile << " <meta content=\"text/html; charset=ISO-8859-1\" " << std::endl;
568  htmlFile << " http-equiv=\"content-type\"> " << std::endl;
569  htmlFile << " <title>Detector Diagnostics LED Monitor</title> " << std::endl;
570  htmlFile << "</head> " << std::endl;
571  htmlFile << "<style type=\"text/css\"> td { font-weight: bold } </style>" << std::endl;
572  htmlFile << "<style type=\"text/css\">"<< std::endl;
573  htmlFile << " td.s0 { font-family: arial, arial ce, helvetica; font-weight: bold; background-color: #FF7700; text-align: center;}"<< std::endl;
574  htmlFile << " td.s1 { font-family: arial, arial ce, helvetica; font-weight: bold; background-color: #FFC169; text-align: center;}"<< std::endl;
575  htmlFile << " td.s2 { font-family: arial, arial ce, helvetica; background-color: red; }"<< std::endl;
576  htmlFile << " td.s3 { font-family: arial, arial ce, helvetica; background-color: yellow; }"<< std::endl;
577  htmlFile << " td.s4 { font-family: arial, arial ce, helvetica; background-color: green; }"<< std::endl;
578  htmlFile << " td.s5 { font-family: arial, arial ce, helvetica; background-color: silver; }"<< std::endl;
579  std::string state[4]={"<td class=\"s2\" align=\"center\">",
580  "<td class=\"s3\" align=\"center\">",
581  "<td class=\"s4\" align=\"center\">",
582  "<td class=\"s5\" align=\"center\">"};
583  htmlFile << "</style>"<< std::endl;
584  htmlFile << "<body> " << std::endl;
585  htmlFile << "<br> " << std::endl;
586  htmlFile << "<h2>Run:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;" << std::endl;
587  htmlFile << "&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <span " << std::endl;
588  htmlFile << " style=\"color: rgb(0, 0, 153);\">" << runNo << "</span></h2>" << std::endl;
589  htmlFile << "<h2>Monitoring task:&nbsp;&nbsp;&nbsp;&nbsp; <span " << std::endl;
590  htmlFile << " style=\"color: rgb(0, 0, 153);\">Detector Diagnostics LED Monitor</span></h2> " << std::endl;
591  htmlFile << "<h2>Events processed:&nbsp;&nbsp;&nbsp;&nbsp;<span " << std::endl;
592  htmlFile << " style=\"color: rgb(0, 0, 153);\">" << ievt_ << "</span></h2>" << std::endl;
593  htmlFile << "<hr>" << std::endl;
595  htmlFile << "<table width=100% border=1>" << std::endl;
596  htmlFile << "<tr>" << std::endl;
597  htmlFile << "<td class=\"s0\" width=15% align=\"center\">SebDet</td>" << std::endl;
598  htmlFile << "<td class=\"s0\" width=17% align=\"center\">Missing</td>" << std::endl;
599  htmlFile << "<td class=\"s0\" width=17% align=\"center\">Unstable</td>" << std::endl;
600  htmlFile << "<td class=\"s0\" width=17% align=\"center\">low/no LED signal</td>" << std::endl;
601  htmlFile << "<td class=\"s0\" width=17% align=\"center\">Bad Timing</td>" << std::endl;
602  htmlFile << "<td class=\"s0\" width=17% align=\"center\">Bad LED signal</td>" << std::endl;
603  htmlFile << "</tr><tr>" << std::endl;
604  int ind1=0,ind2=0,ind3=0,ind4=0,ind5=0;
605  htmlFile << "<td class=\"s1\" align=\"center\">HB+</td>" << std::endl;
606  ind1=3; if(newHBP[0]==0) ind1=2; if(newHBP[0]>0 && newHBP[0]<=12) ind1=1; if(newHBP[0]>=12 && HBP[0]<1296) ind1=0;
607  ind2=3; if(newHBP[1]==0) ind2=2; if(newHBP[1]>0) ind2=1; if(newHBP[1]>21) ind2=0;
608  ind3=3; if(newHBP[2]==0) ind3=2; if(newHBP[2]>0) ind3=1; if(newHBP[2]>21) ind3=0;
609  ind4=3; if((newHBP[3]+newHBP[4])==0) ind4=2; if((newHBP[3]+newHBP[4])>0) ind4=1; if((newHBP[3]+newHBP[4])>21) ind4=0;
610  ind5=3; if((newHBP[5]+newHBP[6])==0) ind5=2; if((newHBP[5]+newHBP[6])>0) ind5=1; if((newHBP[5]+newHBP[6])>21) ind5=0;
611  if(ind1==3) ind2=ind3=ind4=ind5=3;
612  if(ind1==0 || ind2==0 || ind3==0 || ind4==0) status|=2; else if(ind1==1 || ind2==1 || ind3==1 || ind4==1) status|=1;
613  htmlFile << state[ind1] << HBP[0] <<" (1296)</td>" << std::endl;
614  htmlFile << state[ind2] << HBP[1] <<"</td>" << std::endl;
615  htmlFile << state[ind3] << HBP[2] <<"</td>" << std::endl;
616  htmlFile << state[ind5] << HBP[5]+HBP[6] <<"</td>" << std::endl;
617  htmlFile << state[ind4] << HBP[3]+HBP[4] <<"</td>" << std::endl;
618 
619  htmlFile << "</tr><tr>" << std::endl;
620  htmlFile << "<td class=\"s1\" align=\"center\">HB-</td>" << std::endl;
621  ind1=3; if(newHBM[0]==0) ind1=2; if(newHBM[0]>0 && newHBM[0]<=12) ind1=1; if(newHBM[0]>=12 && HBM[0]<1296) ind1=0;
622  ind2=3; if(newHBM[1]==0) ind2=2; if(newHBM[1]>0) ind2=1; if(newHBM[1]>21) ind2=0;
623  ind3=3; if(newHBM[2]==0) ind3=2; if(newHBM[2]>0) ind3=1; if(newHBM[2]>21) ind3=0;
624  ind4=3; if((newHBM[3]+newHBM[4])==0) ind4=2; if((newHBM[3]+newHBM[4])>0) ind4=1; if((newHBM[3]+newHBM[4])>21) ind4=0;
625  ind5=3; if((newHBM[5]+newHBM[6])==0) ind5=2; if((newHBM[5]+newHBM[6])>0) ind5=1; if((newHBM[5]+newHBM[6])>21) ind5=0;
626  if(ind1==3) ind2=ind3=ind4=ind5=3;
627  if(ind1==0 || ind2==0 || ind3==0 || ind4==0) status|=2; else if(ind1==1 || ind2==1 || ind3==1 || ind4==1)status|=1;
628  htmlFile << state[ind1] << HBM[0] <<" (1296)</td>" << std::endl;
629  htmlFile << state[ind2] << HBM[1] <<"</td>" << std::endl;
630  htmlFile << state[ind3] << HBM[2] <<"</td>" << std::endl;
631  htmlFile << state[ind5] << HBM[5]+HBM[6] <<"</td>" << std::endl;
632  htmlFile << state[ind4] << HBM[3]+HBM[4] <<"</td>" << std::endl;
633 
634  htmlFile << "</tr><tr>" << std::endl;
635  htmlFile << "<td class=\"s1\" align=\"center\">HE+</td>" << std::endl;
636  ind1=3; if(newHEP[0]==0) ind1=2; if(newHEP[0]>0 && newHEP[0]<=12) ind1=1; if(newHEP[0]>=12 && HEP[0]<1296) ind1=0;
637  ind2=3; if(newHEP[1]==0) ind2=2; if(newHEP[1]>0) ind2=1; if(newHEP[1]>21) ind2=0;
638  ind3=3; if(newHEP[2]==0) ind3=2; if(newHEP[2]>0) ind3=1; if(newHEP[2]>21) ind3=0;
639  ind4=3; if((newHEP[3]+newHEP[4])==0) ind4=2; if((newHEP[3]+newHEP[4])>0) ind4=1; if((newHEP[3]+newHEP[4])>21) ind4=0;
640  ind5=3; if((newHEP[5]+newHEP[6])==0) ind5=2; if((newHEP[5]+newHEP[6])>0) ind5=1; if((newHEP[5]+newHEP[6])>21) ind5=0;
641  if(ind1==3) ind2=ind3=ind4=ind5=3;
642  if(ind1==0 || ind2==0 || ind3==0 || ind4==0) status|=2; else if(ind1==1 || ind2==1 || ind3==1 || ind4==1)status|=1;
643  htmlFile << state[ind1] << HEP[0] <<" (1296)</td>" << std::endl;
644  htmlFile << state[ind2] << HEP[1] <<"</td>" << std::endl;
645  htmlFile << state[ind3] << HEP[2] <<"</td>" << std::endl;
646  htmlFile << state[ind5] << HEP[5]+HEP[6] <<"</td>" << std::endl;
647  htmlFile << state[ind4] << HEP[3]+HEP[4] <<"</td>" << std::endl;
648 
649  htmlFile << "</tr><tr>" << std::endl;
650  htmlFile << "<td class=\"s1\" align=\"center\">HE-</td>" << std::endl;
651  ind1=3; if(newHEM[0]==0) ind1=2; if(newHEM[0]>0 && newHEM[0]<=12) ind1=1; if(newHEM[0]>=12 && HEM[0]<1296) ind1=0;
652  ind2=3; if(newHEM[1]==0) ind2=2; if(newHEM[1]>0) ind2=1; if(newHEM[1]>21) ind2=0;
653  ind3=3; if(newHEM[2]==0) ind3=2; if(newHEM[2]>0) ind3=1; if(newHEM[2]>21) ind3=0;
654  ind4=3; if((newHEM[3]+newHEM[4])==0) ind4=2; if((newHEM[3]+newHEM[4])>0) ind4=1; if((newHEM[3]+newHEM[4])>21) ind4=0;
655  ind5=3; if((newHEM[5]+newHEM[6])==0) ind5=2; if((newHEM[5]+newHEM[6])>0) ind5=1; if((newHEM[5]+newHEM[6])>21) ind5=0;
656  if(ind1==3) ind2=ind3=ind4=ind5=3;
657  if(ind1==0 || ind2==0 || ind3==0 || ind4==0) status|=2; else if(ind1==1 || ind2==1 || ind3==1 || ind4==1)status|=1;
658  htmlFile << state[ind1] << HEM[0] <<" (1296)</td>" << std::endl;
659  htmlFile << state[ind2] << HEM[1] <<"</td>" << std::endl;
660  htmlFile << state[ind3] << HEM[2] <<"</td>" << std::endl;
661  htmlFile << state[ind5] << HEM[5]+HEM[6] <<"</td>" << std::endl;
662  htmlFile << state[ind4] << HEM[3]+HEM[4] <<"</td>" << std::endl;
663 
664  htmlFile << "</tr><tr>" << std::endl;
665  htmlFile << "<td class=\"s1\" align=\"center\">HF+</td>" << std::endl;
666  ind1=3; if(newHFP[0]==0) ind1=2; if(newHFP[0]>0 && newHFP[0]<=12) ind1=1; if(newHFP[0]>=12 && HFP[0]<864) ind1=0;
667  ind2=3; if(newHFP[1]==0) ind2=2; if(newHFP[1]>0) ind2=1; if(newHFP[1]>21) ind2=0;
668  ind3=3; if(newHFP[2]==0) ind3=2; if(newHFP[2]>0) ind3=1; if(newHFP[2]>21) ind3=0;
669  ind4=3; if((newHFP[3]+newHFP[4])==0) ind4=2; if((newHFP[3]+newHFP[4])>0) ind4=1; if((newHFP[3]+newHFP[4])>21) ind4=0;
670  ind5=3; if((newHFP[5]+newHFP[6])==0) ind5=2; if((newHFP[5]+newHFP[6])>0) ind5=1; if((newHFP[5]+newHFP[6])>21) ind5=0;
671  if(ind1==3) ind2=ind3=ind4=ind5=3;
672  if(ind1==0 || ind2==0 || ind3==0 || ind4==0) status|=2; else if(ind1==1 || ind2==1 || ind3==1 || ind4==1)status|=1;
673  htmlFile << state[ind1] << HFP[0] <<" (864)</td>" << std::endl;
674  htmlFile << state[ind2] << HFP[1] <<"</td>" << std::endl;
675  htmlFile << state[ind3] << HFP[2] <<"</td>" << std::endl;
676  htmlFile << state[ind5] << HFP[5]+HFP[6] <<"</td>" << std::endl;
677  htmlFile << state[ind4] << HFP[3]+HFP[4] <<"</td>" << std::endl;
678 
679  htmlFile << "</tr><tr>" << std::endl;
680  htmlFile << "<td class=\"s1\" align=\"center\">HF-</td>" << std::endl;
681  ind1=3; if(newHFM[0]==0) ind1=2; if(newHFM[0]>0 && newHFM[0]<=12) ind1=1; if(newHFM[0]>=12 && HFM[0]<864) ind1=0;
682  ind2=3; if(newHFM[1]==0) ind2=2; if(newHFM[1]>0) ind2=1; if(newHFM[1]>21) ind2=0;
683  ind3=3; if(newHFM[2]==0) ind3=2; if(newHFM[2]>0) ind3=1; if(newHFM[2]>21) ind3=0;
684  ind4=3; if((HFM[3]+HFM[4])==0) ind4=2; if((HFM[3]+HFM[4])>0) ind4=1; if((HFM[3]+HFM[4])>21) ind4=0;
685  ind5=3; if((HFM[5]+HFM[6])==0) ind5=2; if((HFM[5]+HFM[6])>0) ind5=1; if((HFM[5]+HFM[6])>21) ind5=0;
686  if(ind1==3) ind2=ind3=ind4=ind5=3;
687  if(ind1==0 || ind2==0 || ind3==0 || ind4==0) status|=2; else if(ind1==1 || ind2==1 || ind3==1 || ind4==1)status|=1;
688  htmlFile << state[ind1] << HFM[0] <<" (864)</td>" << std::endl;
689  htmlFile << state[ind2] << HFM[1] <<"</td>" << std::endl;
690  htmlFile << state[ind3] << HFM[2] <<"</td>" << std::endl;
691  htmlFile << state[ind5] << HFM[5]+HFM[6] <<"</td>" << std::endl;
692  htmlFile << state[ind4] << HFM[3]+HFM[4] <<"</td>" << std::endl;
693 
694  htmlFile << "</tr><tr>" << std::endl;
695  htmlFile << "<td class=\"s1\" align=\"center\">HO</td>" << std::endl;
696  ind1=3; if(newHO[0]==0) ind1=2; if(newHO[0]>0 && newHO[0]<=12) ind1=1; if(newHO[0]>=12 && HO[0]<2160) ind1=0;
697  ind2=3; if(newHO[1]==0) ind2=2; if(newHO[1]>0) ind2=1; if(newHO[1]>21) ind2=0;
698  ind3=3; if(newHO[2]==0) ind3=2; if(newHO[2]>0) ind3=1; if(newHO[2]>21) ind3=0;
699  ind4=3; if((newHO[3]+newHO[4])==0) ind4=2; if((newHO[3]+newHO[4])>0) ind4=1; if((newHO[3]+newHO[4])>21) ind4=0;
700  ind5=3; if((newHO[5]+newHO[6])==0) ind5=2; if((newHO[5]+newHO[6])>0) ind5=1; if((newHO[5]+newHO[6])>21) ind5=0;
701  if(ind1==3) ind2=ind3=ind4=ind5=3;
702  if(ind1==0 || ind2==0 || ind3==0 || ind4==0) status|=2; else if(ind1==1 || ind2==1 || ind3==1 || ind4==1)status|=1;
703 
704  htmlFile << state[ind1] << HO[0] <<" (2160)</td>" << std::endl;
705  htmlFile << state[ind2] << HO[1] <<"</td>" << std::endl;
706  htmlFile << state[ind3] << HO[2] <<"</td>" << std::endl;
707  htmlFile << state[ind5] << HO[5]+HO[6] <<"</td>" << std::endl;
708  htmlFile << state[ind4] << HO[3]+HO[4] <<"</td>" << std::endl;
709 
710  htmlFile << "</tr></table>" << std::endl;
711  htmlFile << "<hr>" << std::endl;
713  if((MissingCnt+UnstableCnt+UnstableLEDCnt+BadTimingCnt)>0){
714  htmlFile << "<table width=100% border=1><tr>" << std::endl;
715  if(MissingCnt>0) htmlFile << "<td><a href=\"" << "Missing.html" <<"\">list of missing channels</a></td>";
716  if(UnstableCnt>0) htmlFile << "<td><a href=\"" << "Unstable.html" <<"\">list of unstable channels</a></td>";
717  if(UnstableLEDCnt>0)htmlFile << "<td><a href=\"" << "UnstableLED.html" <<"\">list of low LED signal channels</a></td>";
718  if(BadTimingCnt>0)htmlFile << "<td><a href=\"" << "BadTiming.html" <<"\">list of Bad Timing channels</a></td>";
719  htmlFile << "</tr></table>" << std::endl;
720  }
722 
724  htmlFile << "<h2 align=\"center\">Summary LED plots</h2>" << std::endl;
725  htmlFile << "<table width=100% border=0><tr>" << std::endl;
726  htmlFile << "<tr align=\"left\">" << std::endl;
727  Time2Dhbhehf->SetMaximum(6);
728  Time2Dho->SetMaximum(6);
729  Time2Dhbhehf->SetNdivisions(36,"Y");
730  Time2Dho->SetNdivisions(36,"Y");
731  Time2Dhbhehf->SetStats(0);
732  Time2Dho->SetStats(0);
733  Time2Dhbhehf->Draw("COLZ"); can->SaveAs((htmlDir + "led_timing_hbhehf.gif").c_str());
734  Time2Dho->Draw("COLZ"); can->SaveAs((htmlDir + "led_timing_ho.gif").c_str());
735  htmlFile << "<td align=\"center\"><img src=\"led_timing_hbhehf.gif\" alt=\"led timing distribution\"> </td>" << std::endl;
736  htmlFile << "<td align=\"center\"><img src=\"led_timing_ho.gif\" alt=\"led timing distribution\"> </td>" << std::endl;
737  htmlFile << "</tr>" << std::endl;
738 
739  htmlFile << "<tr align=\"left\">" << std::endl;
740  Energy2Dhbhehf->SetStats(0);
741  Energy2Dho->SetStats(0);
742  Energy2Dhbhehf->SetNdivisions(36,"Y");
743  Energy2Dho->SetNdivisions(36,"Y");
744  Energy2Dhbhehf->Draw("COLZ"); can->SaveAs((htmlDir + "led_energy_hbhehf.gif").c_str());
745  Energy2Dho->Draw("COLZ"); can->SaveAs((htmlDir + "led_energy_ho.gif").c_str());
746  htmlFile << "<td align=\"center\"><img src=\"led_energy_hbhehf.gif\" alt=\"led energy distribution\"> </td>" << std::endl;
747  htmlFile << "<td align=\"center\"><img src=\"led_energy_ho.gif\" alt=\"led energy distribution\"> </td>" << std::endl;
748  htmlFile << "</tr>" << std::endl;
749 
750  can->SetGridy(false);
751  htmlFile << "<tr align=\"left\">" << std::endl;
752  Energy->Draw(); can->SaveAs((htmlDir + "led_energy_distribution.gif").c_str());
753  EnergyRMS->Draw(); can->SaveAs((htmlDir + "led_energy_rms_distribution.gif").c_str());
754  htmlFile << "<td align=\"center\"><img src=\"led_energy_distribution.gif\" alt=\"led energy distribution\"> </td>" << std::endl;
755  htmlFile << "<td align=\"center\"><img src=\"led_energy_rms_distribution.gif\" alt=\"led energy rms distribution\"> </td>" << std::endl;
756  htmlFile << "</tr>" << std::endl;
757  htmlFile << "<tr align=\"left\">" << std::endl;
758  Timing->Draw(); can->SaveAs((htmlDir + "led_timing_distribution.gif").c_str());
759  TimingRMS->Draw(); can->SaveAs((htmlDir + "led_timing_rms_distribution.gif").c_str());
760  htmlFile << "<td align=\"center\"><img src=\"led_timing_distribution.gif\" alt=\"led timing distribution\"> </td>" << std::endl;
761  htmlFile << "<td align=\"center\"><img src=\"led_timing_rms_distribution.gif\" alt=\"led timing rms distribution\"> </td>" << std::endl;
762  htmlFile << "</tr>" << std::endl;
763  htmlFile << "<tr align=\"left\">" << std::endl;
764  EnergyHF->Draw(); can->SaveAs((htmlDir + "led_energyhf_distribution.gif").c_str());
765  EnergyRMSHF->Draw(); can->SaveAs((htmlDir + "led_energyhf_rms_distribution.gif").c_str());
766  htmlFile << "<td align=\"center\"><img src=\"led_energyhf_distribution.gif\" alt=\"hf led energy distribution\"> </td>" << std::endl;
767  htmlFile << "<td align=\"center\"><img src=\"led_energyhf_rms_distribution.gif\" alt=\"hf led energy rms distribution\"> </td>" << std::endl;
768  htmlFile << "</tr>" << std::endl;
769  htmlFile << "<tr align=\"left\">" << std::endl;
770  TimingHF->Draw(); can->SaveAs((htmlDir + "led_timinghf_distribution.gif").c_str());
771  TimingRMSHF->Draw(); can->SaveAs((htmlDir + "led_timinghf_rms_distribution.gif").c_str());
772  htmlFile << "<td align=\"center\"><img src=\"led_timinghf_distribution.gif\" alt=\"hf led timing distribution\"> </td>" << std::endl;
773  htmlFile << "<td align=\"center\"><img src=\"led_timinghf_rms_distribution.gif\" alt=\"hf led timing rms distribution\"> </td>" << std::endl;
774  htmlFile << "</tr>" << std::endl;
775  htmlFile << "</table>" << std::endl;
776 
778  htmlFile << "<h2 align=\"center\">Stability LED plots (Reference run "<<ref_run<<")</h2>" << std::endl;
779  htmlFile << "<table width=100% border=0><tr>" << std::endl;
780  htmlFile << "<tr align=\"left\">" << std::endl;
781  HBPphi->GetXaxis()->SetNdivisions(418,kFALSE);
782  HBMphi->GetXaxis()->SetNdivisions(418,kFALSE);
783  HEPphi->GetXaxis()->SetNdivisions(418,kFALSE);
784  HEMphi->GetXaxis()->SetNdivisions(418,kFALSE);
785 
786  HBPphi->SetMarkerColor(kRed);
787  HBPphi->SetMarkerStyle(23);
788  HBPphi->SetXTitle("HPD Index = RBX*4+RM");
789  HBMphi->SetMarkerColor(kRed);
790  HBMphi->SetMarkerStyle(23);
791  HBMphi->SetXTitle("HPD Index = RBX*4+RM");
792  HBPphi->Draw(); can->SaveAs((htmlDir + "led_hbp_distribution.gif").c_str());
793  HBMphi->Draw(); can->SaveAs((htmlDir + "led_hbm_distribution.gif").c_str());
794  htmlFile << "<td align=\"center\"><img src=\"led_hbp_distribution.gif\" alt=\"led hbp distribution\"> </td>" << std::endl;
795  htmlFile << "<td align=\"center\"><img src=\"led_hbm_distribution.gif\" alt=\"led hbm distribution\"> </td>" << std::endl;
796  htmlFile << "</tr>" << std::endl;
797 
798  htmlFile << "<tr align=\"left\">" << std::endl;
799  HEPphi->SetMarkerColor(kRed);
800  HEPphi->SetMarkerStyle(23);
801  HEPphi->SetXTitle("HPD Index = RBX*4+RM");
802  HEMphi->SetMarkerColor(kRed);
803  HEMphi->SetMarkerStyle(23);
804  HEMphi->SetXTitle("HPD Index = RBX*4+RM");
805  HEPphi->Draw(); can->SaveAs((htmlDir + "led_hep_distribution.gif").c_str());
806  HEMphi->Draw(); can->SaveAs((htmlDir + "led_hem_distribution.gif").c_str());
807  htmlFile << "<td align=\"center\"><img src=\"led_hep_distribution.gif\" alt=\"led hep distribution\"> </td>" << std::endl;
808  htmlFile << "<td align=\"center\"><img src=\"led_hem_distribution.gif\" alt=\"led hem distribution\"> </td>" << std::endl;
809  htmlFile << "</tr>" << std::endl;
810 
811  htmlFile << "<tr align=\"left\">" << std::endl;
812  HFPphi->SetMarkerColor(kRed);
813  HFPphi->SetMarkerStyle(23);
814  HFPphi->SetXTitle("RM Index = RoBox*3+RM");
815  HFMphi->SetMarkerColor(kRed);
816  HFMphi->SetMarkerStyle(23);
817  HFPphi->GetXaxis()->SetNdivisions(312,kFALSE);
818  HFMphi->GetXaxis()->SetNdivisions(312,kFALSE);
819 
820  HFMphi->SetXTitle("RM Index = RoBox*3+RM");
821  HFPphi->Draw(); can->SaveAs((htmlDir + "led_hfp_distribution.gif").c_str());
822  HFMphi->Draw(); can->SaveAs((htmlDir + "led_hfm_distribution.gif").c_str());
823  htmlFile << "<td align=\"center\"><img src=\"led_hfp_distribution.gif\" alt=\"led hfp distribution\"> </td>" << std::endl;
824  htmlFile << "<td align=\"center\"><img src=\"led_hfm_distribution.gif\" alt=\"led hfm distribution\"> </td>" << std::endl;
825  htmlFile << "</tr>" << std::endl;
826 
827  htmlFile << "<tr align=\"left\">" << std::endl;
828  HO1Pphi->SetMarkerColor(kRed);
829  HO1Pphi->SetMarkerStyle(23);
830  HO1Pphi->SetXTitle("HPD Index = RBX*4+RM");
831  HO1Mphi->SetMarkerColor(kRed);
832  HO1Mphi->SetMarkerStyle(23);
833  HO1Mphi->GetXaxis()->SetNdivisions(412,kFALSE);
834  HO1Pphi->GetXaxis()->SetNdivisions(412,kFALSE);
835 
836  HO1Mphi->SetXTitle("HPD Index = RBX*4+RM");
837  HO1Pphi->Draw(); can->SaveAs((htmlDir + "led_ho1p_distribution.gif").c_str());
838  HO1Mphi->Draw(); can->SaveAs((htmlDir + "led_ho1m_distribution.gif").c_str());
839  htmlFile << "<td align=\"center\"><img src=\"led_ho1p_distribution.gif\" alt=\"led ho1p distribution\"> </td>" << std::endl;
840  htmlFile << "<td align=\"center\"><img src=\"led_ho1m_distribution.gif\" alt=\"led ho1m distribution\"> </td>" << std::endl;
841  htmlFile << "</tr>" << std::endl;
842 
843  htmlFile << "<tr align=\"left\">" << std::endl;
844  HO2Pphi->SetMarkerColor(kRed);
845  HO2Pphi->SetMarkerStyle(23);
846  HO2Pphi->SetXTitle("HPD Index = RBX*4+RM");
847  HO2Mphi->SetMarkerColor(kRed);
848  HO2Mphi->SetMarkerStyle(23);
849  HO2Mphi->GetXaxis()->SetNdivisions(412,kFALSE);
850  HO2Pphi->GetXaxis()->SetNdivisions(412,kFALSE);
851 
852  HO2Mphi->SetXTitle("HPD Index = RBX*4+RM");
853  HO2Pphi->Draw(); can->SaveAs((htmlDir + "led_ho2p_distribution.gif").c_str());
854  HO2Mphi->Draw(); can->SaveAs((htmlDir + "led_ho2m_distribution.gif").c_str());
855  htmlFile << "<td align=\"center\"><img src=\"led_ho2p_distribution.gif\" alt=\"led ho2p distribution\"> </td>" << std::endl;
856  htmlFile << "<td align=\"center\"><img src=\"led_ho2m_distribution.gif\" alt=\"led ho2m distribution\"> </td>" << std::endl;
857  htmlFile << "</tr>" << std::endl;
858 
859  htmlFile << "<tr align=\"left\">" << std::endl;
860  HO0phi->SetMarkerColor(kRed);
861  HO0phi->SetMarkerStyle(23);
862  HO0phi->SetXTitle("HPD Index = RBX*4+RM");
863  HO0phi->GetXaxis()->SetNdivisions(412,kFALSE);
864  HO0phi->Draw(); can->SaveAs((htmlDir + "led_ho0_distribution.gif").c_str());
865  htmlFile << "<td align=\"center\"><img src=\"led_ho0_distribution.gif\" alt=\"led ho0 distribution\"> </td>" << std::endl;
866  htmlFile << "</tr>" << std::endl;
867 
868  htmlFile << "</table>" << std::endl;
869 
870  htmlFile << "</body> " << std::endl;
871  htmlFile << "</html> " << std::endl;
872  htmlFile.close();
873  can->Close();
874 }
875 
877 {}
int fiberChannel() const
type
Definition: HCALResponse.h:21
dictionary missing
Definition: combine.py:4
T getUntrackedParameter(std::string const &, T const &) const
const std::string & getName(void) const
get name of ME
int i
Definition: DBlmapReader.cc:9
MonitorElement * ProblemCells
int fiberIndex() const
get the fiber index [1-8] (which of eight fibers carried by a spigot) (valid only for non-trigger-cha...
std::string name()
int pixel() const
TH2F * ChannelStatusUnstableChannels[4]
static void printTableTail(ofstream &file)
HcalDetDiagLEDClient()
Constructors.
std::vector< std::string > HcalEtaPhiHistNames()
int htrSlot() const
get the htr slot
void setup(DQMStore *&m_dbe, std::string Name, std::string Units="")
T eta() const
int readoutVMECrateId() const
get the readout VME crate number
std::string rbx() const
void htmlOutput(std::string)
int depth() const
get the tower depth
Definition: HcalDetId.h:42
list outfile
Definition: EdgesToViz.py:91
HcalElectronicsMap generateHcalElectronicsMap()
std::vector< MonitorElement * > depth
int htrTopBottom() const
get the htr top/bottom (1=top/0=bottom)
std::vector< std::string > problemnames_
int rm() const
EtaPhiHists * ProblemCellsByDepth
int ieta() const
get the cell ieta
Definition: HcalDetId.h:38
int dccid() const
get the (Hcal local) DCC id
void updateChannelStatus(std::map< HcalDetId, unsigned int > &myqual)
HcalSubdetector
Definition: HcalAssistant.h:32
double get_energy(std::string subdet, int eta, int phi, int depth, int type)
TH2F * ChannelStatusMissingChannels[4]
MonitorElement * get(const std::string &path) const
get ME from full pathname (e.g. &quot;my/long/dir/my_histo&quot;)
Definition: DQMStore.cc:1473
double get_channel_status(std::string subdet, int eta, int phi, int depth, int type)
int spigot() const
get the spigot (input number on DCC)
std::string valueString(void) const
const HcalFrontEndId getHcalFrontEndId(const DetId &)
int iphi() const
get the cell iphi
Definition: HcalDetId.h:40
double sd
int adc() const
static void printTableLine(ofstream &file, int ind, HcalDetId &detid, HcalFrontEndId &lmap_entry, HcalElectronicsId &emap_entry, std::string comment="")
std::string subdir_
std::map< HcalDetId, unsigned int > badstatusmap
int CalcEtaBin(int subdet, int ieta, int depth)
char state
Definition: procUtils.cc:75
std::string prefixME_
HcalLogicalMap * logicalMap_
static float SD[]
Definition: sicif.h:67
tuple cout
Definition: gather_cfg.py:121
static void printTableHeader(ofstream &file, std::string header)
void showDirStructure(void) const
Definition: DQMStore.cc:2766
int rmFiber() const
~HcalDetDiagLEDClient()
Destructor.
TH2F * ChannelStatusUnstableLEDsignal[4]
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:850
int qieCard() const
Readout chain identification for Hcal [31:26] Unused (so far) [25] Trigger-chain id flag [24:20] Read...
DQMStore * dqmStore_
const DetId lookup(HcalElectronicsId fId) const
lookup the logical detid associated with the given electronics id
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:434
Definition: DDAxes.h:10
#define comment(par)
Definition: vmac.h:162