| 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,ExpExclCL, 3,0, 2,2); |
| 232 |
< |
TGraph * gobsexcl = plotTools ->GetContour(hexcl,Mzero,Mhalf,ObsExclCL, 3,0, 2,1); |
| 233 |
< |
TGraph * gMCMCexpexcl = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,MCMCExpExclusion,3,0, 3,2); |
| 234 |
< |
TGraph * gMCMCobsexcl = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,MCMCObsExclusion,3,0, 3,1); |
| 235 |
< |
TGraph * gCLsExpExclHT = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL, 3,0, 1,2); |
| 236 |
< |
TGraph * gCLsObsExclHT = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,ObsExclCL, 3,0, 1,1); |
| 237 |
< |
TGraph * gCLsExpExclHTm1 = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,ExpExclCLm1sigma,3,0, 5,2); |
| 238 |
< |
TGraph * gCLsExpExclHTp1 = plotToolsHT ->GetContour(hexcl,Mzero,Mhalf,ExpExclCLp1sigma,3,0, 5,2); |
| 239 |
< |
TGraph * gCLsExp1Sigma = MakeBand(gCLsExpExclHTp1, gCLsExpExclHTm1); |
| 240 |
< |
hexcl->Draw("colz"); |
| 241 |
< |
if (gCLsExp1Sigma) gCLsExp1Sigma->Draw("lf"); |
| 242 |
< |
if (gexpexcl) gexpexcl->Draw("l"); |
| 243 |
< |
if (gCLsExpExclHT) gCLsExpExclHT->Draw("l"); |
| 244 |
< |
if (gMCMCexpexcl) gMCMCexpexcl->Draw("l"); |
| 245 |
< |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l"); |
| 246 |
< |
//if (gobsexcl) gobsexcl->Draw("l"); |
| 247 |
< |
//if (gPLobsexcl) gPLobsexcl->Draw("l"); |
| 248 |
< |
//if (gFCobsexcl) gFCobsexcl->Draw("l"); |
| 249 |
< |
//if (gMCMCobsexcl) gMCMCobsexcl->Draw("l"); |
| 250 |
< |
//if (gFCexpexcl) gFCexpexcl->Draw("l"); |
| 251 |
< |
//if (gPLexpexcl) gPLexpexcl->Draw("l"); |
| 252 |
< |
//if (gFCexpexcl) gFCexpexcl->Draw("l"); |
| 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( 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 |
> |
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 |
| 302 |
|
First->GetYaxis()->SetTitle("m_{1/2} (GeV)"); |
| 303 |
|
TSpline3 *sFirst = new TSpline3("sFirst",First); |
| 304 |
|
sFirst->SetLineColor(kBlue); |
| 305 |
< |
sFirst->SetLineWidth(3); |
| 305 |
> |
//sFirst->SetLineWidth(3); |
| 306 |
|
TSpline3 *sSecond = new TSpline3("sSecond",Second); |
| 307 |
|
sSecond->SetLineColor(kBlue); |
| 308 |
|
sSecond->SetLineStyle(2); |
| 309 |
< |
sSecond->SetLineWidth(3); |
| 309 |
> |
//sSecond->SetLineWidth(3); |
| 310 |
|
TSpline3 *sThird = new TSpline3("sThird",Third); |
| 311 |
< |
sThird->SetLineColor(kGreen+2); |
| 311 |
> |
sThird->SetLineColor(kBlue); |
| 312 |
|
sThird->SetLineStyle(4); |
| 313 |
< |
sThird->SetLineWidth(3); |
| 313 |
> |
//sThird->SetLineWidth(3); |
| 314 |
|
TEV_sn_d0_1->Draw("fsame"); |
| 315 |
|
TEV_sn_d0_2->Draw("fsame"); |
| 316 |
|
TEV_sg_d0->Draw("fsame"); |
| 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); |
| 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,"LO Observed (HT, CLs)","l"); |
| 352 |
< |
if (gCLsExpExclHT) leg->AddEntry(gCLsExpExclHT,"LO Expected (HT, CLs)","l"); |
| 353 |
< |
if (gexpexcl) leg->AddEntry(gexpexcl,"LO Expected (MHT, CLs)","l"); |
| 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,"Observed (FC, RooStat)","l"); |
| 358 |
< |
//if (gFCexpexcl) leg->AddEntry(gFCexpexcl,"Expected (FC, 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 (sFirst) leg->AddEntry(sFirst, "NLO RA1 Observed"); |
| 362 |
< |
if (sSecond) leg->AddEntry(sSecond,"NLO RA1 Expected"); |
| 363 |
< |
if (sThird) leg->AddEntry(sThird,"LO RA1 Observed"); |
| 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"); |
| 389 |
|
if (gobsexcl) gobsexcl->Draw("l"); |
| 390 |
|
c1->SaveAs("results_tb50/SignalContamination_MHT_m0_m12_tb50.pdf"); |
| 391 |
|
|
| 392 |
+ |
// Signal JEC Uncertainty MHT |
| 393 |
+ |
c1->SetLogz(0); |
| 394 |
+ |
TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty", |
| 395 |
+ |
40,200,600,26,140,400); |
| 396 |
+ |
plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC); |
| 397 |
+ |
hsig_jec_mht->SetMinimum(0.0); |
| 398 |
+ |
hsig_jec_mht->SetMaximum(0.25); |
| 399 |
+ |
hsig_jec_mht->SetContour(10); |
| 400 |
+ |
hsig_jec_mht->Draw("colz"); |
| 401 |
+ |
if (gobsexcl) gobsexcl->Draw("l"); |
| 402 |
+ |
c1->SaveAs("results_tb50/SigJEC_MHT_m0_m12_tb50.pdf"); |
| 403 |
+ |
|
| 404 |
+ |
// Signal JEC Uncertainty HT |
| 405 |
+ |
c1->SetLogz(0); |
| 406 |
+ |
TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty", |
| 407 |
+ |
40,200,600,26,140,400); |
| 408 |
+ |
plotTools->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC); |
| 409 |
+ |
hsig_jec_ht->SetMinimum(0.0); |
| 410 |
+ |
hsig_jec_ht->SetMaximum(0.25); |
| 411 |
+ |
hsig_jec_ht->SetContour(10); |
| 412 |
+ |
hsig_jec_ht->Draw("colz"); |
| 413 |
+ |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l"); |
| 414 |
+ |
c1->SaveAs("results_tb50/SigJEC_HT_m0_m12_tb50.pdf"); |
| 415 |
+ |
|
| 416 |
+ |
|
| 417 |
+ |
// Signal MuIso Uncertainty MHT |
| 418 |
+ |
c1->SetLogz(0); |
| 419 |
+ |
TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty", |
| 420 |
+ |
40,200,600,26,140,400); |
| 421 |
+ |
plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso); |
| 422 |
+ |
hsig_MuIso_mht->SetMinimum(0.0); |
| 423 |
+ |
hsig_MuIso_mht->SetMaximum(0.05); |
| 424 |
+ |
hsig_MuIso_mht->SetContour(5); |
| 425 |
+ |
hsig_MuIso_mht->Draw("colz"); |
| 426 |
+ |
if (gobsexcl) gobsexcl->Draw("l"); |
| 427 |
+ |
c1->SaveAs("results_tb50/SigMuIso_MHT_m0_m12_tb50.pdf"); |
| 428 |
+ |
|
| 429 |
+ |
// Signal MuIso Uncertainty HT |
| 430 |
+ |
c1->SetLogz(0); |
| 431 |
+ |
TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty", |
| 432 |
+ |
40,200,600,26,140,400); |
| 433 |
+ |
plotTools->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso); |
| 434 |
+ |
hsig_MuIso_ht->SetMinimum(0.0); |
| 435 |
+ |
hsig_MuIso_ht->SetMaximum(0.05); |
| 436 |
+ |
hsig_MuIso_ht->SetContour(5); |
| 437 |
+ |
hsig_MuIso_ht->Draw("colz"); |
| 438 |
+ |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l"); |
| 439 |
+ |
c1->SaveAs("results_tb50/SigMuIso_HT_m0_m12_tb50.pdf"); |
| 440 |
+ |
|
| 441 |
+ |
|
| 442 |
+ |
// Signal kFactor Uncertainty MHT |
| 443 |
+ |
c1->SetLogz(0); |
| 444 |
+ |
TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty", |
| 445 |
+ |
40,200,600,26,140,400); |
| 446 |
+ |
plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalUncertKfactor); |
| 447 |
+ |
hsig_kFactor_mht->SetMinimum(0.1); |
| 448 |
+ |
hsig_kFactor_mht->SetMaximum(0.2); |
| 449 |
+ |
hsig_kFactor_mht->SetContour(10); |
| 450 |
+ |
hsig_kFactor_mht->Draw("colz"); |
| 451 |
+ |
if (gobsexcl) gobsexcl->Draw("l"); |
| 452 |
+ |
c1->SaveAs("results_tb50/SigkFactor_MHT_m0_m12_tb50.pdf"); |
| 453 |
+ |
|
| 454 |
+ |
// Signal kFactor Uncertainty HT |
| 455 |
+ |
c1->SetLogz(0); |
| 456 |
+ |
TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty", |
| 457 |
+ |
40,200,600,26,140,400); |
| 458 |
+ |
plotTools->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalUncertKfactor); |
| 459 |
+ |
hsig_kFactor_ht->SetMinimum(0.1); |
| 460 |
+ |
hsig_kFactor_ht->SetMaximum(0.2); |
| 461 |
+ |
hsig_kFactor_ht->SetContour(10); |
| 462 |
+ |
hsig_kFactor_ht->Draw("colz"); |
| 463 |
+ |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l"); |
| 464 |
+ |
c1->SaveAs("results_tb50/SigkFactor_HT_m0_m12_tb50.pdf"); |
| 465 |
+ |
|
| 466 |
|
|
| 467 |
|
//----------------------------------------------------------------------------------- |
| 468 |
|
c1->SetLogz(1); |
| 587 |
|
// Exclusion in squark - gluino mass ---------------------------------------- |
| 588 |
|
TH2F*hexcl_qg = new TH2F("hexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion", |
| 589 |
|
40,200,900,30,200,900); |
| 590 |
< |
TGraph * gexpexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL, 3,0, 2,2); |
| 591 |
< |
TGraph * gobsexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ObsExclCL, 3,0, 2,1); |
| 590 |
> |
TGraph * gexpexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino, NLOExpExclCL, 3,0, 2,2); |
| 591 |
> |
TGraph * gexpexcl_qg_LO = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino, ExpExclCL, 3,0, 2,4); |
| 592 |
> |
TGraph * gobsexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino, NLOObsExclCL, 3,0, 2,1); |
| 593 |
|
TGraph * gMCMCexpexcl_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCExpExclusion,3,0, 3,2); |
| 594 |
|
TGraph * gMCMCobsexcl_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCObsExclusion,3,0, 3,1); |
| 595 |
< |
TGraph * gCLsExpExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL, 3,0, 1,2); |
| 596 |
< |
TGraph * gCLsObsExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ObsExclCL, 3,0, 1,1); |
| 597 |
< |
TGraph * gCLsExpExclHTm1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCLm1sigma,3,0, 5,2); |
| 598 |
< |
TGraph * gCLsExpExclHTp1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCLp1sigma,3,0, 5,2); |
| 595 |
> |
TGraph * gCLsExpExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL, 3,0, 1,2); |
| 596 |
> |
TGraph * gCLsExpExclHT_qg_LO = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL, 3,0, 1,4); |
| 597 |
> |
TGraph * gCLsObsExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL, 3,0, 1,1); |
| 598 |
> |
TGraph * gCLsExpExclHTm1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLm1sigma,3,0, 5,2); |
| 599 |
> |
TGraph * gCLsExpExclHTp1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLp1sigma,3,0, 5,2); |
| 600 |
|
TGraph * gCLsExp1Sigma_qg = MakeBand(gCLsExpExclHTp1_qg, gCLsExpExclHTm1_qg); |
| 601 |
|
hexcl_qg->Draw(""); |
| 602 |
|
// |
| 604 |
|
if (gobsexcl_qg) gobsexcl_qg->Draw("l"); |
| 605 |
|
if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->Draw("l"); |
| 606 |
|
if (gexpexcl_qg) gexpexcl_qg->Draw("l"); |
| 607 |
+ |
if (gexpexcl_qg_LO) gexpexcl_qg_LO->Draw("l"); |
| 608 |
|
if (gCLsExpExclHT_qg) gCLsExpExclHT_qg->Draw("l"); |
| 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); |
| 613 |
< |
if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"LO Observed (MHT, CLs)","l"); |
| 614 |
< |
if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"LO Expected (MHT, CLs)","l"); |
| 615 |
< |
if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"LO Observed (HT, CLs)","l"); |
| 616 |
< |
if (gCLsExpExclHT_qg) leg_qg->AddEntry(gCLsExpExclHT_qg,"LO Expected (HT, CLs)","l"); |
| 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"); |
| 616 |
> |
if (gexpexcl_qg_LO) leg_qg->AddEntry(gexpexcl_qg_LO,"LO Expected (MHT, CLs)","l"); |
| 617 |
> |
if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"NLO Observed (HT, CLs)","l"); |
| 618 |
> |
if (gCLsExpExclHT_qg) leg_qg->AddEntry(expLeg_qg,"NLO Expected (HT, CLs)","lf"); |
| 619 |
> |
if (gCLsExpExclHT_qg_LO) leg_qg->AddEntry(gCLsExpExclHT_qg_LO,"LO Expected (HT, CLs)","l"); |
| 620 |
|
//if (gPLobsexcl_qg) leg_qg->AddEntry(gPLobsexcl_qg,"Observed (PL, RooStat)","l"); |
| 621 |
|
//if (gPLexpexcl_qg) leg_qg->AddEntry(gPLexpexcl_qg,"Expected (PL, RooStat)","l"); |
| 622 |
|
leg_qg->Draw(); |
