ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/claudioc/OSNote2010/plotTrgEff/plotTrgEff.C
Revision: 1.2
Committed: Tue Nov 9 18:50:45 2010 UTC (14 years, 5 months ago) by claudioc
Content type: text/plain
Branch: MAIN
CVS Tags: FR1, v4, v3, HEAD
Changes since 1.1: +30 -4 lines
Error occurred while calculating annotation data.
Log Message: 
trigger model changes

File Contents

# Content
1 {
2 gROOT->LoadMacro("triggerSuperModel.cc+");
3 gStyle->SetOptStat(0);
4 gStyle->SetPaintTextFormat(".3f");
5
6
7 //----------------------
8 // plot mu mu first
9 //----------------------
10 TCanvas* cnv = new TCanvas();
11 cnv->Divide(1,2);
12 double mubin[5]={10.,11.,15.,20.,30.};
13 int nmubin = 4;
14 TH2F* mu21 = new TH2F("mu21","mumu - both muons with abs(eta)<2.1", nmubin, mubin, nmubin, mubin);
15 TH2F* mu24 = new TH2F("mu24","mumu - only one muon with abs(eta)<2.1", nmubin, mubin, nmubin, mubin);
16 TH2F* mu24b = new TH2F("mu24b","mumu - both muons with abs(eta)>2.1", nmubin, mubin, nmubin, mubin);
17
18
19
20 int id1=13;
21 int id2=13;
22 int hyp = 0;
23 float eta1=0.;
24 float eta2=0.;
25 for (int i=0; i<=nmubin; i++) {
26 float pt1 = mubin[i]+0.1;
27 for (int j=0; j<=nmubin; j++) {
28 float pt2 = mubin[j]+0.1;
29 if (TMath::Max(pt1,pt2) > 20.) {
30 float eff21 = triggerSuperModelEffic(hyp, pt1, eta1, id1, pt2, eta2, id2);
31 float eff24 = triggerSuperModelEffic(hyp, pt1, eta1, id1, pt2, 2.5, id2);
32 float eff24b = triggerSuperModelEffic(hyp, pt1, 2.5, id1, pt2, 2.5, id2);
33 mu21->Fill(pt1,pt2,eff21);
34 mu24->Fill(pt1,pt2,eff24);
35 mu24b->Fill(pt1,pt2,eff24b);
36 }
37 }
38 }
39 TLine * l1 = new TLine(11.,10.,11.,30.);
40 TLine * l2 = new TLine(15.,10.,15.,30.);
41 TLine * l3 = new TLine(20.,10.,20.,30.);
42 TLine * l4 = new TLine(10.,11.,30.,11.);
43 TLine * l5 = new TLine(10.,15.,30.,15.);
44 TLine * l6 = new TLine(10.,20.,30.,20.);
45 cnv->cd(1);
46 mu21->SetMarkerSize(2.5);
47 mu21->Draw("text");
48 mu21->GetYaxis()->SetTickLength(0);
49 mu21->GetXaxis()->SetTickLength(0);
50 mu21->GetXaxis()->SetTitle("Pt muon 1");
51 mu21->GetYaxis()->SetTitle("Pt muon 2");
52 mu21->GetXaxis()->SetTitleSize(0.05);
53 mu21->GetYaxis()->SetTitleSize(0.05);
54 l1->Draw();
55 l2->Draw();
56 l3->Draw();
57 l4->Draw();
58 l5->Draw();
59 l6->Draw();
60
61 cnv->cd(2);
62 mu24->SetMarkerSize(2.5);
63 mu24->Draw("text");
64 mu24->GetYaxis()->SetTickLength(0);
65 mu24->GetXaxis()->SetTickLength(0);
66 mu24->GetXaxis()->SetTitle("Pt muon of abs(eta)<2.1");
67 mu24->GetYaxis()->SetTitle("Pt muon of abs(eta)>2.1");
68 mu24->GetXaxis()->SetTitleSize(0.05);
69 mu24->GetYaxis()->SetTitleSize(0.05);
70 mu24->Draw("text");
71 mu24->GetYaxis()->SetTickLength(0);
72 mu24->GetXaxis()->SetTickLength(0);
73 l1->Draw();
74 l2->Draw();
75 l3->Draw();
76 l4->Draw();
77 l5->Draw();
78 l6->Draw();
79
80 TCanvas* cnv3 = new TCanvas();
81 cnv3->Divide(1,2);
82 cnv3->cd(1);
83 mu24b->SetMarkerSize(2.5);
84 mu24b->Draw("text");
85 mu24b->GetYaxis()->SetTickLength(0);
86 mu24b->GetXaxis()->SetTickLength(0);
87 mu24b->GetXaxis()->SetTitle("Pt muon");
88 mu24b->GetYaxis()->SetTitle("Pt muon");
89 mu24b->GetXaxis()->SetTitleSize(0.05);
90 mu24b->GetYaxis()->SetTitleSize(0.05);
91 mu24b->Draw("text");
92 mu24b->GetYaxis()->SetTickLength(0);
93 mu24b->GetXaxis()->SetTickLength(0);
94 l1->Draw();
95 l2->Draw();
96 l3->Draw();
97 l4->Draw();
98 l5->Draw();
99 l6->Draw();
100
101
102 //---------------
103 // Now plot ee
104 //---------------
105 TCanvas* cnv2 = new TCanvas();
106 cnv2->Divide(1,2);
107 double ebin[8]={10.,15.,17.,20., 22.,27., 32., 40.};
108 int nebin = 7;
109 TH2F* e21 = new TH2F("e21","emu - muon with abs(eta)<2.1", nebin, ebin, nmubin, mubin);
110 TH2F* e24 = new TH2F("e24","emu - muon with abs(eta)>2.1", nebin, ebin, nmubin, mubin);
111
112
113 int id1=11;
114 int id2=13;
115 int hyp = 1;
116 float eta1=0.;
117 float eta2=0.;
118 for (int i=0; i<=nebin; i++) {
119 float pt1 = ebin[i]+0.1;
120 for (int j=0; j<=nmubin; j++) {
121 float pt2 = mubin[j]+0.1;
122 if (TMath::Max(pt1,pt2) > 20.) {
123 float eff21 = triggerSuperModelEffic(hyp, pt1, eta1, id1, pt2, eta2, id2);
124 float eff24 = triggerSuperModelEffic(hyp, pt1, eta1, id1, pt2, 2.5, id2);
125 e21->Fill(pt1,pt2,eff21);
126 e24->Fill(pt1,pt2,eff24);
127 }
128 }
129 }
130
131
132 TLine * s1 = new TLine(15.,10.,15.,30.);
133 TLine * s2 = new TLine(17.,10.,17.,30.);
134 TLine * s3 = new TLine(20.,10.,20.,30.);
135 TLine * s4 = new TLine(22.,10.,22.,30.);
136 TLine * s5 = new TLine(27.,10.,27.,30.);
137 TLine * s6 = new TLine(32.,10.,32.,30.);
138 TLine * t1 = new TLine(10.,11.,40.,11.);
139 TLine * t2 = new TLine(10.,15.,40.,15.);
140 TLine * t3 = new TLine(10.,20.,40.,20.);
141
142 cnv2->cd(1);
143 e21->SetMarkerSize(2.5);
144 e21->Draw("text");
145 e21->GetYaxis()->SetTickLength(0);
146 e21->GetXaxis()->SetTickLength(0);
147 e21->GetXaxis()->SetTitle("Pt electron");
148 e21->GetYaxis()->SetTitle("Pt muon");
149 e21->GetXaxis()->SetTitleSize(0.05);
150 e21->GetYaxis()->SetTitleSize(0.05);
151 s1->Draw();
152 s2->Draw();
153 s3->Draw();
154 s4->Draw();
155 s5->Draw();
156 s6->Draw();
157 t1->Draw();
158 t2->Draw();
159 t3->Draw();
160
161
162 cnv2->cd(2);
163 e24->SetMarkerSize(2.5);
164 e24->Draw("text");
165 e24->GetYaxis()->SetTickLength(0);
166 e24->GetXaxis()->SetTickLength(0);
167 e24->GetXaxis()->SetTitle("Pt electron");
168 e24->GetYaxis()->SetTitle("Pt muon");
169 e24->GetXaxis()->SetTitleSize(0.05);
170 e24->GetYaxis()->SetTitleSize(0.05);
171 s1->Draw();
172 s2->Draw();
173 s3->Draw();
174 s4->Draw();
175 s5->Draw();
176 s6->Draw();
177 t1->Draw();
178 t2->Draw();
179 t3->Draw();
180
181
182 }
183