|
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
|
void CSCGainsStudy_MakePlots(int debug){
/** Macro to plot histograms produced by CSCGainsStudy.cc
*
* Author: Dominique Fortin - UCR
*/
TFile *file = TFile::Open("csc_strip_gains.root");
// set suffixps to desired file type: e.g. .eps, .jpg, ...
TString suffixps = ".jpg";
for (int i = 1; i < 37; i++ ) {
int j = 0;
if ( j == i ) TString chamber = "All_CSC";
if (++j == i ) TString chamber = "ME_11_27";
if (++j == i ) TString chamber = "ME_11_28";
if (++j == i ) TString chamber = "ME_11_29";
if (++j == i ) TString chamber = "ME_11_30";
if (++j == i ) TString chamber = "ME_11_31";
if (++j == i ) TString chamber = "ME_11_32";
if (++j == i ) TString chamber = "ME_12_27";
if (++j == i ) TString chamber = "ME_12_28";
if (++j == i ) TString chamber = "ME_12_29";
if (++j == i ) TString chamber = "ME_12_30";
if (++j == i ) TString chamber = "ME_12_31";
if (++j == i ) TString chamber = "ME_12_32";
if (++j == i ) TString chamber = "ME_13_27";
if (++j == i ) TString chamber = "ME_13_28";
if (++j == i ) TString chamber = "ME_13_29";
if (++j == i ) TString chamber = "ME_13_30";
if (++j == i ) TString chamber = "ME_13_31";
if (++j == i ) TString chamber = "ME_13_32";
if (++j == i ) TString chamber = "ME_21_14";
if (++j == i ) TString chamber = "ME_21_15";
if (++j == i ) TString chamber = "ME_21_16";
if (++j == i ) TString chamber = "ME_22_27";
if (++j == i ) TString chamber = "ME_22_28";
if (++j == i ) TString chamber = "ME_22_29";
if (++j == i ) TString chamber = "ME_22_30";
if (++j == i ) TString chamber = "ME_22_31";
if (++j == i ) TString chamber = "ME_22_32";
if (++j == i ) TString chamber = "ME_31_14";
if (++j == i ) TString chamber = "ME_31_15";
if (++j == i ) TString chamber = "ME_31_16";
if (++j == i ) TString chamber = "ME_32_27";
if (++j == i ) TString chamber = "ME_32_28";
if (++j == i ) TString chamber = "ME_32_29";
if (++j == i ) TString chamber = "ME_32_30";
if (++j == i ) TString chamber = "ME_32_31";
if (++j == i ) TString chamber = "ME_32_32";
// Set pointers to histograms
hGains = (TH1F *) file->Get("hGain_"+chamber);
hGaindiff = (TH1F *) file->Get("hGaindiff_"+chamber);
hGainvsch = (TH2F *) file->Get("hGainvsch_"+chamber);
gStyle->SetOptFit(0111);
// 1) weight
TString plot1 = "strip_weight_"+chamber+suffixps;
gStyle->SetOptStat(kTRUE);
TCanvas *c1 = new TCanvas("c1","");
c1->SetFillColor(10);
c1->SetLogy(1);
hGains->Draw();
c1->Print(plot1);
// 2) weight difference
TString plot2 = "delta_strip_weight_"+chamber+suffixps;
gStyle->SetOptStat(kTRUE);
TCanvas *c1 = new TCanvas("c1","");
c1->SetLogy(1);
c1->SetFillColor(10);
hGaindiff->Draw();
hGaindiff->Fit("gaus");
c1->Print(plot2);
if (i > 0) {
// 3) strip weight vs channel #
TString plot3 = "strip_weight_vs_channel_"+chamber+suffixps;
gStyle->SetOptStat(kTRUE);
TCanvas *c1 = new TCanvas("c1","");
c1->SetLogy(0);
c1->SetFillColor(10);
hGainvsch->Draw("BOX");
c1->Print(plot3);
}
}
// gROOT->ProcessLine(".q");
}
|