CMS 3D CMS Logo
Functions
zcorrtsnvsampltest.cc File Reference
#include <TChain.h>
#include <TH1D.h>
#include <TH2D.h>
#include <TProfile.h>
#include <TFitResult.h>
#include <TFitResultPtr.h>
#include <TPaveStats.h>
#include <vector>
#include <string>
#include <iomanip>
#include <iostream>
#include <fstream>
#include "TH1.h"
#include "TH2.h"
#include "TCanvas.h"
#include "TROOT.h"
#include <TMath.h>
#include "TStyle.h"
#include "TSystem.h"
#include "TLegend.h"
#include "TText.h"
#include "TAxis.h"
#include "TFile.h"
#include "TLine.h"
#include "TGraph.h"
#include <TPaveText.h>
#include <TPostScript.h>
#include <TClass.h>

Go to the source code of this file.

Functions

int main (int argc, char *argv[])
 

Function Documentation

◆ main()

int main ( int  argc,
char *  argv[] 
)

===============================================================================================================================================================================================

variant2: for each run define phi-averaged A for normalization channel (Dref,16) and then, divide Rijk on it, i.e. get RRijk




















































































eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=29 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=28

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:







Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Sticking ADC counts in neughbouring TS for HB:

Sticking ADC counts in neughbouring TS for HE:

RBX:

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:





Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

errA with average Amplitudes

Summed Amplitude Plots:





Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

Prepare maps of good/bad channels:

Definition at line 38 of file zcorrtsnvsampltest.cc.

References alignmentValidation::c1, gather_cfg::cout, DeadROC_duringRun::dir, mps_fire::i, dqmiolumiharvest::j, AlCaHLTBitMon_QueryRunRegistry::string, and x.

38  {
39  std::string dirnm = "Analyzer";
40  //======================================================================
41  printf("reco: gROOT Reset \n");
42  gROOT->Reset();
43  gROOT->SetStyle("Plain");
44  gStyle->SetOptStat(0); // no statistics _or_
45  // gStyle->SetOptStat(11111111);
46  //gStyle->SetOptStat(1101);// name mean and rms
47  // gStyle->SetOptStat(0101);// name and entries
48  // gStyle->SetOptStat(1100);// mean and rms only !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
49  // gStyle->SetOptStat(1110000);// und over, integral !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
50  // gStyle->SetOptStat(101110); // entries, mean, rms, overflow !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
51  // gStyle->SetOptStat(100000);// over !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
52  //
53  //gStyle->SetOptFit(00010);// constant, mean and sigma only !!
54  // gStyle->SetOptFit(00001);// hi2/nu, constant, mean and sigma only !!
55  gStyle->SetOptFit(0010); // constant, mean and sigma only !!
56  // gStyle->SetOptFit(00011);// constant, mean and sigma only !!
57  // gStyle->SetOptFit(1101);
58  // gStyle->SetOptFit(1011);
59  //
60  //gStyle->SetStatX(0.98);
61  //gStyle->SetStatY(0.99);
62  //gStyle->SetStatW(0.30);
63  //gStyle->SetStatH(0.25);
64  gStyle->SetStatX(0.91);
65  gStyle->SetStatY(0.75);
66  gStyle->SetStatW(0.20);
67  gStyle->SetStatH(0.10);
68  // gStyle->SetStatH(0.35);
69  //
70 
71  // Float_t LeftOffset = 0.12;
72  //Float_t TopOffset = 0.12;
73  Float_t LeftOffset = 0.12;
74  Float_t TopOffset = 0.22;
75 
76  gStyle->SetLineWidth(1);
77  gStyle->SetErrorX(0);
78 
79  //---=[ Titles,Labels ]=-----------
80  gStyle->SetOptTitle(0); // title on/off
81  // gStyle->SetTitleColor(0); // title color
82  gStyle->SetTitleColor(1); // title color
83  // gStyle->SetTitleX(0.35); // title x-position
84  gStyle->SetTitleX(0.15); // title x-position
85  gStyle->SetTitleH(0.15); // title height
86  // gStyle->SetTitleW(0.53); // title width
87  gStyle->SetTitleW(0.60); // title width
88  gStyle->SetTitleFont(42); // title font
89  gStyle->SetTitleFontSize(0.07); // title font size
90 
91  gStyle->SetPalette(1);
92  //---=[ Histogram style ]=----------
93  // gStyle->SetHistFillColor(38);
94  // gStyle->SetFrameFillColor(41);// jelto-kori4nev-svetl
95  // gStyle->SetFrameFillColor(5);// jeltyi
96  // gStyle->SetFrameFillColor(17);// seryi
97  // gStyle->SetFrameFillColor(18);// svetlo seryi
98  // gStyle->SetFrameFillColor(20);// svetlo kori4nev
99  // gStyle->SetFrameFillColor(33);// sine-seryi
100  // gStyle->SetFrameFillColor(40);// fiolet-seryi
101  // gStyle->SetFrameFillColor(23);// sv.kor
102 
103  //---=[ Pad style ]=----------------
104  gStyle->SetPadTopMargin(TopOffset);
105  gStyle->SetPadBottomMargin(LeftOffset);
106  gStyle->SetPadRightMargin(TopOffset);
107  gStyle->SetPadLeftMargin(LeftOffset);
108  //---=[ SetCanvasDef ]=----------------
109  //======================================================================
110  //
111  // Connect the input file and get the 2-d histogram in memory
112  //======================================================================
113  // TBrowser *b = new TBrowser
114 
115  TFile *hfile1 = new TFile("Global_362596.root", "READ");
116  //TFile *hfile1 = new TFile("Global_362365.root", "READ");
117 
118  TPostScript psfile("zcorrtsnvsampltest.ps", 111);
119 
120  //
121 
122  TCanvas *c1 = new TCanvas("c1", "Hcal4test", 200, 10, 700, 900);
123 
124  hfile1->ls();
125  TDirectory *dir = (TDirectory *)hfile1->FindObjectAny(dirnm.c_str());
126 
127  //=============================================================================================== 1
128  //======================================================================
129  //======================================================================
130  //================
131  //======================================================================
132  c1->Clear();
133  c1->Divide(2, 2);
134 
135  c1->cd(1);
136  TH2F *two11 = (TH2F *)dir->FindObjectAny("h2_TSnVsAyear2023_HB");
137  gPad->SetGridy();
138  gPad->SetGridx();
139  two11->SetMarkerStyle(20);
140  two11->SetMarkerSize(0.4);
141  two11->SetYTitle("timing HB \b");
142  two11->SetXTitle("Q HB\b");
143  two11->SetMarkerColor(1);
144  two11->SetLineColor(1);
145  // gPad->SetLogx();
146  // gPad->SetLogy();
147  two11->Draw("BOX");
148  // two11->Draw("");
149  // two11->Draw("COLZ");
150 
151  c1->cd(2);
152  TH2F *two12 = (TH2F *)dir->FindObjectAny("h2_TSnVsAyear2023_HE");
153  gPad->SetGridy();
154  gPad->SetGridx();
155  two12->SetMarkerStyle(20);
156  two12->SetMarkerSize(0.4);
157  two12->SetYTitle("timing HE \b");
158  two12->SetXTitle("Q HE\b");
159  two12->SetMarkerColor(1);
160  two12->SetLineColor(1);
161  // gPad->SetLogy();
162  two12->Draw("BOX");
163  // two12->Draw("SCAT");
164 
165  c1->cd(3);
166  TH2F *two22 = (TH2F *)dir->FindObjectAny("h2_TSnVsAyear2023_HF");
167  gPad->SetGridy();
168  gPad->SetGridx();
169  two22->SetMarkerStyle(20);
170  two22->SetMarkerSize(0.4);
171  two22->SetYTitle("timing HF \b");
172  two22->SetXTitle("Q HF\b");
173  two22->SetMarkerColor(1);
174  two22->SetLineColor(1);
175  two22->Draw("BOX");
176  // two22->Draw("ARR");
177 
178  c1->cd(4);
179  TH2F *two23 = (TH2F *)dir->FindObjectAny("h2_TSnVsAyear2023_HO");
180  gPad->SetGridy();
181  gPad->SetGridx();
182  two23->SetMarkerStyle(20);
183  two23->SetMarkerSize(0.4);
184  two23->SetYTitle("timing HO \b");
185  two23->SetXTitle("Q HO\b");
186  two23->SetMarkerColor(1);
187  two23->SetLineColor(1);
188  two23->Draw("BOX");
189 
190  c1->Update();
191 
192  //========================================================================================== 2
193  //======================================================================
194  //======================================================================
195  //================
196  //======================================================================
197  c1->Clear();
198  c1->Divide(1, 4);
199 
200  c1->cd(1);
201  TH1F *TSNvsQ_HB = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HB");
202  gPad->SetGridy();
203  gPad->SetGridx();
204  gPad->SetLogy();
205  TSNvsQ_HB->SetMarkerStyle(20);
206  TSNvsQ_HB->SetMarkerSize(0.6);
207  TSNvsQ_HB->GetYaxis()->SetLabelSize(0.04);
208  TSNvsQ_HB->SetXTitle("Q,fb HB \b");
209  TSNvsQ_HB->SetYTitle("iev*ieta*iphi*idepth \b");
210  TSNvsQ_HB->SetMarkerColor(4);
211  TSNvsQ_HB->SetLineColor(0);
212  TSNvsQ_HB->SetMinimum(0.8);
213  //TSNvsQ_HB->GetXaxis()->SetRangeUser(0, MaxLumDanila + 5.);
214  TSNvsQ_HB->Draw("E");
215 
216  c1->cd(2);
217  TH1F *TSNvsQ_HE = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HE");
218  gPad->SetGridy();
219  gPad->SetGridx();
220  gPad->SetLogy();
221  TSNvsQ_HE->SetMarkerStyle(20);
222  TSNvsQ_HE->SetMarkerSize(0.6);
223  TSNvsQ_HE->GetYaxis()->SetLabelSize(0.04);
224  TSNvsQ_HE->SetXTitle("Q,fb HE \b");
225  TSNvsQ_HE->SetYTitle("iev*ieta*iphi*idepth \b");
226  TSNvsQ_HE->SetMarkerColor(4);
227  TSNvsQ_HE->SetLineColor(0);
228  TSNvsQ_HE->SetMinimum(0.8);
229  //TSNvsQ_HE->GetXaxis()->SetRangeUser(0, MaxLumDanila + 5.);
230  TSNvsQ_HE->Draw("E");
231 
232  c1->cd(3);
233  TH1F *TSNvsQ_HF = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HF");
234  gPad->SetGridy();
235  gPad->SetGridx();
236  gPad->SetLogy();
237  TSNvsQ_HF->SetMarkerStyle(20);
238  TSNvsQ_HF->SetMarkerSize(0.6);
239  TSNvsQ_HF->GetYaxis()->SetLabelSize(0.04);
240  TSNvsQ_HF->SetXTitle("Q,fb HF \b");
241  TSNvsQ_HF->SetYTitle("iev*ieta*iphi*idepth \b");
242  TSNvsQ_HF->SetMarkerColor(4);
243  TSNvsQ_HF->SetLineColor(0);
244  TSNvsQ_HF->SetMinimum(0.8);
245  //TSNvsQ_HF->GetXaxis()->SetRangeUser(0, MaxLumDanila + 5.);
246  TSNvsQ_HF->Draw("E");
247 
248  c1->cd(4);
249  TH1F *TSNvsQ_HO = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HO");
250  gPad->SetGridy();
251  gPad->SetGridx();
252  gPad->SetLogy();
253  TSNvsQ_HO->SetMarkerStyle(20);
254  TSNvsQ_HO->SetMarkerSize(0.6);
255  TSNvsQ_HO->GetYaxis()->SetLabelSize(0.04);
256  TSNvsQ_HO->SetXTitle("Q,fb HO \b");
257  TSNvsQ_HO->SetYTitle("iev*ieta*iphi*idepth \b");
258  TSNvsQ_HO->SetMarkerColor(4);
259  TSNvsQ_HO->SetLineColor(0);
260  TSNvsQ_HO->SetMinimum(0.8);
261  //TSNvsQ_HO->GetXaxis()->SetRangeUser(0, MaxLumDanila + 5.);
262  TSNvsQ_HO->Draw("E");
263 
264  c1->Update();
265 
266  //========================================================================================== 3
267  //======================================================================
268  //======================================================================
269  //================
270  //======================================================================
271  c1->Clear();
272  c1->Divide(1, 4);
273 
274  c1->cd(1);
275  TH1F *twod1_HB = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear2023_HB");
276  TH1F *twod0_HB = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HB");
277  twod1_HB->Sumw2();
278  twod0_HB->Sumw2();
279  TH1F *Ceff_HB = (TH1F *)twod1_HB->Clone("Ceff_HB");
280  // Ceff_HB->Sumw2();
281  for (int x = 1; x <= twod1_HB->GetXaxis()->GetNbins(); x++) {
282  twod1_HB->SetBinError(float(x), 0.001);
283  } //end x
284  Ceff_HB->Divide(twod1_HB, twod0_HB, 1, 1, "B");
285  gPad->SetGridy();
286  gPad->SetGridx();
287  Ceff_HB->SetMarkerStyle(20);
288  Ceff_HB->SetMarkerSize(0.4);
289  Ceff_HB->SetXTitle("Q, fc \b");
290  Ceff_HB->SetYTitle("<timing>HB \b");
291  Ceff_HB->SetMarkerColor(2);
292  Ceff_HB->SetLineColor(2);
293  Ceff_HB->SetMaximum(140.);
294  Ceff_HB->SetMinimum(30.);
295  Ceff_HB->Draw("E");
296  //Ceff_HB->Draw("COLZ");
297 
298  c1->cd(2);
299  TH1F *twod1_HE = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear2023_HE");
300  TH1F *twod0_HE = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HE");
301  // twod1_HE->Sumw2();
302  // twod0_HE->Sumw2();
303  TH1F *Ceff_HE = (TH1F *)twod1_HE->Clone("Ceff_HE");
304  // Ceff_HE->Sumw2();
305  for (int x = 1; x <= twod1_HE->GetXaxis()->GetNbins(); x++) {
306  twod1_HE->SetBinError(float(x), 0.001);
307  } //end x
308  Ceff_HE->Divide(twod1_HE, twod0_HE, 1, 1, "B");
309  gPad->SetGridy();
310  gPad->SetGridx();
311  Ceff_HE->SetMarkerStyle(20);
312  Ceff_HE->SetMarkerSize(0.4);
313  Ceff_HE->SetXTitle("Q, fc \b");
314  Ceff_HE->SetYTitle("<timing>HE \b");
315  Ceff_HE->SetMarkerColor(2);
316  Ceff_HE->SetLineColor(2);
317  Ceff_HE->SetMaximum(150.);
318  Ceff_HE->SetMinimum(25.);
319  Ceff_HE->Draw("E");
320  //Ceff_HE->Draw("COLZ");
321 
322  c1->cd(3);
323  TH1F *twod1_HF = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear2023_HF");
324  TH1F *twod0_HF = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HF");
325  // twod1_HF->Sumw2();
326  // twod0_HF->Sumw2();
327  TH1F *Ceff_HF = (TH1F *)twod1_HF->Clone("Ceff_HF");
328  // Ceff_HF->Sumw2();
329  for (int x = 1; x <= twod1_HF->GetXaxis()->GetNbins(); x++) {
330  twod1_HF->SetBinError(float(x), 0.001);
331  } //end x
332  Ceff_HF->Divide(twod1_HF, twod0_HF, 1, 1, "B");
333  gPad->SetGridy();
334  gPad->SetGridx();
335  Ceff_HF->SetMarkerStyle(20);
336  Ceff_HF->SetMarkerSize(0.4);
337  Ceff_HF->SetXTitle("Q, fc \b");
338  Ceff_HF->SetYTitle("<timing>HF \b");
339  Ceff_HF->SetMarkerColor(2);
340  Ceff_HF->SetLineColor(2);
341  Ceff_HF->SetMaximum(50.);
342  Ceff_HF->SetMinimum(0.);
343  Ceff_HF->Draw("E");
344  //Ceff_HF->Draw("COLZ");
345 
346  c1->cd(4);
347  TH1F *twod1_HO = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear2023_HO");
348  TH1F *twod0_HO = (TH1F *)dir->FindObjectAny("h1_TSnVsAyear20230_HO");
349  // twod1_HO->Sumw2();
350  // twod0_HO->Sumw2();
351  TH1F *Ceff_HO = (TH1F *)twod1_HO->Clone("Ceff_HO");
352  // Ceff_HO->Sumw2();
353  for (int x = 1; x <= twod1_HO->GetXaxis()->GetNbins(); x++) {
354  twod1_HO->SetBinError(float(x), 0.001);
355  } //end x
356  Ceff_HO->Divide(twod1_HO, twod0_HO, 1, 1, "B");
357  gPad->SetGridy();
358  gPad->SetGridx();
359  Ceff_HO->SetMarkerStyle(20);
360  Ceff_HO->SetMarkerSize(0.4);
361  Ceff_HO->SetXTitle("Q, fc \b");
362  Ceff_HO->SetYTitle("<timing>HO \b");
363  Ceff_HO->SetMarkerColor(2);
364  Ceff_HO->SetLineColor(2);
365  Ceff_HO->SetMaximum(150.);
366  Ceff_HO->SetMinimum(70.);
367  Ceff_HO->Draw("E");
368  //Ceff_HO->Draw("COLZ");
369 
370  c1->Update();
371 
372  //========================================================================================== 4
373  //======================================================================
374  //======================================================================
375  //================
376  //======================================================================
377  c1->Clear();
378  c1->Divide(2, 2);
379 
380  c1->cd(1);
381  TH2F *dva1_HBDepth1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA_HB");
382  TH2F *dva0_HBDepth1 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HB");
383  //dva1_HBDepth1->Sumw2();
384  //dva0_HBDepth1->Sumw2();
385  TH2F *Seff_HBDepth1 = (TH2F *)dva1_HBDepth1->Clone("Seff_HBDepth1");
386  Seff_HBDepth1->Divide(dva1_HBDepth1, dva0_HBDepth1, 25., 1., "B");
387  gPad->SetGridy();
388  gPad->SetGridx();
389  //gPad->SetLogz();
390  Seff_HBDepth1->SetMarkerStyle(20);
391  Seff_HBDepth1->SetMarkerSize(0.4);
392  Seff_HBDepth1->SetXTitle("#eta \b");
393  Seff_HBDepth1->SetYTitle("#phi \b");
394  Seff_HBDepth1->SetZTitle("<timing> HB Depth1 \b");
395  Seff_HBDepth1->SetMarkerColor(2);
396  Seff_HBDepth1->SetLineColor(2);
397  Seff_HBDepth1->SetMaximum(90.);
398  Seff_HBDepth1->SetMinimum(85.);
399  Seff_HBDepth1->Draw("COLZ");
400 
401  c1->cd(2);
402  TH2F *dva1_HBDepth2 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA_HB");
403  TH2F *dva0_HBDepth2 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HB");
404  //dva1_HBDepth2->Sumw2();
405  //dva0_HBDepth2->Sumw2();
406  TH2F *Seff_HBDepth2 = (TH2F *)dva1_HBDepth2->Clone("Seff_HBDepth2");
407  Seff_HBDepth2->Divide(dva1_HBDepth2, dva0_HBDepth2, 25., 1., "B");
408  gPad->SetGridy();
409  gPad->SetGridx();
410  //gPad->SetLogz();
411  Seff_HBDepth2->SetMarkerStyle(20);
412  Seff_HBDepth2->SetMarkerSize(0.4);
413  Seff_HBDepth2->SetXTitle("#eta \b");
414  Seff_HBDepth2->SetYTitle("#phi \b");
415  Seff_HBDepth2->SetZTitle("<timing> HB Depth2 \b");
416  Seff_HBDepth2->SetMarkerColor(2);
417  Seff_HBDepth2->SetLineColor(2);
418  Seff_HBDepth2->SetMaximum(90.);
419  Seff_HBDepth2->SetMinimum(85.);
420  Seff_HBDepth2->Draw("COLZ");
421 
422  c1->cd(3);
423  TH2F *dva1_HBDepth3 = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmeanA_HB");
424  TH2F *dva0_HBDepth3 = (TH2F *)dir->FindObjectAny("h_mapDepth3_HB");
425  //dva1_HBDepth3->Sumw2();
426  //dva0_HBDepth3->Sumw2();
427  TH2F *Seff_HBDepth3 = (TH2F *)dva1_HBDepth3->Clone("Seff_HBDepth3");
428  Seff_HBDepth3->Divide(dva1_HBDepth3, dva0_HBDepth3, 25., 1., "B");
429  gPad->SetGridy();
430  gPad->SetGridx();
431  //gPad->SetLogz();
432  Seff_HBDepth3->SetMarkerStyle(20);
433  Seff_HBDepth3->SetMarkerSize(0.4);
434  Seff_HBDepth3->SetXTitle("#eta \b");
435  Seff_HBDepth3->SetYTitle("#phi \b");
436  Seff_HBDepth3->SetZTitle("<timing> HB Depth3 \b");
437  Seff_HBDepth3->SetMarkerColor(2);
438  Seff_HBDepth3->SetLineColor(2);
439  Seff_HBDepth3->SetMaximum(90.);
440  Seff_HBDepth3->SetMinimum(85.);
441  Seff_HBDepth3->Draw("COLZ");
442 
443  c1->cd(4);
444  TH2F *dva1_HBDepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA_HB");
445  TH2F *dva0_HBDepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4_HB");
446  //dva1_HBDepth4->Sumw2();
447  //dva0_HBDepth4->Sumw2();
448  TH2F *Seff_HBDepth4 = (TH2F *)dva1_HBDepth4->Clone("Seff_HBDepth4");
449  Seff_HBDepth4->Divide(dva1_HBDepth4, dva0_HBDepth4, 25., 1., "B");
450  gPad->SetGridy();
451  gPad->SetGridx();
452  //gPad->SetLogz();
453  Seff_HBDepth4->SetMarkerStyle(20);
454  Seff_HBDepth4->SetMarkerSize(0.4);
455  Seff_HBDepth4->SetXTitle("#eta \b");
456  Seff_HBDepth4->SetYTitle("#phi \b");
457  Seff_HBDepth4->SetZTitle("<timing> HB Depth4 \b");
458  Seff_HBDepth4->SetMarkerColor(2);
459  Seff_HBDepth4->SetLineColor(2);
460  Seff_HBDepth4->SetMaximum(90.);
461  Seff_HBDepth4->SetMinimum(85.);
462  Seff_HBDepth4->Draw("COLZ");
463 
464  c1->Update();
465 
466  //========================================================================================== 5
467  //======================================================================
468  //======================================================================
469  //================
470  //======================================================================
471  c1->Clear();
472  c1->Divide(2, 4);
473 
474  c1->cd(1);
475  TH2F *dva1_HEDepth1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA_HE");
476  TH2F *dva0_HEDepth1 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HE");
477  //dva1_HEDepth1->Sumw2();
478  //dva0_HEDepth1->Sumw2();
479  TH2F *Seff_HEDepth1 = (TH2F *)dva1_HEDepth1->Clone("Seff_HEDepth1");
480  Seff_HEDepth1->Divide(dva1_HEDepth1, dva0_HEDepth1, 25., 1., "B");
481  gPad->SetGridy();
482  gPad->SetGridx();
483  //gPad->SetLogz();
484  Seff_HEDepth1->SetMarkerStyle(20);
485  Seff_HEDepth1->SetMarkerSize(0.4);
486  Seff_HEDepth1->SetXTitle("#eta \b");
487  Seff_HEDepth1->SetYTitle("#phi \b");
488  Seff_HEDepth1->SetZTitle("<timing> HE Depth1 \b");
489  Seff_HEDepth1->SetMarkerColor(2);
490  Seff_HEDepth1->SetLineColor(2);
491  Seff_HEDepth1->SetMaximum(90.);
492  Seff_HEDepth1->SetMinimum(85.);
493  Seff_HEDepth1->Draw("COLZ");
494 
495  c1->cd(2);
496  TH2F *dva1_HEDepth2 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA_HE");
497  TH2F *dva0_HEDepth2 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HE");
498  //dva1_HEDepth2->Sumw2();
499  //dva0_HEDepth2->Sumw2();
500  TH2F *Seff_HEDepth2 = (TH2F *)dva1_HEDepth2->Clone("Seff_HEDepth2");
501  Seff_HEDepth2->Divide(dva1_HEDepth2, dva0_HEDepth2, 25., 1., "B");
502  gPad->SetGridy();
503  gPad->SetGridx();
504  //gPad->SetLogz();
505  Seff_HEDepth2->SetMarkerStyle(20);
506  Seff_HEDepth2->SetMarkerSize(0.4);
507  Seff_HEDepth2->SetXTitle("#eta \b");
508  Seff_HEDepth2->SetYTitle("#phi \b");
509  Seff_HEDepth2->SetZTitle("<timing> HE Depth2 \b");
510  Seff_HEDepth2->SetMarkerColor(2);
511  Seff_HEDepth2->SetLineColor(2);
512  Seff_HEDepth2->SetMaximum(90.);
513  Seff_HEDepth2->SetMinimum(85.);
514  Seff_HEDepth2->Draw("COLZ");
515 
516  c1->cd(3);
517  TH2F *dva1_HEDepth3 = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmeanA_HE");
518  TH2F *dva0_HEDepth3 = (TH2F *)dir->FindObjectAny("h_mapDepth3_HE");
519  //dva1_HEDepth3->Sumw2();
520  //dva0_HEDepth3->Sumw2();
521  TH2F *Seff_HEDepth3 = (TH2F *)dva1_HEDepth3->Clone("Seff_HEDepth3");
522  Seff_HEDepth3->Divide(dva1_HEDepth3, dva0_HEDepth3, 25., 1., "B");
523  gPad->SetGridy();
524  gPad->SetGridx();
525  //gPad->SetLogz();
526  Seff_HEDepth3->SetMarkerStyle(20);
527  Seff_HEDepth3->SetMarkerSize(0.4);
528  Seff_HEDepth3->SetXTitle("#eta \b");
529  Seff_HEDepth3->SetYTitle("#phi \b");
530  Seff_HEDepth3->SetZTitle("<timing> HE Depth3 \b");
531  Seff_HEDepth3->SetMarkerColor(2);
532  Seff_HEDepth3->SetLineColor(2);
533  Seff_HEDepth3->SetMaximum(90.);
534  Seff_HEDepth3->SetMinimum(85.);
535  Seff_HEDepth3->Draw("COLZ");
536 
537  c1->cd(4);
538  TH2F *dva1_HEDepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA_HE");
539  TH2F *dva0_HEDepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4_HE");
540  //dva1_HEDepth4->Sumw2();
541  //dva0_HEDepth4->Sumw2();
542  TH2F *Seff_HEDepth4 = (TH2F *)dva1_HEDepth4->Clone("Seff_HEDepth4");
543  Seff_HEDepth4->Divide(dva1_HEDepth4, dva0_HEDepth4, 25., 1., "B");
544  gPad->SetGridy();
545  gPad->SetGridx();
546  //gPad->SetLogz();
547  Seff_HEDepth4->SetMarkerStyle(20);
548  Seff_HEDepth4->SetMarkerSize(0.4);
549  Seff_HEDepth4->SetXTitle("#eta \b");
550  Seff_HEDepth4->SetYTitle("#phi \b");
551  Seff_HEDepth4->SetZTitle("<timing> HE Depth4 \b");
552  Seff_HEDepth4->SetMarkerColor(2);
553  Seff_HEDepth4->SetLineColor(2);
554  Seff_HEDepth4->SetMaximum(90.);
555  Seff_HEDepth4->SetMinimum(85.);
556  Seff_HEDepth4->Draw("COLZ");
557 
558  c1->cd(5);
559  TH2F *dva1_HEDepth5 = (TH2F *)dir->FindObjectAny("h_mapDepth5TSmeanA_HE");
560  TH2F *dva0_HEDepth5 = (TH2F *)dir->FindObjectAny("h_mapDepth5_HE");
561  //dva1_HEDepth5->Sumw2();
562  //dva0_HEDepth5->Sumw2();
563  TH2F *Seff_HEDepth5 = (TH2F *)dva1_HEDepth5->Clone("Seff_HEDepth5");
564  Seff_HEDepth5->Divide(dva1_HEDepth5, dva0_HEDepth5, 25., 1., "B");
565  gPad->SetGridy();
566  gPad->SetGridx();
567  //gPad->SetLogz();
568  Seff_HEDepth5->SetMarkerStyle(20);
569  Seff_HEDepth5->SetMarkerSize(0.4);
570  Seff_HEDepth5->SetXTitle("#eta \b");
571  Seff_HEDepth5->SetYTitle("#phi \b");
572  Seff_HEDepth5->SetZTitle("<timing> HE Depth5 \b");
573  Seff_HEDepth5->SetMarkerColor(2);
574  Seff_HEDepth5->SetLineColor(2);
575  Seff_HEDepth5->SetMaximum(90.);
576  Seff_HEDepth5->SetMinimum(85.);
577  Seff_HEDepth5->Draw("COLZ");
578 
579  c1->cd(6);
580  TH2F *dva1_HEDepth6 = (TH2F *)dir->FindObjectAny("h_mapDepth6TSmeanA_HE");
581  TH2F *dva0_HEDepth6 = (TH2F *)dir->FindObjectAny("h_mapDepth6_HE");
582  //dva1_HEDepth6->Sumw2();
583  //dva0_HEDepth6->Sumw2();
584  TH2F *Seff_HEDepth6 = (TH2F *)dva1_HEDepth6->Clone("Seff_HEDepth6");
585  Seff_HEDepth6->Divide(dva1_HEDepth6, dva0_HEDepth6, 25., 1., "B");
586  gPad->SetGridy();
587  gPad->SetGridx();
588  //gPad->SetLogz();
589  Seff_HEDepth6->SetMarkerStyle(20);
590  Seff_HEDepth6->SetMarkerSize(0.4);
591  Seff_HEDepth6->SetXTitle("#eta \b");
592  Seff_HEDepth6->SetYTitle("#phi \b");
593  Seff_HEDepth6->SetZTitle("<timing> HE Depth6 \b");
594  Seff_HEDepth6->SetMarkerColor(2);
595  Seff_HEDepth6->SetLineColor(2);
596  Seff_HEDepth6->SetMaximum(90.);
597  Seff_HEDepth6->SetMinimum(85.);
598  Seff_HEDepth6->Draw("COLZ");
599 
600  c1->cd(7);
601  TH2F *dva1_HEDepth7 = (TH2F *)dir->FindObjectAny("h_mapDepth7TSmeanA_HE");
602  TH2F *dva0_HEDepth7 = (TH2F *)dir->FindObjectAny("h_mapDepth7_HE");
603  //dva1_HEDepth7->Sumw2();
604  //dva0_HEDepth7->Sumw2();
605  TH2F *Seff_HEDepth7 = (TH2F *)dva1_HEDepth7->Clone("Seff_HEDepth7");
606  Seff_HEDepth7->Divide(dva1_HEDepth7, dva0_HEDepth7, 25., 1., "B");
607  gPad->SetGridy();
608  gPad->SetGridx();
609  //gPad->SetLogz();
610  Seff_HEDepth7->SetMarkerStyle(20);
611  Seff_HEDepth7->SetMarkerSize(0.4);
612  Seff_HEDepth7->SetXTitle("#eta \b");
613  Seff_HEDepth7->SetYTitle("#phi \b");
614  Seff_HEDepth7->SetZTitle("<timing> HE Depth7 \b");
615  Seff_HEDepth7->SetMarkerColor(2);
616  Seff_HEDepth7->SetLineColor(2);
617  Seff_HEDepth7->SetMaximum(90.);
618  Seff_HEDepth7->SetMinimum(85.);
619  Seff_HEDepth7->Draw("COLZ");
620 
621  c1->Update();
622 
623  //========================================================================================== 6
624  //======================================================================
625  //======================================================================
626  //================
627  //======================================================================
628  c1->Clear();
629  c1->Divide(2, 2);
630 
631  c1->cd(1);
632  TH2F *dva1_HFDepth1 = (TH2F *)dir->FindObjectAny("h_mapDepth1TSmeanA_HF");
633  TH2F *dva0_HFDepth1 = (TH2F *)dir->FindObjectAny("h_mapDepth1_HF");
634  //dva1_HFDepth1->Sumw2();
635  //dva0_HFDepth1->Sumw2();
636  TH2F *Seff_HFDepth1 = (TH2F *)dva1_HFDepth1->Clone("Seff_HFDepth1");
637  Seff_HFDepth1->Divide(dva1_HFDepth1, dva0_HFDepth1, 25., 1., "B");
638  gPad->SetGridy();
639  gPad->SetGridx();
640  //gPad->SetLogz();
641  Seff_HFDepth1->SetMarkerStyle(20);
642  Seff_HFDepth1->SetMarkerSize(0.4);
643  Seff_HFDepth1->SetXTitle("#eta \b");
644  Seff_HFDepth1->SetYTitle("#phi \b");
645  Seff_HFDepth1->SetZTitle("<timing> HF Depth1 \b");
646  Seff_HFDepth1->SetMarkerColor(2);
647  Seff_HFDepth1->SetLineColor(2);
648  Seff_HFDepth1->SetMaximum(27.);
649  Seff_HFDepth1->SetMinimum(23.);
650  Seff_HFDepth1->Draw("COLZ");
651 
652  c1->cd(2);
653  TH2F *dva1_HFDepth2 = (TH2F *)dir->FindObjectAny("h_mapDepth2TSmeanA_HF");
654  TH2F *dva0_HFDepth2 = (TH2F *)dir->FindObjectAny("h_mapDepth2_HF");
655  //dva1_HFDepth2->Sumw2();
656  //dva0_HFDepth2->Sumw2();
657  TH2F *Seff_HFDepth2 = (TH2F *)dva1_HFDepth2->Clone("Seff_HFDepth2");
658  Seff_HFDepth2->Divide(dva1_HFDepth2, dva0_HFDepth2, 25., 1., "B");
659  gPad->SetGridy();
660  gPad->SetGridx();
661  //gPad->SetLogz();
662  Seff_HFDepth2->SetMarkerStyle(20);
663  Seff_HFDepth2->SetMarkerSize(0.4);
664  Seff_HFDepth2->SetXTitle("#eta \b");
665  Seff_HFDepth2->SetYTitle("#phi \b");
666  Seff_HFDepth2->SetZTitle("<timing> HF Depth2 \b");
667  Seff_HFDepth2->SetMarkerColor(2);
668  Seff_HFDepth2->SetLineColor(2);
669  Seff_HFDepth2->SetMaximum(27.);
670  Seff_HFDepth2->SetMinimum(23.);
671  Seff_HFDepth2->Draw("COLZ");
672 
673  c1->cd(3);
674  TH2F *dva1_HFDepth3 = (TH2F *)dir->FindObjectAny("h_mapDepth3TSmeanA_HF");
675  TH2F *dva0_HFDepth3 = (TH2F *)dir->FindObjectAny("h_mapDepth3_HF");
676  //dva1_HFDepth3->Sumw2();
677  //dva0_HFDepth3->Sumw2();
678  TH2F *Seff_HFDepth3 = (TH2F *)dva1_HFDepth3->Clone("Seff_HFDepth3");
679  Seff_HFDepth3->Divide(dva1_HFDepth3, dva0_HFDepth3, 25., 1., "B");
680  gPad->SetGridy();
681  gPad->SetGridx();
682  //gPad->SetLogz();
683  Seff_HFDepth3->SetMarkerStyle(20);
684  Seff_HFDepth3->SetMarkerSize(0.4);
685  Seff_HFDepth3->SetXTitle("#eta \b");
686  Seff_HFDepth3->SetYTitle("#phi \b");
687  Seff_HFDepth3->SetZTitle("<timing> HF Depth3 \b");
688  Seff_HFDepth3->SetMarkerColor(2);
689  Seff_HFDepth3->SetLineColor(2);
690  Seff_HFDepth3->SetMaximum(27.);
691  Seff_HFDepth3->SetMinimum(23.);
692  Seff_HFDepth3->Draw("COLZ");
693 
694  c1->cd(4);
695  TH2F *dva1_HFDepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA_HF");
696  TH2F *dva0_HFDepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4_HF");
697  //dva1_HFDepth4->Sumw2();
698  //dva0_HFDepth4->Sumw2();
699  TH2F *Seff_HFDepth4 = (TH2F *)dva1_HFDepth4->Clone("Seff_HFDepth4");
700  Seff_HFDepth4->Divide(dva1_HFDepth4, dva0_HFDepth4, 25., 1., "B");
701  gPad->SetGridy();
702  gPad->SetGridx();
703  //gPad->SetLogz();
704  Seff_HFDepth4->SetMarkerStyle(20);
705  Seff_HFDepth4->SetMarkerSize(0.4);
706  Seff_HFDepth4->SetXTitle("#eta \b");
707  Seff_HFDepth4->SetYTitle("#phi \b");
708  Seff_HFDepth4->SetZTitle("<timing> HF Depth4 \b");
709  Seff_HFDepth4->SetMarkerColor(2);
710  Seff_HFDepth4->SetLineColor(2);
711  Seff_HFDepth4->SetMaximum(27.);
712  Seff_HFDepth4->SetMinimum(23.);
713  Seff_HFDepth4->Draw("COLZ");
714 
715  c1->Update();
716 
717  //========================================================================================== 7
718  //======================================================================
719  //======================================================================
720  //================
721  //======================================================================
722  c1->Clear();
723  c1->Divide(1, 1);
724 
725  c1->cd(1);
726  TH2F *dva1_HODepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4TSmeanA_HO");
727  for (int i = 1; i <= dva1_HODepth4->GetXaxis()->GetNbins(); i++) {
728  for (int j = 1; j <= dva1_HODepth4->GetYaxis()->GetNbins(); j++) {
729  double ccc1 = dva1_HODepth4->GetBinContent(i, j);
730  if (ccc1 > 0.)
731  std::cout << "****** dva1_HODepth4 ************** i = " << i << " j = " << j << " ccc1 = " << ccc1
732  << std::endl;
733  }
734  }
735  TH2F *dva0_HODepth4 = (TH2F *)dir->FindObjectAny("h_mapDepth4_HO");
736  for (int i = 1; i <= dva0_HODepth4->GetXaxis()->GetNbins(); i++) {
737  for (int j = 1; j <= dva0_HODepth4->GetYaxis()->GetNbins(); j++) {
738  double ccc1 = dva0_HODepth4->GetBinContent(i, j);
739  if (ccc1 > 0.)
740  std::cout << "****** dva0_HODepth4 ************** i = " << i << " j = " << j << " ccc1 = " << ccc1
741  << std::endl;
742  }
743  }
744  TH2F *Seff_HODepth4 = (TH2F *)dva1_HODepth4->Clone("Seff_HODepth4");
745  Seff_HODepth4->Divide(dva1_HODepth4, dva0_HODepth4, 25., 1., "B");
746  for (int i = 1; i <= Seff_HODepth4->GetXaxis()->GetNbins(); i++) {
747  for (int j = 1; j <= Seff_HODepth4->GetYaxis()->GetNbins(); j++) {
748  double ccc1 = Seff_HODepth4->GetBinContent(i, j);
749  if (ccc1 > 0.)
750  std::cout << "****** Seff_HODepth4 ************** i = " << i << " j = " << j << " ccc1 = " << ccc1
751  << std::endl;
752  }
753  }
754  gPad->SetGridy();
755  gPad->SetGridx();
756  Seff_HODepth4->SetMarkerStyle(20);
757  Seff_HODepth4->SetMarkerSize(0.4);
758  Seff_HODepth4->SetXTitle("#eta \b");
759  Seff_HODepth4->SetYTitle("#phi \b");
760  Seff_HODepth4->SetZTitle("<timing> HB Depth1 \b");
761  Seff_HODepth4->SetMarkerColor(2);
762  Seff_HODepth4->SetLineColor(2);
763  Seff_HODepth4->SetMaximum(115.);
764  Seff_HODepth4->SetMinimum(110.);
765  Seff_HODepth4->Draw("COLZ");
766 
767  c1->Update();
768  //======================================================================
769  //======================================================================================== end
770  //======================================================================
771  //======================================================================
772  // close and delete all possible things:
773 
774  // psfile->Close();
775  psfile.Close();
776 
777  hfile1->Close();
778  // hfile1->Delete();
779  hfile1->Close();
780  // hfile1->Delete();
781 
782  // Exit Root
783  gSystem->Exit(0);
784  //======================================================================
785 }
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
x
float x
Definition: beamSpotDipStandalone.cc:55
gather_cfg.cout
cout
Definition: gather_cfg.py:144
Generated for CMSSW Reference Manual by doxygen 1.8.14