|
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
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
|
/*
* Produces 1D plots of chamber efficiency starting from the DQM output files
* of the DTChamberEfficiency and DTChamberEfficiencyTest modules.
*
*
* G. Cerminara 2009
*/
#if !defined(__CINT) || defined(__MAKECINT__)
#include "TH2F.h"
#include "TCanvas.h"
#include "TGraphErrors.h"
#include "TStyle.h"
#include "TFile.h"
// #include "macros.C"
#include <vector>
#endif
using namespace std;
vector<TGraphErrors *> plotEffwh(const TH2F *hEffDQM) {
vector<TGraphErrors *> ret;
if(hEffDQM != 0) {
TCanvas *c = newCanvas(hEffDQM->GetName(), hEffDQM->GetName());
hEffDQM->Draw("COLZ");
double sector[12];
double sectErr[12];
double effiency[4][12];
double efficiencyErr[4][12];
cout << "D-1" << endl;
for(int station = 1; station != 5; ++station) { //loop over stations
for(int sect = 1; sect != 13; ++sect) { // loop over sectors
sector[sect-1] = sect;
effiency[station-1][sect-1] = hEffDQM->GetBinContent(sect,station);
efficiencyErr[station-1][sect-1] = hEffDQM->GetBinError(sect,station);
sectErr[sect-1] = 0.5;
}
}
cout << "D-0.5" << endl;
ret.push_back(new TGraphErrors(12, sector, effiency[0], sectErr, efficiencyErr[0]));
ret.push_back(new TGraphErrors(12, sector, effiency[1], sectErr, efficiencyErr[1]));
ret.push_back(new TGraphErrors(12, sector, effiency[2], sectErr, efficiencyErr[2]));
ret.push_back(new TGraphErrors(12, sector, effiency[3], sectErr, efficiencyErr[3]));
cout << "D0" << endl;
cout << " ret.size() : " << ret.size() << endl;
return ret;
} else return ret;
}
void plotEfficiency(){
// Load needed macros and files
gROOT->LoadMacro("macros.C"); // Load service macros
// Get the style
TStyle * style = getStyle("tdr");
style->cd();
// retrieve the last open file
TFile *file = (TFile *)(gROOT->GetListOfFiles()->Last());
if(file == 0) {
cout << "No file loaded!Exiting." << endl;
return;
}
cout << "Loading file: " << file->GetName() << endl;
// retrieve the original histos from the file
TH2F *histoEff_Whm2 = (TH2F *) file->Get("DQMData/Run 109459/DT/Run summary/05-ChamberEff/EfficiencyMap_All_W-2");
// if(histoEff_Whm2 != 0) {
// TCanvas *c = newCanvas(histoEff_Whm2->GetName(), histoEff_Whm2->GetName());
// histoEff_Whm2->Draw("COLZ");
// double sector[12];
// double sectErr[12];
// double effiency[4][12];
// double efficiencyErr[4][12];
// for(int station = 1; station != 5; ++station) { //loop over stations
// for(int sect = 1; sect != 13; ++sect) { // loop over sectors
// sector[sect-1] = sect;
// effiency[station-1][sect-1] = histoEff_Whm2->GetBinContent(sect,station);
// efficiencyErr[station-1][sect-1] = histoEff_Whm2->GetBinError(sect,station);
// sectErr[sect-1] = 0.5;
// }
// }
// TGraphErrors *mb1Eff_Whm2 = new TGraphErrors(12, sector, effiency[0], sectErr, efficiencyErr[0]);
// TGraphErrors *mb2Eff_Whm2 = new TGraphErrors(12, sector, effiency[1], sectErr, efficiencyErr[1]);
// TGraphErrors *mb3Eff_Whm2 = new TGraphErrors(12, sector, effiency[2], sectErr, efficiencyErr[2]);
// TGraphErrors *mb4Eff_Whm2 = new TGraphErrors(12, sector, effiency[3], sectErr, efficiencyErr[3]);
vector<TGraphErrors *> graphsWhM2 = plotEffwh(histoEff_Whm2);
cout << "DD" << endl;
cout << " size: " << graphsWhM2.size() << endl;
cout << " " << (int) graphsWhM2[0] << " "
<< (int) graphsWhM2[1] << " "
<< (int) graphsWhM2[2] << " "
<< (int) graphsWhM2[3] << " " << endl;
TGraphErrors *mb1Eff_Whm2 = graphsWhM2[0];
TGraphErrors *mb2Eff_Whm2 = graphsWhM2[1];
TGraphErrors *mb3Eff_Whm2 = graphsWhM2[2];
TGraphErrors *mb4Eff_Whm2 = graphsWhM2[3];
cout << "D1" << endl;
TCanvas *c1 = newCanvas("Efficiency_Whm2", "Efficiency_Whm2");
mb1Eff_Whm2->Draw("AP");
mb1Eff_Whm2->SetMarkerColor(1);
mb2Eff_Whm2->Draw("P");
mb2Eff_Whm2->SetMarkerColor(2);
mb3Eff_Whm2->Draw("P");
mb3Eff_Whm2->SetMarkerColor(3);
mb4Eff_Whm2->Draw("P");
mb4Eff_Whm2->SetMarkerColor(4);
// }
}
|