| 4 |
|
#include "PlotTools.h" |
| 5 |
|
#include "TheLimits.h" |
| 6 |
|
#include "GlobalFunctions.h" |
| 7 |
+ |
#include "StyleSettings.h" |
| 8 |
|
|
| 9 |
|
#include "TRint.h" |
| 10 |
|
#include "TROOT.h" |
| 38 |
|
return 1; |
| 39 |
|
} |
| 40 |
|
|
| 41 |
+ |
/* |
| 42 |
|
//Style stuff |
| 43 |
|
gStyle->SetHistFillColor(0); |
| 44 |
|
gStyle->SetPalette(1); |
| 70 |
|
gStyle->SetLabelSize(0.03,"X"); |
| 71 |
|
gStyle->SetLabelSize(0.03,"Y"); |
| 72 |
|
gStyle->SetLabelSize(0.03,"Z"); |
| 73 |
+ |
*/ |
| 74 |
+ |
util::StyleSettings::paper(); |
| 75 |
+ |
gStyle->SetPadBottomMargin(0.18); |
| 76 |
|
|
| 77 |
|
//gROOT->SetStyle("MyStyle"); |
| 78 |
|
TCanvas * c1 = new TCanvas("c1","c1",600,600); |
| 79 |
< |
c1->SetFillStyle ( 4000 ); |
| 80 |
< |
c1->SetLeftMargin ( 0.15 ); |
| 81 |
< |
c1->SetRightMargin ( 0.15 ); |
| 82 |
< |
c1->SetBottomMargin( 0.10 ); |
| 79 |
> |
//c1->SetFillStyle ( 4000 ); |
| 80 |
> |
//c1->SetLeftMargin ( 0.15 ); |
| 81 |
> |
//c1->SetRightMargin ( 0.15 ); |
| 82 |
> |
//c1->SetBottomMargin( 0.10 ); |
| 83 |
|
c1->cd(); |
| 84 |
|
|
| 85 |
|
//Get limits and generator masses --------------------------------------------------- |
| 86 |
|
TheLimits * genpoints = new TheLimits(); |
| 87 |
|
//genpoints->Fill(argc, argv); |
| 88 |
< |
genpoints->Fill("limits_MHT_tb50/filelist.txt"); |
| 88 |
> |
genpoints->Fill("limits_MHT_tb50-20110301/filelist.txt"); |
| 89 |
|
|
| 90 |
|
TheLimits * genpointsHT = new TheLimits(); |
| 91 |
< |
genpointsHT->Fill("limits_HT_tb50/filelist.txt"); |
| 91 |
> |
genpointsHT->Fill("limits_HT_tb50-20110301/filelist.txt"); |
| 92 |
|
|
| 93 |
|
|
| 94 |
|
//Replace read limits with specific numbers |
| 206 |
|
TH2F*texcl = new TH2F("texcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion", |
| 207 |
|
40,200,600,26,140,400); |
| 208 |
|
TH2F*texpexcl=(TH2F*)texcl->Clone(); |
| 209 |
< |
plotTools->Area(texpexcl, Mzero, Mhalf, MCMCExpExclusion); |
| 205 |
< |
TH2F*tobsexcl=(TH2F*)texcl->Clone(); |
| 206 |
< |
plotTools->Area(tobsexcl, Mzero, Mhalf, MCMCObsExclusion); |
| 209 |
> |
plotToolsHT->Area(texpexcl, Mzero, Mhalf, NLOExpExclCLp1sigma); |
| 210 |
|
std::vector<TGraph*> contours = plotTools->GetContours(texpexcl,3); |
| 208 |
– |
//std::vector<TGraph*> contours = plotTools->GetContours(tobsexcl,3); |
| 209 |
– |
//hPLexpexcl |
| 210 |
– |
//hexcl->Draw("colz"); |
| 211 |
– |
//hexpexcl->Draw("colz"); |
| 212 |
– |
//hobsexcl->Draw("colz"); |
| 211 |
|
texpexcl->Draw("colz"); |
| 212 |
|
int col=kBlue-10; |
| 213 |
|
for (std::vector<TGraph*>::iterator cont=contours.begin(); cont!=contours.end(); ++cont){ |
| 227 |
|
|
| 228 |
|
// Exclusion in M0 - M1/2 |
| 229 |
|
TH2F*hexcl = new TH2F("hexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion", |
| 230 |
< |
40,200,600,26,140,400); |
| 231 |
< |
TGraph * gexpexcl = plotTools ->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCL, 3,0, 2,2); |
| 232 |
< |
TGraph * gexpexcl_LO = plotTools ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL, 3,0, 2,4); |
| 233 |
< |
TGraph * gobsexcl = plotTools ->GetContour(hexcl,Mzero,Mhalf,NLOObsExclCL, 3,0, 2,1); |
| 234 |
< |
TGraph * gMCMCexpexcl_LO = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,MCMCExpExclusion,3,0, 3,4); |
| 235 |
< |
TGraph * gMCMCobsexcl_LO = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,MCMCObsExclusion,3,0, 3,1); |
| 236 |
< |
TGraph * gCLsExpExclHT_LO = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL, 3,0, 1,4); |
| 237 |
< |
TGraph * gCLsExpExclHT = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCL, 3,0, 1,2); |
| 238 |
< |
TGraph * gCLsObsExclHT = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,NLOObsExclCL, 3,0, 1,1); |
| 239 |
< |
TGraph * gCLsExpExclHTm1 = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCLm1sigma,3,0, 5,2); |
| 240 |
< |
TGraph * gCLsExpExclHTp1 = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCLp1sigma,3,0, 5,2); |
| 241 |
< |
Smooth( gCLsExpExclHT_LO, 33 ); gCLsExpExclHT_LO->SetLineWidth( 3 ); |
| 230 |
> |
40,200,600,26,140,450); |
| 231 |
> |
TH2F*hs = new TH2F("hs","",40,200,600,26,140,400); |
| 232 |
> |
TGraph * gexpexcl = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCL, 3,0, 2,2); |
| 233 |
> |
TGraph * gexpexcl_LO = plotTools ->GetContour(hs,Mzero,Mhalf,ExpExclCL, 3,0, 2,4); |
| 234 |
> |
TGraph * gobsexcl = plotTools ->GetContour(hs,Mzero,Mhalf,NLOObsExclCL, 3,0, 2,1); |
| 235 |
> |
TGraph * gMCMCexpexcl = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCExpExclusion, 3,0, 3,2); |
| 236 |
> |
TGraph * gMCMCobsexcl = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCObsExclusion, 3,0, 3,1); |
| 237 |
> |
TGraph * gFCexpexclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion, 3,0, 3,2); |
| 238 |
> |
TGraph * gFCexpexclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,FCExpExclusion, 3,0, 3,4); |
| 239 |
> |
TGraph * gFCexpexclMHT = plotTools->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion, 3,0, 3,2); |
| 240 |
> |
TGraph * gFCexpexclMHT_LO = plotTools->GetContour(hs,Mzero,Mhalf,FCExpExclusion, 3,0, 3,4); |
| 241 |
> |
TGraph * gFCobsexclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion, 3,0, 3,1); |
| 242 |
> |
TGraph * gFCobsexclMHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion, 3,0, 3,1); |
| 243 |
> |
TGraph * gCLsExpExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCL, 3,0, 1,2); |
| 244 |
> |
TGraph * gCLsObsExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOObsExclCL, 3,0, 1,1); |
| 245 |
> |
TGraph * gCLsExpExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ExpExclCL, 3,0, 1,4); |
| 246 |
> |
TGraph * gCLsObsExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ObsExclCL, 3,0, 1,1); |
| 247 |
> |
TGraph * gCLsExpExclHTm1 = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma,3,0, 5,2); |
| 248 |
> |
TGraph * gCLsExpExclHTp1 = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma,3,0, 5,2); |
| 249 |
> |
TGraph * gCLsExpExclMHTm1 = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma, 3,0, 5,2); |
| 250 |
> |
TGraph * gCLsExpExclMHTp1 = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma, 3,0, 5,2); |
| 251 |
> |
TGraph * gCheck1 = (TGraph*)gobsexcl->Clone(); |
| 252 |
> |
TGraph * gCheck2 = (TGraph*)gexpexcl->Clone(); |
| 253 |
|
Smooth( gCLsObsExclHT, 33); gCLsObsExclHT->SetLineWidth( 3 ); |
| 254 |
|
Smooth( gCLsExpExclHT, 33 ); gCLsExpExclHT->SetLineWidth( 3 ); |
| 255 |
+ |
Smooth( gCLsObsExclHT_LO, 33); gCLsObsExclHT_LO->SetLineWidth( 3 ); |
| 256 |
+ |
Smooth( gCLsExpExclHT_LO, 33 ); gCLsExpExclHT_LO->SetLineWidth( 3 ); |
| 257 |
|
Smooth( gexpexcl, 33 ); |
| 258 |
|
Smooth( gobsexcl, 33 ); |
| 259 |
< |
Smooth( gMCMCexpexcl_LO, 33 ); |
| 259 |
> |
Smooth( gobsexcl, 33 ); |
| 260 |
> |
Smooth( gexpexcl_LO, 33 ); |
| 261 |
> |
Smooth( gFCexpexclHT, 33 ); |
| 262 |
> |
Smooth( gFCexpexclHT_LO, 33 );gFCexpexclHT_LO->SetLineWidth( 3 ); |
| 263 |
> |
Smooth( gFCexpexclMHT, 33 ); |
| 264 |
> |
Smooth( gFCexpexclMHT_LO, 33 ); |
| 265 |
> |
Smooth( gFCobsexclHT, 33 );gFCobsexclHT->SetLineWidth( 3 ); |
| 266 |
> |
Smooth( gFCobsexclMHT, 33 ); |
| 267 |
|
Smooth( gCLsExpExclHTm1, 33 ); |
| 268 |
|
Smooth( gCLsExpExclHTp1, 33 ); |
| 269 |
< |
TGraph * gCLsExp1Sigma = MakeBand(gCLsExpExclHTp1, gCLsExpExclHTm1); |
| 270 |
< |
hexcl->Draw("colz"); |
| 271 |
< |
gCLsExp1Sigma->Draw("f"); |
| 269 |
> |
Smooth( gCLsExpExclMHTm1, 33 ); |
| 270 |
> |
Smooth( gCLsExpExclMHTp1, 33 ); |
| 271 |
> |
TGraph * CLsObsNLO = plotToolsHT->ChooseBest(gCLsObsExclHT,gobsexcl, gCLsObsExclHT,gobsexcl); |
| 272 |
> |
TGraph * CLsExpNLO = plotToolsHT->ChooseBest(gCLsExpExclHT,gexpexcl, gCLsExpExclHT,gexpexcl); |
| 273 |
> |
TGraph * CLsExpLO = plotToolsHT->ChooseBest(gCLsExpExclHT_LO,gexpexcl_LO, gCLsExpExclHT_LO,gexpexcl_LO); |
| 274 |
> |
TGraph * FCExpLO = plotToolsHT->ChooseBest(gFCexpexclHT_LO,gFCexpexclMHT_LO,gFCexpexclHT_LO,gFCexpexclMHT_LO); |
| 275 |
> |
TGraph * FCObsNLO = plotToolsHT->ChooseBest(gFCobsexclHT,gFCobsexclMHT,gFCobsexclHT,gFCobsexclMHT); |
| 276 |
> |
TGraph * gCLsExpExclp1 = plotToolsHT->ChooseBest(gCLsExpExclHTp1,gCLsExpExclMHTp1, gCLsExpExclHTp1,gCLsExpExclMHTp1); |
| 277 |
> |
TGraph * gCLsExpExclm1 = plotToolsHT->ChooseBest(gCLsExpExclHTm1,gCLsExpExclMHTm1, gCLsExpExclHTm1,gCLsExpExclMHTm1); |
| 278 |
> |
TGraph * gCLsExp1Sigma = MakeBand(gCLsExpExclp1, gCLsExpExclm1); |
| 279 |
> |
hexcl->GetYaxis()->SetTitleOffset(1.3); |
| 280 |
> |
hexcl->Draw("colz"); |
| 281 |
> |
//gCLsExpExclHTp1->Draw("l"); |
| 282 |
> |
//gCLsExpExclMHTp1->Draw("l"); |
| 283 |
> |
//gCLsExpExclHTm1->Draw("l"); |
| 284 |
> |
//gCLsExpExclMHTm1->Draw("l"); |
| 285 |
> |
gCLsExp1Sigma->Draw("f"); |
| 286 |
|
//set old exclusion Limits |
| 287 |
|
TGraph* LEP_ch = set_lep_ch(50); |
| 288 |
|
TGraph* LEP_sl = set_lep_sl(50);//slepton curve |
| 321 |
|
sFirst->Draw("same"); |
| 322 |
|
sSecond->Draw("same"); |
| 323 |
|
sThird->Draw("same"); |
| 292 |
– |
if (gexpexcl) gexpexcl->Draw("l"); |
| 293 |
– |
if (gCLsExpExclHT) gCLsExpExclHT->Draw("l"); |
| 294 |
– |
if (gCLsExpExclHT_LO) gCLsExpExclHT_LO->Draw("l"); |
| 295 |
– |
if (gMCMCexpexcl_LO) gMCMCexpexcl_LO->Draw("l"); |
| 296 |
– |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l"); |
| 324 |
|
TLatex b; b.SetTextSize(0.02); b.SetTextColor(1); |
| 325 |
|
//b.DrawLatex( 10,420,"#tilde{#tau} LSP"); |
| 326 |
|
b.DrawLatex(230,145,"LEP2"); |
| 329 |
|
//b.DrawLatex(100,140,"CDF"); |
| 330 |
|
//b.DrawLatex(230,145,"D0 #tilde{g}, #tilde{q}"); |
| 331 |
|
//b.DrawLatex( 80,180,"D0 #tilde{#nu}"); |
| 332 |
< |
TLegend * leg = new TLegend(0.45,0.7,0.85,0.89); |
| 333 |
< |
leg->SetBorderSize(0);leg->SetFillColor(0); |
| 334 |
< |
TGraph * expLeg = (TGraph*)gCLsExpExclHT->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor()); |
| 332 |
> |
CLsObsNLO->Draw("l"); |
| 333 |
> |
CLsExpNLO->Draw("l"); |
| 334 |
> |
CLsExpLO->Draw("l"); |
| 335 |
> |
FCExpLO->Draw("l"); |
| 336 |
> |
//FCObsNLO->Draw("l"); |
| 337 |
> |
|
| 338 |
> |
//gl500->Draw("c"); gl600->Draw("c"); gl700->Draw("c"); gl800->Draw("c"); gl900->Draw("c"); gl1000->Draw("c"); |
| 339 |
> |
//sq500->Draw("c"); sq600->Draw("c"); sq700->Draw("c"); sq800->Draw("c"); sq900->Draw("c"); sq1000->Draw("c"); |
| 340 |
> |
//TLegend * leg = new TLegend(0.45,0.65,0.84,0.89); |
| 341 |
> |
//leg->SetBorderSize(0);leg->SetFillColor(0); |
| 342 |
> |
TLegend * leg = new TLegend(0.6,0.66,0.96,0.90); |
| 343 |
> |
leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42); |
| 344 |
> |
TGraph * expLeg = (TGraph*)CLsExpNLO->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor()); |
| 345 |
> |
leg->SetHeader("CMS preliminary, 95%CL"); |
| 346 |
> |
leg->AddEntry(expLeg, "Expected (CLs, NLO)","l"); |
| 347 |
> |
leg->AddEntry(CLsObsNLO,"Observed (CLs, NLO)","l"); |
| 348 |
> |
leg->AddEntry(CLsExpLO, "Expected (CLs, LO)","l"); |
| 349 |
> |
leg->AddEntry(FCExpLO, "Expected (FC, LO)","l"); |
| 350 |
|
//if (gobsexcl) leg->AddEntry(gobsexcl,"Observed (MHT, CLs, TLimit)","l"); |
| 351 |
< |
if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT, "NLO Observed (HT, CLs)","l"); |
| 352 |
< |
if (gCLsExpExclHT) leg->AddEntry(expLeg , "NLO Expected (HT, CLs)","l"); |
| 353 |
< |
if (gCLsExpExclHT_LO) leg->AddEntry(gCLsExpExclHT_LO,"LO Expected (HT, CLs)","l"); |
| 354 |
< |
if (gMCMCexpexcl_LO) leg->AddEntry(gMCMCexpexcl_LO, "LO Expected (HT, MCMC, RooStat)","l"); |
| 355 |
< |
if (gexpexcl) leg->AddEntry(gexpexcl, "NLO Expected (MHT, CLs)","l"); |
| 356 |
< |
if (sSecond) leg->AddEntry(sSecond, "NLO #alpha_{T} Expected"); |
| 357 |
< |
if (sFirst) leg->AddEntry(sFirst, "NLO #alpha_{T} Observed"); |
| 358 |
< |
if (sThird) leg->AddEntry(sThird, "LO #alpha_{T} Expected"); |
| 351 |
> |
//if (gCLsExpExclHT) leg->AddEntry(expLeg ,"NLO Expected (HT, CLs)","l"); |
| 352 |
> |
//if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,"NLO Observed (HT, CLs)","l"); |
| 353 |
> |
//if (gCLsObsExclHT_LO) leg->AddEntry(gCLsObsExclHT_LO,"LO Observed (HT, CLs)","l"); |
| 354 |
> |
//if (gCLsExpExclHT_LO) leg->AddEntry(gCLsExpExclHT_LO,"LO Expected (HT, CLs)","l"); |
| 355 |
> |
//if (gPLobsexcl) leg->AddEntry(gPLobsexcl,"Observed (PL, RooStat)","l"); |
| 356 |
> |
//if (gPLexpexcl) leg->AddEntry(gPLexpexcl,"Expected (PL, RooStat)","l"); |
| 357 |
> |
//if (gFCobsexcl) leg->AddEntry(gFCobsexcl,"NLO Observed (HT, FC, RooStat)","l"); |
| 358 |
> |
//if (gFCexpexcl) leg->AddEntry(gFCexpexcl,"NLO Expected (HT, FC, RooStat)","l"); |
| 359 |
> |
//if (gMCMCobsexcl) leg->AddEntry(gMCMCobsexcl,"Observed (MHT, MCMC, RooStat)","l"); |
| 360 |
> |
//if (gMCMCexpexcl) leg->AddEntry(gMCMCexpexcl,"LO Expected (HT, MCMC, RooStat)","l"); |
| 361 |
> |
//if (gexpexcl) leg->AddEntry(gexpexcl,"NLO Expected (MHT, CLs)","l"); |
| 362 |
> |
//if (gobsexcl) leg->AddEntry(gobsexcl,"NLO Observed (MHT, CLs)","l"); |
| 363 |
> |
if (sSecond) leg->AddEntry(sSecond,"#alpha_{T} Expected (FC, NLO)"); |
| 364 |
> |
if (sFirst) leg->AddEntry(sFirst, "#alpha_{T} Observed (FC, NLO)"); |
| 365 |
> |
if (sThird) leg->AddEntry(sThird, "#alpha_{T} Expected (FC, LO)"); |
| 366 |
|
leg->Draw(); |
| 367 |
|
gPad->RedrawAxis(); |
| 368 |
|
c1->SaveAs("results_tb50/Exclusion_m0_m12_tb50.pdf"); |
| 609 |
|
if (gCLsExpExclHT_qg_LO) gCLsExpExclHT_qg_LO->Draw("l"); |
| 610 |
|
// |
| 611 |
|
TLegend * leg_qg = new TLegend(0.45,0.78,0.85,0.89); |
| 612 |
< |
leg_qg->SetBorderSize(0);leg_qg->SetFillColor(0); |
| 612 |
> |
leg_qg->SetBorderSize(0);leg_qg->SetFillColor(0);leg->SetTextFont(42); |
| 613 |
|
TGraph * expLeg_qg = (TGraph*)gCLsExpExclHT_qg->Clone();expLeg_qg->SetFillStyle(gCLsExp1Sigma_qg->GetFillStyle());expLeg_qg->SetFillColor(gCLsExp1Sigma_qg->GetFillColor()); |
| 614 |
|
if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"NLO Observed (MHT, CLs)","l"); |
| 615 |
|
if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"NLO Expected (MHT, CLs)","l"); |
