|
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
|
#include "TH2F.h"
#include "TLine.h"
inline void DrawEB(TH2F* ebmap, float min, float max) {
ebmap->SetXTitle("i#phi");
ebmap->SetYTitle("i#eta");
ebmap->GetXaxis()->SetNdivisions(-418, kFALSE);
ebmap->GetYaxis()->SetNdivisions(-1702, kFALSE);
ebmap->GetXaxis()->SetLabelSize(0.03);
ebmap->GetYaxis()->SetLabelSize(0.03);
ebmap->GetXaxis()->SetTickLength(0.01);
ebmap->GetYaxis()->SetTickLength(0.01);
ebmap->SetMaximum(max);
ebmap->SetMinimum(min);
ebmap->Draw("colz");
TLine* l = new TLine;
l->SetLineWidth(1);
for (int i = 0; i < 17; i++) {
Double_t x = 20. + (i * 20);
l = new TLine(x, -85., x, 85.);
l->Draw();
}
l = new TLine(0., 85., 360., 85.);
l->Draw();
l = new TLine(0., 0., 360., 0.);
l->Draw();
} // DrawEB method
inline void DrawEE(TH2F* endc, float min, float max) {
int ixSectorsEE[202] = {
62, 62, 61, 61, 60, 60, 59, 59, 58, 58, 56, 56, 46, 46, 44, 44, 43, 43, 42, 42, 41, 41, 40, 40, 41, 41,
42, 42, 43, 43, 44, 44, 46, 46, 56, 56, 58, 58, 59, 59, 60, 60, 61, 61, 62, 62, 0, 101, 101, 98, 98, 96,
96, 93, 93, 88, 88, 86, 86, 81, 81, 76, 76, 66, 66, 61, 61, 41, 41, 36, 36, 26, 26, 21, 21, 16, 16, 14,
14, 9, 9, 6, 6, 4, 4, 1, 1, 4, 4, 6, 6, 9, 9, 14, 14, 16, 16, 21, 21, 26, 26, 36, 36, 41,
41, 61, 61, 66, 66, 76, 76, 81, 81, 86, 86, 88, 88, 93, 93, 96, 96, 98, 98, 101, 101, 0, 62, 66, 66, 71,
71, 81, 81, 91, 91, 93, 0, 62, 66, 66, 91, 91, 98, 0, 58, 61, 61, 66, 66, 71, 71, 76, 76, 81, 81, 0,
51, 51, 0, 44, 41, 41, 36, 36, 31, 31, 26, 26, 21, 21, 0, 40, 36, 36, 11, 11, 4, 0, 40, 36, 36, 31,
31, 21, 21, 11, 11, 9, 0, 46, 46, 41, 41, 36, 36, 0, 56, 56, 61, 61, 66, 66};
int iySectorsEE[202] = {51, 56, 56, 58, 58, 59, 59, 60, 60, 61, 61, 62, 62, 61, 61, 60, 60, 59, 59, 58, 58, 56, 56,
46, 46, 44, 44, 43, 43, 42, 42, 41, 41, 40, 40, 41, 41, 42, 42, 43, 43, 44, 44, 46, 46, 51,
0, 51, 61, 61, 66, 66, 76, 76, 81, 81, 86, 86, 88, 88, 93, 93, 96, 96, 98, 98, 101, 101, 98,
98, 96, 96, 93, 93, 88, 88, 86, 86, 81, 81, 76, 76, 66, 66, 61, 61, 41, 41, 36, 36, 26, 26,
21, 21, 16, 16, 14, 14, 9, 9, 6, 6, 4, 4, 1, 1, 4, 4, 6, 6, 9, 9, 14, 14, 16,
16, 21, 21, 26, 26, 36, 36, 41, 41, 51, 0, 46, 46, 41, 41, 36, 36, 31, 31, 26, 26, 0, 51,
51, 56, 56, 61, 61, 0, 61, 61, 66, 66, 71, 71, 76, 76, 86, 86, 88, 0, 62, 101, 0, 61, 61,
66, 66, 71, 71, 76, 76, 86, 86, 88, 0, 51, 51, 56, 56, 61, 61, 0, 46, 46, 41, 41, 36, 36,
31, 31, 26, 26, 0, 40, 31, 31, 16, 16, 6, 0, 40, 31, 31, 16, 16, 6};
TLine* l = new TLine;
l->SetLineWidth(1);
endc->SetXTitle("ix");
endc->SetYTitle("iy");
endc->SetMaximum(max);
endc->SetMinimum(min);
endc->Draw("colz1");
for (int i = 0; i < 201; i = i + 1) {
if ((ixSectorsEE[i] != 0 || iySectorsEE[i] != 0) && (ixSectorsEE[i + 1] != 0 || iySectorsEE[i + 1] != 0)) {
l->DrawLine(ixSectorsEE[i], iySectorsEE[i], ixSectorsEE[i + 1], iySectorsEE[i + 1]);
}
}
} // DrawEE method
inline void DrawEE_Tower(TH2F* endc, TLine* l, double minScale, double maxScale) {
endc->SetStats(false);
endc->SetMinimum(minScale);
endc->SetMaximum(maxScale);
endc->Draw("colz");
int ixSectorsEE[136] = {8, 14, 14, 17, 17, 18, 18, 19, 19, 20, 20, 21, 21, 20, 20, 19, 19, 18, 18, 17, 17, 14, 14,
8, 8, 5, 5, 4, 4, 3, 3, 2, 2, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 8, 8, 8,
9, 9, 10, 10, 12, 12, 13, 13, 12, 12, 10, 10, 9, 9, 10, 10, 0, 11, 11, 0, 10, 9, 9,
8, 8, 7, 7, 6, 6, 5, 5, 0, 12, 13, 13, 14, 14, 15, 15, 16, 16, 17, 17, 0, 9, 8,
8, 3, 3, 1, 0, 13, 14, 14, 19, 19, 21, 0, 9, 8, 8, 7, 7, 5, 5, 3, 3, 2, 0,
13, 14, 14, 15, 15, 17, 17, 19, 19, 20, 0, 14, 14, 13, 13, 12, 12, 0};
int iySectorsEE[136] = {1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 8, 8, 14, 14, 17, 17, 18, 18, 19, 19, 20, 20, 21,
21, 20, 20, 19, 19, 18, 18, 17, 17, 14, 14, 8, 8, 5, 5, 4, 4, 3, 3, 2, 2, 1, 4,
4, 7, 7, 9, 9, 10, 10, 12, 12, 13, 13, 12, 12, 10, 10, 9, 0, 13, 21, 0, 13, 13, 14,
14, 15, 15, 16, 16, 18, 18, 19, 0, 13, 13, 14, 14, 15, 15, 16, 16, 18, 18, 19, 0, 11, 11,
12, 12, 13, 13, 0, 11, 11, 12, 12, 13, 13, 0, 10, 10, 9, 9, 8, 8, 7, 7, 6, 6, 0,
10, 10, 9, 9, 8, 8, 7, 7, 6, 6, 0, 2, 4, 4, 7, 7, 9, 0};
for (int i = 0; i < 136; i = i + 1)
if ((ixSectorsEE[i] != 0 || iySectorsEE[i] != 0) && (ixSectorsEE[i + 1] != 0 || iySectorsEE[i + 1] != 0))
l->DrawLine(ixSectorsEE[i], iySectorsEE[i], ixSectorsEE[i + 1], iySectorsEE[i + 1]);
} //draw EE in case of a tower
inline void drawTable(int nbRows, int nbColumns) {
TLine* l = new TLine;
l->SetLineWidth(1);
for (int i = 1; i < nbRows; i++) {
double y = (double)i;
l = new TLine(0., y, nbColumns, y);
l->Draw();
}
for (int i = 1; i < nbColumns; i++) {
double x = (double)i;
double y = (double)nbRows;
l = new TLine(x, 0., x, y);
l->Draw();
}
}
|