CMS 3D CMS Logo
RemoteMonitoringMAP_Global.cc
Go to the documentation of this file.
1 // How to run:
2 //root -b -q -l RemoteMonitoringMAP.C+
3 //root -b -q -l 'RemoteMonitoringMAP.C+("/afs/cern.ch/cms/CAF/CMSALCA/ALCA_HCALCALIB/HCALMONITORING/RDMweb/histos/LED_214513.root","/afs/cern.ch/cms/CAF/CMSALCA/ALCA_HCALCALIB/HCALMONITORING/RDMweb/histos/LED_214512.root")'
4 //root -b -q -l 'RemoteMonitoringMAP.C+(" /afs/cern.ch/work/d/dtlisov/private/Monitoring/histos/LED_211659.root","/afs/cern.ch/cms/CAF/CMSALCA/ALCA_HCALCALIB/HCALMONITORING/RDMweb/histos/LED_214512.root")'
5 //
6 //
7 //
8 #include "LogEleMapdb.h"
9 
10 #include <iostream>
11 #include <fstream>
12 
13 #include "TH1.h"
14 #include "TH2.h"
15 #include "TCanvas.h"
16 #include "TROOT.h"
17 #include <TMath.h>
18 #include "TStyle.h"
19 #include "TSystem.h"
20 #include "TLegend.h"
21 #include "TText.h"
22 #include "TAxis.h"
23 #include "TFile.h"
24 #include "TLine.h"
25 #include "TGraph.h"
26 #include <TPaveText.h>
27 //#####
28 
29 #include <TChain.h>
30 #include <TH1D.h>
31 #include <TH2D.h>
32 #include <TProfile.h>
33 #include <TFitResult.h>
34 #include <TFitResultPtr.h>
35 #include <TPaveStats.h>
36 #include <vector>
37 #include <string>
38 #include <iomanip>
39 //#####
40 #include <TClass.h>
41 
42 //
43 // https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT
44 // https://cms-cpt-software.web.cern.ch/cms-cpt-software/General/Validation/SVSuite/HcalRemoteMonitoring/GlobalRMT
45 
46 using namespace std;
47 //inline void HERE(const char *msg) { std::cout << msg << std::endl; }
48 
49 int main(int argc, char *argv[]) {
50  std::string dirnm = "Analyzer";
51  gROOT->Reset();
52  gROOT->SetStyle("Plain");
53  gStyle->SetOptStat(0);
54  gStyle->SetOptTitle(1);
55 
56  gStyle->SetStatX(0.91);
57  gStyle->SetStatY(0.75);
58  gStyle->SetStatW(0.20);
59  gStyle->SetStatH(0.10);
60  // gStyle->SetStatH(0.35);
61  //
62 
63  // Float_t LeftOffset = 0.12;
64  //Float_t TopOffset = 0.12;
65  Float_t LeftOffset = 0.12;
66  Float_t TopOffset = 0.22;
67 
68  gStyle->SetLineWidth(1);
69  gStyle->SetErrorX(0);
70 
71  //---=[ Titles,Labels ]=-----------
72  gStyle->SetOptTitle(0); // title on/off
73  // gStyle->SetTitleColor(0); // title color
74  gStyle->SetTitleColor(1); // title color
75  // gStyle->SetTitleX(0.35); // title x-position
76  gStyle->SetTitleX(0.15); // title x-position
77  gStyle->SetTitleH(0.15); // title height
78  // gStyle->SetTitleW(0.53); // title width
79  gStyle->SetTitleW(0.60); // title width
80  gStyle->SetTitleFont(42); // title font
81  gStyle->SetTitleFontSize(0.07); // title font size
82 
83  gStyle->SetPalette(1);
84  //---=[ Pad style ]=----------------
85  gStyle->SetPadTopMargin(TopOffset);
86  gStyle->SetPadBottomMargin(LeftOffset);
87  gStyle->SetPadRightMargin(TopOffset);
88  gStyle->SetPadLeftMargin(LeftOffset);
89 
90  if (argc < 2)
91  return 1;
92  char fname[300];
93  char refname[300];
94  sprintf(fname, "%s", argv[1]);
95  sprintf(refname, "%s", argv[2]);
96 
97  cout << fname << " " << refname << std::endl;
98 
99  //
100 
101  //======================================================================
102  // Connect the input files, parameters and get the 2-d histogram in memory
103  // TFile *hfile= new TFile("GlobalHist.root", "READ");
104  string promt = (string)fname;
105  string runnumber = "";
106  for (unsigned int i = promt.size() - 11; i < promt.size() - 5; i++)
107  runnumber += fname[i];
108  string refrunnumber = "";
109  promt = (string)refname;
110  for (unsigned int i = promt.size() - 11; i < promt.size() - 5; i++)
111  refrunnumber += refname[i];
112 
113  TFile *hfile = new TFile(fname, "READ");
114  hfile->ls();
115  TDirectory *dir = (TDirectory *)hfile->FindObjectAny(dirnm.c_str());
116  //TFile *hreffile = new TFile(refname, "READ");
117  //megatile channels
118  //CUTS: [test][subdetector] CapID(Test=1; ADC amplitude Am(Test= 2); Width for Wm(Test=3); Ratio cut for Rm(Test=4); TS mean for TNm(test=5); TS max for TXm(Test=6);
119  double MIN_M[7][5] = {{0., 0., 0., 0., 0.},
120  {0., 0., 0., 0., 0.},
121  {0, 0., 0., 0., 0.},
122  {0, 1.0, 1.0, 0.2, 0.1},
123  {0, 0.10, 0.10, 0.18, 0.30},
124  {0, 0.8, 0.8, 0.8, 0.1},
125  {0, -0.5, -0.5, -0.5, -0.5}};
126  double MAX_M[7][5] = {{0., 0., 0., 0., 0.},
127  {0., 0., 0., 0., 0.},
128  {0, 900, 900, 9000, 3000},
129  {0, 3.9, 3.9, 4.4, 2.0},
130  {0, 0.95, 0.98, 0.96, 1.04},
131  {0, 8.0, 8.0, 8.0, 2.8},
132  {0, 6.5, 6.5, 6.5, 3.5}};
133 
134  // calibration channels:
135  double MIN_C[7][5] = {{0., 0., 0., 0., 0.},
136  {0., 0., 0., 0., 0.},
137  {0, 120., 120., 120., 60.},
138  {0, 1.0, 1.0, 0.50, 0.2},
139  {0, 0.6, 0.64, 0.25, 0.25},
140  {0, 1.0, 1.0, 1.0, 1.0},
141  {0, 0.5, 0.5, 0.5, 0.5}};
142  double MAX_C[7][5] = {{0., 0., 0., 0., 0.},
143  {0., 0., 0., 0., 0.},
144  {0, 1E20, 1E20, 1E20, 1E20},
145  {0, 2.3, 2.3, 3.0, 2.3},
146  {0, 1., 1., 1., 1.00},
147  {0, 5.5, 5.5, 3.5, 5.2},
148  {0, 8.5, 8.5, 8.5, 9.5}};
149  double porog[5] = {0., 2., 2., 5., 1.}; // Cut for GS test in pro cents
150  // double porog[5] = {0., 200., 200., 100., 100.}; // Cut for GS test in pro cents
151  double Pedest[2][5] = {{0., 0.2, 0.9, 0.1, 0.2}, {0., 0.2, 0.2, 0.1, 0.16}}; //Cuts for Pedestal and pedestal Width
152  //======================================================================
153  // with TfileService implementation, change everywhere below: hfile->Get to dir->FindObjectAny
154  //======================================================================
155  // Prepare histograms and plot them to .png files
156 
157  //TCanvas *cHB = new TCanvas("cHB","cHB",1000,500);
158  TCanvas *cHB = new TCanvas("cHB", "cHB", 1000, 1000);
159  //TCanvas *cHE = new TCanvas("cHE","cHE",1500,500);
160  TCanvas *cHE = new TCanvas("cHE", "cHE", 1500, 1500);
161  // TCanvas *cONE = new TCanvas("cONE","cONE",500,500);
162  TCanvas *cONE = new TCanvas("cONE", "cONE", 1500, 500);
163  TCanvas *cHO = new TCanvas("cHO", "cHO", 500, 500);
164  TCanvas *cPED = new TCanvas("cPED", "cPED", 1000, 500);
165  //TCanvas *cHF = new TCanvas("cHF","cHF",1000,1000);
166  TCanvas *cHF = new TCanvas("cHF", "cHF", 1000, 1000);
167 
168  char *str = (char *)alloca(10000);
169 
170  // before upgrade 2017:
171  // depth: HB depth1,2; HE depth1,2,3; HO depth4; HF depth1,2
172  // 5 depthes: 0(empty), 1,2,3,4
173 
174  // upgrade 2017:
175  // depth: HB depth1,2; HE depth1,2,3,4,5,6,7; HO depth4; HF depth1,2,3,4
176  // 8 depthes: 0(empty), 1,2,3,4,5,6,7
177 
178  // upgrade 2019:
179  // depth: HB depth1,2,3,4; HE depth1,2,3,4,5,6,7; HO depth4; HF depth1,2,3,4
180  // 10 depthes: 0(empty), 1,2,3,4,5,6,7,8,9
181 
182  // Int_t ALLDEPTH = 5;
183  // Int_t ALLDEPTH = 8;
184  Int_t ALLDEPTH = 10;
185  //massive_indx=1 2 3 4 5
186  // 0, HB,HE,HO,HF
187  int k_min[5] = {0, 1, 1, 4, 1}; // minimum depth for each subdet
188 
189  //int k_max[5]={0,2,3,4,2}; // maximum depth for each subdet
190  //int k_max[5]={0,2,7,4,4}; // maximum depth for each subdet
191  int k_max[5] = {0, 4, 7, 4, 4}; // maximum depth for each subdet
192 
193  TH2F *Map_Ampl[33][5][ALLDEPTH]; // 2D histogramm for test,subdet,depth
194  //AZ2023 TH2F *Map_SUBGOOD[5][ALLDEPTH]; // 2d histogramm for subdet, depth
195  TH2F *Map_SUB[5][ALLDEPTH]; // 2d histogramm for subdet, depth
196  TH1F *HistAmplDepth[22][5][ALLDEPTH]; // 1d histogramm for test,subdet, depth
197  TH1F *HistAmpl[22][5]; // 1d histogramm for test,subdet
198  TH2F *Map_SUBTS[5][ALLDEPTH]; // 2d histogramm for subdet, depth in different TSs
199 
200  TH1F *HistPed[3][5][4]; // 1d histogramm for test,subdet, CapID
201  TH2F *Map_Ped[3][5]; // 2d histogramm for test,subdet -> test 33
202  TH1F *hist_GoodTSshape[5]; // 1d histogramm for TS shape subdet -> test 41
203  TH1F *hist_GoodTSshape0[5]; // 1d histogramm for TS shape subdet -> test 41
204  TH1F *hist_BadTSshape[5]; // 1d histogramm for TS shape subdet -> test 41
205  TH1F *hist_BadTSshape0[5]; // 1d histogramm for TS shape subdet -> test 41
206  TH1F *hist_ADC_All[5]; // 1d histogramm for TS shape subdet -> test 42
207  TH1F *hist_ADC_DS[5][ALLDEPTH]; // 1d histogramm for TS shape subdet, depth -> test 42
208  TH1F *hist_SumADC[5][ALLDEPTH]; // 1d histogramm for TS shape subdet, depth -> test 43
209  TH1F *hist_SumADC0[5][ALLDEPTH]; // 1d histogramm for TS shape subdet, depth -> test 43
210  TH1F *hist_SumADC1[5][ALLDEPTH]; // 1d histogramm for TS shape subdet, depth -> test 43
211 
212  Map_SUB[1][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1_HB");
213  Map_SUB[1][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2_HB");
214  Map_SUB[1][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3_HB");
215  Map_SUB[1][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4_HB");
216  Map_SUB[2][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1_HE");
217  Map_SUB[2][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2_HE");
218  Map_SUB[2][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3_HE");
219  Map_SUB[2][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4_HE");
220  Map_SUB[2][5] = (TH2F *)dir->FindObjectAny("h_mapDepth5_HE");
221  Map_SUB[2][6] = (TH2F *)dir->FindObjectAny("h_mapDepth6_HE");
222  Map_SUB[2][7] = (TH2F *)dir->FindObjectAny("h_mapDepth7_HE");
223  Map_SUB[3][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4_HO");
224  Map_SUB[4][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1_HF");
225  Map_SUB[4][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2_HF");
226  Map_SUB[4][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3_HF");
227  Map_SUB[4][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4_HF");
228  //AZ2023:
229  /*
230  Map_SUBGOOD[1][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1_HB");
231  Map_SUBGOOD[1][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2_HB");
232  Map_SUBGOOD[1][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3_HB");
233  Map_SUBGOOD[1][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4_HB");
234  Map_SUBGOOD[2][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1_HE");
235  Map_SUBGOOD[2][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2_HE");
236  Map_SUBGOOD[2][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3_HE");
237  Map_SUBGOOD[2][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4_HE");
238  Map_SUBGOOD[2][5] = (TH2F *)dir->FindObjectAny("h_mapDepth5_HE");
239  Map_SUBGOOD[2][6] = (TH2F *)dir->FindObjectAny("h_mapDepth6_HE");
240  Map_SUBGOOD[2][7] = (TH2F *)dir->FindObjectAny("h_mapDepth7_HE");
241  Map_SUBGOOD[3][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4_HO");
242  Map_SUBGOOD[4][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1_HF");
243  Map_SUBGOOD[4][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2_HF");
244  Map_SUBGOOD[4][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3_HF");
245  Map_SUBGOOD[4][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4_HF");
246 */
247  //+++++++++++++++++++++++++++++
248  //Test 0 Entries
249  //+++++++++++++++++++++++++++++
250 
251  for (int sub = 1; sub <= 4; sub++) { //Subdetector: 1-HB, 2-HE, 3-HF, 4-HO
252  // if (sub==1) cHB->Divide(2,1);
253  if (sub == 1)
254  cHB->Divide(2, 2);
255  // if (sub==2) cHE->Divide(3,1);
256  if (sub == 2)
257  cHE->Divide(3, 3);
258  if (sub == 3)
259  cONE->Divide(1, 1);
260  // if (sub==4) cHF->Divide(2,1);
261  if (sub == 4)
262  cHF->Divide(2, 2);
263  // int k_min[5]={0,1,1,4,1}; // minimum depth for each subdet
264  // int k_max[5]={0,2,3,4,2}; // maximum depth for each subdet
265  // int k_max[5]={0,2,7,4,4}; // maximum depth for each subdet
266  for (int k = k_min[sub]; k <= k_max[sub]; k++) { //Depth
267  if (sub == 1)
268  cHB->cd(k);
269  if (sub == 2)
270  cHE->cd(k);
271  if (sub == 3)
272  cONE->cd(k - 3);
273  if (sub == 4)
274  cHF->cd(k);
275  gPad->SetGridy();
276  gPad->SetGridx();
277  gPad->SetLogz();
278  if (sub == 1)
279  sprintf(str, "HB, Depth%d \b", k);
280  if (sub == 2)
281  sprintf(str, "HE, Depth%d \b", k);
282  if (sub == 3)
283  sprintf(str, "HO, Depth%d \b", k);
284  if (sub == 4)
285  sprintf(str, "HF, Depth%d \b", k);
286  Map_SUB[sub][k]->SetTitle(str);
287  Map_SUB[sub][k]->SetXTitle("#eta \b");
288  Map_SUB[sub][k]->SetYTitle("#phi \b");
289  Map_SUB[sub][k]->SetZTitle("Number of events \b");
290  if (sub == 3)
291  Map_SUB[sub][k]->SetTitleOffset(0.8, "Z");
292  Map_SUB[sub][k]->Draw("COLZ");
293  Map_SUB[sub][k]->GetYaxis()->SetRangeUser(0, 72.);
294  // Map_SUB[sub][k]->GetZaxis()->SetRangeUser(0.0001, 1.);
295  if (sub == 1) {
296  cHB->Modified();
297  cHB->Update();
298  }
299  if (sub == 2) {
300  cHE->Modified();
301  cHE->Update();
302  }
303  if (sub == 3) {
304  cONE->Modified();
305  cONE->Update();
306  }
307  if (sub == 4) {
308  cHF->Modified();
309  cHF->Update();
310  }
311  } //end depth
312 
313  if (sub == 1) {
314  cHB->Print("MapRateEntryHB.png");
315  cHB->Clear();
316  }
317  if (sub == 2) {
318  cHE->Print("MapRateEntryHE.png");
319  cHE->Clear();
320  }
321  if (sub == 3) {
322  cONE->Print("MapRateEntryHO.png");
323  cONE->Clear();
324  }
325  if (sub == 4) {
326  cHF->Print("MapRateEntryHF.png");
327  cHF->Clear();
328  }
329  } // end sub
330 
331  //+++++++++++++++++++++++++++++
332  //Test 1 (Cm) Rate of Cap ID errors
333  //+++++++++++++++++++++++++++++
334 
335  Map_Ampl[1][1][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1Error_HB");
336  Map_Ampl[1][1][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2Error_HB");
337  Map_Ampl[1][2][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1Error_HE");
338  Map_Ampl[1][2][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2Error_HE");
339  Map_Ampl[1][2][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3Error_HE");
340  Map_Ampl[1][3][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4Error_HO");
341  Map_Ampl[1][4][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1Error_HF");
342  Map_Ampl[1][4][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2Error_HF");
343 
344  Map_Ampl[1][2][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4Error_HE");
345  Map_Ampl[1][2][5] = (TH2F *)dir->FindObjectAny("h_mapDepth5Error_HE");
346  Map_Ampl[1][2][6] = (TH2F *)dir->FindObjectAny("h_mapDepth6Error_HE");
347  Map_Ampl[1][2][7] = (TH2F *)dir->FindObjectAny("h_mapDepth7Error_HE");
348  Map_Ampl[1][4][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3Error_HF");
349  Map_Ampl[1][4][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4Error_HF");
350 
351  Map_Ampl[1][1][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3Error_HB");
352  Map_Ampl[1][1][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4Error_HB");
353 
354  for (int sub = 1; sub <= 4; sub++) { //Subdetector: 1-HB, 2-HE, 3-HF, 4-HO
355  // if (sub==1) cHB->Divide(2,1);
356  if (sub == 1)
357  cHB->Divide(2, 2);
358  // if (sub==2) cHE->Divide(3,1);
359  if (sub == 2)
360  cHE->Divide(3, 3);
361  if (sub == 3)
362  cONE->Divide(1, 1);
363  // if (sub==4) cHF->Divide(2,1);
364  if (sub == 4)
365  cHF->Divide(2, 2);
366  // int k_min[5]={0,1,1,4,1}; // minimum depth for each subdet
367  // int k_max[5]={0,2,3,4,2}; // maximum depth for each subdet
368  // int k_max[5]={0,2,7,4,4}; // maximum depth for each subdet
369  for (int k = k_min[sub]; k <= k_max[sub]; k++) { //Depth
370  if (sub == 1)
371  cHB->cd(k);
372  if (sub == 2)
373  cHE->cd(k);
374  if (sub == 3)
375  cONE->cd(k - 3);
376  if (sub == 4)
377  cHF->cd(k);
378  Map_Ampl[1][sub][k]->Divide(Map_Ampl[1][sub][k], Map_SUB[sub][k], 1, 1, "B");
379  gPad->SetGridy();
380  gPad->SetGridx();
381  gPad->SetLogz();
382  if (sub == 1)
383  sprintf(str, "HB, Depth%d \b", k);
384  if (sub == 2)
385  sprintf(str, "HE, Depth%d \b", k);
386  if (sub == 3)
387  sprintf(str, "HO, Depth%d \b", k);
388  if (sub == 4)
389  sprintf(str, "HF, Depth%d \b", k);
390  Map_Ampl[1][sub][k]->SetTitle(str);
391  Map_Ampl[1][sub][k]->SetXTitle("#eta \b");
392  Map_Ampl[1][sub][k]->SetYTitle("#phi \b");
393  Map_Ampl[1][sub][k]->SetZTitle("Rate \b");
394  if (sub == 3)
395  Map_Ampl[1][sub][k]->SetTitleOffset(0.8, "Z");
396  Map_Ampl[1][sub][k]->Draw("COLZ");
397  Map_Ampl[1][sub][k]->GetYaxis()->SetRangeUser(0, 72.);
398  Map_Ampl[1][sub][k]->GetZaxis()->SetRangeUser(0.0001, 1.);
399  if (sub == 1) {
400  cHB->Modified();
401  cHB->Update();
402  }
403  if (sub == 2) {
404  cHE->Modified();
405  cHE->Update();
406  }
407  if (sub == 3) {
408  cONE->Modified();
409  cONE->Update();
410  }
411  if (sub == 4) {
412  cHF->Modified();
413  cHF->Update();
414  }
415  } //end depth
416 
417  if (sub == 1) {
418  cHB->Print("MapRateCapIDHB.png");
419  cHB->Clear();
420  }
421  if (sub == 2) {
422  cHE->Print("MapRateCapIDHE.png");
423  cHE->Clear();
424  }
425  if (sub == 3) {
426  cONE->Print("MapRateCapIDHO.png");
427  cONE->Clear();
428  }
429  if (sub == 4) {
430  cHF->Print("MapRateCapIDHF.png");
431  cHF->Clear();
432  }
433  } // end sub
434 
435  //+++++++++++++++++++++++++++++
436  //Test 2 (Am) ADC amplitude
437  //+++++++++++++++++++++++++++++
438 
439  Map_Ampl[2][1][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1ADCAmpl225_HB");
440  Map_Ampl[2][1][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2ADCAmpl225_HB");
441  Map_Ampl[2][2][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1ADCAmpl225_HE");
442  Map_Ampl[2][2][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2ADCAmpl225_HE");
443  Map_Ampl[2][2][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3ADCAmpl225_HE");
444  Map_Ampl[2][3][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4ADCAmpl225_HO");
445  Map_Ampl[2][4][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1ADCAmpl225_HF");
446  Map_Ampl[2][4][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2ADCAmpl225_HF");
447 
448  Map_Ampl[2][2][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4ADCAmpl225_HE");
449  Map_Ampl[2][2][5] = (TH2F *)dir->FindObjectAny("h_mapDepth5ADCAmpl225_HE");
450  Map_Ampl[2][2][6] = (TH2F *)dir->FindObjectAny("h_mapDepth6ADCAmpl225_HE");
451  Map_Ampl[2][2][7] = (TH2F *)dir->FindObjectAny("h_mapDepth7ADCAmpl225_HE");
452  Map_Ampl[2][4][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3ADCAmpl225_HF");
453  Map_Ampl[2][4][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4ADCAmpl225_HF");
454 
455  Map_Ampl[2][1][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3ADCAmpl225_HB");
456  Map_Ampl[2][1][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4ADCAmpl225_HB");
457 
458  HistAmpl[2][1] = (TH1F *)dir->FindObjectAny("h_ADCAmpl_HB");
459  HistAmpl[2][2] = (TH1F *)dir->FindObjectAny("h_ADCAmpl_HE");
460  HistAmpl[2][3] = (TH1F *)dir->FindObjectAny("h_ADCAmpl_HO");
461  HistAmpl[2][4] = (TH1F *)dir->FindObjectAny("h_ADCAmpl_HF");
462 
463  //+++++++++++++++++++++++++++++
464  //Test 3 (Wm) Rate of RMS
465  //+++++++++++++++++++++++++++++
466 
467  Map_Ampl[3][1][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1Amplitude225_HB");
468  Map_Ampl[3][1][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2Amplitude225_HB");
469  Map_Ampl[3][2][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1Amplitude225_HE");
470  Map_Ampl[3][2][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2Amplitude225_HE");
471  Map_Ampl[3][2][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3Amplitude225_HE");
472  Map_Ampl[3][3][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4Amplitude225_HO");
473  Map_Ampl[3][4][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1Amplitude225_HF");
474  Map_Ampl[3][4][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2Amplitude225_HF");
475 
476  Map_Ampl[3][2][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4Amplitude225_HE");
477  Map_Ampl[3][2][5] = (TH2F *)dir->FindObjectAny("h_mapDepth5Amplitude225_HE");
478  Map_Ampl[3][2][6] = (TH2F *)dir->FindObjectAny("h_mapDepth6Amplitude225_HE");
479  Map_Ampl[3][2][7] = (TH2F *)dir->FindObjectAny("h_mapDepth7Amplitude225_HE");
480  Map_Ampl[3][4][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3Amplitude225_HF");
481  Map_Ampl[3][4][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4Amplitude225_HF");
482 
483  Map_Ampl[3][1][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3Amplitude225_HB");
484  Map_Ampl[3][1][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4Amplitude225_HB");
485 
486  HistAmpl[3][1] = (TH1F *)dir->FindObjectAny("h_Amplitude_HB");
487  HistAmpl[3][2] = (TH1F *)dir->FindObjectAny("h_Amplitude_HE");
488  HistAmpl[3][3] = (TH1F *)dir->FindObjectAny("h_Amplitude_HO");
489  HistAmpl[3][4] = (TH1F *)dir->FindObjectAny("h_Amplitude_HF");
490 
491  //+++++++++++++++++++++++++++++
492  //Test 4 (Rm) Rate of ratio 4 near max TS/ All TS
493  //+++++++++++++++++++++++++++++
494 
495  Map_Ampl[4][1][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1Ampl047_HB");
496  Map_Ampl[4][1][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2Ampl047_HB");
497  Map_Ampl[4][2][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1Ampl047_HE");
498  Map_Ampl[4][2][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2Ampl047_HE");
499  Map_Ampl[4][2][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3Ampl047_HE");
500  Map_Ampl[4][3][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4Ampl047_HO");
501  Map_Ampl[4][4][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1Ampl047_HF");
502  Map_Ampl[4][4][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2Ampl047_HF");
503 
504  Map_Ampl[4][2][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4Ampl047_HE");
505  Map_Ampl[4][2][5] = (TH2F *)dir->FindObjectAny("h_mapDepth5Ampl047_HE");
506  Map_Ampl[4][2][6] = (TH2F *)dir->FindObjectAny("h_mapDepth6Ampl047_HE");
507  Map_Ampl[4][2][7] = (TH2F *)dir->FindObjectAny("h_mapDepth7Ampl047_HE");
508  Map_Ampl[4][4][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3Ampl047_HF");
509  Map_Ampl[4][4][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4Ampl047_HF");
510 
511  Map_Ampl[4][1][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3Ampl047_HB");
512  Map_Ampl[4][1][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4Ampl047_HB");
513 
514  HistAmpl[4][1] = (TH1F *)dir->FindObjectAny("h_Ampl_HB");
515  HistAmpl[4][2] = (TH1F *)dir->FindObjectAny("h_Ampl_HE");
516  HistAmpl[4][3] = (TH1F *)dir->FindObjectAny("h_Ampl_HO");
517  HistAmpl[4][4] = (TH1F *)dir->FindObjectAny("h_Ampl_HF");
518 
519  //+++++++++++++++++++++++++++++
520  //Test 5 (TNm) Mean position in 1-8 TS range
521  //+++++++++++++++++++++++++++++
522 
523  Map_Ampl[5][1][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA225_HB");
524  Map_Ampl[5][1][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA225_HB");
525  Map_Ampl[5][2][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA225_HE");
526  Map_Ampl[5][2][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA225_HE");
527  Map_Ampl[5][2][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmeanA225_HE");
528  Map_Ampl[5][3][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA225_HO");
529  Map_Ampl[5][4][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA225_HF");
530  Map_Ampl[5][4][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA225_HF");
531 
532  Map_Ampl[5][2][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA225_HE");
533  Map_Ampl[5][2][5] = (TH2F *)dir->FindObjectAny("h_mapDepth5TSmeanA225_HE");
534  Map_Ampl[5][2][6] = (TH2F *)dir->FindObjectAny("h_mapDepth6TSmeanA225_HE");
535  Map_Ampl[5][2][7] = (TH2F *)dir->FindObjectAny("h_mapDepth7TSmeanA225_HE");
536  Map_Ampl[5][4][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmeanA225_HF");
537  Map_Ampl[5][4][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA225_HF");
538 
539  Map_Ampl[5][1][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmeanA225_HB");
540  Map_Ampl[5][1][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA225_HB");
541 
542  HistAmpl[5][1] = (TH1F *)dir->FindObjectAny("h_TSmeanA_HB");
543  HistAmpl[5][2] = (TH1F *)dir->FindObjectAny("h_TSmeanA_HE");
544  HistAmpl[5][3] = (TH1F *)dir->FindObjectAny("h_TSmeanA_HO");
545  HistAmpl[5][4] = (TH1F *)dir->FindObjectAny("h_TSmeanA_HF");
546 
547  //+++++++++++++++++++++++++++++
548  //Test 6 (TXm) Maximum position in 1-8 TS range
549  //+++++++++++++++++++++++++++++
550 
551  Map_Ampl[6][1][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmaxA225_HB");
552  Map_Ampl[6][1][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmaxA225_HB");
553  Map_Ampl[6][2][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmaxA225_HE");
554  Map_Ampl[6][2][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmaxA225_HE");
555  Map_Ampl[6][2][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmaxA225_HE");
556  Map_Ampl[6][3][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmaxA225_HO");
557  Map_Ampl[6][4][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmaxA225_HF");
558  Map_Ampl[6][4][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmaxA225_HF");
559 
560  Map_Ampl[6][2][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmaxA225_HE");
561  Map_Ampl[6][2][5] = (TH2F *)dir->FindObjectAny("h_mapDepth5TSmaxA225_HE");
562  Map_Ampl[6][2][6] = (TH2F *)dir->FindObjectAny("h_mapDepth6TSmaxA225_HE");
563  Map_Ampl[6][2][7] = (TH2F *)dir->FindObjectAny("h_mapDepth7TSmaxA225_HE");
564  Map_Ampl[6][4][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmaxA225_HF");
565  Map_Ampl[6][4][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmaxA225_HF");
566 
567  Map_Ampl[6][1][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmaxA225_HB");
568  Map_Ampl[6][1][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmaxA225_HB");
569 
570  HistAmpl[6][1] = (TH1F *)dir->FindObjectAny("h_TSmaxA_HB");
571  HistAmpl[6][2] = (TH1F *)dir->FindObjectAny("h_TSmaxA_HE");
572  HistAmpl[6][3] = (TH1F *)dir->FindObjectAny("h_TSmaxA_HO");
573  HistAmpl[6][4] = (TH1F *)dir->FindObjectAny("h_TSmaxA_HF");
574 
575  for (int test = 2; test <= 6; test++) { //Test: 2-Am, 3-Wm, 4-Rm, 5-TNm, 6-TXm,
576  for (int sub = 1; sub <= 4; sub++) { //Subdetector: 1-HB, 2-HE, 3-HF, 4-HO
577  // if (sub==1) cHB->Divide(2,1);
578  if (sub == 1)
579  cHB->Divide(2, 2);
580  // if (sub==2) cHE->Divide(3,1);
581  if (sub == 2)
582  cHE->Divide(3, 3);
583  if (sub == 3)
584  cONE->Divide(1, 1);
585  // if (sub==4) cHF->Divide(2,1);
586  if (sub == 4)
587  cHF->Divide(2, 2);
588  // int k_min[5]={0,1,1,4,1}; // minimum depth for each subdet
589  // int k_max[5]={0,2,3,4,2}; // maximum depth for each subdet
590  // int k_max[5]={0,2,7,4,4}; // maximum depth for each subdet
591  for (int k = k_min[sub]; k <= k_max[sub]; k++) { //Depth
592  if (sub == 1)
593  cHB->cd(k);
594  if (sub == 2)
595  cHE->cd(k);
596  if (sub == 3)
597  cONE->cd(k - 3);
598  if (sub == 4)
599  cHF->cd(k);
600  Map_Ampl[test][sub][k]->Divide(Map_Ampl[test][sub][k], Map_SUB[sub][k], 1, 1, "B");
601  gPad->SetGridy();
602  gPad->SetGridx();
603  gPad->SetLogz();
604  if (sub == 1)
605  sprintf(str, "HB, Depth%d \b", k);
606  if (sub == 2)
607  sprintf(str, "HE, Depth%d \b", k);
608  if (sub == 3)
609  sprintf(str, "HO, Depth%d \b", k);
610  if (sub == 4)
611  sprintf(str, "HF, Depth%d \b", k);
612  Map_Ampl[test][sub][k]->SetTitle(str);
613  Map_Ampl[test][sub][k]->SetXTitle("#eta \b");
614  Map_Ampl[test][sub][k]->SetYTitle("#phi \b");
615  Map_Ampl[test][sub][k]->SetZTitle("Rate \b");
616  if (sub == 3)
617  Map_Ampl[test][sub][k]->SetTitleOffset(0.8, "Z");
618  Map_Ampl[test][sub][k]->Draw("COLZ");
619  Map_Ampl[test][sub][k]->GetYaxis()->SetRangeUser(0, 72.);
620  Map_Ampl[test][sub][k]->GetZaxis()->SetRangeUser(0.0001, 1.);
621  if (sub == 1) {
622  cHB->Modified();
623  cHB->Update();
624  }
625  if (sub == 2) {
626  cHE->Modified();
627  cHE->Update();
628  }
629  if (sub == 3) {
630  cONE->Modified();
631  cONE->Update();
632  }
633  if (sub == 4) {
634  cHF->Modified();
635  cHF->Update();
636  }
637  } //end depth
638  if (test == 2) {
639  if (sub == 1) {
640  cHB->Print("MapRateAmplHB.png");
641  cHB->Clear();
642  }
643  if (sub == 2) {
644  cHE->Print("MapRateAmplHE.png");
645  cHE->Clear();
646  }
647  if (sub == 3) {
648  cONE->Print("MapRateAmplHO.png");
649  cONE->Clear();
650  }
651  if (sub == 4) {
652  cHF->Print("MapRateAmplHF.png");
653  cHF->Clear();
654  }
655  }
656  if (test == 3) {
657  if (sub == 1) {
658  cHB->Print("MapRateRMSHB.png");
659  cHB->Clear();
660  }
661  if (sub == 2) {
662  cHE->Print("MapRateRMSHE.png");
663  cHE->Clear();
664  }
665  if (sub == 3) {
666  cONE->Print("MapRateRMSHO.png");
667  cONE->Clear();
668  }
669  if (sub == 4) {
670  cHF->Print("MapRateRMSHF.png");
671  cHF->Clear();
672  }
673  }
674  if (test == 4) {
675  if (sub == 1) {
676  cHB->Print("MapRate43TStoAllTSHB.png");
677  cHB->Clear();
678  }
679  if (sub == 2) {
680  cHE->Print("MapRate43TStoAllTSHE.png");
681  cHE->Clear();
682  }
683  if (sub == 3) {
684  cONE->Print("MapRate43TStoAllTSHO.png");
685  cONE->Clear();
686  }
687  if (sub == 4) {
688  cHF->Print("MapRate43TStoAllTSHF.png");
689  cHF->Clear();
690  }
691  }
692  if (test == 5) {
693  if (sub == 1) {
694  cHB->Print("MapRateMeanPosHB.png");
695  cHB->Clear();
696  }
697  if (sub == 2) {
698  cHE->Print("MapRateMeanPosHE.png");
699  cHE->Clear();
700  }
701  if (sub == 3) {
702  cONE->Print("MapRateMeanPosHO.png");
703  cONE->Clear();
704  }
705  if (sub == 4) {
706  cHF->Print("MapRateMeanPosHF.png");
707  cHF->Clear();
708  }
709  }
710  if (test == 6) {
711  if (sub == 1) {
712  cHB->Print("MapRateMaxPosHB.png");
713  cHB->Clear();
714  }
715  if (sub == 2) {
716  cHE->Print("MapRateMaxPosHE.png");
717  cHE->Clear();
718  }
719  if (sub == 3) {
720  cONE->Print("MapRateMaxPosHO.png");
721  cONE->Clear();
722  }
723  if (sub == 4) {
724  cHF->Print("MapRateMaxPosHF.png");
725  cHF->Clear();
726  }
727  }
728 
729  // cONE->Divide(1,1);
730  /*
731  cONE->Divide(2,1);
732  if(test == 2 && sub == 2 ) {
733  cONE->cd(2);
734  TH1F *kjkjkhj2= (TH1F*)dir->FindObjectAny("h_ADCAmpl_HE");kjkjkhj2->Draw("");kjkjkhj2->SetTitle("HE, All Depth: shunt6");
735  }
736  if(test == 2 && sub == 1 ) {
737  cONE->cd(2);
738  TH1F *kjkjkhj1= (TH1F*)dir->FindObjectAny("h_ADCAmpl_HB");kjkjkhj1->Draw("");kjkjkhj1->SetTitle("HB, All Depth: shunt6");
739  }
740 */
741 
742  cONE->Divide(3, 1);
743  if (test == 2 && sub == 2) {
744  cONE->cd(2);
745  TH1F *kjkjkhj2 = (TH1F *)dir->FindObjectAny("h_AmplitudeHEtest1");
746  gPad->SetGridy();
747  gPad->SetGridx();
748  gPad->SetLogy();
749  kjkjkhj2->Draw("");
750  kjkjkhj2->SetXTitle("HE, All Depth: shunt1");
751  cONE->cd(3);
752  TH1F *kjkjkhj3 = (TH1F *)dir->FindObjectAny("h_AmplitudeHEtest6");
753  gPad->SetGridy();
754  gPad->SetGridx();
755  gPad->SetLogy();
756  kjkjkhj3->Draw("");
757  kjkjkhj3->SetXTitle("HE, All Depth: shunt6");
758  }
759  if (test == 2 && sub == 1) {
760  cONE->cd(2);
761  TH1F *kjkjkhb2 = (TH1F *)dir->FindObjectAny("h_AmplitudeHBtest1");
762  gPad->SetGridy();
763  gPad->SetGridx();
764  gPad->SetLogy();
765  kjkjkhb2->Draw("");
766  kjkjkhb2->SetXTitle("HB, All Depth: shunt1");
767  cONE->cd(3);
768  TH1F *kjkjkhb3 = (TH1F *)dir->FindObjectAny("h_AmplitudeHBtest6");
769  gPad->SetGridy();
770  gPad->SetGridx();
771  gPad->SetLogy();
772  kjkjkhb3->Draw("");
773  kjkjkhb3->SetXTitle("HB, All Depth: shunt6");
774  }
775 
776  cONE->cd(1);
777  gPad->SetGridy();
778  gPad->SetGridx();
779  gPad->SetLogy();
780  if (sub == 1)
781  HistAmpl[test][sub]->SetTitle("HB, All Depth: shunt6");
782  if (sub == 2)
783  HistAmpl[test][sub]->SetTitle("HE, All Depth: shunt6");
784  if (sub == 3)
785  HistAmpl[test][sub]->SetTitle("HO, All Depth");
786  if (sub == 4)
787  HistAmpl[test][sub]->SetTitle("HF, All Depth");
788  if (test == 2)
789  HistAmpl[test][sub]->SetXTitle("ADC Amlitude in each event & cell \b");
790  if (test == 3)
791  HistAmpl[test][sub]->SetXTitle("RMS in each event & cell \b");
792  if (test == 4)
793  HistAmpl[test][sub]->SetXTitle("Ratio in each event & cell \b");
794  if (test == 5)
795  HistAmpl[test][sub]->SetXTitle("Mean TS position in each event & cell \b");
796  if (test == 6)
797  HistAmpl[test][sub]->SetXTitle("Max TS position in each event & cell \b");
798  HistAmpl[test][sub]->SetYTitle("Number of cell-events \b");
799  HistAmpl[test][sub]->SetLineColor(4);
800  HistAmpl[test][sub]->SetLineWidth(2);
801  HistAmpl[test][sub]->SetTitleOffset(1.4, "Y");
802  HistAmpl[test][sub]->Draw("");
803  // // HistAmpl[test][sub]->GetYaxis()->SetRangeUser(1., 100.);
804  // if (test==2) {gPad->SetLogx(); HistAmpl[test][sub]->GetXaxis()->SetRangeUser(1., 10000.);}
805  if (test == 2) {
806  gPad->SetLogx();
807  }
808  if (test == 3)
809  HistAmpl[test][sub]->GetXaxis()->SetRangeUser(0., 5.); // width
810  if (test == 4)
811  HistAmpl[test][sub]->GetXaxis()->SetRangeUser(0., 1.); // R
812  if (test == 5)
813  HistAmpl[test][sub]->GetXaxis()->SetRangeUser(0., 9.); // Tn
814  if (test == 6)
815  HistAmpl[test][sub]->GetXaxis()->SetRangeUser(0., 9.); // Tx
816  cONE->Modified();
817  cONE->Update();
818  double min_x[] = {MIN_M[test][sub], MIN_M[test][sub]};
819  double min_y[] = {0., 100000000.};
820  TGraph *MIN = new TGraph(2, min_x, min_y);
821  MIN->SetLineStyle(2);
822  MIN->SetLineColor(2);
823  MIN->SetLineWidth(2 + 100 * 100);
824  MIN->SetFillStyle(3005);
825  MIN->SetFillColor(2);
826  MIN->Draw("L");
827  double max_x[] = {MAX_M[test][sub], MAX_M[test][sub]};
828  double max_y[] = {0., 100000000.};
829  TGraph *MAX = new TGraph(2, max_x, max_y);
830  MAX->SetLineStyle(2);
831  MAX->SetLineColor(2);
832  MAX->SetLineWidth(-2 - 100 * 100);
833  MAX->SetFillStyle(3004);
834  MAX->SetFillColor(2);
835  MAX->Draw("L");
836  if (test == 2) {
837  if (sub == 1) {
838  cONE->Print("HistAmplHB.png");
839  cONE->Clear();
840  }
841  if (sub == 2) {
842  cONE->Print("HistAmplHE.png");
843  cONE->Clear();
844  }
845  if (sub == 3) {
846  cONE->Print("HistAmplHO.png");
847  cONE->Clear();
848  }
849  if (sub == 4) {
850  cONE->Print("HistAmplHF.png");
851  cONE->Clear();
852  }
853  }
854  if (test == 3) {
855  if (sub == 1) {
856  cONE->Print("HistRMSHB.png");
857  cONE->Clear();
858  }
859  if (sub == 2) {
860  cONE->Print("HistRMSHE.png");
861  cONE->Clear();
862  }
863  if (sub == 3) {
864  cONE->Print("HistRMSHO.png");
865  cONE->Clear();
866  }
867  if (sub == 4) {
868  cONE->Print("HistRMSHF.png");
869  cONE->Clear();
870  }
871  }
872  if (test == 4) {
873  if (sub == 1) {
874  cONE->Print("Hist43TStoAllTSHB.png");
875  cONE->Clear();
876  }
877  if (sub == 2) {
878  cONE->Print("Hist43TStoAllTSHE.png");
879  cONE->Clear();
880  }
881  if (sub == 3) {
882  cONE->Print("Hist43TStoAllTSHO.png");
883  cONE->Clear();
884  }
885  if (sub == 4) {
886  cONE->Print("Hist43TStoAllTSHF.png");
887  cONE->Clear();
888  }
889  }
890  if (test == 5) {
891  if (sub == 1) {
892  cONE->Print("HistMeanPosHB.png");
893  cONE->Clear();
894  }
895  if (sub == 2) {
896  cONE->Print("HistMeanPosHE.png");
897  cONE->Clear();
898  }
899  if (sub == 3) {
900  cONE->Print("HistMeanPosHO.png");
901  cONE->Clear();
902  }
903  if (sub == 4) {
904  cONE->Print("HistMeanPosHF.png");
905  cONE->Clear();
906  }
907  }
908  if (test == 6) {
909  if (sub == 1) {
910  cONE->Print("HistMaxPosHB.png");
911  cONE->Clear();
912  }
913  if (sub == 2) {
914  cONE->Print("HistMaxPosHE.png");
915  cONE->Clear();
916  }
917  if (sub == 3) {
918  cONE->Print("HistMaxPosHO.png");
919  cONE->Clear();
920  }
921  if (sub == 4) {
922  cONE->Print("HistMaxPosHF.png");
923  cONE->Clear();
924  }
925  }
926  } // end sub
927  } //end test
928 
929  //+++++++++++++++++++++++++++++++++++
930  //Test 31, 32 Pedestal, pedestalWidths
931  //++++++++++++++++++++++++++++++++++++
932 
933  Map_Ampl[31][1][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1pedestal_HB");
934  Map_Ampl[31][1][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2pedestal_HB");
935  Map_Ampl[31][1][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3pedestal_HB");
936  Map_Ampl[31][1][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4pedestal_HB");
937  Map_Ampl[31][2][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1pedestal_HE");
938  Map_Ampl[31][2][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2pedestal_HE");
939  Map_Ampl[31][2][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3pedestal_HE");
940  Map_Ampl[31][2][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4pedestal_HE");
941  Map_Ampl[31][2][5] = (TH2F *)dir->FindObjectAny("h_mapDepth5pedestal_HE");
942  Map_Ampl[31][2][6] = (TH2F *)dir->FindObjectAny("h_mapDepth6pedestal_HE");
943  Map_Ampl[31][2][7] = (TH2F *)dir->FindObjectAny("h_mapDepth7pedestal_HE");
944  Map_Ampl[31][3][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4pedestal_HO");
945  Map_Ampl[31][4][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1pedestal_HF");
946  Map_Ampl[31][4][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2pedestal_HF");
947  Map_Ampl[31][4][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3pedestal_HF");
948  Map_Ampl[31][4][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4pedestal_HF");
949 
950  Map_Ampl[32][1][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1pedestalw_HB");
951  Map_Ampl[32][1][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2pedestalw_HB");
952  Map_Ampl[32][1][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3pedestalw_HB");
953  Map_Ampl[32][1][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4pedestalw_HB");
954  Map_Ampl[32][2][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1pedestalw_HE");
955  Map_Ampl[32][2][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2pedestalw_HE");
956  Map_Ampl[32][2][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3pedestalw_HE");
957  Map_Ampl[32][2][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4pedestalw_HE");
958  Map_Ampl[32][2][5] = (TH2F *)dir->FindObjectAny("h_mapDepth5pedestalw_HE");
959  Map_Ampl[32][2][6] = (TH2F *)dir->FindObjectAny("h_mapDepth6pedestalw_HE");
960  Map_Ampl[32][2][7] = (TH2F *)dir->FindObjectAny("h_mapDepth7pedestalw_HE");
961  Map_Ampl[32][3][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4pedestalw_HO");
962  Map_Ampl[32][4][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1pedestalw_HF");
963  Map_Ampl[32][4][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2pedestalw_HF");
964  Map_Ampl[32][4][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3pedestalw_HF");
965  Map_Ampl[32][4][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4pedestalw_HF");
966 
967  HistPed[1][1][0] = (TH1F *)dir->FindObjectAny("h_pedestal0_HB");
968  HistPed[1][1][1] = (TH1F *)dir->FindObjectAny("h_pedestal1_HB");
969  HistPed[1][1][2] = (TH1F *)dir->FindObjectAny("h_pedestal2_HB");
970  HistPed[1][1][3] = (TH1F *)dir->FindObjectAny("h_pedestal3_HB");
971  HistPed[2][1][0] = (TH1F *)dir->FindObjectAny("h_pedestalw0_HB");
972  HistPed[2][1][1] = (TH1F *)dir->FindObjectAny("h_pedestalw1_HB");
973  HistPed[2][1][2] = (TH1F *)dir->FindObjectAny("h_pedestalw2_HB");
974  HistPed[2][1][3] = (TH1F *)dir->FindObjectAny("h_pedestalw3_HB");
975 
976  HistPed[1][2][0] = (TH1F *)dir->FindObjectAny("h_pedestal0_HE");
977  HistPed[1][2][1] = (TH1F *)dir->FindObjectAny("h_pedestal1_HE");
978  HistPed[1][2][2] = (TH1F *)dir->FindObjectAny("h_pedestal2_HE");
979  HistPed[1][2][3] = (TH1F *)dir->FindObjectAny("h_pedestal3_HE");
980  HistPed[2][2][0] = (TH1F *)dir->FindObjectAny("h_pedestalw0_HE");
981  HistPed[2][2][1] = (TH1F *)dir->FindObjectAny("h_pedestalw1_HE");
982  HistPed[2][2][2] = (TH1F *)dir->FindObjectAny("h_pedestalw2_HE");
983  HistPed[2][2][3] = (TH1F *)dir->FindObjectAny("h_pedestalw3_HE");
984 
985  HistPed[1][3][0] = (TH1F *)dir->FindObjectAny("h_pedestal0_HO");
986  HistPed[1][3][1] = (TH1F *)dir->FindObjectAny("h_pedestal1_HO");
987  HistPed[1][3][2] = (TH1F *)dir->FindObjectAny("h_pedestal2_HO");
988  HistPed[1][3][3] = (TH1F *)dir->FindObjectAny("h_pedestal3_HO");
989  HistPed[2][3][0] = (TH1F *)dir->FindObjectAny("h_pedestalw0_HO");
990  HistPed[2][3][1] = (TH1F *)dir->FindObjectAny("h_pedestalw1_HO");
991  HistPed[2][3][2] = (TH1F *)dir->FindObjectAny("h_pedestalw2_HO");
992  HistPed[2][3][3] = (TH1F *)dir->FindObjectAny("h_pedestalw3_HO");
993 
994  HistPed[1][4][0] = (TH1F *)dir->FindObjectAny("h_pedestal0_HF");
995  HistPed[1][4][1] = (TH1F *)dir->FindObjectAny("h_pedestal1_HF");
996  HistPed[1][4][2] = (TH1F *)dir->FindObjectAny("h_pedestal2_HF");
997  HistPed[1][4][3] = (TH1F *)dir->FindObjectAny("h_pedestal3_HF");
998  HistPed[2][4][0] = (TH1F *)dir->FindObjectAny("h_pedestalw0_HF");
999  HistPed[2][4][1] = (TH1F *)dir->FindObjectAny("h_pedestalw1_HF");
1000  HistPed[2][4][2] = (TH1F *)dir->FindObjectAny("h_pedestalw2_HF");
1001  HistPed[2][4][3] = (TH1F *)dir->FindObjectAny("h_pedestalw3_HF");
1002 
1003  for (int test = 31; test <= 32; test++) { //Test: 31-Pedestals, 32-pedestal Widths,
1004  for (int sub = 1; sub <= 4; sub++) { //Subdetector: 1-HB, 2-HE, 3-HO, 4-HF
1005  // if (sub==1) cHB->Divide(2,1);
1006  if (sub == 1)
1007  cHB->Divide(2, 2);
1008  // if (sub==2) cHE->Divide(3,1);
1009  if (sub == 2)
1010  cHE->Divide(3, 3);
1011  if (sub == 3)
1012  cONE->Divide(1, 1);
1013  // if (sub==4) cHF->Divide(2,1);
1014  if (sub == 4)
1015  cHF->Divide(2, 2);
1016  // int k_min[5]={0,1,1,4,1}; // minimum depth for each subdet
1017  // int k_max[5]={0,2,3,4,2}; // maximum depth for each subdet
1018  // int k_max[5]={0,2,7,4,4}; // maximum depth for each subdet
1019  for (int k = k_min[sub]; k <= k_max[sub]; k++) { //Depths
1020  if (sub == 1)
1021  cHB->cd(k);
1022  if (sub == 2)
1023  cHE->cd(k);
1024  if (sub == 3)
1025  cONE->cd(k - 3);
1026  if (sub == 4)
1027  cHF->cd(k);
1028  Map_Ampl[test][sub][k]->Divide(Map_Ampl[test][sub][k], Map_SUB[sub][k], 1, 1, "B");
1029  gPad->SetGridy();
1030  gPad->SetGridx();
1031  gPad->SetLogz();
1032  if (sub == 1)
1033  sprintf(str, "HB, Depth%d \b", k);
1034  if (sub == 2)
1035  sprintf(str, "HE, Depth%d \b", k);
1036  if (sub == 3)
1037  sprintf(str, "HO, Depth%d \b", k);
1038  if (sub == 4)
1039  sprintf(str, "HF, Depth%d \b", k);
1040  Map_Ampl[test][sub][k]->SetTitle(str);
1041  Map_Ampl[test][sub][k]->SetXTitle("#eta \b");
1042  Map_Ampl[test][sub][k]->SetYTitle("#phi \b");
1043  Map_Ampl[test][sub][k]->SetZTitle("Rate \b");
1044  if (sub == 3)
1045  Map_Ampl[test][sub][k]->SetTitleOffset(0.8, "Z");
1046  Map_Ampl[test][sub][k]->Draw("COLZ");
1047  Map_Ampl[test][sub][k]->GetYaxis()->SetRangeUser(0, 72.);
1048  Map_Ampl[test][sub][k]->GetZaxis()->SetRangeUser(0.0001, 1.);
1049  if (sub == 1) {
1050  cHB->Modified();
1051  cHB->Update();
1052  }
1053  if (sub == 2) {
1054  cHE->Modified();
1055  cHE->Update();
1056  }
1057  if (sub == 3) {
1058  cONE->Modified();
1059  cONE->Update();
1060  }
1061  if (sub == 4) {
1062  cHF->Modified();
1063  cHF->Update();
1064  }
1065  } //end depth
1066  if (test == 31) {
1067  if (sub == 1) {
1068  cHB->Print("MapRatePedHB.png");
1069  cHB->Clear();
1070  }
1071  if (sub == 2) {
1072  cHE->Print("MapRatePedHE.png");
1073  cHE->Clear();
1074  }
1075  if (sub == 3) {
1076  cONE->Print("MapRatePedHO.png");
1077  cONE->Clear();
1078  }
1079  if (sub == 4) {
1080  cHF->Print("MapRatePedHF.png");
1081  cHF->Clear();
1082  }
1083  }
1084  if (test == 32) {
1085  if (sub == 1) {
1086  cHB->Print("MapRatePedWidthsHB.png");
1087  cHB->Clear();
1088  }
1089  if (sub == 2) {
1090  cHE->Print("MapRatePedWidthsHE.png");
1091  cHE->Clear();
1092  }
1093  if (sub == 3) {
1094  cONE->Print("MapRatePedWidthsHO.png");
1095  cONE->Clear();
1096  }
1097  if (sub == 4) {
1098  cHF->Print("MapRatePedWidthsHF.png");
1099  cHF->Clear();
1100  }
1101  }
1102 
1104 
1105  cPED->Divide(2, 2);
1106  for (int cap = 0; cap <= 3; cap++) {
1107  cPED->cd(cap + 1);
1108  gPad->SetGridy();
1109  gPad->SetGridx();
1110  gPad->SetLogy();
1111 
1112  if (sub == 1)
1113  sprintf(str, "HB, Cap%d, all depth\b", cap);
1114  if (sub == 2)
1115  sprintf(str, "HE, Cap%d, all depth\b", cap);
1116  if (sub == 3)
1117  sprintf(str, "HO, Cap%d, all depth\b", cap);
1118  if (sub == 4)
1119  sprintf(str, "HF, Cap%d, all depth\b", cap);
1120 
1121  HistPed[test - 30][sub][cap]->SetTitle(str);
1122 
1123  if (test == 31)
1124  HistPed[test - 30][sub][cap]->SetXTitle("Pedestals in each event & cell \b");
1125  if (test == 32)
1126  HistPed[test - 30][sub][cap]->SetXTitle("Pedestal Widths in each event & cell \b");
1127 
1128  HistPed[test - 30][sub][cap]->SetYTitle("Number of channel-events \b");
1129  HistPed[test - 30][sub][cap]->SetLineColor(4);
1130  HistPed[test - 30][sub][cap]->SetLineWidth(2);
1131  HistPed[test - 30][sub][cap]->SetTitleOffset(1.4, "Y");
1132  HistPed[test - 30][sub][cap]->Draw("");
1133  // HistPed[test-30][sub][cap]->GetYaxis()->SetRangeUser(1., 100.);
1134  // if (test==31) {gPad->SetLogx(); HistPed[test-30][sub][cap]->GetXaxis()->SetRangeUser(1., 10000.);}
1135  // if (test==32) HistPed[test-30][sub][cap]->GetXaxis()->SetRangeUser(0., 5.);
1136 
1137  cPED->Modified();
1138  cPED->Update();
1139  double min_x[] = {Pedest[test - 31][sub], Pedest[test - 31][sub]};
1140  double min_y[] = {0., 100000000.};
1141  TGraph *MIN = new TGraph(2, min_x, min_y);
1142  MIN->SetLineStyle(2);
1143  MIN->SetLineColor(2);
1144  MIN->SetLineWidth(2 + 100 * 100);
1145  MIN->SetFillStyle(3005);
1146  MIN->SetFillColor(2);
1147  MIN->Draw("L");
1148  }
1149  if (test == 31) {
1150  if (sub == 1) {
1151  cPED->Print("HistPedestalsHB.png");
1152  cPED->Clear();
1153  }
1154  if (sub == 2) {
1155  cPED->Print("HistPedestalsHE.png");
1156  cPED->Clear();
1157  }
1158  if (sub == 3) {
1159  cPED->Print("HistPedestalsHO.png");
1160  cPED->Clear();
1161  }
1162  if (sub == 4) {
1163  cPED->Print("HistPedestalsHF.png");
1164  cPED->Clear();
1165  }
1166  }
1167  if (test == 32) {
1168  if (sub == 1) {
1169  cPED->Print("HistPedestalWidthsHB.png");
1170  cPED->Clear();
1171  }
1172  if (sub == 2) {
1173  cPED->Print("HistPedestalWidthsHE.png");
1174  cPED->Clear();
1175  }
1176  if (sub == 3) {
1177  cPED->Print("HistPedestalWidthsHO.png");
1178  cPED->Clear();
1179  }
1180  if (sub == 4) {
1181  cPED->Print("HistPedestalWidthsHF.png");
1182  cPED->Clear();
1183  }
1184  }
1185  } // end sub
1186  } //end test 31,32
1187 
1188  //+++++++++++++++++++++++++++++++++++
1189  //Test 33 Correlation of Pedestal, pedestalWidths Vs fullAmplitude
1190  //++++++++++++++++++++++++++++++++++++
1191 
1192  cPED->Clear();
1193  Map_Ped[1][1] = (TH2F *)dir->FindObjectAny("h2_pedvsampl_HB");
1194  Map_Ped[1][2] = (TH2F *)dir->FindObjectAny("h2_pedvsampl_HE");
1195  Map_Ped[1][3] = (TH2F *)dir->FindObjectAny("h2_pedvsampl_HO");
1196  Map_Ped[1][4] = (TH2F *)dir->FindObjectAny("h2_pedvsampl_HF");
1197  Map_Ped[2][1] = (TH2F *)dir->FindObjectAny("h2_pedwvsampl_HB");
1198  Map_Ped[2][2] = (TH2F *)dir->FindObjectAny("h2_pedwvsampl_HE");
1199  Map_Ped[2][3] = (TH2F *)dir->FindObjectAny("h2_pedwvsampl_HO");
1200  Map_Ped[2][4] = (TH2F *)dir->FindObjectAny("h2_pedwvsampl_HF");
1201  for (int sub = 1; sub <= 4; sub++) { //Subdetector: 1-HB, 2-HE, 3-HO, 4-HF
1202  cPED->Divide(2, 1);
1203  for (int test = 1; test <= 2; test++) {
1204  cPED->cd(test);
1205  gPad->SetGridy();
1206  gPad->SetGridx();
1207  gPad->SetLogz();
1208  if (test == 1)
1209  Map_Ped[test][sub]->SetXTitle("Pedestal, fC \b");
1210  if (test == 2)
1211  Map_Ped[test][sub]->SetXTitle("pedestal Width, fC \b");
1212  Map_Ped[test][sub]->SetYTitle("Amplitude, fC \b");
1213  Map_Ped[test][sub]->SetZTitle("entries \b");
1214  if (test == 1)
1215  sprintf(str, "Cap0 Pedestal vs Amplitude \b");
1216  if (test == 2)
1217  sprintf(str, "Cap0 pedestalWidth vs Amplitude \b");
1218  Map_Ped[test][sub]->SetTitle(str);
1219  Map_Ped[test][sub]->Draw("COLZ");
1220  // Map_Ped[test][sub]->GetYaxis()->SetRangeUser(0, 72.);
1221  // Map_Ped[test][sub]->GetZaxis()->SetRangeUser(0.0001, 1.);
1222  cPED->Modified();
1223  cPED->Update();
1224  } // test 1,2
1225  if (sub == 1) {
1226  cPED->Print("CorrelationsMapPedestalVsfullAmplitudeHB.png");
1227  cPED->Clear();
1228  }
1229  if (sub == 2) {
1230  cPED->Print("CorrelationsMapPedestalVsfullAmplitudeHE.png");
1231  cPED->Clear();
1232  }
1233  if (sub == 3) {
1234  cPED->Print("CorrelationsMapPedestalVsfullAmplitudeHO.png");
1235  cPED->Clear();
1236  }
1237  if (sub == 4) {
1238  cPED->Print("CorrelationsMapPedestalVsfullAmplitudeHF.png");
1239  cPED->Clear();
1240  }
1241  } // end sub
1242 
1243  //+++++++++++++++++++++++++++++++++++
1244  //Test 41 Time Slices shape for good and bad channels
1245  //++++++++++++++++++++++++++++++++++++
1246 
1247  cONE->Clear();
1248  hist_GoodTSshape[1] = (TH1F *)dir->FindObjectAny("h_shape_good_channels_HB");
1249  hist_GoodTSshape[2] = (TH1F *)dir->FindObjectAny("h_shape_good_channels_HE");
1250  hist_GoodTSshape[3] = (TH1F *)dir->FindObjectAny("h_shape_good_channels_HO");
1251  hist_GoodTSshape[4] = (TH1F *)dir->FindObjectAny("h_shape_good_channels_HF");
1252 
1253  hist_GoodTSshape0[1] = (TH1F *)dir->FindObjectAny("h_shape0_good_channels_HB");
1254  hist_GoodTSshape0[2] = (TH1F *)dir->FindObjectAny("h_shape0_good_channels_HE");
1255  hist_GoodTSshape0[3] = (TH1F *)dir->FindObjectAny("h_shape0_good_channels_HO");
1256  hist_GoodTSshape0[4] = (TH1F *)dir->FindObjectAny("h_shape0_good_channels_HF");
1257 
1258  hist_BadTSshape[1] = (TH1F *)dir->FindObjectAny("h_shape_bad_channels_HB");
1259  hist_BadTSshape[2] = (TH1F *)dir->FindObjectAny("h_shape_bad_channels_HE");
1260  hist_BadTSshape[3] = (TH1F *)dir->FindObjectAny("h_shape_bad_channels_HO");
1261  hist_BadTSshape[4] = (TH1F *)dir->FindObjectAny("h_shape_bad_channels_HF");
1262 
1263  hist_BadTSshape0[1] = (TH1F *)dir->FindObjectAny("h_shape0_bad_channels_HB");
1264  hist_BadTSshape0[2] = (TH1F *)dir->FindObjectAny("h_shape0_bad_channels_HE");
1265  hist_BadTSshape0[3] = (TH1F *)dir->FindObjectAny("h_shape0_bad_channels_HO");
1266  hist_BadTSshape0[4] = (TH1F *)dir->FindObjectAny("h_shape0_bad_channels_HF");
1267 
1268  cONE->cd(1);
1269 
1270  for (int sub = 1; sub <= 4; sub++) { //Subdetector: 1-HB, 2-HE, 3-HO, 4-HF
1271 
1272  gPad->SetGridy();
1273  gPad->SetGridx();
1274  gPad->SetLogz();
1275  hist_GoodTSshape[sub]->Divide(hist_GoodTSshape[sub], hist_GoodTSshape0[sub], 1, 1, "B");
1276  hist_GoodTSshape[sub]->SetXTitle("Time slice \b");
1277  hist_GoodTSshape[sub]->SetYTitle("ADC counts \b");
1278  sprintf(str, "Mean ADC Shape \b");
1279  hist_GoodTSshape[sub]->SetTitle(str);
1280  hist_GoodTSshape[sub]->Draw("");
1281  // hist_GoodTSshape[sub]->GetYaxis()->SetRangeUser(0, 72.);
1282  // hist_GoodTSshape[sub]->GetZaxis()->SetRangeUser(0.0001, 1.);
1283  cONE->Modified();
1284  cONE->Update();
1285  if (sub == 1) {
1286  cONE->Print("HistGoodTSshapesHB.png");
1287  cONE->Clear();
1288  }
1289  if (sub == 2) {
1290  cONE->Print("HistGoodTSshapesHE.png");
1291  cONE->Clear();
1292  }
1293  if (sub == 3) {
1294  cONE->Print("HistGoodTSshapesHO.png");
1295  cONE->Clear();
1296  }
1297  if (sub == 4) {
1298  cONE->Print("HistGoodTSshapesHF.png");
1299  cONE->Clear();
1300  }
1301  } // end sub
1302 
1303  for (int sub = 1; sub <= 4; sub++) { //Subdetector: 1-HB, 2-HE, 3-HO, 4-HF
1304 
1305  gPad->SetGridy();
1306  gPad->SetGridx();
1307  gPad->SetLogz();
1308  hist_BadTSshape[sub]->Divide(hist_BadTSshape[sub], hist_BadTSshape0[sub], 1, 1, "B");
1309  hist_BadTSshape[sub]->SetXTitle("Time slice \b");
1310  hist_BadTSshape[sub]->SetYTitle("ADC counts \b");
1311  sprintf(str, "Mean ADC Shape \b");
1312  hist_BadTSshape[sub]->SetTitle(str);
1313  hist_BadTSshape[sub]->Draw("");
1314  // hist_BadTSshape[sub]->GetYaxis()->SetRangeUser(0, 72.);
1315  // hist_BadTSshape[sub]->GetZaxis()->SetRangeUser(0.0001, 1.);
1316  cONE->Modified();
1317  cONE->Update();
1318  if (sub == 1) {
1319  cONE->Print("HistBadTSshapesHB.png");
1320  cONE->Clear();
1321  }
1322  if (sub == 2) {
1323  cONE->Print("HistBadTSshapesHE.png");
1324  cONE->Clear();
1325  }
1326  if (sub == 3) {
1327  cONE->Print("HistBadTSshapesHO.png");
1328  cONE->Clear();
1329  }
1330  if (sub == 4) {
1331  cONE->Print("HistBadTSshapesHF.png");
1332  cONE->Clear();
1333  }
1334  } // end sub
1335 
1336  //+++++++++++++++++++++++++++++
1337  //Entries in different TSs:
1338  //+++++++++++++++++++++++++++++
1339  Map_SUBTS[1][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1TS2_HB");
1340  Map_SUBTS[1][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2TS2_HB");
1341  Map_SUBTS[1][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3TS2_HB");
1342  Map_SUBTS[1][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4TS2_HB");
1343 
1344  Map_SUBTS[2][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1TS2_HE");
1345  Map_SUBTS[2][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2TS2_HE");
1346  Map_SUBTS[2][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3TS2_HE");
1347  Map_SUBTS[2][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4TS2_HE");
1348  Map_SUBTS[2][5] = (TH2F *)dir->FindObjectAny("h_mapDepth5TS2_HE");
1349  Map_SUBTS[2][6] = (TH2F *)dir->FindObjectAny("h_mapDepth6TS2_HE");
1350  Map_SUBTS[2][7] = (TH2F *)dir->FindObjectAny("h_mapDepth7TS2_HE");
1351 
1352  Map_SUBTS[3][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4TS012_HO");
1353 
1354  Map_SUBTS[4][1] = (TH2F *)dir->FindObjectAny("h_mapDepth1TS1_HF");
1355  Map_SUBTS[4][2] = (TH2F *)dir->FindObjectAny("h_mapDepth2TS1_HF");
1356  Map_SUBTS[4][3] = (TH2F *)dir->FindObjectAny("h_mapDepth3TS1_HF");
1357  Map_SUBTS[4][4] = (TH2F *)dir->FindObjectAny("h_mapDepth4TS1_HF");
1358 
1359  for (int sub = 1; sub <= 4; sub++) { //Subdetector: 1-HB, 2-HE, 3-HF, 4-HO
1360  // if (sub==1) cHB->Divide(2,1);
1361  if (sub == 1)
1362  cHB->Divide(2, 2);
1363  // if (sub==2) cHE->Divide(3,1);
1364  if (sub == 2)
1365  cHE->Divide(3, 3);
1366  if (sub == 3)
1367  cONE->Divide(1, 1);
1368  // if (sub==4) cHF->Divide(2,1);
1369  if (sub == 4)
1370  cHF->Divide(2, 2);
1371  // int k_min[5]={0,1,1,4,1}; // minimum depth for each subdet
1372  // int k_max[5]={0,2,3,4,2}; // maximum depth for each subdet
1373  // int k_max[5]={0,2,7,4,4}; // maximum depth for each subdet
1374  for (int k = k_min[sub]; k <= k_max[sub]; k++) { //Depth
1375  if (sub == 1)
1376  cHB->cd(k);
1377  if (sub == 2)
1378  cHE->cd(k);
1379  if (sub == 3)
1380  cONE->cd(k - 3);
1381  if (sub == 4)
1382  cHF->cd(k);
1383  Map_SUBTS[sub][k]->Divide(Map_SUBTS[sub][k], Map_SUB[sub][k], 1, 1, "B");
1384  gPad->SetGridy();
1385  gPad->SetGridx();
1386  gPad->SetLogz();
1387  if (sub == 1)
1388  sprintf(str, "HB, Depth%d \b", k);
1389  if (sub == 2)
1390  sprintf(str, "HE, Depth%d \b", k);
1391  if (sub == 3)
1392  sprintf(str, "HO, Depth%d \b", k);
1393  if (sub == 4)
1394  sprintf(str, "HF, Depth%d \b", k);
1395  Map_SUBTS[sub][k]->SetTitle(str);
1396  Map_SUBTS[sub][k]->SetXTitle("#eta \b");
1397  Map_SUBTS[sub][k]->SetYTitle("#phi \b");
1398  Map_SUBTS[sub][k]->SetZTitle("Number of events \b");
1399  if (sub == 3)
1400  Map_SUBTS[sub][k]->SetTitleOffset(0.8, "Z");
1401  Map_SUBTS[sub][k]->Draw("COLZ");
1402  Map_SUBTS[sub][k]->GetYaxis()->SetRangeUser(0, 72.);
1403  // Map_SUBTS[sub][k]->GetZaxis()->SetRangeUser(0.0001, 1.);
1404  if (sub == 1) {
1405  cHB->Modified();
1406  cHB->Update();
1407  }
1408  if (sub == 2) {
1409  cHE->Modified();
1410  cHE->Update();
1411  }
1412  if (sub == 3) {
1413  cONE->Modified();
1414  cONE->Update();
1415  }
1416  if (sub == 4) {
1417  cHF->Modified();
1418  cHF->Update();
1419  }
1420  } //end depth
1421 
1422  if (sub == 1) {
1423  cHB->Print("Hist_mapDepthAllTS2_HB.png");
1424  cHB->Clear();
1425  }
1426  if (sub == 2) {
1427  cHE->Print("Hist_mapDepthAllTS2_HE.png");
1428  cHE->Clear();
1429  }
1430  if (sub == 3) {
1431  cONE->Print("Hist_mapDepthAllTS012_HO.png");
1432  cONE->Clear();
1433  }
1434  if (sub == 4) {
1435  cHF->Print("Hist_mapDepthAllTS1_HF.png");
1436  cHF->Clear();
1437  }
1438  } // end sub
1439 
1440  //======================================================================
1441 
1442  //AZ2023: std::cout << " We are here to print general 2D MAP " << std::endl;
1443 
1444  //======================================================================
1445 
1446  //======================================================================
1448  // i - Eta; j - Phi
1449  //Subdetector: 1-HB, 2-HE, 3-HF, 4-HO
1450  // int k_min[5]={0,1,1,4,1}; // minimum depth for each subdet
1451  // int k_max[5]={0,2,3,4,2}; // maximum depth for each subdet before upgrade
1452  // int k_max[5]={0,2,7,4,4}; // maximum depth for each subdet
1453  //k-Depth
1454 
1455  // TH2F *Map_ALL = new TH2F("Map_All", "Map_all", 82, -41, 40, 72, 0, 71);
1456  //AZ2023:
1457  //AZ2023 TH2F *Map_ALL = new TH2F("Map_All", "Map_all", 82, -41, 41, 72, 0, 72);
1458 
1459  /*
1460  int nx = Map_ALL->GetXaxis()->GetNbins();
1461  int ny = Map_ALL->GetYaxis()->GetNbins();
1462  cout << " nx= " << nx << " ny= " << ny << endl;
1463 */
1464  // int NBad = 0;
1465  // int NWarn = 0;
1466  // int NCalib = 0;
1467  // int NPed = 0;
1468  // // int Eta[3][10000]={0};
1469  // int Eta[4][10000] = {0};
1470  // int Phi[4][10000] = {0};
1471  // int Sub[4][10000] = {0};
1472  // int Depth[4][10000] = {0};
1473  // string Comment[4][10000] = {""};
1474  // string Text[33] = {"", "Cm", "Am", "Wm", "Rm", "TNm", "TXm", "", "", "", "", "Cc", "Ac", "Wc", "Rc", "TNc", "TXc",
1475  // "", "", "", "", "GS", "", "", "", "", "", "", "", "", "", "Pm", "pWm"};
1476  // int flag_W = 0;
1477  // int flag_B = 0;
1478  // int flag_P = 0;
1479  //AZ2023:
1480  /*
1481  int fffffflag = 0;
1482  std::cout << " Map_ALL SUBGOOD update " << std::endl;
1483  for (int sub = 1; sub <= 4; sub++) {
1484  for (int k = k_min[sub]; k <= k_max[sub]; k++) {
1485  for (int i = 1; i <= nx; i++) {
1486  for (int j = 1; j <= ny; j++) {
1487  if (Map_SUB[sub][k]->GetBinContent(i, j) != 0) {
1488  //AZ2023 Map_SUBGOOD[sub][k]->SetBinContent(i, j, 0.5);
1489  //AZ2023 Map_ALL->SetBinContent(i, j, 0.5);
1490  }
1491  }
1492  }
1493  }
1494  }
1495 */
1496  //AZ2023:
1497  /*
1498  std::cout << " Map_ALL SUBGOOD filling............... " << std::endl;
1499  for (int sub = 1; sub <= 4; sub++) {
1500  for (int k = k_min[sub]; k <= k_max[sub]; k++) {
1501  for (int i = 1; i <= nx; i++) {
1502  for (int j = 1; j <= ny; j++) {
1503  // flag_W = 0;
1504  // flag_B = 0;
1505  // flag_P = 0;
1506  // CapID(Test=1; ADC amplitude Am(Test= 2); Width for Wm(Test=3); Ratio cut for Rm(Test=4); TS mean for TNm(test=5); TS max for TXm(Test=6);
1507  for (int test = 3; test <= 6; test++) {
1508  // cout<<" test= "<<test<<" sbd= "<<sub<<" depth= "<<k<<" eta= "<<i<<" , phi= "<<j<<endl;
1509  // cout<<" initial content Map_Ampl[test][sub][k]->GetBinContent(i, j)= "<< Map_Ampl[test][sub][k]->GetBinContent(i, j) <<endl;
1510 
1511  //Bad
1512  //Rate 0.1 for displaying on whole detector map and subdetector map
1513  if (Map_Ampl[test][sub][k]->GetBinContent(i, j) > 0.1) {
1514  //AZ2023 Map_ALL->SetBinContent(i, j, 1.);
1515  //AZ2023 Map_SUBGOOD[sub][k]->SetBinContent(i, j, 1.);
1516  fffffflag = 1;
1517  }
1518 
1519  if ((Map_Ampl[test][sub][k]->GetBinContent(i, j) != 0.) &&
1520  (Map_Ampl[test][sub][k]->GetBinContent(i, j) < 0.001)) {
1521  if (Map_SUBGOOD[sub][k]->GetBinContent(i, j) != 1.)
1522  //AZ2023 Map_SUBGOOD[sub][k]->SetBinContent(i, j, 0.75);
1523  //AZ2023 if (Map_ALL->GetBinContent(i, j) != 1.)Map_ALL->SetBinContent(i, j, 0.75);
1524  fffffflag = 2;
1525  }
1527 
1528  // if(fffffflag != 0) cout<<"Map_Ampl["<<test<<"]["<<sub<<"]["<<k<<"]->GetBinContent("<<i<<","<<j<<")= "<<Map_Ampl[test][sub][k]->GetBinContent(i,j) << "fffffflag = "<< fffffflag <<endl;
1529 
1530  } //end test
1531 
1532  // std::cout << " RUN3 2022 MAPS_SUB: Pedestals......"<< std::endl;
1533  //Pedestals
1534  for (int test = 31; test <= 32; test++) {
1535  // cout<<"Pedestals test= "<<test<<" sbd= "<<sub<<" depth= "<<k<<" eta= "<<i<<" , phi= "<<j<<endl;
1536  if (Map_Ampl[test][sub][k]->GetBinContent(i, j) > 0.9) {
1537  //AZ2023 if (Map_SUBGOOD[sub][k]->GetBinContent(i, j) != 1.0) Map_SUBGOOD[sub][k]->SetBinContent(i, j, 0.15);
1538  //AZ2023 if (Map_ALL->GetBinContent(i, j) != 1.) Map_ALL->SetBinContent(i, j, 0.15);
1539  }
1540  // cout<<"Pedestals Map_Ampl["<<test<<"]["<<sub<<"]["<<k<<"]->GetBinContent("<<i<<","<<j<<")= "<<Map_Ampl[test][sub][k]->GetBinContent(i,j)<<endl;
1541  } //end test
1542  }
1543  }
1544  }
1545  }
1546 */
1547  //AZ2023:
1548  /*
1549  std::cout << " RUN3: 2022 Plots with MAPS_SUB: start ..............................." << std::endl;
1550  // subdet maps
1551  for (int sub = 1; sub <= 4; sub++) { //Subdetector: 1-HB, 2-HE, 3-HF, 4-HO
1552 
1553  std::cout << " RUN3: 2022 MAPS_SUB= " << sub << std::endl;
1554  // if (sub==1) cHB->Divide(2,1);
1555  if (sub == 1)
1556  cHB->Divide(2, 2);
1557  // if (sub==2) cHE->Divide(3,1);
1558  if (sub == 2)
1559  cHE->Divide(3, 3);
1560  if (sub == 3)
1561  cONE->Divide(1, 1);
1562  // if (sub==4) cHB->Divide(2,1);
1563  if (sub == 4)
1564  cHF->Divide(2, 2);
1565  // int k_min[5]={0,1,1,4,1}; // minimum depth for each subdet
1566  // int k_max[5]={0,2,3,4,2}; // maximum depth for each subdet
1567  // int k_max[5]={0,2,7,4,4}; // maximum depth for each subdet
1568  //k = Depth
1569  for (int k = k_min[sub]; k <= k_max[sub]; k++) {
1570  if (sub == 1)
1571  cHB->cd(k);
1572  if (sub == 2)
1573  cHE->cd(k);
1574  if (sub == 3)
1575  cONE->cd(k - 3);
1576  if (sub == 4)
1577  cHF->cd(k);
1578  gPad->SetGridy();
1579  gPad->SetGridx();
1580  gPad->SetLogz();
1581  // gStyle->SetTitleOffset(0.5, "Y");
1582  if (sub == 1)
1583  sprintf(str, "HB, Depth%d \b", k);
1584  if (sub == 2)
1585  sprintf(str, "HE, Depth%d \b", k);
1586  if (sub == 3)
1587  sprintf(str, "HO, Depth%d \b", k);
1588  if (sub == 4)
1589  sprintf(str, "HF, Depth%d \b", k);
1590  Map_SUBGOOD[sub][k]->SetTitle(str);
1591  Map_SUBGOOD[sub][k]->SetXTitle("#eta \b");
1592  Map_SUBGOOD[sub][k]->SetYTitle("#phi \b");
1593  Map_SUBGOOD[sub][k]->Draw("COLZ");
1594  Map_SUBGOOD[sub][k]->GetYaxis()->SetRangeUser(0, 72.);
1595  Map_SUBGOOD[sub][k]->GetZaxis()->SetRangeUser(0., 1.);
1596 
1597  if (sub == 1) {
1598  cHB->Modified();
1599  cHB->Update();
1600  }
1601  if (sub == 2) {
1602  cHE->Modified();
1603  cHE->Update();
1604  }
1605  if (sub == 3) {
1606  cONE->Modified();
1607  cONE->Update();
1608  }
1609  if (sub == 4) {
1610  cHF->Modified();
1611  cHF->Update();
1612  }
1613  } //end depth
1614  if (sub == 1) {
1615  cHB->Print("MAPHB.png");
1616  cHB->Clear();
1617  }
1618  if (sub == 2) {
1619  cHE->Print("MAPHE.png");
1620  cHE->Clear();
1621  }
1622  if (sub == 3) {
1623  cONE->Print("MAPHO.png");
1624  cONE->Clear();
1625  }
1626  if (sub == 4) {
1627  cHF->Print("MAPHF.png");
1628  cHF->Clear();
1629  }
1630  } // end sub
1631 */
1633 
1634  TCanvas *cmain1 = new TCanvas("cmain1", "cmain1", 200, 10, 1400, 1800);
1635  cmain1->Divide(2, 2);
1636 
1637  cmain1->cd(1);
1638  TH1F *JDBEYESJ0 = (TH1F *)dir->FindObjectAny("h_totalAmplitudeHBperEvent");
1639  JDBEYESJ0->SetStats(0);
1640  JDBEYESJ0->SetMarkerStyle(20);
1641  JDBEYESJ0->SetMarkerSize(0.8);
1642  JDBEYESJ0->GetYaxis()->SetLabelSize(0.04);
1643  JDBEYESJ0->SetXTitle("iEvent \b");
1644  JDBEYESJ0->SetYTitle("totalAmplitude perEvent \b");
1645  JDBEYESJ0->SetTitle("HB \b");
1646  JDBEYESJ0->SetMarkerColor(2);
1647  JDBEYESJ0->SetLineColor(1);
1648  JDBEYESJ0->SetMinimum(0.8);
1649  JDBEYESJ0->Draw("HIST same P0");
1650  //JDBEYESJ0->Clear();
1651 
1652  cmain1->cd(2);
1653  TH1F *JDBEYESJ1 = (TH1F *)dir->FindObjectAny("h_totalAmplitudeHEperEvent");
1654  JDBEYESJ1->SetStats(0);
1655  JDBEYESJ1->SetMarkerStyle(20);
1656  JDBEYESJ1->SetMarkerSize(0.8);
1657  JDBEYESJ1->GetYaxis()->SetLabelSize(0.04);
1658  JDBEYESJ1->SetXTitle("iEvent \b");
1659  JDBEYESJ1->SetYTitle("totalAmplitude perEvent \b");
1660  JDBEYESJ1->SetTitle("HE \b");
1661  JDBEYESJ1->SetMarkerColor(2);
1662  JDBEYESJ1->SetLineColor(1);
1663  JDBEYESJ1->SetMinimum(0.8);
1664  JDBEYESJ1->Draw("HIST same P0");
1665  //JDBEYESJ1->Clear();
1666 
1667  cmain1->cd(3);
1668  TH1F *JDBEYESJ2 = (TH1F *)dir->FindObjectAny("h_totalAmplitudeHFperEvent");
1669  JDBEYESJ2->SetStats(0);
1670  JDBEYESJ2->SetMarkerStyle(20);
1671  JDBEYESJ2->SetMarkerSize(0.8);
1672  JDBEYESJ2->GetYaxis()->SetLabelSize(0.04);
1673  JDBEYESJ2->SetXTitle("iEvent \b");
1674  JDBEYESJ2->SetYTitle("totalAmplitude perEvent \b");
1675  JDBEYESJ2->SetTitle("HF \b");
1676  JDBEYESJ2->SetMarkerColor(2);
1677  JDBEYESJ2->SetLineColor(1);
1678  JDBEYESJ2->SetMinimum(0.8);
1679  JDBEYESJ2->Draw("HIST same P0");
1680  //JDBEYESJ2->Clear();
1681 
1682  cmain1->cd(4);
1683  TH1F *JDBEYESJ3 = (TH1F *)dir->FindObjectAny("h_totalAmplitudeHOperEvent");
1684  JDBEYESJ3->SetStats(0);
1685  JDBEYESJ3->SetMarkerStyle(20);
1686  JDBEYESJ3->SetMarkerSize(0.8);
1687  JDBEYESJ3->GetYaxis()->SetLabelSize(0.04);
1688  JDBEYESJ3->SetXTitle("iEvent \b");
1689  JDBEYESJ3->SetYTitle("totalAmplitude perEvent \b");
1690  JDBEYESJ3->SetTitle("HO \b");
1691  JDBEYESJ3->SetMarkerColor(2);
1692  JDBEYESJ3->SetLineColor(1);
1693  JDBEYESJ3->SetMinimum(0.8);
1694  JDBEYESJ3->Draw("HIST same P0");
1695  //JDBEYESJ3->Clear();
1696  cmain1->Modified();
1697  cmain1->Update();
1698  cmain1->Print("EVENTDEPENDENCE.png");
1699  cmain1->Clear();
1700  // std::cout << " EVENTDEPENDENCE " << std::endl;
1702  //AZ2023:
1703  /*
1704  // ALL SubDet
1705  gStyle->SetOptTitle(0);
1706  TCanvas *cmain = new TCanvas("cmain", "MAP", 1000, 1000);
1707  gPad->SetGridy();
1708  gPad->SetGridx();
1709  gPad->SetLogz();
1710  Map_ALL->SetTitleOffset(1.3, "Y");
1711  Map_ALL->SetXTitle("#eta \b");
1712  Map_ALL->SetYTitle("#phi \b");
1713  Map_ALL->Draw("COLZ");
1714  Map_ALL->GetYaxis()->SetRangeUser(0, 72.);
1715  Map_ALL->GetZaxis()->SetRangeUser(0, 1.);
1716  cmain->Modified();
1717  cmain->Update();
1718  cmain->Print("MAP.png");
1719  cmain->Clear();
1720  std::cout << " MAP_ALL " << std::endl;
1721 */
1722 
1724  // 7. Correlation of Charge(=Amplitude) vs timing, fc (=TSN * 25 ns)
1725  // three plots -----------------------------------------------------------------
1726 
1727  TCanvas *corravstsn = new TCanvas("corravstsn", "corravstsn", 200, 10, 1400, 1800);
1728  // three plots for HB:
1729  corravstsn->Divide(2, 2);
1730  corravstsn->cd(1);
1731  TH2F *two11 = (TH2F *)dir->FindObjectAny("h2_TSnVsAyear2023_HB");
1732  gPad->SetGridy();
1733  gPad->SetGridx();
1734  two11->SetMarkerStyle(20);
1735  two11->SetMarkerSize(0.4);
1736  two11->SetYTitle("timing HB \b");
1737  two11->SetXTitle("Q,fc HB\b");
1738  two11->SetMarkerColor(1);
1739  two11->SetLineColor(1);
1740  two11->Draw("BOX");
1741  corravstsn->cd(2);
1742  TH1F *TSNvsQ_HB = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HB");
1743  gPad->SetGridy();
1744  gPad->SetGridx();
1745  gPad->SetLogy();
1746  TSNvsQ_HB->SetMarkerStyle(20);
1747  TSNvsQ_HB->SetMarkerSize(0.6);
1748  TSNvsQ_HB->GetYaxis()->SetLabelSize(0.04);
1749  TSNvsQ_HB->SetXTitle("Q,fc HB \b");
1750  TSNvsQ_HB->SetYTitle("iev*ieta*iphi*idepth \b");
1751  TSNvsQ_HB->SetMarkerColor(4);
1752  TSNvsQ_HB->SetLineColor(0);
1753  TSNvsQ_HB->SetMinimum(0.8);
1754  TSNvsQ_HB->Draw("E");
1755  corravstsn->cd(3);
1756  TH1F *twod1_HB = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear2023_HB");
1757  TH1F *twod0_HB = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HB");
1758  // twod1_HB->Sumw2();
1759  // twod0_HB->Sumw2();
1760  TH1F *Ceff_HB = (TH1F *)twod1_HB->Clone("Ceff_HB");
1761  for (int x = 1; x <= twod1_HB->GetXaxis()->GetNbins(); x++) {
1762  twod1_HB->SetBinError(float(x), 0.001);
1763  } //end x
1764  Ceff_HB->Divide(twod1_HB, twod0_HB, 1, 1, "B");
1765  gPad->SetGridy();
1766  gPad->SetGridx();
1767  Ceff_HB->SetMarkerStyle(20);
1768  Ceff_HB->SetMarkerSize(0.4);
1769  Ceff_HB->SetXTitle("Q,fc \b");
1770  Ceff_HB->SetYTitle("<timing>HB \b");
1771  Ceff_HB->SetMarkerColor(2);
1772  Ceff_HB->SetLineColor(2);
1773  Ceff_HB->SetMaximum(140.);
1774  Ceff_HB->SetMinimum(30.);
1775  Ceff_HB->Draw("E");
1776  corravstsn->Modified();
1777  corravstsn->Update();
1778  corravstsn->Print("corravstsnPLOTSHB.png");
1779  corravstsn->Clear();
1780  // std::cout << " corravstsnPLOTSHB.png created " << std::endl;
1781  // three plots for HE:
1782  corravstsn->Divide(2, 2);
1783  corravstsn->cd(1);
1784  TH2F *twoHE = (TH2F *)dir->FindObjectAny("h2_TSnVsAyear2023_HE");
1785  gPad->SetGridy();
1786  gPad->SetGridx();
1787  twoHE->SetMarkerStyle(20);
1788  twoHE->SetMarkerSize(0.4);
1789  twoHE->SetYTitle("timing HE \b");
1790  twoHE->SetXTitle("Q,fc HE\b");
1791  twoHE->SetMarkerColor(1);
1792  twoHE->SetLineColor(1);
1793  twoHE->Draw("BOX");
1794  corravstsn->cd(2);
1795  TH1F *TSNvsQ_HE = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HE");
1796  gPad->SetGridy();
1797  gPad->SetGridx();
1798  gPad->SetLogy();
1799  TSNvsQ_HE->SetMarkerStyle(20);
1800  TSNvsQ_HE->SetMarkerSize(0.6);
1801  TSNvsQ_HE->GetYaxis()->SetLabelSize(0.04);
1802  TSNvsQ_HE->SetXTitle("Q,fc HE \b");
1803  TSNvsQ_HE->SetYTitle("iev*ieta*iphi*idepth \b");
1804  TSNvsQ_HE->SetMarkerColor(4);
1805  TSNvsQ_HE->SetLineColor(0);
1806  TSNvsQ_HE->SetMinimum(0.8);
1807  TSNvsQ_HE->Draw("E");
1808  corravstsn->cd(3);
1809  TH1F *twod1_HE = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear2023_HE");
1810  TH1F *twod0_HE = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HE");
1811  // twod1_HE->Sumw2();
1812  // twod0_HE->Sumw2();
1813  TH1F *Ceff_HE = (TH1F *)twod1_HE->Clone("Ceff_HE");
1814  for (int x = 1; x <= twod1_HE->GetXaxis()->GetNbins(); x++) {
1815  twod1_HE->SetBinError(float(x), 0.001);
1816  } //end x
1817  Ceff_HE->Divide(twod1_HE, twod0_HE, 1, 1, "B");
1818  gPad->SetGridy();
1819  gPad->SetGridx();
1820  Ceff_HE->SetMarkerStyle(20);
1821  Ceff_HE->SetMarkerSize(0.4);
1822  Ceff_HE->SetXTitle("Q,fc \b");
1823  Ceff_HE->SetYTitle("<timing>HE \b");
1824  Ceff_HE->SetMarkerColor(2);
1825  Ceff_HE->SetLineColor(2);
1826  Ceff_HE->SetMaximum(150.);
1827  Ceff_HE->SetMinimum(25.);
1828  Ceff_HE->Draw("E");
1829  corravstsn->Modified();
1830  corravstsn->Update();
1831  corravstsn->Print("corravstsnPLOTSHE.png");
1832  corravstsn->Clear();
1833  // std::cout << " corravstsnPLOTSHE.png created " << std::endl;
1834  // three plots for HF:
1835  corravstsn->Divide(2, 2);
1836  corravstsn->cd(1);
1837  TH2F *twoHF = (TH2F *)dir->FindObjectAny("h2_TSnVsAyear2023_HF");
1838  gPad->SetGridy();
1839  gPad->SetGridx();
1840  twoHF->SetMarkerStyle(20);
1841  twoHF->SetMarkerSize(0.4);
1842  twoHF->SetYTitle("timing HF \b");
1843  twoHF->SetXTitle("Q,fc HF\b");
1844  twoHF->SetMarkerColor(1);
1845  twoHF->SetLineColor(1);
1846  twoHF->Draw("BOX");
1847  corravstsn->cd(2);
1848  TH1F *TSNvsQ_HF = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HF");
1849  gPad->SetGridy();
1850  gPad->SetGridx();
1851  gPad->SetLogy();
1852  TSNvsQ_HF->SetMarkerStyle(20);
1853  TSNvsQ_HF->SetMarkerSize(0.6);
1854  TSNvsQ_HF->GetYaxis()->SetLabelSize(0.04);
1855  TSNvsQ_HF->SetXTitle("Q,fc HF \b");
1856  TSNvsQ_HF->SetYTitle("iev*ieta*iphi*idepth \b");
1857  TSNvsQ_HF->SetMarkerColor(4);
1858  TSNvsQ_HF->SetLineColor(0);
1859  TSNvsQ_HF->SetMinimum(0.8);
1860  TSNvsQ_HF->Draw("E");
1861  corravstsn->cd(3);
1862  TH1F *twod1_HF = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear2023_HF");
1863  TH1F *twod0_HF = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HF");
1864  // twod1_HF->Sumw2();
1865  // twod0_HF->Sumw2();
1866  TH1F *Ceff_HF = (TH1F *)twod1_HF->Clone("Ceff_HF");
1867  for (int x = 1; x <= twod1_HF->GetXaxis()->GetNbins(); x++) {
1868  twod1_HF->SetBinError(float(x), 0.001);
1869  } //end x
1870  Ceff_HF->Divide(twod1_HF, twod0_HF, 1, 1, "B");
1871  gPad->SetGridy();
1872  gPad->SetGridx();
1873  Ceff_HF->SetMarkerStyle(20);
1874  Ceff_HF->SetMarkerSize(0.4);
1875  Ceff_HF->SetXTitle("Q,fc \b");
1876  Ceff_HF->SetYTitle("<timing>HF \b");
1877  Ceff_HF->SetMarkerColor(2);
1878  Ceff_HF->SetLineColor(2);
1879  Ceff_HF->SetMaximum(50.);
1880  Ceff_HF->SetMinimum(0.);
1881  Ceff_HF->Draw("E");
1882  corravstsn->Modified();
1883  corravstsn->Update();
1884  corravstsn->Print("corravstsnPLOTSHF.png");
1885  corravstsn->Clear();
1886  // std::cout << " corravstsnPLOTSHF.png created " << std::endl;
1887  // three plots for HO:
1888  corravstsn->Divide(2, 2);
1889  corravstsn->cd(1);
1890  TH2F *twoHO = (TH2F *)dir->FindObjectAny("h2_TSnVsAyear2023_HO");
1891  gPad->SetGridy();
1892  gPad->SetGridx();
1893  twoHO->SetMarkerStyle(20);
1894  twoHO->SetMarkerSize(0.4);
1895  twoHO->SetYTitle("timing HO \b");
1896  twoHO->SetXTitle("Q,fc HO\b");
1897  twoHO->SetMarkerColor(1);
1898  twoHO->SetLineColor(1);
1899  twoHO->Draw("BOX");
1900  corravstsn->cd(2);
1901  TH1F *TSNvsQ_HO = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HO");
1902  gPad->SetGridy();
1903  gPad->SetGridx();
1904  gPad->SetLogy();
1905  TSNvsQ_HO->SetMarkerStyle(20);
1906  TSNvsQ_HO->SetMarkerSize(0.6);
1907  TSNvsQ_HO->GetYaxis()->SetLabelSize(0.04);
1908  TSNvsQ_HO->SetXTitle("Q,fc HO \b");
1909  TSNvsQ_HO->SetYTitle("iev*ieta*iphi*idepth \b");
1910  TSNvsQ_HO->SetMarkerColor(4);
1911  TSNvsQ_HO->SetLineColor(0);
1912  TSNvsQ_HO->SetMinimum(0.8);
1913  TSNvsQ_HO->Draw("E");
1914  corravstsn->cd(3);
1915  TH1F *twod1_HO = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear2023_HO");
1916  TH1F *twod0_HO = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HO");
1917  // twod1_HO->Sumw2();
1918  // twod0_HO->Sumw2();
1919  TH1F *Ceff_HO = (TH1F *)twod1_HO->Clone("Ceff_HO");
1920  for (int x = 1; x <= twod1_HO->GetXaxis()->GetNbins(); x++) {
1921  twod1_HO->SetBinError(float(x), 0.001);
1922  } //end x
1923  Ceff_HO->Divide(twod1_HO, twod0_HO, 1, 1, "B");
1924  gPad->SetGridy();
1925  gPad->SetGridx();
1926  Ceff_HO->SetMarkerStyle(20);
1927  Ceff_HO->SetMarkerSize(0.4);
1928  Ceff_HO->SetXTitle("Q,fc \b");
1929  Ceff_HO->SetYTitle("<timing>HO \b");
1930  Ceff_HO->SetMarkerColor(2);
1931  Ceff_HO->SetLineColor(2);
1932  Ceff_HO->SetMaximum(150.);
1933  Ceff_HO->SetMinimum(70.);
1934  Ceff_HO->Draw("E");
1935  corravstsn->Modified();
1936  corravstsn->Update();
1937  corravstsn->Print("corravstsnPLOTSHO.png");
1938  corravstsn->Clear();
1939  // std::cout << " corravstsnPLOTSHO.png created " << std::endl;
1940 
1941  // 2D plots (from 1 to 7) <TSn> -----------------------------------------------------------------
1942  TCanvas *cHBnew = new TCanvas("cHBnew", "cHBnew", 0, 10, 1400, 1800);
1943  // 4 plots for HB:
1944  cHBnew->Clear();
1945  cHBnew->Divide(2, 2);
1946  cHBnew->cd(1);
1947  TH2F *dva1_HBDepth1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA_HB");
1948  TH2F *dva0_HBDepth1 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HB");
1949  // dva1_HBDepth1->Sumw2();
1950  // dva0_HBDepth1->Sumw2();
1951  TH2F *Seff_HBDepth1 = (TH2F *)dva1_HBDepth1->Clone("Seff_HBDepth1");
1952  Seff_HBDepth1->Divide(dva1_HBDepth1, dva0_HBDepth1, 25., 1., "B");
1953  gPad->SetGridy();
1954  gPad->SetGridx();
1955  Seff_HBDepth1->SetMarkerStyle(20);
1956  Seff_HBDepth1->SetMarkerSize(0.4);
1957  Seff_HBDepth1->SetXTitle("#eta \b");
1958  Seff_HBDepth1->SetYTitle("#phi \b");
1959  Seff_HBDepth1->SetZTitle("<timing> HB Depth1 \b");
1960  Seff_HBDepth1->SetMarkerColor(2);
1961  Seff_HBDepth1->SetLineColor(2);
1962  Seff_HBDepth1->SetMaximum(100.);
1963  Seff_HBDepth1->SetMinimum(80.);
1964  Seff_HBDepth1->Draw("COLZ");
1965  cHBnew->cd(2);
1966  TH2F *dva1_HBDepth2 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA_HB");
1967  TH2F *dva0_HBDepth2 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HB");
1968  TH2F *Seff_HBDepth2 = (TH2F *)dva1_HBDepth2->Clone("Seff_HBDepth2");
1969  Seff_HBDepth2->Divide(dva1_HBDepth2, dva0_HBDepth2, 25., 1., "B");
1970  gPad->SetGridy();
1971  gPad->SetGridx();
1972  Seff_HBDepth2->SetMarkerStyle(20);
1973  Seff_HBDepth2->SetMarkerSize(0.4);
1974  Seff_HBDepth2->SetXTitle("#eta \b");
1975  Seff_HBDepth2->SetYTitle("#phi \b");
1976  Seff_HBDepth2->SetZTitle("<timing> HB Depth2 \b");
1977  Seff_HBDepth2->SetMarkerColor(2);
1978  Seff_HBDepth2->SetLineColor(2);
1979  Seff_HBDepth2->SetMaximum(100.);
1980  Seff_HBDepth2->SetMinimum(80.);
1981  Seff_HBDepth2->Draw("COLZ");
1982  cHBnew->cd(3);
1983  TH2F *dva1_HBDepth3 = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmeanA_HB");
1984  TH2F *dva0_HBDepth3 = (TH2F *)dir->FindObjectAny("h_mapDepth3_HB");
1985  TH2F *Seff_HBDepth3 = (TH2F *)dva1_HBDepth3->Clone("Seff_HBDepth3");
1986  Seff_HBDepth3->Divide(dva1_HBDepth3, dva0_HBDepth3, 25., 1., "B");
1987  gPad->SetGridy();
1988  gPad->SetGridx();
1989  Seff_HBDepth3->SetMarkerStyle(20);
1990  Seff_HBDepth3->SetMarkerSize(0.4);
1991  Seff_HBDepth3->SetXTitle("#eta \b");
1992  Seff_HBDepth3->SetYTitle("#phi \b");
1993  Seff_HBDepth3->SetZTitle("<timing> HB Depth3 \b");
1994  Seff_HBDepth3->SetMarkerColor(2);
1995  Seff_HBDepth3->SetLineColor(2);
1996  Seff_HBDepth3->SetMaximum(100.);
1997  Seff_HBDepth3->SetMinimum(80.);
1998  Seff_HBDepth3->Draw("COLZ");
1999  cHBnew->cd(4);
2000  TH2F *dva1_HBDepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA_HB");
2001  TH2F *dva0_HBDepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4_HB");
2002  TH2F *Seff_HBDepth4 = (TH2F *)dva1_HBDepth4->Clone("Seff_HBDepth4");
2003  Seff_HBDepth4->Divide(dva1_HBDepth4, dva0_HBDepth4, 25., 1., "B");
2004  gPad->SetGridy();
2005  gPad->SetGridx();
2006  Seff_HBDepth4->SetMarkerStyle(20);
2007  Seff_HBDepth4->SetMarkerSize(0.4);
2008  Seff_HBDepth4->SetXTitle("#eta \b");
2009  Seff_HBDepth4->SetYTitle("#phi \b");
2010  Seff_HBDepth4->SetZTitle("<timing> HB Depth4 \b");
2011  Seff_HBDepth4->SetMarkerColor(2);
2012  Seff_HBDepth4->SetLineColor(2);
2013  Seff_HBDepth4->SetMaximum(100.);
2014  Seff_HBDepth4->SetMinimum(80.);
2015  Seff_HBDepth4->Draw("COLZ");
2016  cHBnew->Modified();
2017  cHBnew->Update();
2018  cHBnew->Print("2DcorravstsnPLOTSHB.png");
2019  cHBnew->Clear();
2020  // std::cout << " 2DcorravstsnPLOTSHB.png created " << std::endl;
2021  // clean-up
2022  if (dva1_HBDepth1)
2023  delete dva1_HBDepth1;
2024  if (dva0_HBDepth1)
2025  delete dva0_HBDepth1;
2026  if (Seff_HBDepth1)
2027  delete Seff_HBDepth1;
2028  if (dva1_HBDepth2)
2029  delete dva1_HBDepth2;
2030  if (dva0_HBDepth2)
2031  delete dva0_HBDepth2;
2032  if (Seff_HBDepth2)
2033  delete Seff_HBDepth2;
2034  if (dva1_HBDepth3)
2035  delete dva1_HBDepth3;
2036  if (dva0_HBDepth3)
2037  delete dva0_HBDepth3;
2038  if (Seff_HBDepth3)
2039  delete Seff_HBDepth3;
2040  if (dva1_HBDepth4)
2041  delete dva1_HBDepth4;
2042  if (dva0_HBDepth4)
2043  delete dva0_HBDepth4;
2044  if (Seff_HBDepth4)
2045  delete Seff_HBDepth4;
2046 
2047  // 7 plots for HE:
2048  TCanvas *cHEnew = new TCanvas("cHEnew", "cHEnew", 5, 10, 1400, 1800);
2049  cHEnew->Clear();
2050  cHEnew->Divide(2, 4);
2051  cHEnew->cd(1);
2052  TH2F *dva1_HEDepth1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA_HE");
2053  TH2F *dva0_HEDepth1 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HE");
2054  // dva1_HEDepth1->Sumw2();
2055  // dva0_HEDepth1->Sumw2();
2056  TH2F *Seff_HEDepth1 = (TH2F *)dva1_HEDepth1->Clone("Seff_HEDepth1");
2057  Seff_HEDepth1->Divide(dva1_HEDepth1, dva0_HEDepth1, 25., 1., "B");
2058  gPad->SetGridy();
2059  gPad->SetGridx();
2060  Seff_HEDepth1->SetMarkerStyle(20);
2061  Seff_HEDepth1->SetMarkerSize(0.4);
2062  Seff_HEDepth1->SetXTitle("#eta \b");
2063  Seff_HEDepth1->SetYTitle("#phi \b");
2064  Seff_HEDepth1->SetZTitle("<timing> HE Depth1 \b");
2065  Seff_HEDepth1->SetMarkerColor(2);
2066  Seff_HEDepth1->SetLineColor(2);
2067  Seff_HEDepth1->SetMaximum(100.);
2068  Seff_HEDepth1->SetMinimum(80.);
2069  Seff_HEDepth1->Draw("COLZ");
2070  cHEnew->cd(2);
2071  TH2F *dva1_HEDepth2 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA_HE");
2072  TH2F *dva0_HEDepth2 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HE");
2073  TH2F *Seff_HEDepth2 = (TH2F *)dva1_HEDepth2->Clone("Seff_HEDepth2");
2074  Seff_HEDepth2->Divide(dva1_HEDepth2, dva0_HEDepth2, 25., 1., "B");
2075  gPad->SetGridy();
2076  gPad->SetGridx();
2077  Seff_HEDepth2->SetMarkerStyle(20);
2078  Seff_HEDepth2->SetMarkerSize(0.4);
2079  Seff_HEDepth2->SetXTitle("#eta \b");
2080  Seff_HEDepth2->SetYTitle("#phi \b");
2081  Seff_HEDepth2->SetZTitle("<timing> HE Depth2 \b");
2082  Seff_HEDepth2->SetMarkerColor(2);
2083  Seff_HEDepth2->SetLineColor(2);
2084  Seff_HEDepth2->SetMaximum(100.);
2085  Seff_HEDepth2->SetMinimum(80.);
2086  Seff_HEDepth2->Draw("COLZ");
2087  cHEnew->cd(3);
2088  TH2F *dva1_HEDepth3 = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmeanA_HE");
2089  TH2F *dva0_HEDepth3 = (TH2F *)dir->FindObjectAny("h_mapDepth3_HE");
2090  TH2F *Seff_HEDepth3 = (TH2F *)dva1_HEDepth3->Clone("Seff_HEDepth3");
2091  Seff_HEDepth3->Divide(dva1_HEDepth3, dva0_HEDepth3, 25., 1., "B");
2092  gPad->SetGridy();
2093  gPad->SetGridx();
2094  Seff_HEDepth3->SetMarkerStyle(20);
2095  Seff_HEDepth3->SetMarkerSize(0.4);
2096  Seff_HEDepth3->SetXTitle("#eta \b");
2097  Seff_HEDepth3->SetYTitle("#phi \b");
2098  Seff_HEDepth3->SetZTitle("<timing> HE Depth3 \b");
2099  Seff_HEDepth3->SetMarkerColor(2);
2100  Seff_HEDepth3->SetLineColor(2);
2101  Seff_HEDepth3->SetMaximum(100.);
2102  Seff_HEDepth3->SetMinimum(80.);
2103  Seff_HEDepth3->Draw("COLZ");
2104  cHEnew->cd(4);
2105  TH2F *dva1_HEDepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA_HE");
2106  TH2F *dva0_HEDepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4_HE");
2107  TH2F *Seff_HEDepth4 = (TH2F *)dva1_HEDepth4->Clone("Seff_HEDepth4");
2108  Seff_HEDepth4->Divide(dva1_HEDepth4, dva0_HEDepth4, 25., 1., "B");
2109  gPad->SetGridy();
2110  gPad->SetGridx();
2111  Seff_HEDepth4->SetMarkerStyle(20);
2112  Seff_HEDepth4->SetMarkerSize(0.4);
2113  Seff_HEDepth4->SetXTitle("#eta \b");
2114  Seff_HEDepth4->SetYTitle("#phi \b");
2115  Seff_HEDepth4->SetZTitle("<timing> HE Depth4 \b");
2116  Seff_HEDepth4->SetMarkerColor(2);
2117  Seff_HEDepth4->SetLineColor(2);
2118  Seff_HEDepth4->SetMaximum(100.);
2119  Seff_HEDepth4->SetMinimum(80.);
2120  Seff_HEDepth4->Draw("COLZ");
2121  cHEnew->cd(5);
2122  TH2F *dva1_HEDepth5 = (TH2F *)dir->FindObjectAny("h_mapDepth5TSmeanA_HE");
2123  TH2F *dva0_HEDepth5 = (TH2F *)dir->FindObjectAny("h_mapDepth5_HE");
2124  TH2F *Seff_HEDepth5 = (TH2F *)dva1_HEDepth5->Clone("Seff_HEDepth5");
2125  Seff_HEDepth5->Divide(dva1_HEDepth5, dva0_HEDepth5, 25., 1., "B");
2126  gPad->SetGridy();
2127  gPad->SetGridx();
2128  Seff_HEDepth5->SetMarkerStyle(20);
2129  Seff_HEDepth5->SetMarkerSize(0.4);
2130  Seff_HEDepth5->SetXTitle("#eta \b");
2131  Seff_HEDepth5->SetYTitle("#phi \b");
2132  Seff_HEDepth5->SetZTitle("<timing> HE Depth5 \b");
2133  Seff_HEDepth5->SetMarkerColor(2);
2134  Seff_HEDepth5->SetLineColor(2);
2135  Seff_HEDepth5->SetMaximum(100.);
2136  Seff_HEDepth5->SetMinimum(80.);
2137  Seff_HEDepth5->Draw("COLZ");
2138  cHEnew->cd(6);
2139  TH2F *dva1_HEDepth6 = (TH2F *)dir->FindObjectAny("h_mapDepth6TSmeanA_HE");
2140  TH2F *dva0_HEDepth6 = (TH2F *)dir->FindObjectAny("h_mapDepth6_HE");
2141  TH2F *Seff_HEDepth6 = (TH2F *)dva1_HEDepth6->Clone("Seff_HEDepth6");
2142  Seff_HEDepth6->Divide(dva1_HEDepth6, dva0_HEDepth6, 25., 1., "B");
2143  gPad->SetGridy();
2144  gPad->SetGridx();
2145  Seff_HEDepth6->SetMarkerStyle(20);
2146  Seff_HEDepth6->SetMarkerSize(0.4);
2147  Seff_HEDepth6->SetXTitle("#eta \b");
2148  Seff_HEDepth6->SetYTitle("#phi \b");
2149  Seff_HEDepth6->SetZTitle("<timing> HE Depth6 \b");
2150  Seff_HEDepth6->SetMarkerColor(2);
2151  Seff_HEDepth6->SetLineColor(2);
2152  Seff_HEDepth6->SetMaximum(100.);
2153  Seff_HEDepth6->SetMinimum(80.);
2154  Seff_HEDepth6->Draw("COLZ");
2155  cHEnew->cd(7);
2156  TH2F *dva1_HEDepth7 = (TH2F *)dir->FindObjectAny("h_mapDepth7TSmeanA_HE");
2157  TH2F *dva0_HEDepth7 = (TH2F *)dir->FindObjectAny("h_mapDepth7_HE");
2158  TH2F *Seff_HEDepth7 = (TH2F *)dva1_HEDepth7->Clone("Seff_HEDepth7");
2159  Seff_HEDepth7->Divide(dva1_HEDepth7, dva0_HEDepth7, 25., 1., "B");
2160  gPad->SetGridy();
2161  gPad->SetGridx();
2162  Seff_HEDepth7->SetMarkerStyle(20);
2163  Seff_HEDepth7->SetMarkerSize(0.4);
2164  Seff_HEDepth7->SetXTitle("#eta \b");
2165  Seff_HEDepth7->SetYTitle("#phi \b");
2166  Seff_HEDepth7->SetZTitle("<timing> HE Depth7 \b");
2167  Seff_HEDepth7->SetMarkerColor(2);
2168  Seff_HEDepth7->SetLineColor(2);
2169  Seff_HEDepth7->SetMaximum(100.);
2170  Seff_HEDepth7->SetMinimum(80.);
2171  Seff_HEDepth7->Draw("COLZ");
2172  cHEnew->Modified();
2173  cHEnew->Update();
2174  cHEnew->Print("2DcorravstsnPLOTSHE.png");
2175  cHEnew->Clear();
2176  // std::cout << " 2DcorravstsnPLOTSHE.png created " << std::endl;
2177  // clean-up
2178  if (dva1_HEDepth1)
2179  delete dva1_HEDepth1;
2180  if (dva0_HEDepth1)
2181  delete dva0_HEDepth1;
2182  if (Seff_HEDepth1)
2183  delete Seff_HEDepth1;
2184  if (dva1_HEDepth2)
2185  delete dva1_HEDepth2;
2186  if (dva0_HEDepth2)
2187  delete dva0_HEDepth2;
2188  if (Seff_HEDepth2)
2189  delete Seff_HEDepth2;
2190  if (dva1_HEDepth3)
2191  delete dva1_HEDepth3;
2192  if (dva0_HEDepth3)
2193  delete dva0_HEDepth3;
2194  if (Seff_HEDepth3)
2195  delete Seff_HEDepth3;
2196  if (dva1_HEDepth4)
2197  delete dva1_HEDepth4;
2198  if (dva0_HEDepth4)
2199  delete dva0_HEDepth4;
2200  if (Seff_HEDepth4)
2201  delete Seff_HEDepth4;
2202  if (dva1_HEDepth5)
2203  delete dva1_HEDepth5;
2204  if (dva0_HEDepth5)
2205  delete dva0_HEDepth5;
2206  if (Seff_HEDepth5)
2207  delete Seff_HEDepth5;
2208  if (dva1_HEDepth6)
2209  delete dva1_HEDepth6;
2210  if (dva0_HEDepth6)
2211  delete dva0_HEDepth6;
2212  if (Seff_HEDepth6)
2213  delete Seff_HEDepth6;
2214  if (dva1_HEDepth7)
2215  delete dva1_HEDepth7;
2216  if (dva0_HEDepth7)
2217  delete dva0_HEDepth7;
2218  if (Seff_HEDepth7)
2219  delete Seff_HEDepth7;
2220  // 4 plots for HF:
2221  TCanvas *cHFnew = new TCanvas("cHFnew", "cHFnew", 200, 10, 1400, 1800);
2222  cHFnew->Clear();
2223  cHFnew->Divide(2, 2);
2224  cHFnew->cd(1);
2225  TH2F *dva1_HFDepth1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA_HF");
2226  TH2F *dva0_HFDepth1 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HF");
2227  // dva1_HFDepth1->Sumw2();
2228  // dva0_HFDepth1->Sumw2();
2229  TH2F *Seff_HFDepth1 = (TH2F *)dva1_HFDepth1->Clone("Seff_HFDepth1");
2230  Seff_HFDepth1->Divide(dva1_HFDepth1, dva0_HFDepth1, 25., 1., "B");
2231  /*
2232  for (int i=1;i<=Seff_HFDepth1->GetXaxis()->GetNbins();i++) {
2233  for (int j=1;j<=Seff_HFDepth1->GetYaxis()->GetNbins();j++) {
2234  double ccc1 = Seff_HFDepth1->GetBinContent(i,j) ;
2235  // if(ccc1 > 0.) std::cout << "******************** i = " << i << " j = " << j << " ccc1 = " << ccc1 << std::endl;
2236 
2237  Seff_HFDepth1->SetBinContent(i,j,0.);
2238  if(ccc1 > 0.) Seff_HFDepth1->SetBinContent(i,j,ccc1);
2239  }
2240  }
2241 */
2242  gPad->SetGridy();
2243  gPad->SetGridx();
2244  Seff_HFDepth1->SetMarkerStyle(20);
2245  Seff_HFDepth1->SetMarkerSize(0.4);
2246  Seff_HFDepth1->SetXTitle("#eta \b");
2247  Seff_HFDepth1->SetYTitle("#phi \b");
2248  Seff_HFDepth1->SetZTitle("<timing> HF Depth1 \b");
2249  Seff_HFDepth1->SetMarkerColor(2);
2250  Seff_HFDepth1->SetLineColor(2);
2251  Seff_HFDepth1->SetMaximum(50.);
2252  Seff_HFDepth1->SetMinimum(20.);
2253  Seff_HFDepth1->Draw("COLZ");
2254  cHFnew->cd(2);
2255  TH2F *dva1_HFDepth2 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA_HF");
2256  TH2F *dva0_HFDepth2 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HF");
2257  TH2F *Seff_HFDepth2 = (TH2F *)dva1_HFDepth2->Clone("Seff_HFDepth2");
2258  Seff_HFDepth2->Divide(dva1_HFDepth2, dva0_HFDepth2, 25., 1., "B");
2259  /*
2260  for (int i=1;i<=Seff_HFDepth2->GetXaxis()->GetNbins();i++) {
2261  for (int j=1;j<=Seff_HFDepth2->GetYaxis()->GetNbins();j++) {
2262  double ccc1 = Seff_HFDepth2->GetBinContent(i,j) ;
2263  Seff_HFDepth2->SetBinContent(i,j,0.);
2264  if(ccc1 > 0. ) Seff_HFDepth2->SetBinContent(i,j,ccc1);
2265  }
2266  }
2267  */
2268  gPad->SetGridy();
2269  gPad->SetGridx();
2270  Seff_HFDepth2->SetMarkerStyle(20);
2271  Seff_HFDepth2->SetMarkerSize(0.4);
2272  Seff_HFDepth2->SetXTitle("#eta \b");
2273  Seff_HFDepth2->SetYTitle("#phi \b");
2274  Seff_HFDepth2->SetZTitle("<timing> HF Depth2 \b");
2275  Seff_HFDepth2->SetMarkerColor(2);
2276  Seff_HFDepth2->SetLineColor(2);
2277  Seff_HFDepth2->SetMaximum(50.);
2278  Seff_HFDepth2->SetMinimum(20.);
2279  Seff_HFDepth2->Draw("COLZ");
2280  cHFnew->cd(3);
2281  TH2F *dva1_HFDepth3 = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmeanA_HF");
2282  TH2F *dva0_HFDepth3 = (TH2F *)dir->FindObjectAny("h_mapDepth3_HF");
2283  TH2F *Seff_HFDepth3 = (TH2F *)dva1_HFDepth3->Clone("Seff_HFDepth3");
2284  Seff_HFDepth3->Divide(dva1_HFDepth3, dva0_HFDepth3, 25., 1., "B");
2285  gPad->SetGridy();
2286  gPad->SetGridx();
2287  Seff_HFDepth3->SetMarkerStyle(20);
2288  Seff_HFDepth3->SetMarkerSize(0.4);
2289  Seff_HFDepth3->SetXTitle("#eta \b");
2290  Seff_HFDepth3->SetYTitle("#phi \b");
2291  Seff_HFDepth3->SetZTitle("<timing> HF Depth3 \b");
2292  Seff_HFDepth3->SetMarkerColor(2);
2293  Seff_HFDepth3->SetLineColor(2);
2294  Seff_HFDepth3->SetMaximum(50.);
2295  Seff_HFDepth3->SetMinimum(20.);
2296  Seff_HFDepth3->Draw("COLZ");
2297  cHFnew->cd(4);
2298  TH2F *dva1_HFDepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA_HF");
2299  TH2F *dva0_HFDepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4_HF");
2300  TH2F *Seff_HFDepth4 = (TH2F *)dva1_HFDepth4->Clone("Seff_HFDepth4");
2301  Seff_HFDepth4->Divide(dva1_HFDepth4, dva0_HFDepth4, 25., 1., "B");
2302  gPad->SetGridy();
2303  gPad->SetGridx();
2304  Seff_HFDepth4->SetMarkerStyle(20);
2305  Seff_HFDepth4->SetMarkerSize(0.4);
2306  Seff_HFDepth4->SetXTitle("#eta \b");
2307  Seff_HFDepth4->SetYTitle("#phi \b");
2308  Seff_HFDepth4->SetZTitle("<timing> HF Depth4 \b");
2309  Seff_HFDepth4->SetMarkerColor(2);
2310  Seff_HFDepth4->SetLineColor(2);
2311  Seff_HFDepth4->SetMaximum(50.);
2312  Seff_HFDepth4->SetMinimum(20.);
2313  Seff_HFDepth4->Draw("COLZ");
2314  cHFnew->Modified();
2315  cHFnew->Update();
2316  cHFnew->Print("2DcorravstsnPLOTSHF.png");
2317  cHFnew->Clear();
2318  // std::cout << " 2DcorravstsnPLOTSHF.png created " << std::endl;
2319  if (dva1_HFDepth1)
2320  delete dva1_HFDepth1;
2321  if (dva0_HFDepth1)
2322  delete dva0_HFDepth1;
2323  if (Seff_HFDepth1)
2324  delete Seff_HFDepth1;
2325  if (dva1_HFDepth2)
2326  delete dva1_HFDepth2;
2327  if (dva0_HFDepth2)
2328  delete dva0_HFDepth2;
2329  if (Seff_HFDepth2)
2330  delete Seff_HFDepth2;
2331  if (dva1_HFDepth3)
2332  delete dva1_HFDepth3;
2333  if (dva0_HFDepth3)
2334  delete dva0_HFDepth3;
2335  if (Seff_HFDepth3)
2336  delete Seff_HFDepth3;
2337  if (dva1_HFDepth4)
2338  delete dva1_HFDepth4;
2339  if (dva0_HFDepth4)
2340  delete dva0_HFDepth4;
2341  if (Seff_HFDepth4)
2342  delete Seff_HFDepth4;
2343  // 1 plot for HO:
2344  TCanvas *cHOnew = new TCanvas("cHOnew", "cHOnew", 200, 10, 1400, 1800);
2345  // TCanvas *cHOnew = new TCanvas("cHOnew", "cHOnew", 1500, 500);
2346  cHOnew->Clear();
2347  cHOnew->Divide(1, 1);
2348  cHOnew->cd(1);
2349  TH2F *dva1_HODepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA_HO");
2350  /*
2351  for (int i=1;i<=dva1_HODepth4->GetXaxis()->GetNbins();i++) {
2352  for (int j=1;j<=dva1_HODepth4->GetYaxis()->GetNbins();j++) {
2353  double ccc1 = dva1_HODepth4->GetBinContent(i,j) ;
2354  if(ccc1 > 0.) std::cout << "****** dva1_HODepth4 ************** i = " << i << " j = " << j << " ccc1 = " << ccc1 << std::endl;
2355  }
2356  }
2357 */
2358  TH2F *dva0_HODepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4_HO");
2359  /*
2360  for (int i=1;i<=dva0_HODepth4->GetXaxis()->GetNbins();i++) {
2361  for (int j=1;j<=dva0_HODepth4->GetYaxis()->GetNbins();j++) {
2362  double ccc1 = dva0_HODepth4->GetBinContent(i,j) ;
2363  if(ccc1 > 0.) std::cout << "****** dva0_HODepth4 ************** i = " << i << " j = " << j << " ccc1 = " << ccc1 << std::endl;
2364  }
2365  }
2366  */
2367  TH2F *Seff_HODepth4 = (TH2F *)dva1_HODepth4->Clone("Seff_HODepth4");
2368  /*
2369  for (int x = 1; x <= Seff_HODepth4->GetXaxis()->GetNbins(); x++) {
2370  for (int y = 1; y <= Seff_HODepth4->GetYaxis()->GetNbins(); y++) {
2371  // dva1_HODepth4->SetBinError(float(x), float(y), 0.001);
2372  // Seff_HODepth4->SetBinContent(float(x), float(y), 0.0);
2373  } //end x
2374  } //end y
2375 */
2376  Seff_HODepth4->Divide(dva1_HODepth4, dva0_HODepth4, 25., 1., "B");
2377  /*
2378  for (int i=1;i<=Seff_HODepth4->GetXaxis()->GetNbins();i++) {
2379  for (int j=1;j<=Seff_HODepth4->GetYaxis()->GetNbins();j++) {
2380  double ccc1 = Seff_HODepth4->GetBinContent(i,j);
2381  if(ccc1 > 0.) std::cout << "****** Seff_HODepth4 ************** i = " << i << " j = " << j << " ccc1 = " << ccc1 << std::endl;
2382  }
2383  }
2384 */
2385  gPad->SetGridy();
2386  gPad->SetGridx();
2387  Seff_HODepth4->SetMarkerStyle(20);
2388  Seff_HODepth4->SetMarkerSize(0.4);
2389  Seff_HODepth4->SetXTitle("#eta \b");
2390  Seff_HODepth4->SetYTitle("#phi \b");
2391  Seff_HODepth4->SetZTitle("<timing> HO Depth4 \b");
2392  Seff_HODepth4->SetMarkerColor(2);
2393  Seff_HODepth4->SetLineColor(2);
2394  Seff_HODepth4->SetMaximum(130.);
2395  Seff_HODepth4->SetMinimum(70.);
2396  Seff_HODepth4->Draw("COLZ");
2397  cHOnew->Modified();
2398  cHOnew->Update();
2399  cHOnew->Print("2DcorravstsnPLOTSHO.png");
2400  cHOnew->Clear();
2401  // std::cout << " 2DcorravstsnPLOTSHO.png created " << std::endl;
2402  if (dva1_HODepth4)
2403  delete dva1_HODepth4;
2404  if (dva0_HODepth4)
2405  delete dva0_HODepth4;
2406  if (Seff_HODepth4)
2407  delete Seff_HODepth4;
2408 
2409  // std::cout << " END OF 2023 " << std::endl;
2410  //
2412  //======================================================================
2413 
2414  //====================================================================== html pages CREATING:
2415  std::cout << " html pages CREATING: " << std::endl;
2416  //======================================================================
2417  // Creating each test kind for each subdet html pages:
2418  std::string raw_class, raw_class1, raw_class2, raw_class3;
2419  int ind = 0;
2420 
2421  for (int sub = 1; sub <= 4; sub++) { //Subdetector: 1-HB, 2-HE, 3-HO, 4-HF
2422  ofstream htmlFileT, htmlFileC, htmlFileD, htmlFileP, htmlFileS;
2423  if (sub == 1) {
2424  htmlFileT.open("HB_Tile.html");
2425  htmlFileC.open("HB_Calib.html");
2426  htmlFileD.open("HB_Drift.html");
2427  htmlFileP.open("HB_Pedestals.html");
2428  htmlFileS.open("HB_Shapes.html");
2429  }
2430  if (sub == 2) {
2431  htmlFileT.open("HE_Tile.html");
2432  htmlFileC.open("HE_Calib.html");
2433  htmlFileD.open("HE_Drift.html");
2434  htmlFileP.open("HE_Pedestals.html");
2435  htmlFileS.open("HE_Shapes.html");
2436  }
2437  if (sub == 3) {
2438  htmlFileT.open("HO_Tile.html");
2439  htmlFileC.open("HO_Calib.html");
2440  htmlFileD.open("HO_Drift.html");
2441  htmlFileP.open("HO_Pedestals.html");
2442  htmlFileS.open("HO_Shapes.html");
2443  }
2444  if (sub == 4) {
2445  htmlFileT.open("HF_Tile.html");
2446  htmlFileC.open("HF_Calib.html");
2447  htmlFileD.open("HF_Drift.html");
2448  htmlFileP.open("HF_Pedestals.html");
2449  htmlFileS.open("HF_Shapes.html");
2450  }
2451 
2452  // Megatile channels
2453  htmlFileT << "</html><html xmlns=\"http://www.w3.org/1999/xhtml\">" << std::endl;
2454  htmlFileT << "<head>" << std::endl;
2455  htmlFileT << "<meta http-equiv=\"Content-Type\" content=\"text/html\"/>" << std::endl;
2456  htmlFileT << "<title> Remote Monitoring Tool Global</title>" << std::endl;
2457  htmlFileT << "<style type=\"text/css\">" << std::endl;
2458  htmlFileT << " body,td{ background-color: #FFFFCC; font-family: arial, arial ce, helvetica; font-size: 12px; }"
2459  << std::endl;
2460  htmlFileT << " td.s0 { font-family: arial, arial ce, helvetica; }" << std::endl;
2461  htmlFileT << " td.s1 { font-family: arial, arial ce, helvetica; font-weight: bold; background-color: #FFC169; "
2462  "text-align: center;}"
2463  << std::endl;
2464  htmlFileT << " td.s2 { font-family: arial, arial ce, helvetica; background-color: #eeeeee; }" << std::endl;
2465  htmlFileT << " td.s3 { font-family: arial, arial ce, helvetica; background-color: #d0d0d0; }" << std::endl;
2466  htmlFileT << " td.s4 { font-family: arial, arial ce, helvetica; background-color: #FFC169; }" << std::endl;
2467  htmlFileT << "</style>" << std::endl;
2468  htmlFileT << "<body>" << std::endl;
2469 
2470  if (sub == 1)
2471  htmlFileT << "<h1> Criteria for megatile channels for HB, RUN = " << runnumber << " </h1>" << std::endl;
2472  if (sub == 2)
2473  htmlFileT << "<h1> Criteria for megatile channels for HE, RUN = " << runnumber << " </h1>" << std::endl;
2474  if (sub == 3)
2475  htmlFileT << "<h1> Criteria for megatile channels for HO, RUN = " << runnumber << " </h1>" << std::endl;
2476  if (sub == 4)
2477  htmlFileT << "<h1> Criteria for megatile channels for HF, RUN = " << runnumber << " </h1>" << std::endl;
2478  htmlFileT << "<br>" << std::endl;
2479 
2480  // Test Entries
2481 
2482  htmlFileT << "<h2> 0. Entries for each channel.</h3>" << std::endl;
2483  htmlFileT << "<h3> 0.A. Entries in each channel for each depth.</h3>" << std::endl;
2484  htmlFileT << "<h4> Channel legend: color is number of hits in digi collection </h4>" << std::endl;
2485  if (sub == 1)
2486  htmlFileT << " <img src=\"MapRateEntryHB.png\" />" << std::endl;
2487  if (sub == 2)
2488  htmlFileT << " <img src=\"MapRateEntryHE.png\" />" << std::endl;
2489  if (sub == 3)
2490  htmlFileT << " <img src=\"MapRateEntryHO.png\" />" << std::endl;
2491  if (sub == 4)
2492  htmlFileT << " <img src=\"MapRateEntryHF.png\" />" << std::endl;
2493  htmlFileT << "<br>" << std::endl;
2494 
2495  // Test Cm
2496  htmlFileT << "<h2> 1. Cm criterion: CapID errors for each channel.</h3>" << std::endl;
2497  htmlFileT << "<h3> 1.A. Rate of CapId failures in each channel for each depth.</h3>" << std::endl;
2498  htmlFileT << "<h4> Channel legend: white - good, other colour - bad. </h4>" << std::endl;
2499  if (sub == 1)
2500  htmlFileT << " <img src=\"MapRateCapIDHB.png\" />" << std::endl;
2501  if (sub == 2)
2502  htmlFileT << " <img src=\"MapRateCapIDHE.png\" />" << std::endl;
2503  if (sub == 3)
2504  htmlFileT << " <img src=\"MapRateCapIDHO.png\" />" << std::endl;
2505  if (sub == 4)
2506  htmlFileT << " <img src=\"MapRateCapIDHF.png\" />" << std::endl;
2507  htmlFileT << "<br>" << std::endl;
2508 
2509  // Am
2510  htmlFileT << "<h2> 2. Am criterion: ADC amplitude collected over all TSs(Full Amplitude) for each channel. </h3>"
2511  << std::endl;
2512  htmlFileT << "<h3> 2.A. Full ADC amplitude distribution over all events, channels and depths.</h3>" << std::endl;
2513  htmlFileT << "<h4> Legend: Bins less " << MIN_M[2][sub] << " correpond to bad ADC amplitude </h4>" << std::endl;
2514  if (sub == 1)
2515  htmlFileT << " <img src=\"HistAmplHB.png\" />" << std::endl;
2516  if (sub == 2)
2517  htmlFileT << " <img src=\"HistAmplHE.png\" />" << std::endl;
2518  if (sub == 3)
2519  htmlFileT << " <img src=\"HistAmplHO.png\" />" << std::endl;
2520  if (sub == 4)
2521  htmlFileT << " <img src=\"HistAmplHF.png\" />" << std::endl;
2522  htmlFileT << "<br>" << std::endl;
2523  htmlFileT << "<h3> 2.B. Rate of bad ADC amplitude (<" << MIN_M[2][sub] << ") in each channel for each depth. </h3>"
2524  << std::endl;
2525  htmlFileT << "<h4> Channel legend: white - good, other colours - bad. </h4>" << std::endl;
2526  if (sub == 1)
2527  htmlFileT << " <img src=\"MapRateAmplHB.png\" />" << std::endl;
2528  if (sub == 2)
2529  htmlFileT << " <img src=\"MapRateAmplHE.png\" />" << std::endl;
2530  if (sub == 3)
2531  htmlFileT << " <img src=\"MapRateAmplHO.png\" />" << std::endl;
2532  if (sub == 4)
2533  htmlFileT << " <img src=\"MapRateAmplHF.png\" />" << std::endl;
2534  htmlFileT << "<br>" << std::endl;
2535 
2536  // Test Wm
2537  htmlFileT << "<h2> 3. Wm criterion: RMS (width) of ADC amplutude for each channel.</h3>" << std::endl;
2538  htmlFileT << "<h3> 3.A. RMS distribution over all events, channel and depth.</h3>" << std::endl;
2539  htmlFileT << "<h4> Legend: Bins less " << MIN_M[3][sub] << " and more " << MAX_M[3][sub]
2540  << " correpond to bad RMS </h4>" << std::endl;
2541  if (sub == 1)
2542  htmlFileT << " <img src=\"HistRMSHB.png\" />" << std::endl;
2543  if (sub == 2)
2544  htmlFileT << " <img src=\"HistRMSHE.png\" />" << std::endl;
2545  if (sub == 3)
2546  htmlFileT << " <img src=\"HistRMSHO.png\" />" << std::endl;
2547  if (sub == 4)
2548  htmlFileT << " <img src=\"HistRMSHF.png\" />" << std::endl;
2549  htmlFileT << "<br>" << std::endl;
2550  htmlFileT << "<h3> 3.B. Rate of bad RMS (<" << MIN_M[3][sub] << ",>" << MAX_M[3][sub]
2551  << ") in each channel for each depth.</h3>" << std::endl;
2552  htmlFileT << "<h4> Channel legend: white - good, other colour - bad. </h4>" << std::endl;
2553  if (sub == 1)
2554  htmlFileT << " <img src=\"MapRateRMSHB.png\" />" << std::endl;
2555  if (sub == 2)
2556  htmlFileT << " <img src=\"MapRateRMSHE.png\" />" << std::endl;
2557  if (sub == 3)
2558  htmlFileT << " <img src=\"MapRateRMSHO.png\" />" << std::endl;
2559  if (sub == 4)
2560  htmlFileT << " <img src=\"MapRateRMSHF.png\" />" << std::endl;
2561  htmlFileT << "<br>" << std::endl;
2562 
2563  // Rm
2564  htmlFileT << "<h2> 4. Rm criterion: Ratio ADC value sum over four near maximum (-2, -1, max, +1) TS to ADC value "
2565  "sum over all TS for each channel. </h3>"
2566  << std::endl;
2567  htmlFileT << "<h3> 4.A. Ratio distribution over all events, channels and depths.</h3>" << std::endl;
2568  htmlFileT << "<h4> Legend: Bins less " << MIN_M[4][sub] << " and more " << MAX_M[4][sub]
2569  << " correpond to bad ratio </h4>" << std::endl;
2570  if (sub == 1)
2571  htmlFileT << " <img src=\"Hist43TStoAllTSHB.png\" />" << std::endl;
2572  if (sub == 2)
2573  htmlFileT << " <img src=\"Hist43TStoAllTSHE.png\" />" << std::endl;
2574  if (sub == 3)
2575  htmlFileT << " <img src=\"Hist43TStoAllTSHO.png\" />" << std::endl;
2576  if (sub == 4)
2577  htmlFileT << " <img src=\"Hist43TStoAllTSHF.png\" />" << std::endl;
2578  htmlFileT << "<br>" << std::endl;
2579  htmlFileT << "<h3> 4.B. Rate of bad ratio (<" << MIN_M[4][sub] << ", >" << MAX_M[4][sub]
2580  << ") in each channel for each depth.</h3>" << std::endl;
2581  htmlFileT << "<h4> Channel legend: white - good, other colour - bad. </h4>" << std::endl;
2582  if (sub == 1)
2583  htmlFileT << " <img src=\"MapRate43TStoAllTSHB.png\" />" << std::endl;
2584  if (sub == 2)
2585  htmlFileT << " <img src=\"MapRate43TStoAllTSHE.png\" />" << std::endl;
2586  if (sub == 3)
2587  htmlFileT << " <img src=\"MapRate43TStoAllTSHO.png\" />" << std::endl;
2588  if (sub == 4)
2589  htmlFileT << " <img src=\"MapRate43TStoAllTSHF.png\" />" << std::endl;
2590  htmlFileT << "<br>" << std::endl;
2591 
2592  // TNm
2593  htmlFileT << "<h2> 5. TNm criterion: Mean TS position for each channel.</h3>" << std::endl;
2594  htmlFileT << "<h3> 5.A. TN position distribution over all events, channels and depths.</h3>" << std::endl;
2595  htmlFileT << "<h4> Legend: Bins less " << MIN_M[5][sub] << " and more " << MAX_M[5][sub]
2596  << " correpond to bad mean position </h4>" << std::endl;
2597  if (sub == 1)
2598  htmlFileT << " <img src=\"HistMeanPosHB.png\" />" << std::endl;
2599  if (sub == 2)
2600  htmlFileT << " <img src=\"HistMeanPosHE.png\" />" << std::endl;
2601  if (sub == 3)
2602  htmlFileT << " <img src=\"HistMeanPosHO.png\" />" << std::endl;
2603  if (sub == 4)
2604  htmlFileT << " <img src=\"HistMeanPosHF.png\" />" << std::endl;
2605  htmlFileT << "<br>" << std::endl;
2606  htmlFileT << "<h3> 5.B. Rate of bad TN position (<" << MIN_M[5][sub] << ", >" << MAX_M[5][sub]
2607  << ") in each channel for each depth. </h3>" << std::endl;
2608  htmlFileT << "<h4> Channel legend: white - good, other colour - bad. </h4>" << std::endl;
2609  if (sub == 1)
2610  htmlFileT << " <img src=\"MapRateMeanPosHB.png\" />" << std::endl;
2611  if (sub == 2)
2612  htmlFileT << " <img src=\"MapRateMeanPosHE.png\" />" << std::endl;
2613  if (sub == 3)
2614  htmlFileT << " <img src=\"MapRateMeanPosHO.png\" />" << std::endl;
2615  if (sub == 4)
2616  htmlFileT << " <img src=\"MapRateMeanPosHF.png\" />" << std::endl;
2617  htmlFileT << "<br>" << std::endl;
2618 
2619  // TXm
2620  htmlFileT << "<h2> 6.TXm criterion: Maximum TS position for each channel.</h3>" << std::endl;
2621  htmlFileT << "<h3> 6.A. TX position distribution over all events, channel and depth.</h3>" << std::endl;
2622  htmlFileT << "<h4> Legend: Bins less " << MIN_M[6][sub] << " and more " << MAX_M[6][sub]
2623  << " correpond to bad position </h4>" << std::endl;
2624  if (sub == 1)
2625  htmlFileT << " <img src=\"HistMaxPosHB.png\" />" << std::endl;
2626  if (sub == 2)
2627  htmlFileT << " <img src=\"HistMaxPosHE.png\" />" << std::endl;
2628  if (sub == 3)
2629  htmlFileT << " <img src=\"HistMaxPosHO.png\" />" << std::endl;
2630  if (sub == 4)
2631  htmlFileT << " <img src=\"HistMaxPosHF.png\" />" << std::endl;
2632  htmlFileT << "<br>" << std::endl;
2633  htmlFileT << "<h3> 6.B. Rate of bad TX position (<" << MIN_M[6][sub] << ", >" << MAX_M[6][sub]
2634  << ") in each channel for each depth. </h3>" << std::endl;
2635  htmlFileT << "<h4> Channel legend: white - good, other colour - bad. </h4>" << std::endl;
2636  if (sub == 1)
2637  htmlFileT << " <img src=\"MapRateMaxPosHB.png\" />" << std::endl;
2638  if (sub == 2)
2639  htmlFileT << " <img src=\"MapRateMaxPosHE.png\" />" << std::endl;
2640  if (sub == 3)
2641  htmlFileT << " <img src=\"MapRateMaxPosHO.png\" />" << std::endl;
2642  if (sub == 4)
2643  htmlFileT << " <img src=\"MapRateMaxPosHF.png\" />" << std::endl;
2644  htmlFileT << "<br>" << std::endl;
2645 
2646  // Correlation of A vs TSn done in 2023 for Run3 and for GlobalRMT only
2647  htmlFileT << "<h2> 7....... Correlation of A(=Q) vs timing(=25ns*MeanTSposition) </h3>" << std::endl;
2648 
2649  htmlFileT << "<h3> 7.A..... 1)2D-correlation of timing vs Q,fc;.......... 2)Q,fc;................ 3)mean timing vs "
2650  "Q,fc ....... </h3>"
2651  << std::endl;
2652  if (sub == 1)
2653  htmlFileT << " <img src=\"corravstsnPLOTSHB.png\" />" << std::endl;
2654  if (sub == 2)
2655  htmlFileT << " <img src=\"corravstsnPLOTSHE.png\" />" << std::endl;
2656  if (sub == 3)
2657  htmlFileT << " <img src=\"corravstsnPLOTSHO.png\" />" << std::endl;
2658  if (sub == 4)
2659  htmlFileT << " <img src=\"corravstsnPLOTSHF.png\" />" << std::endl;
2660  htmlFileT << "<br>" << std::endl;
2661 
2662  htmlFileT << "<h3> 7.B....... Mean timing in 2D space of eta-phi for different Depthes........ </h3>" << std::endl;
2663  if (sub == 1)
2664  htmlFileT << " <img src=\"2DcorravstsnPLOTSHB.png\" />" << std::endl;
2665  if (sub == 2)
2666  htmlFileT << " <img src=\"2DcorravstsnPLOTSHE.png\" />" << std::endl;
2667  if (sub == 3)
2668  htmlFileT << " <img src=\"2DcorravstsnPLOTSHO.png\" />" << std::endl;
2669  if (sub == 4)
2670  htmlFileT << " <img src=\"2DcorravstsnPLOTSHF.png\" />" << std::endl;
2671  htmlFileT << "<br>" << std::endl;
2672 
2673  htmlFileT << "</body> " << std::endl;
2674  htmlFileT << "</html> " << std::endl;
2675  htmlFileT.close();
2676 
2677  // Pedestals
2678  htmlFileP << "</html><html xmlns=\"http://www.w3.org/1999/xhtml\">" << std::endl;
2679  htmlFileP << "<head>" << std::endl;
2680  htmlFileP << "<meta http-equiv=\"Content-Type\" content=\"text/html\"/>" << std::endl;
2681  htmlFileP << "<title> Remote Monitoring Tool Global</title>" << std::endl;
2682  htmlFileP << "<style type=\"text/css\">" << std::endl;
2683  htmlFileP << " body,td{ background-color: #FFFFCC; font-family: arial, arial ce, helvetica; font-size: 12px; }"
2684  << std::endl;
2685  htmlFileP << " td.s0 { font-family: arial, arial ce, helvetica; }" << std::endl;
2686  htmlFileP << " td.s1 { font-family: arial, arial ce, helvetica; font-weight: bold; background-color: #FFC169; "
2687  "text-align: center;}"
2688  << std::endl;
2689  htmlFileP << " td.s2 { font-family: arial, arial ce, helvetica; background-color: #eeeeee; }" << std::endl;
2690  htmlFileP << " td.s3 { font-family: arial, arial ce, helvetica; background-color: #d0d0d0; }" << std::endl;
2691  htmlFileP << " td.s4 { font-family: arial, arial ce, helvetica; background-color: #FFC169; }" << std::endl;
2692  htmlFileP << "</style>" << std::endl;
2693  htmlFileP << "<body>" << std::endl;
2694 
2695  if (sub == 1)
2696  htmlFileP << "<h1> Pedestals for HB, RUN = " << runnumber << " </h1>" << std::endl;
2697  if (sub == 2)
2698  htmlFileP << "<h1> Pedestals for HE, RUN = " << runnumber << " </h1>" << std::endl;
2699  if (sub == 3)
2700  htmlFileP << "<h1> Pedestals for HO, RUN = " << runnumber << " </h1>" << std::endl;
2701  if (sub == 4)
2702  htmlFileP << "<h1> Pedestals for HF, RUN = " << runnumber << " </h1>" << std::endl;
2703  htmlFileP << "<br>" << std::endl;
2704 
2705  // Pedestal:
2706  htmlFileP << "<h2> 1.Pm criterion: Pedestals for each CapID .</h3>" << std::endl;
2707  htmlFileP << "<h3> 1.A. Pedestal distribution over all events, channels for each CapID and all depths.</h3>"
2708  << std::endl;
2709  htmlFileP << "<h4> Legend: Bins less " << Pedest[0][sub] << " correpond to bad Pedestals </h4>" << std::endl;
2710  if (sub == 1)
2711  htmlFileP << " <img src=\"HistPedestalsHB.png\" />" << std::endl;
2712  if (sub == 2)
2713  htmlFileP << " <img src=\"HistPedestalsHE.png\" />" << std::endl;
2714  if (sub == 3)
2715  htmlFileP << " <img src=\"HistPedestalsHO.png\" />" << std::endl;
2716  if (sub == 4)
2717  htmlFileP << " <img src=\"HistPedestalsHF.png\" />" << std::endl;
2718  htmlFileP << "<br>" << std::endl;
2719  htmlFileP << "<h3> 1.B. Rate of channels at very low Pedestals at least in one CapID for each depth.</h3>"
2720  << std::endl;
2721  htmlFileP << "<h4> Channel legend: white - good, other colour - bad. </h4>" << std::endl;
2722  if (sub == 1)
2723  htmlFileP << " <img src=\"MapRatePedHB.png\" />" << std::endl;
2724  if (sub == 2)
2725  htmlFileP << " <img src=\"MapRatePedHE.png\" />" << std::endl;
2726  if (sub == 3)
2727  htmlFileP << " <img src=\"MapRatePedHO.png\" />" << std::endl;
2728  if (sub == 4)
2729  htmlFileP << " <img src=\"MapRatePedHF.png\" />" << std::endl;
2730 
2731  // PedestalWidth:
2732  htmlFileP << "<h2> 2.pWm criterion: Pedestal Widths for each CapID .</h3>" << std::endl;
2733  htmlFileP << "<h3> 2.A. Pedestal Widths distribution over all events, channels for each CapID and all depths.</h3>"
2734  << std::endl;
2735  htmlFileP << "<h4> Legend: Bins less " << Pedest[1][sub] << " correpond to bad Pedestal Widths </h4>" << std::endl;
2736  if (sub == 1)
2737  htmlFileP << " <img src=\"HistPedestalWidthsHB.png\" />" << std::endl;
2738  if (sub == 2)
2739  htmlFileP << " <img src=\"HistPedestalWidthsHE.png\" />" << std::endl;
2740  if (sub == 3)
2741  htmlFileP << " <img src=\"HistPedestalWidthsHO.png\" />" << std::endl;
2742  if (sub == 4)
2743  htmlFileP << " <img src=\"HistPedestalWidthsHF.png\" />" << std::endl;
2744  htmlFileP << "<br>" << std::endl;
2745  htmlFileP << "<h3> 2.B. Rate of channels at very low Pedestal Widths at least in one CapID for each depth.</h3>"
2746  << std::endl;
2747  htmlFileP << "<h4> Channel legend: white - good, other colour - bad. </h4>" << std::endl;
2748  if (sub == 1)
2749  htmlFileP << " <img src=\"MapRatePedWidthsHB.png\" />" << std::endl;
2750  if (sub == 2)
2751  htmlFileP << " <img src=\"MapRatePedWidthsHE.png\" />" << std::endl;
2752  if (sub == 3)
2753  htmlFileP << " <img src=\"MapRatePedWidthsHO.png\" />" << std::endl;
2754  if (sub == 4)
2755  htmlFileP << " <img src=\"MapRatePedWidthsHF.png\" />" << std::endl;
2756 
2757  // Correlations of Pedestal(Width) and fullAmplitude:
2758  htmlFileP << "<h2> 3.Pedestal and pedestalWidths vs Amplitude .</h3>" << std::endl;
2759  htmlFileP << "<h3> 3.A. Correlation of Pedestal(pedestalWidths) and Amplitude over all channels and events .</h3>"
2760  << std::endl;
2761  htmlFileP << "<h4> Legend: colour - entries </h4>" << std::endl;
2762  if (sub == 1)
2763  htmlFileP << "<img src=\"CorrelationsMapPedestalVsfullAmplitudeHB.png\" />" << std::endl;
2764  if (sub == 2)
2765  htmlFileP << "<img src=\"CorrelationsMapPedestalVsfullAmplitudeHE.png\" />" << std::endl;
2766  if (sub == 3)
2767  htmlFileP << "<img src=\"CorrelationsMapPedestalVsfullAmplitudeHO.png\" />" << std::endl;
2768  if (sub == 4)
2769  htmlFileP << "<img src=\"CorrelationsMapPedestalVsfullAmplitudeHF.png\" />" << std::endl;
2770  htmlFileP << "<br>" << std::endl;
2771 
2772  // TSs Shapes:
2773 
2774  htmlFileS << "</html><html xmlns=\"http://www.w3.org/1999/xhtml\">" << std::endl;
2775  htmlFileS << "<head>" << std::endl;
2776  htmlFileS << "<meta http-equiv=\"Content-Type\" content=\"text/html\"/>" << std::endl;
2777  htmlFileS << "<title> Remote Monitoring Tool Global</title>" << std::endl;
2778  htmlFileS << "<style type=\"text/css\">" << std::endl;
2779  htmlFileS << " body,td{ background-color: #FFFFCC; font-family: arial, arial ce, helvetica; font-size: 12px; }"
2780  << std::endl;
2781  htmlFileS << " td.s0 { font-family: arial, arial ce, helvetica; }" << std::endl;
2782  htmlFileS << " td.s1 { font-family: arial, arial ce, helvetica; font-weight: bold; background-color: #FFC169; "
2783  "text-align: center;}"
2784  << std::endl;
2785  htmlFileS << " td.s2 { font-family: arial, arial ce, helvetica; background-color: #eeeeee; }" << std::endl;
2786  htmlFileS << " td.s3 { font-family: arial, arial ce, helvetica; background-color: #d0d0d0; }" << std::endl;
2787  htmlFileS << " td.s4 { font-family: arial, arial ce, helvetica; background-color: #FFC169; }" << std::endl;
2788  htmlFileS << "</style>" << std::endl;
2789  htmlFileS << "<body>" << std::endl;
2790 
2791  if (sub == 1)
2792  htmlFileS << "<h1> ADC Shape for HB, RUN = " << runnumber << " </h1>" << std::endl;
2793  if (sub == 2)
2794  htmlFileS << "<h1> ADC Shape for HE, RUN = " << runnumber << " </h1>" << std::endl;
2795  if (sub == 3)
2796  htmlFileS << "<h1> ADC Shape for HO, RUN = " << runnumber << " </h1>" << std::endl;
2797  if (sub == 4)
2798  htmlFileS << "<h1> ADC Shape for HF, RUN = " << runnumber << " </h1>" << std::endl;
2799  htmlFileP << "<br>" << std::endl;
2800 
2801  htmlFileS << "<h2> 1.Mean ADC Shape </h3>" << std::endl;
2802  htmlFileS << "<h3> 1.A. ADC shape averaged over all good channels, depth and events.</h3>" << std::endl;
2803  // htmlFileS << "<h4> Legend: Bins less "<<Pedest[0][sub]<<" correpond to bad Pedestals </h4>"<< std::endl;
2804  if (sub == 1)
2805  htmlFileS << " <img src=\"HistGoodTSshapesHB.png\" />" << std::endl;
2806  if (sub == 2)
2807  htmlFileS << " <img src=\"HistGoodTSshapesHE.png\" />" << std::endl;
2808  if (sub == 3)
2809  htmlFileS << " <img src=\"HistGoodTSshapesHO.png\" />" << std::endl;
2810  if (sub == 4)
2811  htmlFileS << " <img src=\"HistGoodTSshapesHF.png\" />" << std::endl;
2812  htmlFileS << "<br>" << std::endl;
2813  htmlFileS << "<h3> 1.B. ADC shape averaged over all bad channels, depth and events. Bad channels are selected by 5 "
2814  "criteria: CapId, A, W, P, Pw</h3>"
2815  << std::endl;
2816  // htmlFileS << "<h4> Channel legend: white - good, other colour - bad. </h4>"<< std::endl;
2817  if (sub == 1)
2818  htmlFileS << " <img src=\"HistBadTSshapesHB.png\" />" << std::endl;
2819  if (sub == 2)
2820  htmlFileS << " <img src=\"HistBadTSshapesHE.png\" />" << std::endl;
2821  if (sub == 3)
2822  htmlFileS << " <img src=\"HistBadTSshapesHO.png\" />" << std::endl;
2823  if (sub == 4)
2824  htmlFileS << " <img src=\"HistBadTSshapesHF.png\" />" << std::endl;
2825 
2826  htmlFileS << "<h2> 2. Pattern of channels for Sub-Detector over depth,eta.phi </h3>" << std::endl;
2827  htmlFileS << "<h3> 2.A. reminder:.......................... for HBHE, TS=2;...................................... "
2828  "for HF, TS=1;..................................... for HO, TS=0,1,2 </h3>"
2829  << std::endl;
2830  if (sub == 1)
2831  htmlFileS << " <img src=\"Hist_mapDepthAllTS2_HB.png\" />" << std::endl;
2832  if (sub == 2)
2833  htmlFileS << " <img src=\"Hist_mapDepthAllTS2_HE.png\" />" << std::endl;
2834  if (sub == 3)
2835  htmlFileS << " <img src=\"Hist_mapDepthAllTS012_HO.png\" />" << std::endl;
2836  if (sub == 4)
2837  htmlFileS << " <img src=\"Hist_mapDepthAllTS1_HF.png\" />" << std::endl;
2838  htmlFileS << "<br>" << std::endl;
2839 
2840  htmlFileS.close();
2841  } // end sub
2842 
2843  //======================================================================
2844 
2845  //======================================================================
2846  // Creating subdet html pages:
2847 
2848  for (int sub = 1; sub <= 4; sub++) { //Subdetector: 1-HB, 2-HE, 3-HO, 4-HF
2849  ofstream htmlFile;
2850  if (sub == 1)
2851  htmlFile.open("HB.html");
2852  if (sub == 2)
2853  htmlFile.open("HE.html");
2854  if (sub == 3)
2855  htmlFile.open("HO.html");
2856  if (sub == 4)
2857  htmlFile.open("HF.html");
2858 
2859  htmlFile << "</html><html xmlns=\"http://www.w3.org/1999/xhtml\">" << std::endl;
2860  htmlFile << "<head>" << std::endl;
2861  htmlFile << "<meta http-equiv=\"Content-Type\" content=\"text/html\"/>" << std::endl;
2862  htmlFile << "<title> Remote Monitoring Tool </title>" << std::endl;
2863  htmlFile << "<style type=\"text/css\">" << std::endl;
2864  htmlFile << " body,td{ background-color: #FFFFCC; font-family: arial, arial ce, helvetica; font-size: 12px; }"
2865  << std::endl;
2866  htmlFile << " td.s0 { font-family: arial, arial ce, helvetica; }" << std::endl;
2867  htmlFile << " td.s1 { font-family: arial, arial ce, helvetica; font-weight: bold; background-color: #FFC169; "
2868  "text-align: center;}"
2869  << std::endl;
2870  htmlFile << " td.s2 { font-family: arial, arial ce, helvetica; background-color: #eeeeee; }" << std::endl;
2871  htmlFile << " td.s3 { font-family: arial, arial ce, helvetica; background-color: #d0d0d0; }" << std::endl;
2872  htmlFile << " td.s4 { font-family: arial, arial ce, helvetica; background-color: #FFC169; }" << std::endl;
2873  htmlFile << " td.s5 { font-family: arial, arial ce, helvetica; background-color: #FF00FF; }" << std::endl;
2874  htmlFile << " td.s6 { font-family: arial, arial ce, helvetica; background-color: #9ACD32; }" << std::endl;
2875  htmlFile << " td.s7 { font-family: arial, arial ce, helvetica; background-color: #32CD32; }" << std::endl;
2876  htmlFile << " td.s8 { font-family: arial, arial ce, helvetica; background-color: #00FFFF; }" << std::endl;
2877  htmlFile << " td.s9 { font-family: arial, arial ce, helvetica; background-color: #FFE4E1; }" << std::endl;
2878  htmlFile << " td.s10 { font-family: arial, arial ce, helvetica; background-color: #A0522D; }" << std::endl;
2879  htmlFile << " td.s11 { font-family: arial, arial ce, helvetica; background-color: #1E90FF; }" << std::endl;
2880  htmlFile << " td.s12 { font-family: arial, arial ce, helvetica; background-color: #00BFFF; }" << std::endl;
2881  htmlFile << " td.s13 { font-family: arial, arial ce, helvetica; background-color: #FFFF00; }" << std::endl;
2882  htmlFile << " td.s14 { font-family: arial, arial ce, helvetica; background-color: #B8860B; }" << std::endl;
2883  htmlFile << "</style>" << std::endl;
2884  htmlFile << "<body>" << std::endl;
2885  if (sub == 1)
2886  htmlFile << "<h1> HCAL BARREL, RUN = " << runnumber << " </h1>" << std::endl;
2887  if (sub == 2)
2888  htmlFile << "<h1> HCAL ENDCAP, RUN = " << runnumber << " </h1>" << std::endl;
2889  if (sub == 3)
2890  htmlFile << "<h1> HCAL OUTER, RUN = " << runnumber << " </h1>" << std::endl;
2891  if (sub == 4)
2892  htmlFile << "<h1> HCAL FORWARD, RUN = " << runnumber << " </h1>" << std::endl;
2893  htmlFile << "<br>" << std::endl;
2894 
2895  htmlFile << "<a name=\"Top\"></a>\n";
2896  htmlFile << "<b>Contents:<br>\n";
2897  htmlFile << "1. <a href=\"#AnalysisResults\">Analysis results</a><br>\n";
2898  htmlFile << "2. <a href=\"#Status\">Status</a><br>\n";
2899  htmlFile << "2A. <a href=\"#ChannelMap\">Channel map</a><br>\n";
2900  // htmlFile << "2B. <a href=\"#BadChannels\">List of bad channels</a><br>\n";
2901  // htmlFile << "2C. <a href=\"#BadPedestals\">List of channels with bad pedestals</a><br>\n";
2902 
2903  htmlFile << "<a name=\"AnalysisResults\"></a>\n";
2904  if (sub == 1)
2905  htmlFile << "<h2> 1. Analysis results for HB</h2>" << std::endl;
2906  if (sub == 2)
2907  htmlFile << "<h2> 1. Analysis results for HE</h2>" << std::endl;
2908  if (sub == 3)
2909  htmlFile << "<h2> 1. Analysis results for HO</h2>" << std::endl;
2910  if (sub == 4)
2911  htmlFile << "<h2> 1. Analysis results for HF</h2>" << std::endl;
2912  htmlFile << "<table width=\"400\">" << std::endl;
2913  htmlFile << "<tr>" << std::endl;
2914 
2915  if (sub == 1) {
2916  // AZ 12.03.2019
2917  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
2918  << runnumber << "/HB_Tile.html\">Megatile Channels</a></td>" << std::endl;
2919  // htmlFile << " <td><a href=\"HB_Tile.html\">Megatile Channels</a></td>"<< std::endl;
2920 
2921  // htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"<<runnumber<<"/HB_Calib.html\">Calibration Channels</a></td>"<< std::endl;
2922  // htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"<<runnumber<<"/HB_Drift.html\">Gain Stability</a></td>"<< std::endl;
2923  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
2924  << runnumber << "/HB_Pedestals.html\">Pedestals</a></td>" << std::endl;
2925  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
2926  << runnumber << "/HB_Shapes.html\">ADC Shapes</a></td>" << std::endl;
2927  }
2928  if (sub == 2) {
2929  // AZ 12.03.2019
2930  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
2931  << runnumber << "/HE_Tile.html\">Megatile Channels</a></td>" << std::endl;
2932  // htmlFile << " <td><a href=\"HE_Tile.html\">Megatile Channels</a></td>"<< std::endl;
2933 
2934  // htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"<<runnumber<<"/HE_Calib.html\">Calibration Channels</a></td>"<< std::endl;
2935  // htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"<<runnumber<<"/HE_Drift.html\">Gain Stability</a></td>"<< std::endl;
2936  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
2937  << runnumber << "/HE_Pedestals.html\">Pedestals</a></td>" << std::endl;
2938  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
2939  << runnumber << "/HE_Shapes.html\">ADC Shapes</a></td>" << std::endl;
2940  }
2941  if (sub == 3) {
2942  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
2943  << runnumber << "/HO_Tile.html\">Megatile Channels</a></td>" << std::endl;
2944  // htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"<<runnumber<<"/HO_Calib.html\">Calibration Channels</a></td>"<< std::endl;
2945  // htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"<<runnumber<<"/HO_Drift.html\">Gain Stability</a></td>"<< std::endl;
2946  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
2947  << runnumber << "/HO_Pedestals.html\">Pedestals</a></td>" << std::endl;
2948  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
2949  << runnumber << "/HO_Shapes.html\">ADC Shapes</a></td>" << std::endl;
2950  }
2951  if (sub == 4) {
2952  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
2953  << runnumber << "/HF_Tile.html\">Megatile Channels</a></td>" << std::endl;
2954  // htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"<<runnumber<<"/HF_Calib.html\">Calibration Channels</a></td>"<< std::endl;
2955  // htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"<<runnumber<<"/HF_Drift.html\">Gain Stability</a></td>"<< std::endl;
2956  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
2957  << runnumber << "/HF_Pedestals.html\">Pedestals</a></td>" << std::endl;
2958  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
2959  << runnumber << "/HF_Shapes.html\">ADC Shapes</a></td>" << std::endl;
2960  }
2961 
2962  htmlFile << "</tr>" << std::endl;
2963  htmlFile << "</table>" << std::endl;
2964  htmlFile << "<br>" << std::endl;
2965 
2966  //AZ2023:
2967  /*
2968  htmlFile << "<a name=\"Status\"></a>\n";
2969  if (sub == 1)
2970  htmlFile << "<h2> 2.Status HB over all criteria </h2>" << std::endl;
2971  if (sub == 2)
2972  htmlFile << "<h2> 2.Status HE over all criteria </h2>" << std::endl;
2973  if (sub == 3)
2974  htmlFile << "<h2> 2.Status HO over all criteria </h2>" << std::endl;
2975  if (sub == 4)
2976  htmlFile << "<h2> 2.Status HF over all criteria </h2>" << std::endl;
2977 
2978  htmlFile << "<a name=\"ChannelMap\"></a>\n";
2979  htmlFile << "<h3> 2.A.Channel map for each Depth </h3>" << std::endl;
2980  htmlFile << "<h4> Channel legend: yellow - good, white - "
2981  "not applicable or out of range </h4>"
2982  << std::endl;
2983  if (sub == 1)
2984  htmlFile << " <img src=\"MAPHB.png\" />" << std::endl;
2985  if (sub == 2)
2986  htmlFile << " <img src=\"MAPHE.png\" />" << std::endl;
2987  if (sub == 3)
2988  htmlFile << " <img src=\"MAPHO.png\" />" << std::endl;
2989  if (sub == 4)
2990  htmlFile << " <img src=\"MAPHF.png\" />" << std::endl;
2991  htmlFile << "<br>" << std::endl;
2992  htmlFile << "<a href=\"#Top\">to top</a><br>\n";
2993 */
2995  /*
2996 // htmlFile << "<h3> 2.B.List of Bad channels (rate > 0.1) and its rates for each RMT criteria (for GS - %) </h3>"<< std::endl;
2997 
2998  htmlFile << "<a name=\"BadChannels\"></a>\n";
2999  htmlFile << "<h3> 2.B.List of Bad channels (rate > 0.1) and its rates for each RMT criteria </h3>"<< std::endl;
3000 
3001  //htmlFile << " <td><a href=\"HELP.html\"> Description of criteria for bad channel selection</a></td>"<< std::endl;
3002  // htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"<<runnumber<<"/HELP.html\"> Description of criteria for bad channel selection</a></td>"<< std::endl;
3003  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/HELP.html\"> Description of criteria for bad channel selection</a></td>"<< std::endl;
3004 
3005  htmlFile << "<table>"<< std::endl;
3006  htmlFile << "<tr>";
3007  htmlFile << "<td class=\"s4\" align=\"center\">#</td>" << std::endl;
3008  htmlFile << "<td class=\"s1\" align=\"center\">ETA</td>" << std::endl;
3009  htmlFile << "<td class=\"s1\" align=\"center\">PHI</td>" << std::endl;
3010  htmlFile << "<td class=\"s1\" align=\"center\">DEPTH</td>"<< std::endl;
3011  htmlFile << "<td class=\"s1\" align=\"center\">RBX</td>" << std::endl;
3012  htmlFile << "<td class=\"s1\" align=\"center\">RM</td>" << std::endl;
3013  htmlFile << "<td class=\"s1\" align=\"center\">PIXEL</td>" << std::endl;
3014  htmlFile << "<td class=\"s1\" align=\"center\">RM_FIBER</td>" << std::endl;
3015  htmlFile << "<td class=\"s1\" align=\"center\">FIBER_CH</td>" << std::endl;
3016  htmlFile << "<td class=\"s1\" align=\"center\">QIE</td>" << std::endl;
3017  htmlFile << "<td class=\"s1\" align=\"center\">ADC</td>" << std::endl;
3018  htmlFile << "<td class=\"s1\" align=\"center\">CRATE</td>" << std::endl;
3019  htmlFile << "<td class=\"s1\" align=\"center\">DCC</td>" << std::endl;
3020  htmlFile << "<td class=\"s1\" align=\"center\">SPIGOT</td>" << std::endl;
3021  htmlFile << "<td class=\"s1\" align=\"center\">HTR_FIBER</td>" << std::endl;
3022  htmlFile << "<td class=\"s1\" align=\"center\">HTR_SLOT</td>" << std::endl;
3023  htmlFile << "<td class=\"s1\" align=\"center\">HTR_FPGA</td>" << std::endl;
3024  htmlFile << "<td class=\"s5\" align=\"center\">Cm</td>" << std::endl;
3025  htmlFile << "<td class=\"s5\" align=\"center\">Am</td>" << std::endl;
3026  htmlFile << "<td class=\"s5\" align=\"center\">Wm</td>" << std::endl;
3027  htmlFile << "<td class=\"s5\" align=\"center\">Rm</td>" << std::endl;
3028  htmlFile << "<td class=\"s5\" align=\"center\">TNm</td>" << std::endl;
3029  htmlFile << "<td class=\"s5\" align=\"center\">TXm</td>" << std::endl;
3030 // htmlFile << "<td class=\"s8\" align=\"center\">Cc</td>" << std::endl;
3031 // htmlFile << "<td class=\"s8\" align=\"center\">Ac</td>" << std::endl;
3032 // htmlFile << "<td class=\"s8\" align=\"center\">Wc</td>" << std::endl;
3033 // htmlFile << "<td class=\"s8\" align=\"center\">Rc</td>" << std::endl;
3034 // htmlFile << "<td class=\"s8\" align=\"center\">TNc</td>" << std::endl;
3035 // htmlFile << "<td class=\"s8\" align=\"center\">TXc</td>" << std::endl;
3036 // htmlFile << "<td class=\"s9\" align=\"center\">GS (%)</td>" << std::endl;
3037  htmlFile << "<td class=\"s4\" align=\"center\">Pm</td>" << std::endl;
3038  htmlFile << "<td class=\"s4\" align=\"center\">pWm</td>" << std::endl;
3039  htmlFile << "</tr>" << std::endl;
3040 
3041  for (int i=1;i<=NBad;i++) {
3042  if((ind%2)==1){
3043  raw_class="<td class=\"s2\" align=\"center\">";
3044  raw_class1="<td class=\"s6\" align=\"center\">";
3045  raw_class2="<td class=\"s11\" align=\"center\">";
3046  raw_class3="<td class=\"s13\" align=\"center\">";
3047 
3048  }else{
3049  raw_class="<td class=\"s3\" align=\"center\">";
3050  raw_class1="<td class=\"s7\" align=\"center\">";
3051  raw_class2="<td class=\"s12\" align=\"center\">";
3052  raw_class3="<td class=\"s14\" align=\"center\">";
3053  }
3054  const CellDB db;
3055  CellDB ce;
3056  if ((ce.size()>=1)&&(Sub[2][i]==sub)) {
3057  // AZ 19.03.2018
3058 
3059 // AZ 19 if (Sub[2][i]==1) {
3060 // AZ 19 ce = db.find("subdet", "HB").find("Eta", Eta[2][i]).find("Phi", Phi[2][i]).find("Depth", Depth[2][i]);
3061 // AZ 19 if (ce.size()==0) {cout<<"Error: No such HB, Eta="<< Eta[2][i] <<", Phi="<< Phi[2][i] <<", Depth="<< Depth[2][i] <<" in database"<<endl; continue;}
3062 // AZ 19 else if (ce.size()>1) {cout<<"Warning: More than one line correspond to such HB, Eta="<< Eta[2][i] <<", Phi="<< Phi[2][i] <<", Depth="<< Depth[2][i] <<" in database"<<endl;}
3063 // AZ 19 }
3064 // AZ 19 if (Sub[2][i]==2) {
3065 // AZ 19 ce = db.find("subdet", "HE").find("Eta", Eta[2][i]).find("Phi", Phi[2][i]).find("Depth", Depth[2][i]);
3066 // AZ 19 if (ce.size()==0) {cout<<"Error: No such HE, Eta="<< Eta[2][i] <<", Phi="<< Phi[2][i] <<", Depth="<< Depth[2][i] <<" in database"<<endl;continue;}
3067 // AZ 19 else if (ce.size()>1) {cout<<"Warning: More than one line correspond to such HE, Eta="<< Eta[2][i] <<", Phi="<< Phi[2][i] <<", Depth="<< Depth[2][i] <<" in database"<<endl;}
3068 // AZ 19 }
3069 // AZ 19 if (Sub[2][i]==3) {
3070 // AZ 19 ce = db.find("subdet", "HO").find("Eta", Eta[2][i]).find("Phi", Phi[2][i]).find("Depth", Depth[2][i]);
3071 // AZ 19 if (ce.size()==0) {cout<<"Error: No such HO, Eta="<< Eta[2][i] <<", Phi="<< Phi[2][i] <<", Depth="<< Depth[2][i] <<" in database"<<endl;continue;}
3072 // AZ 19 else if (ce.size()>1) {cout<<"Warning: More than one line correspond to such HO, Eta="<< Eta[2][i] <<", Phi="<< Phi[2][i] <<", Depth="<< Depth[2][i] <<" in database"<<endl;}
3073 // AZ 19 }
3074 // AZ 19 if (Sub[2][i]==4) {
3075 // AZ 19 ce = db.find("subdet", "HF").find("Eta", Eta[2][i]).find("Phi", Phi[2][i]).find("Depth", Depth[2][i]);
3076 // AZ 19 if (ce.size()==0) {cout<<"Error: No such HF, Eta="<< Eta[2][i] <<", Phi="<< Phi[2][i] <<", Depth="<< Depth[2][i] <<" in database"<<endl;continue;}
3077 // AZ 19 else if (ce.size()>1) {cout<<"Warning: More than one line correspond to such HF, Eta="<< Eta[2][i] <<", Phi="<< Phi[2][i] <<", Depth="<< Depth[2][i] <<" in database"<<endl;}
3078 // AZ 19 }
3079 
3080  htmlFile << "<tr>"<< std::endl;
3081  htmlFile << "<td class=\"s4\" align=\"center\">" << ind+1 <<"</td>"<< std::endl;
3082  htmlFile << raw_class<< Eta[2][i]<<"</td>"<< std::endl;
3083  htmlFile << raw_class<< Phi[2][i]<<"</td>"<< std::endl;
3084  htmlFile << raw_class<< Depth[2][i] <<"</td>"<< std::endl;
3085  htmlFile << raw_class<< ce[0].RBX <<"</td>"<< std::endl;
3086  htmlFile << raw_class<< ce[0].RM <<"</td>"<< std::endl;
3087  htmlFile << raw_class<< ce[0].Pixel <<"</td>"<< std::endl;
3088  htmlFile << raw_class<< ce[0].RMfiber <<"</td>"<< std::endl;
3089  htmlFile << raw_class<< ce[0].FiberCh <<"</td>"<< std::endl;
3090  htmlFile << raw_class<< ce[0].QIE <<"</td>"<< std::endl;
3091  htmlFile << raw_class<< ce[0].ADC<<"</td>"<< std::endl;
3092  htmlFile << raw_class<< ce[0].VMECardID <<"</td>"<< std::endl;
3093  htmlFile << raw_class<< ce[0].dccID <<"</td>"<< std::endl;
3094  htmlFile << raw_class<< ce[0].Spigot <<"</td>"<< std::endl;
3095  htmlFile << raw_class<< ce[0].FiberIndex <<"</td>"<< std::endl;
3096  htmlFile << raw_class<< ce[0].HtrSlot <<"</td>"<< std::endl;
3097  htmlFile << raw_class<< ce[0].HtrTB <<"</td>"<< std::endl;
3098  htmlFile << raw_class1<< Map_Ampl[1][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3099  htmlFile << raw_class1<< Map_Ampl[2][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3100  htmlFile << raw_class1<< Map_Ampl[3][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3101  htmlFile << raw_class1<< Map_Ampl[4][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3102  htmlFile << raw_class1<< Map_Ampl[5][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3103  htmlFile << raw_class1<< Map_Ampl[6][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3104 // htmlFile << raw_class2<< Map_Ampl[11][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3105 // htmlFile << raw_class2<< Map_Ampl[12][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3106 // htmlFile << raw_class2<< Map_Ampl[13][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3107 // htmlFile << raw_class2<< Map_Ampl[14][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3108 // htmlFile << raw_class2<< Map_Ampl[15][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3109 // htmlFile << raw_class2<< Map_Ampl[16][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3110 // htmlFile << raw_class3<< Map_Ampl[21][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3111  htmlFile << raw_class<< Map_Ampl[31][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3112  htmlFile << raw_class<< Map_Ampl[32][Sub[2][i]][Depth[2][i]]->GetBinContent(Eta[2][i]+41,Phi[2][i]+1)<<"</td>"<< std::endl;
3113  htmlFile << "</tr>" << std::endl;
3114 
3115  ind+=1;
3116  }
3117  }
3118  htmlFile << "</table>" << std::endl;
3119  htmlFile << "<br>"<< std::endl;
3120  htmlFile << "<a href=\"#Top\">to top</a><br>\n";
3121 
3122 
3123  htmlFile << "<h3> 2.C.List of Gain unstable channels and its value in % (for other criterias - rate)</h3>"<< std::endl;
3124  htmlFile << "<table>"<< std::endl;
3125  htmlFile << "<tr>";
3126  htmlFile << "<td class=\"s4\" align=\"center\">#</td>" << std::endl;
3127  htmlFile << "<td class=\"s1\" align=\"center\">ETA</td>" << std::endl;
3128  htmlFile << "<td class=\"s1\" align=\"center\">PHI</td>" << std::endl;
3129  htmlFile << "<td class=\"s1\" align=\"center\">DEPTH</td>"<< std::endl;
3130  htmlFile << "<td class=\"s1\" align=\"center\">RBX</td>" << std::endl;
3131  htmlFile << "<td class=\"s1\" align=\"center\">RM</td>" << std::endl;
3132  htmlFile << "<td class=\"s1\" align=\"center\">PIXEL</td>" << std::endl;
3133  htmlFile << "<td class=\"s1\" align=\"center\">RM_FIBER</td>" << std::endl;
3134  htmlFile << "<td class=\"s1\" align=\"center\">FIBER_CH</td>" << std::endl;
3135  htmlFile << "<td class=\"s1\" align=\"center\">QIE</td>" << std::endl;
3136  htmlFile << "<td class=\"s1\" align=\"center\">ADC</td>" << std::endl;
3137  htmlFile << "<td class=\"s1\" align=\"center\">CRATE</td>" << std::endl;
3138  htmlFile << "<td class=\"s1\" align=\"center\">DCC</td>" << std::endl;
3139  htmlFile << "<td class=\"s1\" align=\"center\">SPIGOT</td>" << std::endl;
3140  htmlFile << "<td class=\"s1\" align=\"center\">HTR_FIBER</td>" << std::endl;
3141  htmlFile << "<td class=\"s1\" align=\"center\">HTR_SLOT</td>" << std::endl;
3142  htmlFile << "<td class=\"s1\" align=\"center\">HTR_FPGA</td>" << std::endl;
3143  htmlFile << "<td class=\"s5\" align=\"center\">Cm</td>" << std::endl;
3144  htmlFile << "<td class=\"s5\" align=\"center\">Am</td>" << std::endl;
3145  htmlFile << "<td class=\"s5\" align=\"center\">Wm</td>" << std::endl;
3146  htmlFile << "<td class=\"s5\" align=\"center\">Rm</td>" << std::endl;
3147  htmlFile << "<td class=\"s5\" align=\"center\">TNm</td>" << std::endl;
3148  htmlFile << "<td class=\"s5\" align=\"center\">TXm</td>" << std::endl;
3149 // htmlFile << "<td class=\"s8\" align=\"center\">Cc</td>" << std::endl;
3150 // htmlFile << "<td class=\"s8\" align=\"center\">Ac</td>" << std::endl;
3151 // htmlFile << "<td class=\"s8\" align=\"center\">Wc</td>" << std::endl;
3152 // htmlFile << "<td class=\"s8\" align=\"center\">Rc</td>" << std::endl;
3153 // htmlFile << "<td class=\"s8\" align=\"center\">TNc</td>" << std::endl;
3154 // htmlFile << "<td class=\"s8\" align=\"center\">TXc</td>" << std::endl;
3155 // htmlFile << "<td class=\"s9\" align=\"center\">GS(%)</td>" << std::endl;
3156  htmlFile << "<td class=\"s4\" align=\"center\">Pm</td>" << std::endl;
3157  htmlFile << "<td class=\"s4\" align=\"center\">pWm</td>" << std::endl;
3158  htmlFile << "</tr>" << std::endl;
3159 
3160  for (int i=1;i<=NWarn;i++) {
3161  if((ind%2)==1){
3162  raw_class="<td class=\"s2\" align=\"center\">";
3163  raw_class1="<td class=\"s6\" align=\"center\">";
3164  raw_class2="<td class=\"s11\" align=\"center\">";
3165  raw_class3="<td class=\"s13\" align=\"center\">";
3166 
3167  }else{
3168  raw_class="<td class=\"s3\" align=\"center\">";
3169  raw_class1="<td class=\"s7\" align=\"center\">";
3170  raw_class2="<td class=\"s12\" align=\"center\">";
3171  raw_class3="<td class=\"s14\" align=\"center\">";
3172  }
3173  const CellDB db;
3174  CellDB ce;
3175  if ((ce.size()>=1)&&(Sub[1][i]==sub)) {
3176  if (Sub[1][i]==1) {
3177  ce = db.find("subdet", "HB").find("Eta", Eta[1][i]).find("Phi", Phi[1][i]).find("Depth", Depth[1][i]);
3178  if (ce.size()==0) {cout<<"Error: No such HB, Eta="<< Eta[1][i] <<", Phi="<< Phi[1][i] <<", Depth="<< Depth[1][i] <<" in database"<<endl;}
3179  else if (ce.size()>1) {cout<<"Warning: More than one line correspond to such HB, Eta="<< Eta[1][i] <<", Phi="<< Phi[1][i] <<", Depth="<< Depth[1][i] <<" in database"<<endl;}
3180  }
3181  if (Sub[1][i]==2) {
3182  ce = db.find("subdet", "HE").find("Eta", Eta[1][i]).find("Phi", Phi[1][i]).find("Depth", Depth[1][i]);
3183  if (ce.size()==0) {cout<<"Error: No such HE, Eta="<< Eta[1][i] <<", Phi="<< Phi[1][i] <<", Depth="<< Depth[1][i] <<" in database"<<endl;}
3184  else if (ce.size()>1) {cout<<"Warning: More than one line correspond to such HE, Eta="<< Eta[1][i] <<", Phi="<< Phi[1][i] <<", Depth="<< Depth[1][i] <<" in database"<<endl;}
3185  }
3186  if (Sub[1][i]==3) {
3187  ce = db.find("subdet", "HO").find("Eta", Eta[1][i]).find("Phi", Phi[1][i]).find("Depth", Depth[1][i]);
3188  if (ce.size()==0) {cout<<"Error: No such HO, Eta="<< Eta[1][i] <<", Phi="<< Phi[1][i] <<", Depth="<< Depth[1][i] <<" in database"<<endl;}
3189  else if (ce.size()>1) {cout<<"Warning: More than one line correspond to such HO, Eta="<< Eta[1][i] <<", Phi="<< Phi[1][i] <<", Depth="<< Depth[1][i] <<" in database"<<endl;}
3190  }
3191  if (Sub[1][i]==4) {
3192  ce = db.find("subdet", "HF").find("Eta", Eta[1][i]).find("Phi", Phi[1][i]).find("Depth", Depth[1][i]);
3193  if (ce.size()==0) {cout<<"Error: No such HF, Eta="<< Eta[1][i] <<", Phi="<< Phi[1][i] <<", Depth="<< Depth[1][i] <<" in database"<<endl;}
3194  else if (ce.size()>1) {cout<<"Warning: More than one line correspond to such HF, Eta="<< Eta[1][i] <<", Phi="<< Phi[1][i] <<", Depth="<< Depth[1][i] <<" in database"<<endl;}
3195  }
3196  htmlFile << "<td class=\"s4\" align=\"center\">" << ind+1 <<"</td>"<< std::endl;
3197  htmlFile << raw_class<< Eta[1][i]<<"</td>"<< std::endl;
3198  htmlFile << raw_class<< Phi[1][i]<<"</td>"<< std::endl;
3199  htmlFile << raw_class<< Depth[1][i] <<"</td>"<< std::endl;
3200  htmlFile << raw_class<< ce[0].RBX <<"</td>"<< std::endl;
3201  htmlFile << raw_class<< ce[0].RM <<"</td>"<< std::endl;
3202  htmlFile << raw_class<< ce[0].Pixel <<"</td>"<< std::endl;
3203  htmlFile << raw_class<< ce[0].RMfiber <<"</td>"<< std::endl;
3204  htmlFile << raw_class<< ce[0].FiberCh <<"</td>"<< std::endl;
3205  htmlFile << raw_class<< ce[0].QIE <<"</td>"<< std::endl;
3206  htmlFile << raw_class<< ce[0].ADC<<"</td>"<< std::endl;
3207  htmlFile << raw_class<< ce[0].VMECardID <<"</td>"<< std::endl;
3208  htmlFile << raw_class<< ce[0].dccID <<"</td>"<< std::endl;
3209  htmlFile << raw_class<< ce[0].Spigot <<"</td>"<< std::endl;
3210  htmlFile << raw_class<< ce[0].FiberIndex <<"</td>"<< std::endl;
3211  htmlFile << raw_class<< ce[0].HtrSlot <<"</td>"<< std::endl;
3212  htmlFile << raw_class<< ce[0].HtrTB <<"</td>"<< std::endl;
3213  htmlFile << raw_class1<< Map_Ampl[1][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3214  htmlFile << raw_class1<< Map_Ampl[2][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3215  htmlFile << raw_class1<< Map_Ampl[3][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3216  htmlFile << raw_class1<< Map_Ampl[4][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3217  htmlFile << raw_class1<< Map_Ampl[5][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3218  htmlFile << raw_class1<< Map_Ampl[6][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3219 // htmlFile << raw_class2<< Map_Ampl[11][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3220 // htmlFile << raw_class2<< Map_Ampl[12][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3221 // htmlFile << raw_class2<< Map_Ampl[13][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3222 // htmlFile << raw_class2<< Map_Ampl[14][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3223 // htmlFile << raw_class2<< Map_Ampl[15][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3224 // htmlFile << raw_class2<< Map_Ampl[16][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3225 // htmlFile << raw_class3<< Map_Ampl[21][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3226  htmlFile << raw_class<< Map_Ampl[31][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3227  htmlFile << raw_class<< Map_Ampl[32][Sub[1][i]][Depth[1][i]]->GetBinContent(Eta[1][i]+41,Phi[1][i]+1)<<"</td>"<< std::endl;
3228  htmlFile << "</tr>" << std::endl;
3229  htmlFile << "</tr>" << std::endl;
3230  ind+=1;
3231  }
3232  }
3233  htmlFile << "</table>" << std::endl;
3234  htmlFile << "<br>"<< std::endl;
3235 
3236 
3237 // htmlFile << "<h3> 2.D.List of channels with Bad Pedestals (rate > 0.1) and its rates (for GS - %)</h3>"<< std::endl;
3238  htmlFile << "<a name=\"BadPedestals\"></a>\n";
3239  htmlFile << "<h3> 2.C.List of channels with Bad Pedestals (rate > 0.1) and its rates </h3>"<< std::endl;
3240  htmlFile << "<table>"<< std::endl;
3241  htmlFile << "<tr>";
3242  htmlFile << "<td class=\"s4\" align=\"center\">#</td>" << std::endl;
3243  htmlFile << "<td class=\"s1\" align=\"center\">ETA</td>" << std::endl;
3244  htmlFile << "<td class=\"s1\" align=\"center\">PHI</td>" << std::endl;
3245  htmlFile << "<td class=\"s1\" align=\"center\">DEPTH</td>"<< std::endl;
3246  htmlFile << "<td class=\"s1\" align=\"center\">RBX</td>" << std::endl;
3247  htmlFile << "<td class=\"s1\" align=\"center\">RM</td>" << std::endl;
3248  htmlFile << "<td class=\"s1\" align=\"center\">PIXEL</td>" << std::endl;
3249  htmlFile << "<td class=\"s1\" align=\"center\">RM_FIBER</td>" << std::endl;
3250  htmlFile << "<td class=\"s1\" align=\"center\">FIBER_CH</td>" << std::endl;
3251  htmlFile << "<td class=\"s1\" align=\"center\">QIE</td>" << std::endl;
3252  htmlFile << "<td class=\"s1\" align=\"center\">ADC</td>" << std::endl;
3253  htmlFile << "<td class=\"s1\" align=\"center\">CRATE</td>" << std::endl;
3254  htmlFile << "<td class=\"s1\" align=\"center\">DCC</td>" << std::endl;
3255  htmlFile << "<td class=\"s1\" align=\"center\">SPIGOT</td>" << std::endl;
3256  htmlFile << "<td class=\"s1\" align=\"center\">HTR_FIBER</td>" << std::endl;
3257  htmlFile << "<td class=\"s1\" align=\"center\">HTR_SLOT</td>" << std::endl;
3258  htmlFile << "<td class=\"s1\" align=\"center\">HTR_FPGA</td>" << std::endl;
3259  htmlFile << "<td class=\"s5\" align=\"center\">Cm</td>" << std::endl;
3260  htmlFile << "<td class=\"s5\" align=\"center\">Am</td>" << std::endl;
3261  htmlFile << "<td class=\"s5\" align=\"center\">Wm</td>" << std::endl;
3262  htmlFile << "<td class=\"s5\" align=\"center\">Rm</td>" << std::endl;
3263  htmlFile << "<td class=\"s5\" align=\"center\">TNm</td>" << std::endl;
3264  htmlFile << "<td class=\"s5\" align=\"center\">TXm</td>" << std::endl;
3265 // htmlFile << "<td class=\"s8\" align=\"center\">Cc</td>" << std::endl;
3266 // htmlFile << "<td class=\"s8\" align=\"center\">Ac</td>" << std::endl;
3267 // htmlFile << "<td class=\"s8\" align=\"center\">Wc</td>" << std::endl;
3268 // htmlFile << "<td class=\"s8\" align=\"center\">Rc</td>" << std::endl;
3269 // htmlFile << "<td class=\"s8\" align=\"center\">TNc</td>" << std::endl;
3270 // htmlFile << "<td class=\"s8\" align=\"center\">TXc</td>" << std::endl;
3271 // htmlFile << "<td class=\"s9\" align=\"center\">GS(%)</td>" << std::endl;
3272  htmlFile << "<td class=\"s4\" align=\"center\">Pm</td>" << std::endl;
3273  htmlFile << "<td class=\"s4\" align=\"center\">pWm</td>" << std::endl;
3274  htmlFile << "</tr>" << std::endl;
3275 
3276  for (int i=1;i<=NPed;i++) {
3277  if((ind%2)==1){
3278  raw_class="<td class=\"s2\" align=\"center\">";
3279  raw_class1="<td class=\"s6\" align=\"center\">";
3280  raw_class2="<td class=\"s11\" align=\"center\">";
3281  raw_class3="<td class=\"s13\" align=\"center\">";
3282 
3283  }else{
3284  raw_class="<td class=\"s3\" align=\"center\">";
3285  raw_class1="<td class=\"s7\" align=\"center\">";
3286  raw_class2="<td class=\"s12\" align=\"center\">";
3287  raw_class3="<td class=\"s14\" align=\"center\">";
3288  }
3289  const CellDB db;
3290  CellDB ce;
3291  if ((ce.size()>=1)&&(Sub[3][i]==sub)) {
3292 
3293  // AZ 19.03.2018
3294 // AZ 19 if (Sub[3][i]==1) {
3295 // AZ 19 ce = db.find("subdet", "HB").find("Eta", Eta[3][i]).find("Phi", Phi[3][i]).find("Depth", Depth[3][i]);
3296 // AZ 19 if (ce.size()==0) {cout<<"Error: No such HB, Eta="<< Eta[3][i] <<", Phi="<< Phi[3][i] <<", Depth="<< Depth[3][i] <<" in database"<<endl;continue;}
3297 // AZ 19 else if (ce.size()>1) {cout<<"Warning: More than one line correspond to such HB, Eta="<< Eta[3][i] <<", Phi="<< Phi[3][i] <<", Depth="<< Depth[3][i] <<" in database"<<endl;}
3298 // AZ 19 }
3299 // AZ 19 if (Sub[3][i]==2) {
3300 // AZ 19 ce = db.find("subdet", "HE").find("Eta", Eta[3][i]).find("Phi", Phi[3][i]).find("Depth", Depth[3][i]);
3301 // AZ 19 if (ce.size()==0) {cout<<"Error: No such HE, Eta="<< Eta[3][i] <<", Phi="<< Phi[3][i] <<", Depth="<< Depth[3][i] <<" in database"<<endl;continue;}
3302 // AZ 19 else if (ce.size()>1) {cout<<"Warning: More than one line correspond to such HE, Eta="<< Eta[3][i] <<", Phi="<< Phi[3][i] <<", Depth="<< Depth[3][i] <<" in database"<<endl;}
3303 // AZ 19 }
3304 // AZ 19 if (Sub[3][i]==3) {
3305 // AZ 19 ce = db.find("subdet", "HO").find("Eta", Eta[3][i]).find("Phi", Phi[3][i]).find("Depth", Depth[3][i]);
3306 // AZ 19 if (ce.size()==0) {cout<<"Error: No such HO, Eta="<< Eta[3][i] <<", Phi="<< Phi[3][i] <<", Depth="<< Depth[3][i] <<" in database"<<endl;continue;}
3307 // AZ 19 else if (ce.size()>1) {cout<<"Warning: More than one line correspond to such HO, Eta="<< Eta[3][i] <<", Phi="<< Phi[3][i] <<", Depth="<< Depth[3][i] <<" in database"<<endl;}
3308 // AZ 19 }
3309 // AZ 19 if (Sub[3][i]==4) {
3310 // AZ 19 ce = db.find("subdet", "HF").find("Eta", Eta[3][i]).find("Phi", Phi[3][i]).find("Depth", Depth[3][i]);
3311 // AZ 19 if (ce.size()==0) {cout<<"Error: No such HF, Eta="<< Eta[3][i] <<", Phi="<< Phi[3][i] <<", Depth="<< Depth[3][i] <<" in database"<<endl;continue;}
3312 // AZ 19 else if (ce.size()>1) {cout<<"Warning: More than one line correspond to such HF, Eta="<< Eta[3][i] <<", Phi="<< Phi[3][i] <<", Depth="<< Depth[3][i] <<" in database"<<endl;}
3313 // AZ 19 }
3314 
3315  htmlFile << "<td class=\"s4\" align=\"center\">" << ind+1 <<"</td>"<< std::endl;
3316  htmlFile << raw_class<< Eta[3][i]<<"</td>"<< std::endl;
3317  htmlFile << raw_class<< Phi[3][i]<<"</td>"<< std::endl;
3318  htmlFile << raw_class<< Depth[3][i] <<"</td>"<< std::endl;
3319  htmlFile << raw_class<< ce[0].RBX <<"</td>"<< std::endl;
3320  htmlFile << raw_class<< ce[0].RM <<"</td>"<< std::endl;
3321  htmlFile << raw_class<< ce[0].Pixel <<"</td>"<< std::endl;
3322  htmlFile << raw_class<< ce[0].RMfiber <<"</td>"<< std::endl;
3323  htmlFile << raw_class<< ce[0].FiberCh <<"</td>"<< std::endl;
3324  htmlFile << raw_class<< ce[0].QIE <<"</td>"<< std::endl;
3325  htmlFile << raw_class<< ce[0].ADC<<"</td>"<< std::endl;
3326  htmlFile << raw_class<< ce[0].VMECardID <<"</td>"<< std::endl;
3327  htmlFile << raw_class<< ce[0].dccID <<"</td>"<< std::endl;
3328  htmlFile << raw_class<< ce[0].Spigot <<"</td>"<< std::endl;
3329  htmlFile << raw_class<< ce[0].FiberIndex <<"</td>"<< std::endl;
3330  htmlFile << raw_class<< ce[0].HtrSlot <<"</td>"<< std::endl;
3331  htmlFile << raw_class<< ce[0].HtrTB <<"</td>"<< std::endl;
3332  htmlFile << raw_class1<< Map_Ampl[1][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3333  htmlFile << raw_class1<< Map_Ampl[2][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3334  htmlFile << raw_class1<< Map_Ampl[3][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3335  htmlFile << raw_class1<< Map_Ampl[4][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3336  htmlFile << raw_class1<< Map_Ampl[5][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3337  htmlFile << raw_class1<< Map_Ampl[6][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3338 // htmlFile << raw_class2<< Map_Ampl[11][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3339 // htmlFile << raw_class2<< Map_Ampl[12][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3340 // htmlFile << raw_class2<< Map_Ampl[13][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3341 // htmlFile << raw_class2<< Map_Ampl[14][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3342 // htmlFile << raw_class2<< Map_Ampl[15][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3343 // htmlFile << raw_class2<< Map_Ampl[16][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3344 // htmlFile << raw_class3<< Map_Ampl[21][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3345  htmlFile << raw_class<< Map_Ampl[31][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3346  htmlFile << raw_class<< Map_Ampl[32][Sub[3][i]][Depth[3][i]]->GetBinContent(Eta[3][i]+41,Phi[3][i]+1)<<"</td>"<< std::endl;
3347  htmlFile << "</tr>" << std::endl;
3348  htmlFile << "</tr>" << std::endl;
3349  ind+=1;
3350  }
3351  }
3352  htmlFile << "</table><br>" << std::endl;
3353  htmlFile << "<a href=\"#Top\">to top</a><br>\n";
3354 */
3356 
3357  htmlFile << "</body> " << std::endl;
3358  htmlFile << "</html> " << std::endl;
3359  htmlFile.close();
3360  }
3361 
3362  //======================================================================
3363  // Creating description html file:
3364  ofstream htmlFile;
3365 
3366  //======================================================================
3367  /*
3368  htmlFile.open("HELP.html");
3369  htmlFile << "</html><html xmlns=\"http://www.w3.org/1999/xhtml\">"<< std::endl;
3370  htmlFile << "<head>"<< std::endl;
3371  htmlFile << "<meta http-equiv=\"Content-Type\" content=\"text/html\"/>"<< std::endl;
3372  htmlFile << "<title> Remote Monitoring Tool </title>"<< std::endl;
3373  htmlFile << "<style type=\"text/css\">"<< std::endl;
3374  htmlFile << " body,td{ background-color: #FFFFCC; font-family: arial, arial ce, helvetica; font-size: 12px; }"<< std::endl;
3375  htmlFile << " td.s0 { font-family: arial, arial ce, helvetica; }"<< std::endl;
3376  htmlFile << " td.s1 { font-family: arial, arial ce, helvetica; font-weight: bold; background-color: #FFC169; text-align: center;}"<< std::endl;
3377  htmlFile << " td.s2 { font-family: arial, arial ce, helvetica; background-color: #eeeeee; }"<< std::endl;
3378  htmlFile << " td.s3 { font-family: arial, arial ce, helvetica; background-color: #d0d0d0; }"<< std::endl;
3379  htmlFile << " td.s4 { font-family: arial, arial ce, helvetica; background-color: #FFC169; }"<< std::endl;
3380  htmlFile << "</style>"<< std::endl;
3381  htmlFile << "<body>"<< std::endl;
3382  htmlFile << "<h1> Description of Remote Monitoring Tool criteria for bad channel selection</h1>"<< std::endl;
3383  htmlFile << "<br>"<< std::endl;
3384  htmlFile << "<h3> - C means CAPID Errors assuming we inspect CAPID non-rotation,error & validation bits, and for this criterion - no need to apply any cuts to select bcs.</h3> "<< std::endl;
3385  htmlFile << "<br>"<< std::endl;
3386  htmlFile << "<h3> - A means full amplitude, collected over all time slices </h3> "<< std::endl;
3387  htmlFile << "<h3> - R means ratio criterion where we define as a bad, the channels, for which the signal portion in 4 middle TSs(plus one, minus two around TS with maximal amplitude) is out of some range of reasonable values </h3> "<< std::endl;
3388  htmlFile << "<br>"<< std::endl;
3389  htmlFile << "<h3> - W means width of shape distribution. Width is defined as square root from dispersion. </h3> "<< std::endl;
3390  htmlFile << "<br>"<< std::endl;
3391  htmlFile << "<h3> - TN means mean time position of adc signal. </h3> "<< std::endl;
3392  htmlFile << "<br>"<< std::endl;
3393  htmlFile << "<h3> - TX means TS number of maximum signal </h3> "<< std::endl;
3394  htmlFile << "<br>"<< std::endl;
3395  htmlFile << "<h3> - m means megatile channels. For example Am means Amplitude criteria for megatile channels </h3> "<< std::endl;
3396  htmlFile << "<br>"<< std::endl;
3397  htmlFile << "<h3> - c means calibration channels. For example Ac means Amplitude criteria for calibration channels </h3> "<< std::endl;
3398  htmlFile << "<br>"<< std::endl;
3399  htmlFile << "<h3> - Pm means Pedestals. </h3> "<< std::endl;
3400  htmlFile << "<br>"<< std::endl;
3401  htmlFile << "<h3> - pWm means pedestal Width. </h3> "<< std::endl;
3402  htmlFile << "<br>"<< std::endl;
3403  htmlFile << "</body> " << std::endl;
3404  htmlFile << "</html> " << std::endl;
3405  htmlFile.close();
3406 */
3407  //======================================================================
3408 
3409  //======================================================================
3410  // Creating main html file:
3411  htmlFile.open("MAP.html");
3412  htmlFile << "</html><html xmlns=\"http://www.w3.org/1999/xhtml\">" << std::endl;
3413  htmlFile << "<head>" << std::endl;
3414  htmlFile << "<meta http-equiv=\"Content-Type\" content=\"text/html\"/>" << std::endl;
3415  htmlFile << "<title> Remote Monitoring Tool </title>" << std::endl;
3416  htmlFile << "<style type=\"text/css\">" << std::endl;
3417  htmlFile << " body,td{ background-color: #FFFFCC; font-family: arial, arial ce, helvetica; font-size: 12px; }"
3418  << std::endl;
3419  htmlFile << " td.s0 { font-family: arial, arial ce, helvetica; }" << std::endl;
3420  htmlFile << " td.s1 { font-family: arial, arial ce, helvetica; font-weight: bold; background-color: #FFC169; "
3421  "text-align: center;}"
3422  << std::endl;
3423  htmlFile << " td.s2 { font-family: arial, arial ce, helvetica; background-color: #eeeeee; }" << std::endl;
3424  htmlFile << " td.s3 { font-family: arial, arial ce, helvetica; background-color: #d0d0d0; }" << std::endl;
3425  htmlFile << " td.s4 { font-family: arial, arial ce, helvetica; background-color: #FFC169; }" << std::endl;
3426  htmlFile << " td.s5 { font-family: arial, arial ce, helvetica; background-color: #FF00FF; }" << std::endl;
3427  htmlFile << " td.s6 { font-family: arial, arial ce, helvetica; background-color: #9ACD32; }" << std::endl;
3428  htmlFile << " td.s7 { font-family: arial, arial ce, helvetica; background-color: #32CD32; }" << std::endl;
3429  htmlFile << "</style>" << std::endl;
3430  htmlFile << "<body>" << std::endl;
3431 
3432  htmlFile << "<h1> Remote Monitoring Tool, RUN = " << runnumber << ". </h1>" << std::endl;
3433  htmlFile << "<br>" << std::endl;
3434 
3435  htmlFile << "<h2> 1. Analysis results for subdetectors </h2>" << std::endl;
3436  htmlFile << "<table width=\"400\">" << std::endl;
3437  htmlFile << "<tr>" << std::endl;
3438 
3439  // AZ 12.03.2019
3440  /*
3441  htmlFile << " <td><a href=\"HB.html\">HB</a></td>"<< std::endl;
3442  htmlFile << " <td><a href=\"HE.html\">HE</a></td>"<< std::endl;
3443  htmlFile << " <td><a href=\"HO.html\">HO</a></td>"<< std::endl;
3444  htmlFile << " <td><a href=\"HF.html\">HF</a></td>"<< std::endl;
3445 */
3446 
3447  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
3448  << runnumber << "/HB.html\">HB</a></td>" << std::endl;
3449  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
3450  << runnumber << "/HE.html\">HE</a></td>" << std::endl;
3451  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
3452  << runnumber << "/HO.html\">HO</a></td>" << std::endl;
3453  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"
3454  << runnumber << "/HF.html\">HF</a></td>" << std::endl;
3455 
3456  htmlFile << "</tr>" << std::endl;
3457  htmlFile << "</table>" << std::endl;
3458  htmlFile << "<br>" << std::endl;
3459  //AZ2023:
3460  /*
3461  htmlFile << "<h2> 2. HCAL status over all criteria and subdetectors </h2>" << std::endl;
3462  htmlFile << "<h3> 2.A. Channels in detector space </h3>" << std::endl;
3463  htmlFile << "<h4> Legend for channel status: green - good, others - may be a problems, white - not applicable or out "
3464  "of range </h4>"
3465  << std::endl;
3466  htmlFile << " <img src=\"MAP.png\" />" << std::endl;
3467  htmlFile << "<br>" << std::endl;
3468 
3469  htmlFile << "<h3> 2.B. List of Bad channels </h3>" << std::endl;
3470 
3471  //htmlFile << " <td><a href=\"HELP.html\"> Description of criteria for bad channel selection</a></td>"<< std::endl;
3472  // htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"<<runnumber<<"/HELP.html\"> Description of criteria for bad channel selection</a></td>"<< std::endl;
3473  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/HELP.html\"> "
3474  "Description of criteria for bad channel selection</a></td>"
3475  << std::endl;
3476 
3477  htmlFile << "<table>" << std::endl;
3478  htmlFile << "<tr>";
3479  htmlFile << "<td class=\"s4\" align=\"center\">#</td>" << std::endl;
3480  htmlFile << "<td class=\"s1\" align=\"center\">ETA</td>" << std::endl;
3481  htmlFile << "<td class=\"s1\" align=\"center\">PHI</td>" << std::endl;
3482  htmlFile << "<td class=\"s1\" align=\"center\">DEPTH</td>" << std::endl;
3483  htmlFile << "<td class=\"s1\" align=\"center\">RBX</td>" << std::endl;
3484  htmlFile << "<td class=\"s1\" align=\"center\">RM</td>" << std::endl;
3485  htmlFile << "<td class=\"s1\" align=\"center\">PIXEL</td>" << std::endl;
3486  htmlFile << "<td class=\"s1\" align=\"center\">RM_FIBER</td>" << std::endl;
3487  htmlFile << "<td class=\"s1\" align=\"center\">FIBER_CH</td>" << std::endl;
3488  htmlFile << "<td class=\"s1\" align=\"center\">QIE</td>" << std::endl;
3489  htmlFile << "<td class=\"s1\" align=\"center\">ADC</td>" << std::endl;
3490  htmlFile << "<td class=\"s1\" align=\"center\">CRATE</td>" << std::endl;
3491  htmlFile << "<td class=\"s1\" align=\"center\">DCC</td>" << std::endl;
3492  htmlFile << "<td class=\"s1\" align=\"center\">SPIGOT</td>" << std::endl;
3493  htmlFile << "<td class=\"s1\" align=\"center\">HTR_FIBER</td>" << std::endl;
3494  htmlFile << "<td class=\"s1\" align=\"center\">HTR_SLOT</td>" << std::endl;
3495  htmlFile << "<td class=\"s1\" align=\"center\">HTR_FPGA</td>" << std::endl;
3496  htmlFile << "<td class=\"s5\" align=\"center\">RMT-criteria</td>" << std::endl;
3497  htmlFile << "</tr>" << std::endl;
3498 */
3499  //AZ2023 ind = 0;
3500  // AZ 19.03.2018
3501  /*
3502  for (int i=1;i<=NBad;i++) {
3503  if((ind%2)==1){
3504  raw_class="<td class=\"s2\" align=\"center\">";
3505  raw_class1="<td class=\"s6\" align=\"center\">";
3506  }else{
3507  raw_class="<td class=\"s3\" align=\"center\">";
3508  raw_class1="<td class=\"s7\" align=\"center\">";
3509  }
3510  const CellDB db;
3511  const CellDB ce = db.find("Eta", Eta[2][i]).find("Phi", Phi[2][i]).find("Depth", Depth[2][i]);
3512  // if (ce.size()==0) {cout<<"Error: No such Eta="<< Eta[2][i] <<", Phi="<< Phi[2][i] <<", Depth="<< Depth[2][i] <<" in database"<<endl;continue;}
3513 // else if (ce.size()>1) { cout<<"Warning: More than one line correspond to such Eta="<< Eta[2][i] <<", Phi="<< Phi[2][i] <<", Depth="<< Depth[2][i] <<" in database"<<endl;}
3514 
3515  if (ce.size()>=1) {
3516  htmlFile << "<tr>"<< std::endl;
3517  htmlFile << "<td class=\"s1\" align=\"center\">" << ind+1 <<"</td>"<< std::endl;
3518  htmlFile << raw_class<< Eta[2][i]<<"</td>"<< std::endl;
3519  htmlFile << raw_class<< Phi[2][i]<<"</td>"<< std::endl;
3520  htmlFile << raw_class<< Depth[2][i] <<"</td>"<< std::endl;
3521  htmlFile << raw_class<< ce[0].RBX <<"</td>"<< std::endl;
3522  htmlFile << raw_class<< ce[0].RM <<"</td>"<< std::endl;
3523  htmlFile << raw_class<< ce[0].Pixel <<"</td>"<< std::endl;
3524  htmlFile << raw_class<< ce[0].RMfiber <<"</td>"<< std::endl;
3525  htmlFile << raw_class<< ce[0].FiberCh <<"</td>"<< std::endl;
3526  htmlFile << raw_class<< ce[0].QIE <<"</td>"<< std::endl;
3527  htmlFile << raw_class<< ce[0].ADC<<"</td>"<< std::endl;
3528  htmlFile << raw_class<< ce[0].VMECardID <<"</td>"<< std::endl;
3529  htmlFile << raw_class<< ce[0].dccID <<"</td>"<< std::endl;
3530  htmlFile << raw_class<< ce[0].Spigot <<"</td>"<< std::endl;
3531  htmlFile << raw_class<< ce[0].FiberIndex <<"</td>"<< std::endl;
3532  htmlFile << raw_class<< ce[0].HtrSlot <<"</td>"<< std::endl;
3533  htmlFile << raw_class<< ce[0].HtrTB <<"</td>"<< std::endl;
3534  htmlFile << raw_class1<< Comment[2][i]<<"</td>"<< std::endl;
3535  htmlFile << "</tr>" << std::endl;
3536 
3537  ind+=1;
3538  }
3539  } /// end loop
3540 */
3541  //AZ2023 htmlFile << "</table>" << std::endl;
3542  //AZ2023 htmlFile << "<br>" << std::endl;
3543  /*
3544  htmlFile << "<h3> 2.C.List of Gain unstable channels </h3>"<< std::endl;
3545  //htmlFile << " <td><a href=\"HELP.html\"> Description of criteria for bad channel selection</a></td>"<< std::endl;
3546  // htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"<<runnumber<<"/HELP.html\"> Description of criteria for bad channel selection</a></td>"<< std::endl;
3547  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/HELP.html\"> Description of criteria for bad channel selection</a></td>"<< std::endl;
3548 
3549  htmlFile << "<table>"<< std::endl;
3550  htmlFile << "<tr>";
3551  htmlFile << "<td class=\"s4\" align=\"center\">#</td>" << std::endl;
3552  htmlFile << "<td class=\"s1\" align=\"center\">ETA</td>" << std::endl;
3553  htmlFile << "<td class=\"s1\" align=\"center\">PHI</td>" << std::endl;
3554  htmlFile << "<td class=\"s1\" align=\"center\">DEPTH</td>"<< std::endl;
3555  htmlFile << "<td class=\"s1\" align=\"center\">RBX</td>" << std::endl;
3556  htmlFile << "<td class=\"s1\" align=\"center\">RM</td>" << std::endl;
3557  htmlFile << "<td class=\"s1\" align=\"center\">PIXEL</td>" << std::endl;
3558  htmlFile << "<td class=\"s1\" align=\"center\">RM_FIBER</td>" << std::endl;
3559  htmlFile << "<td class=\"s1\" align=\"center\">FIBER_CH</td>" << std::endl;
3560  htmlFile << "<td class=\"s1\" align=\"center\">QIE</td>" << std::endl;
3561  htmlFile << "<td class=\"s1\" align=\"center\">ADC</td>" << std::endl;
3562  htmlFile << "<td class=\"s1\" align=\"center\">CRATE</td>" << std::endl;
3563  htmlFile << "<td class=\"s1\" align=\"center\">DCC</td>" << std::endl;
3564  htmlFile << "<td class=\"s1\" align=\"center\">SPIGOT</td>" << std::endl;
3565  htmlFile << "<td class=\"s1\" align=\"center\">HTR_FIBER</td>" << std::endl;
3566  htmlFile << "<td class=\"s1\" align=\"center\">HTR_SLOT</td>" << std::endl;
3567  htmlFile << "<td class=\"s1\" align=\"center\">HTR_FPGA</td>" << std::endl;
3568  htmlFile << "<td class=\"s5\" align=\"center\">Failed criteria</td>" << std::endl;
3569  htmlFile << "</tr>" << std::endl;
3570 
3571  for (int i=1;i<=NWarn;i++) {
3572  if((ind%2)==1){
3573  raw_class="<td class=\"s2\" align=\"center\">";
3574  raw_class1="<td class=\"s6\" align=\"center\">";
3575  }else{
3576  raw_class="<td class=\"s3\" align=\"center\">";
3577  raw_class1="<td class=\"s7\" align=\"center\">";
3578  }
3579  const CellDB db;
3580  const CellDB ce = db.find("Eta", Eta[1][i]).find("Phi", Phi[1][i]).find("Depth", Depth[1][i]);
3581  // if (ce.size()==0) {cout<<"Error: No such Eta="<< Eta[1][i] <<", Phi="<< Phi[1][i] <<", Depth="<< Depth[1][i] <<" in database"<<endl;continue;}
3582 // else if (ce.size()>1) { cout<<"Warning: More than one line correspond to such Eta="<< Eta[1][i] <<", Phi="<< Phi[1][i] <<", Depth="<< Depth[1][i] <<" in database"<<endl;}
3583 
3584  if (ce.size()>=1) {
3585  htmlFile << "<tr>"<< std::endl;
3586  htmlFile << "<td class=\"s1\" align=\"center\">" << ind+1 <<"</td>"<< std::endl;
3587  htmlFile << raw_class<< Eta[1][i]<<"</td>"<< std::endl;
3588  htmlFile << raw_class<< Phi[1][i]<<"</td>"<< std::endl;
3589  htmlFile << raw_class<< Depth[1][i] <<"</td>"<< std::endl;
3590  htmlFile << raw_class<< ce[0].RBX <<"</td>"<< std::endl;
3591  htmlFile << raw_class<< ce[0].RM <<"</td>"<< std::endl;
3592  htmlFile << raw_class<< ce[0].Pixel <<"</td>"<< std::endl;
3593  htmlFile << raw_class<< ce[0].RMfiber <<"</td>"<< std::endl;
3594  htmlFile << raw_class<< ce[0].FiberCh <<"</td>"<< std::endl;
3595  htmlFile << raw_class<< ce[0].QIE <<"</td>"<< std::endl;
3596  htmlFile << raw_class<< ce[0].ADC<<"</td>"<< std::endl;
3597  htmlFile << raw_class<< ce[0].VMECardID <<"</td>"<< std::endl;
3598  htmlFile << raw_class<< ce[0].dccID <<"</td>"<< std::endl;
3599  htmlFile << raw_class<< ce[0].Spigot <<"</td>"<< std::endl;
3600  htmlFile << raw_class<< ce[0].FiberIndex <<"</td>"<< std::endl;
3601  htmlFile << raw_class<< ce[0].HtrSlot <<"</td>"<< std::endl;
3602  htmlFile << raw_class<< ce[0].HtrTB <<"</td>"<< std::endl;
3603  htmlFile << raw_class1<< Comment[1][i]<<"</td>"<< std::endl;
3604  htmlFile << "</tr>" << std::endl;
3605 
3606  ind+=1;
3607  }
3608  }
3609 
3610 
3611  htmlFile << "</table>" << std::endl;
3612  htmlFile << "<br>"<< std::endl;
3613 
3614 
3615  htmlFile << "<h3> 2.D.List of channels with bad Pedestals </h3>"<< std::endl;
3616  // htmlFile << " <td><a href=\"HELP.html\"> Description of criteria for bad channel selection</a></td>"<< std::endl;
3617  // htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/GLOBAL_"<<runnumber<<"/HELP.html\"> Description of criteria for bad channel selection</a></td>"<< std::endl;
3618  htmlFile << " <td><a href=\"https://cms-conddb.cern.ch/eosweb/hcal/HcalRemoteMonitoring/GlobalRMT/HELP.html\"> Description of criteria for bad channel selection</a></td>"<< std::endl;
3619 
3620  htmlFile << "<table>"<< std::endl;
3621  htmlFile << "<tr>";
3622  htmlFile << "<td class=\"s4\" align=\"center\">#</td>" << std::endl;
3623  htmlFile << "<td class=\"s1\" align=\"center\">ETA</td>" << std::endl;
3624  htmlFile << "<td class=\"s1\" align=\"center\">PHI</td>" << std::endl;
3625  htmlFile << "<td class=\"s1\" align=\"center\">DEPTH</td>"<< std::endl;
3626  htmlFile << "<td class=\"s1\" align=\"center\">RBX</td>" << std::endl;
3627  htmlFile << "<td class=\"s1\" align=\"center\">RM</td>" << std::endl;
3628  htmlFile << "<td class=\"s1\" align=\"center\">PIXEL</td>" << std::endl;
3629  htmlFile << "<td class=\"s1\" align=\"center\">RM_FIBER</td>" << std::endl;
3630  htmlFile << "<td class=\"s1\" align=\"center\">FIBER_CH</td>" << std::endl;
3631  htmlFile << "<td class=\"s1\" align=\"center\">QIE</td>" << std::endl;
3632  htmlFile << "<td class=\"s1\" align=\"center\">ADC</td>" << std::endl;
3633  htmlFile << "<td class=\"s1\" align=\"center\">CRATE</td>" << std::endl;
3634  htmlFile << "<td class=\"s1\" align=\"center\">DCC</td>" << std::endl;
3635  htmlFile << "<td class=\"s1\" align=\"center\">SPIGOT</td>" << std::endl;
3636  htmlFile << "<td class=\"s1\" align=\"center\">HTR_FIBER</td>" << std::endl;
3637  htmlFile << "<td class=\"s1\" align=\"center\">HTR_SLOT</td>" << std::endl;
3638  htmlFile << "<td class=\"s1\" align=\"center\">HTR_FPGA</td>" << std::endl;
3639  htmlFile << "<td class=\"s5\" align=\"center\">Failed criteria</td>" << std::endl;
3640  htmlFile << "</tr>" << std::endl;
3641 
3642  for (int i=1;i<=NPed;i++) {
3643  if((ind%2)==1){
3644  raw_class="<td class=\"s2\" align=\"center\">";
3645  raw_class1="<td class=\"s6\" align=\"center\">";
3646  }else{
3647  raw_class="<td class=\"s3\" align=\"center\">";
3648  raw_class1="<td class=\"s7\" align=\"center\">";
3649  }
3650  const CellDB db;
3651  const CellDB ce = db.find("Eta", Eta[3][i]).find("Phi", Phi[3][i]).find("Depth", Depth[3][i]);
3652  // if (ce.size()==0) {cout<<"Error: No such Eta="<< Eta[3][i] <<", Phi="<< Phi[3][i] <<", Depth="<< Depth[3][i] <<" in database"<<endl;continue;}
3653 // else if (ce.size()>1) { cout<<"Warning: More than one line correspond to such Eta="<< Eta[1][i] <<", Phi="<< Phi[1][i] <<", Depth="<< Depth[1][i] <<" in database"<<endl;}
3654 
3655  if (ce.size()>=1) {
3656  htmlFile << "<tr>"<< std::endl;
3657  htmlFile << "<td class=\"s1\" align=\"center\">" << ind+1 <<"</td>"<< std::endl;
3658  htmlFile << raw_class<< Eta[3][i]<<"</td>"<< std::endl;
3659  htmlFile << raw_class<< Phi[3][i]<<"</td>"<< std::endl;
3660  htmlFile << raw_class<< Depth[3][i] <<"</td>"<< std::endl;
3661  htmlFile << raw_class<< ce[0].RBX <<"</td>"<< std::endl;
3662  htmlFile << raw_class<< ce[0].RM <<"</td>"<< std::endl;
3663  htmlFile << raw_class<< ce[0].Pixel <<"</td>"<< std::endl;
3664  htmlFile << raw_class<< ce[0].RMfiber <<"</td>"<< std::endl;
3665  htmlFile << raw_class<< ce[0].FiberCh <<"</td>"<< std::endl;
3666  htmlFile << raw_class<< ce[0].QIE <<"</td>"<< std::endl;
3667  htmlFile << raw_class<< ce[0].ADC<<"</td>"<< std::endl;
3668  htmlFile << raw_class<< ce[0].VMECardID <<"</td>"<< std::endl;
3669  htmlFile << raw_class<< ce[0].dccID <<"</td>"<< std::endl;
3670  htmlFile << raw_class<< ce[0].Spigot <<"</td>"<< std::endl;
3671  htmlFile << raw_class<< ce[0].FiberIndex <<"</td>"<< std::endl;
3672  htmlFile << raw_class<< ce[0].HtrSlot <<"</td>"<< std::endl;
3673  htmlFile << raw_class<< ce[0].HtrTB <<"</td>"<< std::endl;
3674  htmlFile << raw_class1<< Comment[3][i]<<"</td>"<< std::endl;
3675  htmlFile << "</tr>" << std::endl;
3676 
3677  ind+=1;
3678  }
3679  }
3680 
3681 
3682  htmlFile << "</table>" << std::endl;
3683 */
3684  htmlFile << "</body> " << std::endl;
3685  htmlFile << "</html> " << std::endl;
3686  htmlFile.close();
3687  //======================================================================
3688 
3689  //======================================================================
3690  // Close and delete all possible things:
3691  hfile->Close();
3692  // hfile->Delete();
3693  // Exit Root
3694  gSystem->Exit(0);
3695  //======================================================================
3696 }
mps_fire.i
i
Definition: mps_fire.py:429
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
DeadROC_duringRun.dir
dir
Definition: DeadROC_duringRun.py:23
ctpps_dqm_sourceclient-live_cfg.test
test
Definition: ctpps_dqm_sourceclient-live_cfg.py:7
alignmentValidation.fname
string fname
main script
Definition: alignmentValidation.py:959
x
float x
Definition: beamSpotDipStandalone.cc:55
main
int main(int argc, char *argv[])
Definition: RemoteMonitoringMAP_Global.cc:49
gather_cfg.cout
cout
Definition: gather_cfg.py:144
str
#define str(s)
Definition: TestProcessor.cc:56
HLT_2022v15_cff.cap
cap
Definition: HLT_2022v15_cff.py:4618
Generated for CMSSW Reference Manual by doxygen 1.8.14