221 int UnstableLEDCnt=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};
232 if (
debug_>0)
std::cout <<
"<HcalDetDiagLEDClient::htmlOutput> Preparing html output ..." << std::endl;
236 for(
int i=0;
i<4;++
i){
241 s=
subdir_+
"channel status/"+name[
i]+
" Missing Channels";
242 me=ig.
get(s.c_str());
244 s=
subdir_+
"channel status/"+name[
i]+
" Unstable Channels";
245 me=ig.
get(s.c_str());
247 s=
subdir_+
"channel status/"+name[
i]+
" Unstable LED";
248 me=ig.
get(s.c_str());
250 s=
subdir_+
"channel status/"+name[
i]+
" LED Mean";
251 me=ig.
get(s.c_str());
253 s=
subdir_+
"channel status/"+name[
i]+
" LED RMS";
254 me=ig.
get(s.c_str());
256 s=
subdir_+
"channel status/"+name[
i]+
" Time Mean";
257 me=ig.
get(s.c_str());
259 s=
subdir_+
"channel status/"+name[
i]+
" Time RMS";
260 me=ig.
get(s.c_str());
262 s=
subdir_+
"Summary Plots/"+name[
i]+
" Channel LED Energy";
263 me=ig.
get(s.c_str());
265 s=
subdir_+
"Summary Plots/"+name[
i]+
" Channel LED Energy Reference";
266 me=ig.
get(s.c_str());
272 int feta=0,teta=0,
fdepth=0,tdepth=0;
273 if(
sd==0){ feta=-16; teta=16 ;
fdepth=1; tdepth=2; }
274 if(
sd==1){ feta=-29; teta=29 ;
fdepth=1; tdepth=3; }
275 if(
sd==2){ feta=-15; teta=15 ;
fdepth=4; tdepth=4; }
276 if(
sd==3){ feta=-42; teta=42 ;
fdepth=1; tdepth=2; }
278 if(
sd==3 &&
eta>-29 &&
eta<29)
continue;
279 double problem[7]={0,0,0,0,0,0,0};
280 for(
int i=0;
i<6;
i++){
289 if(
sd==0){
if(
eta>0){
296 if(
sd==1){
if(
eta>0){
307 if(
sd==3){
if(
eta>0){
319 std::ofstream Missing;
320 Missing.open((htmlDir +
"Missing.html").c_str());
323 std::ofstream BadTiming;
324 BadTiming.open((htmlDir +
"BadTiming.html").c_str());
327 std::ofstream Unstable;
328 Unstable.open((htmlDir +
"Unstable.html").c_str());
331 std::ofstream BadLED;
332 BadLED.open((htmlDir +
"UnstableLED.html").c_str());
337 if(
sd==0 && ((HBM[0]+HBP[0])==0 || (HBM[0]+HBP[0])==(1296*2)))
continue;
338 if(
sd==1 && ((HEM[0]+HEP[0])==0 || (HEM[0]+HEP[0])==(1296*2)))
continue;
339 if(
sd==2 && ((HO[0])==0 || HO[0]==2160))
continue;
340 if(
sd==3 && ((HFM[0]+HFP[0])==0 || (HFM[0]+HFP[0])==(864*2)))
continue;
341 Missing <<
"<tr><td align=\"center\"><h3>"<< subdet[
sd] <<
"</h3></td></tr>" << std::endl;
342 int feta=0,teta=0,
fdepth=0,tdepth=0;
343 if(
sd==0){ feta=-16; teta=16 ;
fdepth=1; tdepth=2;
if(HBM[0]==1296) feta=0;
if(HBP[0]==1296) teta=0;}
344 if(
sd==1){ feta=-29; teta=29 ;
fdepth=1; tdepth=3;
if(HEM[0]==1296) feta=0;
if(HEP[0]==1296) teta=0;}
345 if(
sd==2){ feta=-15; teta=15 ;
fdepth=4; tdepth=4;
if(HO[0] ==2160) {feta=0; teta=0; }}
346 if(
sd==3){ feta=-42; teta=42 ;
fdepth=1; tdepth=2;
if(HFM[0]==864) feta=0;
if(HFP[0]==864) teta=0; }
348 if(
sd==3 &&
eta>-29 &&
eta<29)
continue;
360 printTableLine(Missing,cnt++,*detid,lmap_entry,emap_entry,s); MissingCnt++;
368 if(
sd==0 && (HBM[5]+HBP[5])==0)
continue;
369 if(
sd==1 && (HEM[5]+HEP[5])==0)
continue;
370 if(
sd==2 && (HO[5])==0)
continue;
371 if(
sd==3 && (HFM[5]+HFP[5])==0)
continue;
372 BadTiming <<
"<tr><td align=\"center\"><h3>"<< subdet[
sd] <<
"</h3></td></tr>" << std::endl;
373 int feta=0,teta=0,
fdepth=0,tdepth=0;
374 if(
sd==0){ feta=-16; teta=16 ;
fdepth=1; tdepth=2;
if(HBM[0]==1296) feta=0;
if(HBP[0]==1296) teta=0;}
375 if(
sd==1){ feta=-29; teta=29 ;
fdepth=1; tdepth=3;
if(HEM[0]==1296) feta=0;
if(HEP[0]==1296) teta=0;}
376 if(
sd==2){ feta=-15; teta=15 ;
fdepth=4; tdepth=4;
if(HO[0] ==2160) {feta=0; teta=0; }}
377 if(
sd==3){ feta=-42; teta=42 ;
fdepth=1; tdepth=2;
if(HFM[0]==864) feta=0;
if(HFP[0]==864) teta=0; }
379 if(
sd==3 &&
eta>-29 &&
eta<29)
continue;
383 char comment[100]; sprintf(comment,
"Time-mean=%.1f\n",badtiming);
391 printTableLine(BadTiming,cnt++,*detid,lmap_entry,emap_entry,comment); BadTimingCnt++;
400 if(
sd==0 && (HBM[1]+HBP[1])==0)
continue;
401 if(
sd==1 && (HEM[1]+HEP[1])==0)
continue;
402 if(
sd==2 && (HO[1])==0)
continue;
403 if(
sd==3 && (HFM[1]+HFP[1])==0)
continue;
404 Unstable <<
"<tr><td align=\"center\"><h3>"<< subdet[
sd] <<
"</h3></td></tr>" << std::endl;
405 int feta=0,teta=0,
fdepth=0,tdepth=0;
406 if(
sd==0){ feta=-16; teta=16 ;
fdepth=1; tdepth=2;}
407 if(
sd==1){ feta=-29; teta=29 ;
fdepth=1; tdepth=3;}
408 if(
sd==2){ feta=-15; teta=15 ;
fdepth=4; tdepth=4;}
409 if(
sd==3){ feta=-42; teta=42 ;
fdepth=1; tdepth=2;}
411 if(
sd==3 &&
eta>-29 &&
eta<29)
continue;
415 char comment[100]; sprintf(comment,
"%.3f%%\n",(1.0-unstable)*100.0);
423 printTableLine(Unstable,cnt++,*detid,lmap_entry,emap_entry,comment); UnstableCnt++;
432 if(
sd==0 && (HBM[2]+HBP[2])==0)
continue;
433 if(
sd==1 && (HEM[2]+HEP[2])==0)
continue;
434 if(
sd==2 && (HO[2])==0)
continue;
435 if(
sd==3 && (HFM[2]+HFP[2])==0)
continue;
436 BadLED <<
"<tr><td align=\"center\"><h3>"<< subdet[
sd] <<
"</h3></td></tr>" << std::endl;
437 int feta=0,teta=0,
fdepth=0,tdepth=0;
438 if(
sd==0){ feta=-16; teta=16 ;
fdepth=1; tdepth=2;}
439 if(
sd==1){ feta=-29; teta=29 ;
fdepth=1; tdepth=3;}
440 if(
sd==2){ feta=-15; teta=15 ;
fdepth=4; tdepth=4;}
441 if(
sd==3){ feta=-42; teta=42 ;
fdepth=1; tdepth=2;}
443 if(
sd==3 &&
eta>-29 &&
eta<29)
continue;
447 char comment[100]; sprintf(comment,
"%.3f%%\n",(badled)*100.0);
457 printTableLine(BadLED,cnt++,*detid,lmap_entry,emap_entry,s); UnstableLEDCnt++;
472 int ievt_ = -1,runNo=-1;
475 me = ig.
get(s.c_str());
478 sscanf((s.substr(2,s.length()-2)).c_str(),
"%d", &ievt_);
480 s=
subdir_+
"HcalDetDiagLEDMonitor Run Number";
481 me = ig.
get(s.c_str());
484 sscanf((s.substr(2,s.length()-2)).c_str(),
"%d", &runNo);
486 s=
subdir_+
"HcalDetDiagLEDMonitor Reference Run";
487 me = ig.
get(s.c_str());
491 sscanf((s.substr(2,s.length()-2)).c_str(),
"%s", str);
494 TH1F *Energy=0,*
Timing=0,*EnergyHF=0,*TimingHF=0,*EnergyRMS=0,*TimingRMS=0,*EnergyRMSHF=0,*TimingRMSHF=0;
495 TH2F *Time2Dhbhehf=0,*Time2Dho=0,*Energy2Dhbhehf=0,*Energy2Dho=0;
496 TH2F *HBPphi=0,*HBMphi=0,*HEPphi=0,*HEMphi=0,*HFPphi=0,*HFMphi=0,*HO0phi=0,*HO1Pphi=0,*HO2Pphi=0,*HO1Mphi=0,*HO2Mphi=0;
498 s=
subdir_+
"Summary Plots/HBHEHO LED Energy Distribution"; me=ig.
get(s.c_str());
499 if(me!=0) Energy=HcalUtilsClient::getHisto<TH1F*>(
me,
cloneME_, Energy,
debug_);
else return;
500 s=
subdir_+
"Summary Plots/HBHEHO LED Timing Distribution"; me=ig.
get(s.c_str());
502 s=
subdir_+
"Summary Plots/HBHEHO LED Energy RMS_div_Energy Distribution"; me=ig.
get(s.c_str());
503 if(me!=0) EnergyRMS=HcalUtilsClient::getHisto<TH1F*>(
me,
cloneME_, EnergyRMS,
debug_);
else return;
504 s=
subdir_+
"Summary Plots/HBHEHO LED Timing RMS Distribution"; me=ig.
get(s.c_str());
505 if(me!=0) TimingRMS=HcalUtilsClient::getHisto<TH1F*>(
me,
cloneME_, TimingRMS,
debug_);
else return;
506 s=
subdir_+
"Summary Plots/HF LED Energy Distribution"; me=ig.
get(s.c_str());
507 if(me!=0) EnergyHF=HcalUtilsClient::getHisto<TH1F*>(
me,
cloneME_, EnergyHF,
debug_);
else return;
508 s=
subdir_+
"Summary Plots/HF LED Timing Distribution"; me=ig.
get(s.c_str());
509 if(me!=0) TimingHF=HcalUtilsClient::getHisto<TH1F*>(
me,
cloneME_, TimingHF,
debug_);
else return;
510 s=
subdir_+
"Summary Plots/HF LED Energy RMS_div_Energy Distribution"; me=ig.
get(s.c_str());
511 if(me!=0) EnergyRMSHF=HcalUtilsClient::getHisto<TH1F*>(
me,
cloneME_, EnergyRMSHF,
debug_);
else return;
512 s=
subdir_+
"Summary Plots/HF LED Timing RMS Distribution"; me=ig.
get(s.c_str());
513 if(me!=0) TimingRMSHF=HcalUtilsClient::getHisto<TH1F*>(
me,
cloneME_, TimingRMSHF,
debug_);
else return;
515 s=
subdir_+
"Summary Plots/LED Timing HBHEHF"; me=ig.
get(s.c_str());
516 if(me!=0) Time2Dhbhehf=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, Time2Dhbhehf,
debug_);
else return;
517 s=
subdir_+
"Summary Plots/LED Timing HO"; me=ig.
get(s.c_str());
518 if(me!=0) Time2Dho=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, Time2Dho,
debug_);
else return;
519 s=
subdir_+
"Summary Plots/LED Energy HBHEHF"; me=ig.
get(s.c_str());
520 if(me!=0) Energy2Dhbhehf=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, Energy2Dhbhehf,
debug_);
else return;
521 s=
subdir_+
"Summary Plots/LED Energy HO"; me=ig.
get(s.c_str());
522 if(me!=0) Energy2Dho=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, Energy2Dho,
debug_);
else return;
524 s=
subdir_+
"Summary Plots/HBP Average over HPD LED Ref"; me=ig.
get(s.c_str());
525 if(me!=0) HBPphi=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, HBPphi,
debug_);
else return;
526 s=
subdir_+
"Summary Plots/HBM Average over HPD LED Ref"; me=ig.
get(s.c_str());
527 if(me!=0) HBMphi=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, HBMphi,
debug_);
else return;
528 s=
subdir_+
"Summary Plots/HEP Average over HPD LED Ref"; me=ig.
get(s.c_str());
529 if(me!=0) HEPphi=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, HEPphi,
debug_);
else return;
530 s=
subdir_+
"Summary Plots/HEM Average over HPD LED Ref"; me=ig.
get(s.c_str());
531 if(me!=0) HEMphi=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, HEMphi,
debug_);
else return;
532 s=
subdir_+
"Summary Plots/HFP Average over RM LED Ref"; me=ig.
get(s.c_str());
533 if(me!=0) HFPphi=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, HFPphi,
debug_);
else return;
534 s=
subdir_+
"Summary Plots/HFM Average over RM LED Ref"; me=ig.
get(s.c_str());
535 if(me!=0) HFMphi=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, HFMphi,
debug_);
else return;
537 s=
subdir_+
"Summary Plots/HO0 Average over HPD LED Ref"; me=ig.
get(s.c_str());
538 if(me!=0) HO0phi=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, HO0phi,
debug_);
else return;
539 s=
subdir_+
"Summary Plots/HO1P Average over HPD LED Ref"; me=ig.
get(s.c_str());
540 if(me!=0) HO1Pphi=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, HO1Pphi,
debug_);
else return;
541 s=
subdir_+
"Summary Plots/HO2P Average over HPD LED Ref"; me=ig.
get(s.c_str());
542 if(me!=0) HO2Pphi=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, HO2Pphi,
debug_);
else return;
543 s=
subdir_+
"Summary Plots/HO1M Average over HPD LED Ref"; me=ig.
get(s.c_str());
544 if(me!=0) HO1Mphi=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, HO1Mphi,
debug_);
else return;
545 s=
subdir_+
"Summary Plots/HO2M Average over HPD LED Ref"; me=ig.
get(s.c_str());
546 if(me!=0) HO2Mphi=HcalUtilsClient::getHisto<TH2F*>(
me,
cloneME_, HO2Mphi,
debug_);
else return;
548 gROOT->SetBatch(
true);
549 gStyle->SetCanvasColor(0);
550 gStyle->SetPadColor(0);
551 gStyle->SetOptStat(111110);
552 gStyle->SetPalette(1);
554 TCanvas *
can=
new TCanvas(
"HcalDetDiagLEDClient",
"HcalDetDiagLEDClient",0,0,500,350);
559 std::ofstream htmlFile;
561 htmlFile.open(outfile.c_str());
563 htmlFile <<
"<!DOCTYPE html PUBLIC \"-//W3C//DTD HTML 4.01 Transitional//EN\"> " << std::endl;
564 htmlFile <<
"<html> " << std::endl;
565 htmlFile <<
"<head> " << std::endl;
566 htmlFile <<
" <meta content=\"text/html; charset=ISO-8859-1\" " << std::endl;
567 htmlFile <<
" http-equiv=\"content-type\"> " << std::endl;
568 htmlFile <<
" <title>Detector Diagnostics LED Monitor</title> " << std::endl;
569 htmlFile <<
"</head> " << std::endl;
570 htmlFile <<
"<style type=\"text/css\"> td { font-weight: bold } </style>" << std::endl;
571 htmlFile <<
"<style type=\"text/css\">"<< std::endl;
572 htmlFile <<
" td.s0 { font-family: arial, arial ce, helvetica; font-weight: bold; background-color: #FF7700; text-align: center;}"<< std::endl;
573 htmlFile <<
" td.s1 { font-family: arial, arial ce, helvetica; font-weight: bold; background-color: #FFC169; text-align: center;}"<< std::endl;
574 htmlFile <<
" td.s2 { font-family: arial, arial ce, helvetica; background-color: red; }"<< std::endl;
575 htmlFile <<
" td.s3 { font-family: arial, arial ce, helvetica; background-color: yellow; }"<< std::endl;
576 htmlFile <<
" td.s4 { font-family: arial, arial ce, helvetica; background-color: green; }"<< std::endl;
577 htmlFile <<
" td.s5 { font-family: arial, arial ce, helvetica; background-color: silver; }"<< std::endl;
578 std::string state[4]={
"<td class=\"s2\" align=\"center\">",
579 "<td class=\"s3\" align=\"center\">",
580 "<td class=\"s4\" align=\"center\">",
581 "<td class=\"s5\" align=\"center\">"};
582 htmlFile <<
"</style>"<< std::endl;
583 htmlFile <<
"<body> " << std::endl;
584 htmlFile <<
"<br> " << std::endl;
585 htmlFile <<
"<h2>Run: " << std::endl;
586 htmlFile <<
" <span " << std::endl;
587 htmlFile <<
" style=\"color: rgb(0, 0, 153);\">" << runNo <<
"</span></h2>" << std::endl;
588 htmlFile <<
"<h2>Monitoring task: <span " << std::endl;
589 htmlFile <<
" style=\"color: rgb(0, 0, 153);\">Detector Diagnostics LED Monitor</span></h2> " << std::endl;
590 htmlFile <<
"<h2>Events processed: <span " << std::endl;
591 htmlFile <<
" style=\"color: rgb(0, 0, 153);\">" << ievt_ <<
"</span></h2>" << std::endl;
592 htmlFile <<
"<hr>" << std::endl;
594 htmlFile <<
"<table width=100% border=1>" << std::endl;
595 htmlFile <<
"<tr>" << std::endl;
596 htmlFile <<
"<td class=\"s0\" width=15% align=\"center\">SebDet</td>" << std::endl;
597 htmlFile <<
"<td class=\"s0\" width=17% align=\"center\">Missing</td>" << std::endl;
598 htmlFile <<
"<td class=\"s0\" width=17% align=\"center\">Unstable</td>" << std::endl;
599 htmlFile <<
"<td class=\"s0\" width=17% align=\"center\">low/no LED signal</td>" << std::endl;
600 htmlFile <<
"<td class=\"s0\" width=17% align=\"center\">Bad Timing</td>" << std::endl;
601 htmlFile <<
"<td class=\"s0\" width=17% align=\"center\">Bad LED signal</td>" << std::endl;
602 htmlFile <<
"</tr><tr>" << std::endl;
603 int ind1=0,ind2=0,ind3=0,ind4=0,ind5=0;
604 htmlFile <<
"<td class=\"s1\" align=\"center\">HB+</td>" << std::endl;
605 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;
606 ind2=3;
if(newHBP[1]==0) ind2=2;
if(newHBP[1]>0) ind2=1;
if(newHBP[1]>21) ind2=0;
607 ind3=3;
if(newHBP[2]==0) ind3=2;
if(newHBP[2]>0) ind3=1;
if(newHBP[2]>21) ind3=0;
608 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;
609 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;
610 if(ind1==3) ind2=ind3=ind4=ind5=3;
611 if(ind1==0 || ind2==0 || ind3==0 || ind4==0)
status|=2;
else if(ind1==1 || ind2==1 || ind3==1 || ind4==1)
status|=1;
612 htmlFile << state[ind1] << HBP[0] <<
" (1296)</td>" << std::endl;
613 htmlFile << state[ind2] << HBP[1] <<
"</td>" << std::endl;
614 htmlFile << state[ind3] << HBP[2] <<
"</td>" << std::endl;
615 htmlFile << state[ind5] << HBP[5]+HBP[6] <<
"</td>" << std::endl;
616 htmlFile << state[ind4] << HBP[3]+HBP[4] <<
"</td>" << std::endl;
618 htmlFile <<
"</tr><tr>" << std::endl;
619 htmlFile <<
"<td class=\"s1\" align=\"center\">HB-</td>" << std::endl;
620 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;
621 ind2=3;
if(newHBM[1]==0) ind2=2;
if(newHBM[1]>0) ind2=1;
if(newHBM[1]>21) ind2=0;
622 ind3=3;
if(newHBM[2]==0) ind3=2;
if(newHBM[2]>0) ind3=1;
if(newHBM[2]>21) ind3=0;
623 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;
624 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;
625 if(ind1==3) ind2=ind3=ind4=ind5=3;
626 if(ind1==0 || ind2==0 || ind3==0 || ind4==0)
status|=2;
else if(ind1==1 || ind2==1 || ind3==1 || ind4==1)
status|=1;
627 htmlFile << state[ind1] << HBM[0] <<
" (1296)</td>" << std::endl;
628 htmlFile << state[ind2] << HBM[1] <<
"</td>" << std::endl;
629 htmlFile << state[ind3] << HBM[2] <<
"</td>" << std::endl;
630 htmlFile << state[ind5] << HBM[5]+HBM[6] <<
"</td>" << std::endl;
631 htmlFile << state[ind4] << HBM[3]+HBM[4] <<
"</td>" << std::endl;
633 htmlFile <<
"</tr><tr>" << std::endl;
634 htmlFile <<
"<td class=\"s1\" align=\"center\">HE+</td>" << std::endl;
635 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;
636 ind2=3;
if(newHEP[1]==0) ind2=2;
if(newHEP[1]>0) ind2=1;
if(newHEP[1]>21) ind2=0;
637 ind3=3;
if(newHEP[2]==0) ind3=2;
if(newHEP[2]>0) ind3=1;
if(newHEP[2]>21) ind3=0;
638 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;
639 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;
640 if(ind1==3) ind2=ind3=ind4=ind5=3;
641 if(ind1==0 || ind2==0 || ind3==0 || ind4==0)
status|=2;
else if(ind1==1 || ind2==1 || ind3==1 || ind4==1)
status|=1;
642 htmlFile << state[ind1] << HEP[0] <<
" (1296)</td>" << std::endl;
643 htmlFile << state[ind2] << HEP[1] <<
"</td>" << std::endl;
644 htmlFile << state[ind3] << HEP[2] <<
"</td>" << std::endl;
645 htmlFile << state[ind5] << HEP[5]+HEP[6] <<
"</td>" << std::endl;
646 htmlFile << state[ind4] << HEP[3]+HEP[4] <<
"</td>" << std::endl;
648 htmlFile <<
"</tr><tr>" << std::endl;
649 htmlFile <<
"<td class=\"s1\" align=\"center\">HE-</td>" << std::endl;
650 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;
651 ind2=3;
if(newHEM[1]==0) ind2=2;
if(newHEM[1]>0) ind2=1;
if(newHEM[1]>21) ind2=0;
652 ind3=3;
if(newHEM[2]==0) ind3=2;
if(newHEM[2]>0) ind3=1;
if(newHEM[2]>21) ind3=0;
653 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;
654 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;
655 if(ind1==3) ind2=ind3=ind4=ind5=3;
656 if(ind1==0 || ind2==0 || ind3==0 || ind4==0)
status|=2;
else if(ind1==1 || ind2==1 || ind3==1 || ind4==1)
status|=1;
657 htmlFile << state[ind1] << HEM[0] <<
" (1296)</td>" << std::endl;
658 htmlFile << state[ind2] << HEM[1] <<
"</td>" << std::endl;
659 htmlFile << state[ind3] << HEM[2] <<
"</td>" << std::endl;
660 htmlFile << state[ind5] << HEM[5]+HEM[6] <<
"</td>" << std::endl;
661 htmlFile << state[ind4] << HEM[3]+HEM[4] <<
"</td>" << std::endl;
663 htmlFile <<
"</tr><tr>" << std::endl;
664 htmlFile <<
"<td class=\"s1\" align=\"center\">HF+</td>" << std::endl;
665 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;
666 ind2=3;
if(newHFP[1]==0) ind2=2;
if(newHFP[1]>0) ind2=1;
if(newHFP[1]>21) ind2=0;
667 ind3=3;
if(newHFP[2]==0) ind3=2;
if(newHFP[2]>0) ind3=1;
if(newHFP[2]>21) ind3=0;
668 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;
669 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;
670 if(ind1==3) ind2=ind3=ind4=ind5=3;
671 if(ind1==0 || ind2==0 || ind3==0 || ind4==0)
status|=2;
else if(ind1==1 || ind2==1 || ind3==1 || ind4==1)
status|=1;
672 htmlFile << state[ind1] << HFP[0] <<
" (864)</td>" << std::endl;
673 htmlFile << state[ind2] << HFP[1] <<
"</td>" << std::endl;
674 htmlFile << state[ind3] << HFP[2] <<
"</td>" << std::endl;
675 htmlFile << state[ind5] << HFP[5]+HFP[6] <<
"</td>" << std::endl;
676 htmlFile << state[ind4] << HFP[3]+HFP[4] <<
"</td>" << std::endl;
678 htmlFile <<
"</tr><tr>" << std::endl;
679 htmlFile <<
"<td class=\"s1\" align=\"center\">HF-</td>" << std::endl;
680 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;
681 ind2=3;
if(newHFM[1]==0) ind2=2;
if(newHFM[1]>0) ind2=1;
if(newHFM[1]>21) ind2=0;
682 ind3=3;
if(newHFM[2]==0) ind3=2;
if(newHFM[2]>0) ind3=1;
if(newHFM[2]>21) ind3=0;
683 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;
684 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;
685 if(ind1==3) ind2=ind3=ind4=ind5=3;
686 if(ind1==0 || ind2==0 || ind3==0 || ind4==0)
status|=2;
else if(ind1==1 || ind2==1 || ind3==1 || ind4==1)
status|=1;
687 htmlFile << state[ind1] << HFM[0] <<
" (864)</td>" << std::endl;
688 htmlFile << state[ind2] << HFM[1] <<
"</td>" << std::endl;
689 htmlFile << state[ind3] << HFM[2] <<
"</td>" << std::endl;
690 htmlFile << state[ind5] << HFM[5]+HFM[6] <<
"</td>" << std::endl;
691 htmlFile << state[ind4] << HFM[3]+HFM[4] <<
"</td>" << std::endl;
693 htmlFile <<
"</tr><tr>" << std::endl;
694 htmlFile <<
"<td class=\"s1\" align=\"center\">HO</td>" << std::endl;
695 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;
696 ind2=3;
if(newHO[1]==0) ind2=2;
if(newHO[1]>0) ind2=1;
if(newHO[1]>21) ind2=0;
697 ind3=3;
if(newHO[2]==0) ind3=2;
if(newHO[2]>0) ind3=1;
if(newHO[2]>21) ind3=0;
698 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;
699 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;
700 if(ind1==3) ind2=ind3=ind4=ind5=3;
701 if(ind1==0 || ind2==0 || ind3==0 || ind4==0)
status|=2;
else if(ind1==1 || ind2==1 || ind3==1 || ind4==1)
status|=1;
703 htmlFile << state[ind1] << HO[0] <<
" (2160)</td>" << std::endl;
704 htmlFile << state[ind2] << HO[1] <<
"</td>" << std::endl;
705 htmlFile << state[ind3] << HO[2] <<
"</td>" << std::endl;
706 htmlFile << state[ind5] << HO[5]+HO[6] <<
"</td>" << std::endl;
707 htmlFile << state[ind4] << HO[3]+HO[4] <<
"</td>" << std::endl;
709 htmlFile <<
"</tr></table>" << std::endl;
710 htmlFile <<
"<hr>" << std::endl;
712 if((MissingCnt+UnstableCnt+UnstableLEDCnt+BadTimingCnt)>0){
713 htmlFile <<
"<table width=100% border=1><tr>" << std::endl;
714 if(MissingCnt>0) htmlFile <<
"<td><a href=\"" <<
"Missing.html" <<
"\">list of missing channels</a></td>";
715 if(UnstableCnt>0) htmlFile <<
"<td><a href=\"" <<
"Unstable.html" <<
"\">list of unstable channels</a></td>";
716 if(UnstableLEDCnt>0)htmlFile <<
"<td><a href=\"" <<
"UnstableLED.html" <<
"\">list of low LED signal channels</a></td>";
717 if(BadTimingCnt>0)htmlFile <<
"<td><a href=\"" <<
"BadTiming.html" <<
"\">list of Bad Timing channels</a></td>";
718 htmlFile <<
"</tr></table>" << std::endl;
723 htmlFile <<
"<h2 align=\"center\">Summary LED plots</h2>" << std::endl;
724 htmlFile <<
"<table width=100% border=0><tr>" << std::endl;
725 htmlFile <<
"<tr align=\"left\">" << std::endl;
726 Time2Dhbhehf->SetMaximum(6);
727 Time2Dho->SetMaximum(6);
728 Time2Dhbhehf->SetNdivisions(36,
"Y");
729 Time2Dho->SetNdivisions(36,
"Y");
730 Time2Dhbhehf->SetStats(0);
731 Time2Dho->SetStats(0);
732 Time2Dhbhehf->Draw(
"COLZ"); can->SaveAs((htmlDir +
"led_timing_hbhehf.gif").c_str());
733 Time2Dho->Draw(
"COLZ"); can->SaveAs((htmlDir +
"led_timing_ho.gif").c_str());
734 htmlFile <<
"<td align=\"center\"><img src=\"led_timing_hbhehf.gif\" alt=\"led timing distribution\"> </td>" << std::endl;
735 htmlFile <<
"<td align=\"center\"><img src=\"led_timing_ho.gif\" alt=\"led timing distribution\"> </td>" << std::endl;
736 htmlFile <<
"</tr>" << std::endl;
738 htmlFile <<
"<tr align=\"left\">" << std::endl;
739 Energy2Dhbhehf->SetStats(0);
740 Energy2Dho->SetStats(0);
741 Energy2Dhbhehf->SetNdivisions(36,
"Y");
742 Energy2Dho->SetNdivisions(36,
"Y");
743 Energy2Dhbhehf->Draw(
"COLZ"); can->SaveAs((htmlDir +
"led_energy_hbhehf.gif").c_str());
744 Energy2Dho->Draw(
"COLZ"); can->SaveAs((htmlDir +
"led_energy_ho.gif").c_str());
745 htmlFile <<
"<td align=\"center\"><img src=\"led_energy_hbhehf.gif\" alt=\"led energy distribution\"> </td>" << std::endl;
746 htmlFile <<
"<td align=\"center\"><img src=\"led_energy_ho.gif\" alt=\"led energy distribution\"> </td>" << std::endl;
747 htmlFile <<
"</tr>" << std::endl;
749 can->SetGridy(
false);
750 htmlFile <<
"<tr align=\"left\">" << std::endl;
751 Energy->Draw(); can->SaveAs((htmlDir +
"led_energy_distribution.gif").c_str());
752 EnergyRMS->Draw(); can->SaveAs((htmlDir +
"led_energy_rms_distribution.gif").c_str());
753 htmlFile <<
"<td align=\"center\"><img src=\"led_energy_distribution.gif\" alt=\"led energy distribution\"> </td>" << std::endl;
754 htmlFile <<
"<td align=\"center\"><img src=\"led_energy_rms_distribution.gif\" alt=\"led energy rms distribution\"> </td>" << std::endl;
755 htmlFile <<
"</tr>" << std::endl;
756 htmlFile <<
"<tr align=\"left\">" << std::endl;
757 Timing->Draw(); can->SaveAs((htmlDir +
"led_timing_distribution.gif").c_str());
758 TimingRMS->Draw(); can->SaveAs((htmlDir +
"led_timing_rms_distribution.gif").c_str());
759 htmlFile <<
"<td align=\"center\"><img src=\"led_timing_distribution.gif\" alt=\"led timing distribution\"> </td>" << std::endl;
760 htmlFile <<
"<td align=\"center\"><img src=\"led_timing_rms_distribution.gif\" alt=\"led timing rms distribution\"> </td>" << std::endl;
761 htmlFile <<
"</tr>" << std::endl;
762 htmlFile <<
"<tr align=\"left\">" << std::endl;
763 EnergyHF->Draw(); can->SaveAs((htmlDir +
"led_energyhf_distribution.gif").c_str());
764 EnergyRMSHF->Draw(); can->SaveAs((htmlDir +
"led_energyhf_rms_distribution.gif").c_str());
765 htmlFile <<
"<td align=\"center\"><img src=\"led_energyhf_distribution.gif\" alt=\"hf led energy distribution\"> </td>" << std::endl;
766 htmlFile <<
"<td align=\"center\"><img src=\"led_energyhf_rms_distribution.gif\" alt=\"hf led energy rms distribution\"> </td>" << std::endl;
767 htmlFile <<
"</tr>" << std::endl;
768 htmlFile <<
"<tr align=\"left\">" << std::endl;
769 TimingHF->Draw(); can->SaveAs((htmlDir +
"led_timinghf_distribution.gif").c_str());
770 TimingRMSHF->Draw(); can->SaveAs((htmlDir +
"led_timinghf_rms_distribution.gif").c_str());
771 htmlFile <<
"<td align=\"center\"><img src=\"led_timinghf_distribution.gif\" alt=\"hf led timing distribution\"> </td>" << std::endl;
772 htmlFile <<
"<td align=\"center\"><img src=\"led_timinghf_rms_distribution.gif\" alt=\"hf led timing rms distribution\"> </td>" << std::endl;
773 htmlFile <<
"</tr>" << std::endl;
774 htmlFile <<
"</table>" << std::endl;
777 htmlFile <<
"<h2 align=\"center\">Stability LED plots (Reference run "<<ref_run<<
")</h2>" << std::endl;
778 htmlFile <<
"<table width=100% border=0><tr>" << std::endl;
779 htmlFile <<
"<tr align=\"left\">" << std::endl;
780 HBPphi->GetXaxis()->SetNdivisions(418,kFALSE);
781 HBMphi->GetXaxis()->SetNdivisions(418,kFALSE);
782 HEPphi->GetXaxis()->SetNdivisions(418,kFALSE);
783 HEMphi->GetXaxis()->SetNdivisions(418,kFALSE);
785 HBPphi->SetMarkerColor(kRed);
786 HBPphi->SetMarkerStyle(23);
787 HBPphi->SetXTitle(
"HPD Index = RBX*4+RM");
788 HBMphi->SetMarkerColor(kRed);
789 HBMphi->SetMarkerStyle(23);
790 HBMphi->SetXTitle(
"HPD Index = RBX*4+RM");
791 HBPphi->Draw(); can->SaveAs((htmlDir +
"led_hbp_distribution.gif").c_str());
792 HBMphi->Draw(); can->SaveAs((htmlDir +
"led_hbm_distribution.gif").c_str());
793 htmlFile <<
"<td align=\"center\"><img src=\"led_hbp_distribution.gif\" alt=\"led hbp distribution\"> </td>" << std::endl;
794 htmlFile <<
"<td align=\"center\"><img src=\"led_hbm_distribution.gif\" alt=\"led hbm distribution\"> </td>" << std::endl;
795 htmlFile <<
"</tr>" << std::endl;
797 htmlFile <<
"<tr align=\"left\">" << std::endl;
798 HEPphi->SetMarkerColor(kRed);
799 HEPphi->SetMarkerStyle(23);
800 HEPphi->SetXTitle(
"HPD Index = RBX*4+RM");
801 HEMphi->SetMarkerColor(kRed);
802 HEMphi->SetMarkerStyle(23);
803 HEMphi->SetXTitle(
"HPD Index = RBX*4+RM");
804 HEPphi->Draw(); can->SaveAs((htmlDir +
"led_hep_distribution.gif").c_str());
805 HEMphi->Draw(); can->SaveAs((htmlDir +
"led_hem_distribution.gif").c_str());
806 htmlFile <<
"<td align=\"center\"><img src=\"led_hep_distribution.gif\" alt=\"led hep distribution\"> </td>" << std::endl;
807 htmlFile <<
"<td align=\"center\"><img src=\"led_hem_distribution.gif\" alt=\"led hem distribution\"> </td>" << std::endl;
808 htmlFile <<
"</tr>" << std::endl;
810 htmlFile <<
"<tr align=\"left\">" << std::endl;
811 HFPphi->SetMarkerColor(kRed);
812 HFPphi->SetMarkerStyle(23);
813 HFPphi->SetXTitle(
"RM Index = RoBox*3+RM");
814 HFMphi->SetMarkerColor(kRed);
815 HFMphi->SetMarkerStyle(23);
816 HFPphi->GetXaxis()->SetNdivisions(312,kFALSE);
817 HFMphi->GetXaxis()->SetNdivisions(312,kFALSE);
819 HFMphi->SetXTitle(
"RM Index = RoBox*3+RM");
820 HFPphi->Draw(); can->SaveAs((htmlDir +
"led_hfp_distribution.gif").c_str());
821 HFMphi->Draw(); can->SaveAs((htmlDir +
"led_hfm_distribution.gif").c_str());
822 htmlFile <<
"<td align=\"center\"><img src=\"led_hfp_distribution.gif\" alt=\"led hfp distribution\"> </td>" << std::endl;
823 htmlFile <<
"<td align=\"center\"><img src=\"led_hfm_distribution.gif\" alt=\"led hfm distribution\"> </td>" << std::endl;
824 htmlFile <<
"</tr>" << std::endl;
826 htmlFile <<
"<tr align=\"left\">" << std::endl;
827 HO1Pphi->SetMarkerColor(kRed);
828 HO1Pphi->SetMarkerStyle(23);
829 HO1Pphi->SetXTitle(
"HPD Index = RBX*4+RM");
830 HO1Mphi->SetMarkerColor(kRed);
831 HO1Mphi->SetMarkerStyle(23);
832 HO1Mphi->GetXaxis()->SetNdivisions(412,kFALSE);
833 HO1Pphi->GetXaxis()->SetNdivisions(412,kFALSE);
835 HO1Mphi->SetXTitle(
"HPD Index = RBX*4+RM");
836 HO1Pphi->Draw(); can->SaveAs((htmlDir +
"led_ho1p_distribution.gif").c_str());
837 HO1Mphi->Draw(); can->SaveAs((htmlDir +
"led_ho1m_distribution.gif").c_str());
838 htmlFile <<
"<td align=\"center\"><img src=\"led_ho1p_distribution.gif\" alt=\"led ho1p distribution\"> </td>" << std::endl;
839 htmlFile <<
"<td align=\"center\"><img src=\"led_ho1m_distribution.gif\" alt=\"led ho1m distribution\"> </td>" << std::endl;
840 htmlFile <<
"</tr>" << std::endl;
842 htmlFile <<
"<tr align=\"left\">" << std::endl;
843 HO2Pphi->SetMarkerColor(kRed);
844 HO2Pphi->SetMarkerStyle(23);
845 HO2Pphi->SetXTitle(
"HPD Index = RBX*4+RM");
846 HO2Mphi->SetMarkerColor(kRed);
847 HO2Mphi->SetMarkerStyle(23);
848 HO2Mphi->GetXaxis()->SetNdivisions(412,kFALSE);
849 HO2Pphi->GetXaxis()->SetNdivisions(412,kFALSE);
851 HO2Mphi->SetXTitle(
"HPD Index = RBX*4+RM");
852 HO2Pphi->Draw(); can->SaveAs((htmlDir +
"led_ho2p_distribution.gif").c_str());
853 HO2Mphi->Draw(); can->SaveAs((htmlDir +
"led_ho2m_distribution.gif").c_str());
854 htmlFile <<
"<td align=\"center\"><img src=\"led_ho2p_distribution.gif\" alt=\"led ho2p distribution\"> </td>" << std::endl;
855 htmlFile <<
"<td align=\"center\"><img src=\"led_ho2m_distribution.gif\" alt=\"led ho2m distribution\"> </td>" << std::endl;
856 htmlFile <<
"</tr>" << std::endl;
858 htmlFile <<
"<tr align=\"left\">" << std::endl;
859 HO0phi->SetMarkerColor(kRed);
860 HO0phi->SetMarkerStyle(23);
861 HO0phi->SetXTitle(
"HPD Index = RBX*4+RM");
862 HO0phi->GetXaxis()->SetNdivisions(412,kFALSE);
863 HO0phi->Draw(); can->SaveAs((htmlDir +
"led_ho0_distribution.gif").c_str());
864 htmlFile <<
"<td align=\"center\"><img src=\"led_ho0_distribution.gif\" alt=\"led ho0 distribution\"> </td>" << std::endl;
865 htmlFile <<
"</tr>" << std::endl;
867 htmlFile <<
"</table>" << std::endl;
869 htmlFile <<
"</body> " << std::endl;
870 htmlFile <<
"</html> " << std::endl;
TH2F * ChannelStatusUnstableChannels[4]
static void printTableLine(std::ofstream &file, int ind, HcalDetId &detid, HcalFrontEndId &lmap_entry, HcalElectronicsId &emap_entry, std::string comment="")
static void printTableTail(std::ofstream &file)
TH2F * ChannelStatusTimeRMS[4]
MonitorElement * get(const std::string &path)
std::vector< std::string > HcalEtaPhiHistNames()
static void printTableHeader(std::ofstream &file, std::string header)
HcalElectronicsMap generateHcalElectronicsMap()
TH2F * ChannelStatusLEDMean[4]
TH2F * ChannelsLEDEnergy[4]
TH2F * ChannelStatusMissingChannels[4]
TH2F * ChannelsLEDEnergyRef[4]
double get_channel_status(std::string subdet, int eta, int phi, int depth, int type)
std::string valueString(void) const
const HcalFrontEndId getHcalFrontEndId(const DetId &)
std::map< HcalDetId, unsigned int > badstatusmap
HcalLogicalMap * logicalMap_
TH2F * ChannelStatusLEDRMS[4]
TH2F * ChannelStatusTimeMean[4]
TH2F * ChannelStatusUnstableLEDsignal[4]
Readout chain identification for Hcal.
const DetId lookup(HcalElectronicsId fId) const
lookup the logical detid associated with the given electronics id