|
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
|
////////////////////////////
// Script to plot sensor efficiency layer by layer
// The input files are the output from CalibTracker/SiStripHitEfficiency/src/HitRes.cc
// Usage in root: root> .x plotHitEffSummaryQualGoodCRAFT09.C
// Original Author: Keith Ulmer University of Colorado
// keith.ulmer@colorado.edu
////////////////////////////
{
#include <vector>
#include <tdrstyle.C>
gROOT->Reset();
setTDRStyle();
int nLayers = 21;
float RatioNoBad[nLayers+1];
float errRatio[nLayers+1];
float RatioAll[nLayers+1];
float errRatioAll[nLayers+1];
float Mod[nLayers+1];
float ModErr[nLayers+1];
int cont[nLayers+1];
TH1F * found = new TH1F("found","found",nLayers+1,0,nLayers+1);
TH1F * all = new TH1F("all","all",nLayers+1,0,nLayers+1);
float SNoBad[nLayers+1],BNoBad[nLayers+1];
float SAll[nLayers+1],BAll[nLayers+1];
SNoBad[0] = 0;
BNoBad[0] = 6;
SAll[0] = 0;
BAll[0] = 6;
for (Long_t i=1; i< nLayers+1; ++i) {
if (i==10) i++;
TString file = TString("HitEff_V11_69912_layer") + i;
TFile* f = TFile::Open(file+".root");
f->cd("anEff");
bool isTEC = false;
if (i>13) isTEC = true;
TString cut;
if (isTEC) cut = TString("SiStripQualBad==0&&(Id<402674821||Id>470078964||(Id<470048885&Id>402674977))&&abs(TrajLocY)>2");
else cut = TString("SiStripQualBad==0&&(Id<402674821||Id>470078964||(Id<470048885&Id>402674977))");
traj->Draw("ModIsBad",cut,"goff");
//traj->Draw("ModIsBad","SiStripQualBad==0&&(Id<402674821||Id>470078964||(Id<470048885&Id>402674977))&&abs(TrajLocY)>2","goff");
//if (!isTEC)
//traj->Draw("ModIsBad","SiStripQualBad==0&&(Id<402674821||Id>470078964||(Id<470048885&Id>402674977))","goff");
//else
//traj->Draw("ModIsBad","SiStripQualBad==0&&(Id<402674821||Id>470078964||(Id<470048885&Id>402674977))","goff");
//traj->Draw("ModIsBad","SiStripQualBad==0","goff");
//exclude bad modules from TID 2: 402674821-402674977
// and TEC 3: 470048885-470078964
SNoBad[i] = htemp->GetBinContent(1);
if (i==13) SNoBad[i] = htemp->GetBinContent(3); //not sure why the entries are in bin 3 for layer 13?
BNoBad[i] = htemp->GetBinContent(2);
cout << "filling for layer " << i << " s = " << htemp->GetBinContent(1) << " b = " << htemp->GetBinContent(2) << endl;
if ((SNoBad[i]+BNoBad[i]) > 5) {
float ratio = (SNoBad[i]*1. / (SNoBad[i]+BNoBad[i])*1.);
RatioNoBad[i] = ratio;
float deno = (SNoBad[i] + BNoBad[i]) * (SNoBad[i] + BNoBad[i]) * (SNoBad[i] + BNoBad[i]);
errRatio[i] = sqrt( (SNoBad[i]*BNoBad[i]) / deno*1.);
found->SetBinContent(i,SNoBad[i]);
all->SetBinContent(i,SNoBad[i]+BNoBad[i]);
} else {
RatioNoBad[i] = -1.0;
errRatio[i] = 0.0;
found->SetBinContent(i,0);
all->SetBinContent(i,10);
}
cout << i << " SNoBad " << SNoBad[i] << " BNoBad " << BNoBad[i] << " ratio " << ratio*100. << " +o- " << errRatio[i]*100. << endl;
Mod[i] = i;
ModErr[i] = 0.;
}
TCanvas *c7 =new TCanvas("c7"," test ",10,10,800,600);
c7->SetFillColor(0);
c7->SetGrid();
found->Sumw2();
all->Sumw2();
gr = new TGraphAsymmErrors(nLayers+1);
gr->BayesDivide(found,all);
for(int j = 0; j<nLayers+1; j++){
gr->SetPointError(j, 0., 0., gr->GetErrorYlow(j),gr->GetErrorYhigh(j) );
}
gr->GetXaxis()->SetLimits(0,nLayers);
gr->SetMarkerColor(2);
gr->SetMarkerSize(1.2);
gr->SetLineColor(2);
gr->SetLineWidth(4);
gr->SetMarkerStyle(20);
gr->SetMinimum(0.98);
gr->SetMaximum(1.001);
gr->GetYaxis()->SetTitle("Efficiency");
for ( int j=1; j<nLayers+1; j++) {
if (j==10) j++;
TString label;
if (j<5) {
label = TString("TIB ")+j;
} else if (j>4&&j<11) {
label = TString("TOB ")+(j-4);
} else if (j>10&&j<14) {
label = TString("TID ")+(j-10);
} else if (j>13) {
label = TString("TEC ")+(j-13);
}
gr->GetXaxis()->SetBinLabel((j*100)/(nLayers)-2,label);
}
gr->Draw("AP");
gr->GetXaxis()->SetNdivisions(36);
TPaveText *pt = new TPaveText(0.3,0.3,0.5,0.4,"blNDC");
pt->SetBorderSize(0);
pt->SetFillColor(0);
TText *text = pt->AddText("CMS 2008");
pt->Draw("same");
c7->SaveAs("HitEffSummary69912QualGood_TECcut.png");
c7->SaveAs("HitEffSummary69912QualGood_TECcut.eps");
c7->SaveAs("HitEffSummary69912QualGood_TECcut.pdf");
}
|