ScansICHEP2012/code/T1tttt_lipstick.C

{
   TCanvas *c1_n7 = new TCanvas("c1_n7", "c1_n7",1,49,700,500);
   gStyle->SetOptStat(0);
   gStyle->SetOptTitle(0);
   c1_n7->Range(258.2583,-90.03164,1176.798,890.981);
   c1_n7->SetFillColor(0);
   c1_n7->SetBorderMode(0);
   c1_n7->SetBorderSize(2);
   c1_n7->SetTickx(1);
   c1_n7->SetTicky(1);
   c1_n7->SetLeftMargin(0.1407035);
   c1_n7->SetRightMargin(0.07);
   c1_n7->SetTopMargin(0.08);
   c1_n7->SetBottomMargin(0.13);
   c1_n7->SetFrameBorderMode(0);
   c1_n7->SetFrameBorderMode(0);
   
   TH2F *empty__1 = new TH2F("empty__1","T1 model",29,387.5,1112.5,31,37.5,812.5);
   empty__1->SetDirectory(0);
   empty__1->SetStats(0);
   empty__1->GetXaxis()->SetTitle("m(#tilde{g}) GeV");
   empty__1->GetXaxis()->SetTitleSize(0.060);
   empty__1->GetYaxis()->SetTitle("m(#tilde{#chi}_{1}^{0}) GeV");
   empty__1->GetYaxis()->SetTitleSize(0.060);
   empty__1->GetXaxis()->SetLabelSize(0.060);
   empty__1->GetYaxis()->SetLabelSize(0.060);
   empty__1->GetYaxis()->SetTitleOffset(1.1);
   empty__1->GetZaxis()->SetTitleOffset(1.2);
   empty__1->Draw("");
   TLine *line = new TLine(387.5,37.5,1112.5,762.5);
   line->SetLineWidth(3);
   line->Draw();
   TLatex *   tex = new TLatex(460,843.5,"CMS Preliminary, #sqrt{s} = 8 TeV, L_{int} = 3.95 fb^{-1}");
   tex->SetTextSize(0.055);
   tex->SetLineWidth(2);
   tex->Draw();
      tex = new TLatex(420,750.5,"Same Sign dileptons with btag selection");
   tex->SetTextSize(0.055);
   tex->SetLineWidth(2);
   tex->Draw();
      tex = new TLatex(499.75,688.5,"Observed Limit #sigma^{prod} = #sigma^{NLO+NLL} #pm 1 #sigma");
   tex->SetTextSize(0.055);
   tex->SetLineWidth(2);
   tex->Draw();
   
   Double_t *dum = 0;
   TPolyLine *pline = new TPolyLine(8,dum,dum,"");

   Int_t ci;   // for color index setting
   ci = TColor::GetColor("#0000ff");
   pline->SetFillColor(kCyan-10);
   //pline->SetFillStyle(3244);

   ci = TColor::GetColor("#0000ff");
   pline->SetLineColor(ci);
   pline->SetPoint(0,420,685);
   pline->SetPoint(1,492.5,685);
   pline->SetPoint(2,492.5,725);
   pline->SetPoint(3,420,725);
   pline->SetPoint(4,420,685);
   pline->Draw("fl");
   TLine * blah = new TLine(420, 705, 492.5, 705);
   blah->SetLineColor(4);
   blah->SetLineWidth(2);
   blah->Draw();

   tex = new TLatex(499.75,630.5,"Expected Limit #pm 1 stat. #sigma");
   tex->SetTextSize(0.055);
   tex->SetLineWidth(2);
   tex->Draw();
   TLine* bb  = new TLine(420, 645, 492.5, 645);
   TLine* bb2 = new TLine(420, 630, 492.5, 630);
   TLine* bb3 = new TLine(420, 660, 492.5, 660);
   bb->SetLineColor(2);
   bb->SetLineWidth(2);
   bb->SetLineStyle(2);
   bb2->SetLineColor(2);
   bb2->SetLineWidth(2);
   bb2->SetLineStyle(3);
   bb3->SetLineColor(2);
   bb3->SetLineWidth(2);
   bb3->SetLineStyle(3);
   bb->Draw();
   bb2->Draw();
   bb3->Draw();

   
   Double_t *dum = 0;
   pline = new TPolyLine(16,dum,dum,"");

   ci = TColor::GetColor("#0000ff");
   pline->SetFillColor(kCyan-10);
   //pline->SetFillStyle(3244);

   ci = TColor::GetColor("#0000ff");
   pline->SetLineColor(ci);
   int kk=0;
   float mtop=175.;
   pline->SetPoint(kk++,869,37.5);
   pline->SetPoint(kk++,869,100);
   pline->SetPoint(kk++,860,200);
   pline->SetPoint(kk++,860,250);
   pline->SetPoint(kk++,856,300);
   pline->SetPoint(kk++,834,350);
   pline->SetPoint(kk++,820,400);
   pline->SetPoint(kk++,800,800.-2*mtop);
   pline->SetPoint(kk++,820,820-2*mtop);
   pline->SetPoint(kk++,830,450);
   pline->SetPoint(kk++,860,400);
   pline->SetPoint(kk++,895,350);
   pline->SetPoint(kk++,900,300);
   pline->SetPoint(kk++,909,200);
   pline->SetPoint(kk++,916,100);
   pline->SetPoint(kk++,916,37.5);
   pline->SetPoint(kk++,869,37.5);
   pline->Draw("fl");

   // This is the observed line
   Double_t x[7],y[7];
   x[0] = (869+916)/2.;
   y[0] = 37.5;
   x[1] = (869+916)/2.;
   y[1] = 100.;
   x[2] = (909+860)/2.;
   y[2] = 200;
   x[3] = (900+856)/2.;
   y[3] = 300;
   x[4] = (895+834)/2.;
   y[4] = 350;
   x[5] = (860+820)/2.;
   y[5] = 400;
   x[6] = 810.;
   y[6] = 810-2*mtop;
   gr = new TGraph(7,x,y);
   gr->SetLineColor(4);
   gr->SetLineWidth(2);
   gr->Draw();

   // The expected line is about 7 GeV better than the observed line, 
   // except at the top where the SR switches
   for (int i=0; i<4; i++) {
     x[i] = x[i] + 6.;
   }
   x[5] = x[5] - 5.;
   x[6] = x[6] - 5.;
   y[6] = x[6] - 2*mtop;
   gr2 = new TGraph(7,x,y);
   gr2->SetLineColor(2);
   gr2->SetLineWidth(2);
   gr2->SetLineStyle(2);
   gr2->Draw();
   
   // -1 sigma...
   x[0] = x[0] - 40.;
   x[1] = x[1] - 40.;
   x[2] = x[2] - 40.;
   x[3] = x[3] - 50.;
   x[4] = x[4] - 55.;
   x[5] = 780.;
   y[5] = x[5] - 2*mtop;
   gr3 = new TGraph(6,x,y);
   gr3->SetLineColor(2);
   gr3->SetLineWidth(2);
   gr3->SetLineStyle(3);
   gr3->Draw();

   // +1 sigma...
   x[0] = 925.;
   x[1] = 925.;
   x[2] = 921;
   x[3] = 905;
   x[4] = 905;
   x[5] = 890;
   x[6] = 860;
   y[6] = 860 - 2*mtop;
   gr4 = new TGraph(7,x,y);
   gr4->SetLineColor(2);
   gr4->SetLineWidth(2);
   gr4->SetLineStyle(3);
   gr4->Draw();


   TPaveText *pt = new TPaveText(0.01,0.945,0.167931,0.995,"blNDC");
   pt->SetName("title");
   pt->SetBorderSize(1);
   pt->SetFillColor(0);
   text = pt->AddText("T1 model");
   pt->Draw();
   c1_n7->Modified();
   c1_n7->cd();
   c1_n7->SetSelected(c1_n7);
}
Revision:	1.1
Committed:	Mon Jul 2 04:42:05 2012 UTC (12 years, 10 months ago) by claudioc
Content type:	text/plain
Branch:	MAIN
CVS Tags:	HEAD
Log Message:	lipstick
#	User	Rev	Content
1	claudioc	1.1	{
2			TCanvas *c1_n7 = new TCanvas("c1_n7", "c1_n7",1,49,700,500);
3			gStyle->SetOptStat(0);
4			gStyle->SetOptTitle(0);
5			c1_n7->Range(258.2583,-90.03164,1176.798,890.981);
6			c1_n7->SetFillColor(0);
7			c1_n7->SetBorderMode(0);
8			c1_n7->SetBorderSize(2);
9			c1_n7->SetTickx(1);
10			c1_n7->SetTicky(1);
11			c1_n7->SetLeftMargin(0.1407035);
12			c1_n7->SetRightMargin(0.07);
13			c1_n7->SetTopMargin(0.08);
14			c1_n7->SetBottomMargin(0.13);
15			c1_n7->SetFrameBorderMode(0);
16			c1_n7->SetFrameBorderMode(0);
17
18			TH2F *empty__1 = new TH2F("empty__1","T1 model",29,387.5,1112.5,31,37.5,812.5);
19			empty__1->SetDirectory(0);
20			empty__1->SetStats(0);
21			empty__1->GetXaxis()->SetTitle("m(#tilde{g}) GeV");
22			empty__1->GetXaxis()->SetTitleSize(0.060);
23			empty__1->GetYaxis()->SetTitle("m(#tilde{#chi}_{1}^{0}) GeV");
24			empty__1->GetYaxis()->SetTitleSize(0.060);
25			empty__1->GetXaxis()->SetLabelSize(0.060);
26			empty__1->GetYaxis()->SetLabelSize(0.060);
27			empty__1->GetYaxis()->SetTitleOffset(1.1);
28			empty__1->GetZaxis()->SetTitleOffset(1.2);
29			empty__1->Draw("");
30			TLine *line = new TLine(387.5,37.5,1112.5,762.5);
31			line->SetLineWidth(3);
32			line->Draw();
33			TLatex * tex = new TLatex(460,843.5,"CMS Preliminary, #sqrt{s} = 8 TeV, L_{int} = 3.95 fb^{-1}");
34			tex->SetTextSize(0.055);
35			tex->SetLineWidth(2);
36			tex->Draw();
37			tex = new TLatex(420,750.5,"Same Sign dileptons with btag selection");
38			tex->SetTextSize(0.055);
39			tex->SetLineWidth(2);
40			tex->Draw();
41			tex = new TLatex(499.75,688.5,"Observed Limit #sigma^{prod} = #sigma^{NLO+NLL} #pm 1 #sigma");
42			tex->SetTextSize(0.055);
43			tex->SetLineWidth(2);
44			tex->Draw();
45
46			Double_t *dum = 0;
47			TPolyLine *pline = new TPolyLine(8,dum,dum,"");
48
49			Int_t ci; // for color index setting
50			ci = TColor::GetColor("#0000ff");
51			pline->SetFillColor(kCyan-10);
52			//pline->SetFillStyle(3244);
53
54			ci = TColor::GetColor("#0000ff");
55			pline->SetLineColor(ci);
56			pline->SetPoint(0,420,685);
57			pline->SetPoint(1,492.5,685);
58			pline->SetPoint(2,492.5,725);
59			pline->SetPoint(3,420,725);
60			pline->SetPoint(4,420,685);
61			pline->Draw("fl");
62			TLine * blah = new TLine(420, 705, 492.5, 705);
63			blah->SetLineColor(4);
64			blah->SetLineWidth(2);
65			blah->Draw();
66
67			tex = new TLatex(499.75,630.5,"Expected Limit #pm 1 stat. #sigma");
68			tex->SetTextSize(0.055);
69			tex->SetLineWidth(2);
70			tex->Draw();
71			TLine* bb = new TLine(420, 645, 492.5, 645);
72			TLine* bb2 = new TLine(420, 630, 492.5, 630);
73			TLine* bb3 = new TLine(420, 660, 492.5, 660);
74			bb->SetLineColor(2);
75			bb->SetLineWidth(2);
76			bb->SetLineStyle(2);
77			bb2->SetLineColor(2);
78			bb2->SetLineWidth(2);
79			bb2->SetLineStyle(3);
80			bb3->SetLineColor(2);
81			bb3->SetLineWidth(2);
82			bb3->SetLineStyle(3);
83			bb->Draw();
84			bb2->Draw();
85			bb3->Draw();
86
87
88			Double_t *dum = 0;
89			pline = new TPolyLine(16,dum,dum,"");
90
91			ci = TColor::GetColor("#0000ff");
92			pline->SetFillColor(kCyan-10);
93			//pline->SetFillStyle(3244);
94
95			ci = TColor::GetColor("#0000ff");
96			pline->SetLineColor(ci);
97			int kk=0;
98			float mtop=175.;
99			pline->SetPoint(kk++,869,37.5);
100			pline->SetPoint(kk++,869,100);
101			pline->SetPoint(kk++,860,200);
102			pline->SetPoint(kk++,860,250);
103			pline->SetPoint(kk++,856,300);
104			pline->SetPoint(kk++,834,350);
105			pline->SetPoint(kk++,820,400);
106			pline->SetPoint(kk++,800,800.-2*mtop);
107			pline->SetPoint(kk++,820,820-2*mtop);
108			pline->SetPoint(kk++,830,450);
109			pline->SetPoint(kk++,860,400);
110			pline->SetPoint(kk++,895,350);
111			pline->SetPoint(kk++,900,300);
112			pline->SetPoint(kk++,909,200);
113			pline->SetPoint(kk++,916,100);
114			pline->SetPoint(kk++,916,37.5);
115			pline->SetPoint(kk++,869,37.5);
116			pline->Draw("fl");
117
118			// This is the observed line
119			Double_t x[7],y[7];
120			x[0] = (869+916)/2.;
121			y[0] = 37.5;
122			x[1] = (869+916)/2.;
123			y[1] = 100.;
124			x[2] = (909+860)/2.;
125			y[2] = 200;
126			x[3] = (900+856)/2.;
127			y[3] = 300;
128			x[4] = (895+834)/2.;
129			y[4] = 350;
130			x[5] = (860+820)/2.;
131			y[5] = 400;
132			x[6] = 810.;
133			y[6] = 810-2*mtop;
134			gr = new TGraph(7,x,y);
135			gr->SetLineColor(4);
136			gr->SetLineWidth(2);
137			gr->Draw();
138
139			// The expected line is about 7 GeV better than the observed line,
140			// except at the top where the SR switches
141			for (int i=0; i<4; i++) {
142			x[i] = x[i] + 6.;
143			}
144			x[5] = x[5] - 5.;
145			x[6] = x[6] - 5.;
146			y[6] = x[6] - 2*mtop;
147			gr2 = new TGraph(7,x,y);
148			gr2->SetLineColor(2);
149			gr2->SetLineWidth(2);
150			gr2->SetLineStyle(2);
151			gr2->Draw();
152
153			// -1 sigma...
154			x[0] = x[0] - 40.;
155			x[1] = x[1] - 40.;
156			x[2] = x[2] - 40.;
157			x[3] = x[3] - 50.;
158			x[4] = x[4] - 55.;
159			x[5] = 780.;
160			y[5] = x[5] - 2*mtop;
161			gr3 = new TGraph(6,x,y);
162			gr3->SetLineColor(2);
163			gr3->SetLineWidth(2);
164			gr3->SetLineStyle(3);
165			gr3->Draw();
166
167			// +1 sigma...
168			x[0] = 925.;
169			x[1] = 925.;
170			x[2] = 921;
171			x[3] = 905;
172			x[4] = 905;
173			x[5] = 890;
174			x[6] = 860;
175			y[6] = 860 - 2*mtop;
176			gr4 = new TGraph(7,x,y);
177			gr4->SetLineColor(2);
178			gr4->SetLineWidth(2);
179			gr4->SetLineStyle(3);
180			gr4->Draw();
181
182
183			TPaveText *pt = new TPaveText(0.01,0.945,0.167931,0.995,"blNDC");
184			pt->SetName("title");
185			pt->SetBorderSize(1);
186			pt->SetFillColor(0);
187			text = pt->AddText("T1 model");
188			pt->Draw();
189			c1_n7->Modified();
190			c1_n7->cd();
191			c1_n7->SetSelected(c1_n7);
192			}