|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
|
#include "SiStripQualityHistoryPlots.h"
#include <iostream>
#include <algorithm>
#include <vector>
#include <string>
#include <map>
#include "TPad.h"
#include "TFile.h"
#include "TH2F.h"
#include "TH1F.h"
#include "TProfile.h"
#include "TGraph.h"
#include "DPGAnalysis/SiStripTools/interface/CommonAnalyzer.h"
#include "TCanvas.h"
#include "TStyle.h"
TH1D* AverageRunBadChannels(TFile& ff, const char* module, const char* histo, const bool excludeLastBins) {
CommonAnalyzer camult(&ff, "", module);
TH1D* badchannels = new TH1D("badchannels", "Average Number of Bad Channels vs run", 10, 0., 10.);
badchannels->SetCanExtend(TH1::kXaxis);
std::vector<unsigned int> runs = camult.getRunList();
std::sort(runs.begin(), runs.end());
{
for (unsigned int i = 0; i < runs.size(); ++i) {
char runlabel[100];
sprintf(runlabel, "%d", runs[i]);
char runpath[100];
sprintf(runpath, "run_%d", runs[i]);
camult.setPath(runpath);
TProfile* multvstime = nullptr;
if (multvstime == nullptr)
multvstime = (TProfile*)camult.getObject(histo);
if (multvstime) {
// compute mean exlucing the last filled bins
if (excludeLastBins) {
int lastbin = multvstime->GetNbinsX() + 1;
int firstbin = 1;
for (int ibin = multvstime->GetNbinsX() + 1; ibin > 0; --ibin) {
if (multvstime->GetBinEntries(ibin) != 0) {
lastbin = ibin;
break;
}
}
std::cout << "Restricted range: " << firstbin << " " << lastbin << std::endl;
multvstime->GetXaxis()->SetRangeUser(multvstime->GetBinLowEdge(firstbin),
multvstime->GetBinLowEdge(lastbin - 2));
}
// fill the summary
badchannels->Fill(runlabel, multvstime->GetMean(2));
}
}
}
return badchannels;
}
TCanvas* StripCompletePlot(TFile& ff, const char* module, const bool excludeLastBins) {
TCanvas* cc = new TCanvas();
TH1D* cabling = AverageRunBadChannels(ff, module, "badmodrun_Cabling", excludeLastBins);
TH1D* runinfo = AverageRunBadChannels(ff, module, "badmodrun_RunInfo", excludeLastBins);
TH1D* badchannel = AverageRunBadChannels(ff, module, "badmodrun_BadChannel", excludeLastBins);
TH1D* dcs = AverageRunBadChannels(ff, module, "badmodrun_DCS", excludeLastBins);
TH1D* badfiber = AverageRunBadChannels(ff, module, "badmodrun_BadFiber", excludeLastBins);
cabling->SetLineColor(kRed);
runinfo->SetLineColor(kMagenta);
badchannel->SetLineColor(kCyan);
dcs->SetLineColor(kGreen);
badfiber->SetLineColor(kBlue);
badfiber->Draw();
dcs->Draw("same");
badchannel->Draw("same");
runinfo->Draw("same");
cabling->Draw("same");
return cc;
}
|