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
Log Message: 
trigger model changes

File Contents

# User Rev Content
1 claudioc 1.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 claudioc 1.2 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 claudioc 1.1
18 claudioc 1.2
19    
20 claudioc 1.1 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 claudioc 1.2 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 claudioc 1.1 mu21->Fill(pt1,pt2,eff21);
34     mu24->Fill(pt1,pt2,eff24);
35 claudioc 1.2 mu24b->Fill(pt1,pt2,eff24b);
36 claudioc 1.1 }
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 claudioc 1.2 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 claudioc 1.1
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