ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/claudioc/ScansICHEP2012/code/B2_150_lipstick.C
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
Error occurred while calculating annotation data.
Log Message:
lipstick

File Contents

# Content
1 {
2 //=========Macro generated from canvas: c1/c1
3 //========= (Fri Apr 27 20:31:51 2012) by ROOT version5.30/02
4 TCanvas *c1 = new TCanvas("c1", "c1",0,0,700,500);
5 gStyle->SetOptStat(0);
6 gStyle->SetOptTitle(0);
7 c1->Range(210.2564,99.99999,1107.692,1100);
8 c1->SetFillColor(0);
9 c1->SetBorderMode(0);
10 c1->SetBorderSize(2);
11 c1->SetTickx(1);
12 c1->SetTicky(1);
13 // c1->SetRightMargin(0.12);
14
15 c1->SetLeftMargin(0.1407035);
16 c1->SetRightMargin(0.07);
17 c1->SetTopMargin(0.08);
18 c1->SetBottomMargin(0.13);
19 c1->SetFrameBorderMode(0);
20 c1->SetFrameBorderMode(0);
21
22
23 c1->SetFrameBorderMode(0);
24 c1->SetFrameBorderMode(0);
25
26 TH2F *empty = new TH2F("empty","Gluino Sbottom model",100,300,1000,100,200,1000);
27 empty->SetStats(0);
28 empty->GetXaxis()->SetTitle("m(#tilde{g}) GeV");
29 empty->GetXaxis()->SetTitleSize(0.060);
30 empty->GetXaxis()->SetLabelSize(0.060);
31 empty->GetYaxis()->SetTitle("m(#tilde{b}_{1}) GeV");
32 empty->GetYaxis()->SetTitleSize(0.060);
33 empty->GetYaxis()->SetTitleOffset(1.2);
34 empty->GetYaxis()->SetLabelSize(0.060);
35 empty->Draw("");
36
37
38 Double_t *dum = 0;
39 TPolyLine *pline = new TPolyLine(16,dum,dum,"");
40
41 ci = TColor::GetColor("#ff0000");
42 pline->SetFillColor(kCyan-10);
43 // pline->SetFillStyle(3244);
44 // pline->SetFillStyle(3305);
45
46 ci = TColor::GetColor("#ff0000");
47 int kk=0;
48 float mb=5.;
49 pline->SetLineColor(kCyan-10);
50 pline->SetPoint(kk++,873,325);
51 pline->SetPoint(kk++,888,400);
52 pline->SetPoint(kk++,913,450);
53 pline->SetPoint(kk++,933,500);
54 // pline->SetPoint(kk++,940,650);
55 pline->SetPoint(kk++,948,700);
56 pline->SetPoint(kk++,858,858-mb);
57 pline->SetPoint(kk++,828,828-mb);
58 pline->SetPoint(kk++,896,700);
59 pline->SetPoint(kk++,883,500);
60 pline->SetPoint(kk++,863,450);
61 // pline->SetPoint(kk++,865,400);
62 pline->SetPoint(kk++,828,325);
63 pline->SetPoint(kk++,873,325);
64 pline->Draw("fl");
65 pline->Draw();
66
67 // a line down the middle
68 Float_t x[10],y[10];
69 x[0] = (873+828)/2;
70 x[1] = (863+913)/2;
71 x[2] = (883+913)/2;
72 x[3] = (896+940)/2;
73 x[4] = (828+858)/2;
74 y[0] = 325;
75 y[1] = 450;
76 y[2] = 500;
77 y[3] = 700;
78 y[4] = x[4] - mb;
79 TGraph* gr00 = new TGraph(5,x,y);
80 gr00->SetLineColor(4);
81 gr00->SetLineWidth(2);
82 gr00->Draw();
83
84 // the expected line is shifted by about 7 gev
85 for (int i=0; i<5; i++) {
86 x[i] = x[i] + 7;
87 }
88 y[4] = x[4] - mb;
89 TGraph* gr01 = new TGraph(5,x,y);
90 gr01->SetLineColor(2);
91 gr01->SetLineWidth(2);
92 gr01->SetLineStyle(2);
93 gr01->Draw();
94
95 // the +1 sigma is 40 gev off the expected
96 for (int i=0; i<5; i++) {
97 x[i] = x[i] - 40;
98 }
99 y[4] = x[4] - mb;
100 TGraph* gr02 = new TGraph(5,x,y);
101 gr02->SetLineColor(2);
102 gr02->SetLineWidth(2);
103 gr02->SetLineStyle(3);
104 gr02->Draw();
105
106 // the -1 sigma is 65 GeV is higher
107 for (int i=0; i<5; i++) {
108 x[i] = x[i] + 65;
109 }
110 y[4] = x[4] - mb;
111 TGraph* gr03 = new TGraph(5,x,y);
112 gr03->SetLineColor(2);
113 gr03->SetLineWidth(2);
114 gr03->SetLineStyle(3);
115 gr03->Draw();
116
117
118 Double_t *dum = 0;
119 pline = new TPolyLine(16,dum,dum,"");
120
121 ci = TColor::GetColor("#0000ff");
122 pline->SetFillColor(ci);
123 // pline->SetFillStyle(3344);
124 pline->SetFillStyle(3395);
125
126 Double_t *dum2 = 0;
127 TPolyLine *pline2 = new TPolyLine(16,dum2,dum2,"");
128 pline2->SetPoint(0, 325, 876);
129 pline2->SetPoint(1, 375, 876);
130 pline2->SetPoint(2, 375, 916);
131 pline2->SetPoint(3, 325, 916);
132 pline2->SetPoint(0, 325, 876);
133 pline2->SetFillColor(kCyan-10);
134 pline2->Draw("fl");
135 TLine * blah = new TLine(325,896,375,896);
136 blah->SetLineColor(4);
137 blah->SetLineWidth(2);
138 blah->Draw();
139
140 TLine * blah2 = new TLine(325,825,375,825);
141 blah2->SetLineColor(2);
142 blah2->SetLineStyle(2);
143 blah2->SetLineWidth(2);
144 blah2->Draw();
145
146 TLine * blah3 = new TLine(325,845,375,845);
147 blah3->SetLineColor(2);
148 blah3->SetLineStyle(3);
149 blah3->SetLineWidth(2);
150 blah3->Draw();
151
152 TLine * blah4 = new TLine(325,805,375,805);
153 blah4->SetLineColor(2);
154 blah4->SetLineStyle(3);
155 blah4->SetLineWidth(2);
156 blah4->Draw();
157
158
159 tex = new TLatex(380,876,"Observed Limit #sigma^{prod} = #sigma^{NLO+NLL} #pm 1 #sigma");
160 tex->SetTextSize(0.055);
161 tex->SetLineWidth(2);
162 tex->Draw();
163
164 tex = new TLatex(380,810,"Expected Limit #pm 1 stat. #sigma");
165 tex->SetTextSize(0.055);
166 tex->SetLineWidth(2);
167 tex->Draw();
168
169
170
171 TLatex * tex = new TLatex(370,1032,"CMS Preliminary, #sqrt{s} = 8 TeV, L_{int} = 3.95 fb^{-1}");
172 tex->SetTextSize(0.055);
173 tex->SetLineWidth(2);
174 tex->Draw();
175 tex = new TLatex(330,936,"Same Sign dileptons with btag selection");
176 tex->SetTextSize(0.055);
177 tex->SetLineWidth(2);
178 tex->Draw();
179 tex = new TLatex(330,750,"m(#tilde{#chi}^{0}_{1}) = 50 GeV");
180 tex->SetTextSize(0.055);
181 tex->SetLineWidth(2);
182 tex->Draw();
183
184 tex = new TLatex(330,660,"#color[1]{m(#tilde{#chi}_{1}^{+}) = 150 GeV}");
185 tex->SetTextSize(0.055);
186 tex->SetLineWidth(2);
187 tex->Draw();
188 tex = new TLatex(540,495,"#color[4]{m(#tilde{#chi}_{1}^{+}) = 300 GeV}");
189 tex->SetTextSize(0.055);
190 tex->SetLineWidth(2);
191 // tex->Draw();
192
193
194 line = new TLine(330,325,1000,325);
195 line->SetLineStyle(1);
196 line->SetLineWidth(2);
197
198 ci = TColor::GetColor("#ff0000");
199 line->SetLineColor(1);
200 line->Draw();
201 line = new TLine(330,325,1000,995);
202 line->SetLineStyle(1);
203
204 ci = TColor::GetColor("#ff0000");
205 line->SetLineColor(1);
206 line->SetLineWidth(2);
207 line->Draw();
208
209
210
211 c1->Modified();
212 c1->cd();
213 c1->SetSelected(c1);
214 }