ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/auterman/SusyScan/PlotScript/plot50tb.cc
(Generate patch)

Comparing UserCode/auterman/SusyScan/PlotScript/plot50tb.cc (file contents):
Revision 1.4 by auterman, Sun Feb 13 14:39:33 2011 UTC vs.
Revision 1.5 by auterman, Tue Feb 22 13:59:19 2011 UTC

# Line 230 | Line 230 | int plot(int argc, char** argv)
230     // Exclusion in M0 - M1/2
231     TH2F*hexcl = new TH2F("hexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
232                       40,200,600,26,140,400);
233 <   TGraph * gexpexcl         = plotTools     ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL,       3,0, 2,2);
234 <   TGraph * gobsexcl         = plotTools     ->GetContour(hexcl,Mzero,Mhalf,ObsExclCL,       3,0, 2,1);
235 <   TGraph * gMCMCexpexcl     = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,MCMCExpExclusion,3,0, 3,2);
236 <   TGraph * gMCMCobsexcl     = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,MCMCObsExclusion,3,0, 3,1);
237 <   TGraph * gCLsExpExclHT    = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL,       3,0, 1,2);
238 <   TGraph * gCLsObsExclHT    = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ObsExclCL,       3,0, 1,1);
239 <   TGraph * gCLsExpExclHTm1  = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ExpExclCLm1sigma,3,0, 5,2);
240 <   TGraph * gCLsExpExclHTp1  = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ExpExclCLp1sigma,3,0, 5,2);
233 >   TGraph * gexpexcl         = plotTools     ->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCL,       3,0, 2,2);
234 >   TGraph * gexpexcl_LO      = plotTools     ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL,       3,0, 2,4);
235 >   TGraph * gobsexcl         = plotTools     ->GetContour(hexcl,Mzero,Mhalf,NLOObsExclCL,       3,0, 2,1);
236 >   TGraph * gMCMCexpexcl_LO  = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,MCMCExpExclusion,3,0, 3,4);
237 >   TGraph * gMCMCobsexcl_LO  = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,MCMCObsExclusion,3,0, 3,1);
238 >   TGraph * gCLsExpExclHT_LO = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL,       3,0, 1,4);
239 >   TGraph * gCLsExpExclHT    = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCL,       3,0, 1,2);
240 >   TGraph * gCLsObsExclHT    = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,NLOObsExclCL,       3,0, 1,1);
241 >   TGraph * gCLsExpExclHTm1  = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCLm1sigma,3,0, 5,2);
242 >   TGraph * gCLsExpExclHTp1  = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCLp1sigma,3,0, 5,2);
243 >   Smooth( gCLsExpExclHT_LO, 33 ); gCLsExpExclHT_LO->SetLineWidth( 3 );
244 >   Smooth( gCLsObsExclHT, 33); gCLsObsExclHT->SetLineWidth( 3 );
245 >   Smooth( gCLsExpExclHT, 33 ); gCLsExpExclHT->SetLineWidth( 3 );
246 >   Smooth( gexpexcl, 33 );
247 >   Smooth( gobsexcl, 33 );
248 >   Smooth( gMCMCexpexcl_LO, 33 );
249 >   Smooth( gCLsExpExclHTm1, 33 );
250 >   Smooth( gCLsExpExclHTp1, 33 );
251     TGraph * gCLsExp1Sigma    = MakeBand(gCLsExpExclHTp1, gCLsExpExclHTm1);
252 <   hexcl->Draw("colz");  
253 <   if (gCLsExp1Sigma)    gCLsExp1Sigma->Draw("lf");
244 <   if (gexpexcl)         gexpexcl->Draw("l");
245 <   if (gCLsExpExclHT)    gCLsExpExclHT->Draw("l");
246 <   if (gMCMCexpexcl)     gMCMCexpexcl->Draw("l");
247 <   if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
248 <   //if (gobsexcl)       gobsexcl->Draw("l");
249 <   //if (gPLobsexcl)     gPLobsexcl->Draw("l");
250 <   //if (gFCobsexcl)     gFCobsexcl->Draw("l");
251 <   //if (gMCMCobsexcl)   gMCMCobsexcl->Draw("l");
252 <   //if (gFCexpexcl)     gFCexpexcl->Draw("l");
253 <   //if (gPLexpexcl)     gPLexpexcl->Draw("l");
254 <   //if (gFCexpexcl)     gFCexpexcl->Draw("l");
252 >   hexcl->Draw("colz");
253 >   gCLsExp1Sigma->Draw("f");  
254     //set old exclusion Limits
255     TGraph* LEP_ch = set_lep_ch(50);
256     TGraph* LEP_sl = set_lep_sl(50);//slepton curve
# Line 271 | Line 270 | int plot(int argc, char** argv)
270     First->GetYaxis()->SetTitle("m_{1/2} (GeV)");
271     TSpline3 *sFirst = new TSpline3("sFirst",First);
272     sFirst->SetLineColor(kBlue);
273 <   sFirst->SetLineWidth(3);
273 >   //sFirst->SetLineWidth(3);
274     TSpline3 *sSecond = new TSpline3("sSecond",Second);
275     sSecond->SetLineColor(kBlue);
276     sSecond->SetLineStyle(2);
277 <   sSecond->SetLineWidth(3);
277 >   //sSecond->SetLineWidth(3);
278     TSpline3 *sThird = new TSpline3("sThird",Third);
279 <   sThird->SetLineColor(kGreen+2);
279 >   sThird->SetLineColor(kBlue);
280     sThird->SetLineStyle(4);
281 <   sThird->SetLineWidth(3);
281 >   //sThird->SetLineWidth(3);
282     TEV_sn_d0_1->Draw("fsame");
283     TEV_sn_d0_2->Draw("fsame");
284     TEV_sg_d0->Draw("fsame");
# Line 290 | Line 289 | int plot(int argc, char** argv)
289     sFirst->Draw("same");
290     sSecond->Draw("same");
291     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");
297     TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
298     //b.DrawLatex( 10,420,"#tilde{#tau} LSP");
299     b.DrawLatex(230,145,"LEP2");
# Line 300 | Line 304 | int plot(int argc, char** argv)
304     //b.DrawLatex( 80,180,"D0 #tilde{#nu}");
305     TLegend * leg = new TLegend(0.45,0.7,0.85,0.89);
306     leg->SetBorderSize(0);leg->SetFillColor(0);
307 +   TGraph * expLeg = (TGraph*)gCLsExpExclHT->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor());
308     //if (gobsexcl) leg->AddEntry(gobsexcl,"Observed (MHT, CLs, TLimit)","l");
309 <   if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,"LO Observed (HT, CLs)","l");
310 <   if (gCLsExpExclHT) leg->AddEntry(gCLsExpExclHT,"LO Expected (HT, CLs)","l");
311 <   if (gexpexcl) leg->AddEntry(gexpexcl,"LO Expected (MHT, CLs)","l");
312 <   //if (gPLobsexcl) leg->AddEntry(gPLobsexcl,"Observed (PL, RooStat)","l");
313 <   //if (gPLexpexcl) leg->AddEntry(gPLexpexcl,"Expected (PL, RooStat)","l");
314 <   //if (gFCobsexcl) leg->AddEntry(gFCobsexcl,"Observed (FC, RooStat)","l");
315 <   //if (gFCexpexcl) leg->AddEntry(gFCexpexcl,"Expected (FC, RooStat)","l");
316 <   //if (gMCMCobsexcl) leg->AddEntry(gMCMCobsexcl,"Observed (MHT, MCMC, RooStat)","l");
312 <   if (gMCMCexpexcl) leg->AddEntry(gMCMCexpexcl,"LO Expected (HT, MCMC, RooStat)","l");
313 <   if (sFirst)  leg->AddEntry(sFirst, "NLO RA1 Observed");
314 <   if (sSecond) leg->AddEntry(sSecond,"NLO RA1 Expected");
315 <   if (sThird) leg->AddEntry(sThird,"LO RA1 Observed");
309 >   if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,      "NLO Observed (HT, CLs)","l");
310 >   if (gCLsExpExclHT) leg->AddEntry(expLeg       ,      "NLO Expected (HT, CLs)","l");
311 >   if (gCLsExpExclHT_LO) leg->AddEntry(gCLsExpExclHT_LO,"LO Expected (HT, CLs)","l");
312 >   if (gMCMCexpexcl_LO) leg->AddEntry(gMCMCexpexcl_LO,  "LO Expected (HT, MCMC, RooStat)","l");
313 >   if (gexpexcl) leg->AddEntry(gexpexcl,                "NLO Expected (MHT, CLs)","l");
314 >   if (sSecond) leg->AddEntry(sSecond,                  "NLO #alpha_{T} Expected");
315 >   if (sFirst)  leg->AddEntry(sFirst,                   "NLO #alpha_{T} Observed");
316 >   if (sThird) leg->AddEntry(sThird,                    "LO #alpha_{T} Expected");
317     leg->Draw();
318     gPad->RedrawAxis();
319     c1->SaveAs("results_tb50/Exclusion_m0_m12_tb50.pdf");
# Line 339 | Line 340 | int plot(int argc, char** argv)
340     if (gobsexcl)    gobsexcl->Draw("l");
341     c1->SaveAs("results_tb50/SignalContamination_MHT_m0_m12_tb50.pdf");
342    
343 +   // Signal JEC Uncertainty  MHT  
344 +   c1->SetLogz(0);
345 +   TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
346 +                     40,200,600,26,140,400);
347 +   plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
348 +   hsig_jec_mht->SetMinimum(0.0);
349 +   hsig_jec_mht->SetMaximum(0.25);
350 +   hsig_jec_mht->SetContour(10);
351 +   hsig_jec_mht->Draw("colz");
352 +   if (gobsexcl)    gobsexcl->Draw("l");
353 +   c1->SaveAs("results_tb50/SigJEC_MHT_m0_m12_tb50.pdf");
354 +  
355 +   // Signal JEC Uncertainty  HT  
356 +   c1->SetLogz(0);
357 +   TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
358 +                     40,200,600,26,140,400);
359 +   plotTools->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
360 +   hsig_jec_ht->SetMinimum(0.0);
361 +   hsig_jec_ht->SetMaximum(0.25);
362 +   hsig_jec_ht->SetContour(10);
363 +   hsig_jec_ht->Draw("colz");
364 +   if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
365 +   c1->SaveAs("results_tb50/SigJEC_HT_m0_m12_tb50.pdf");
366 +  
367 +
368 +   // Signal MuIso Uncertainty  MHT  
369 +   c1->SetLogz(0);
370 +   TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
371 +                     40,200,600,26,140,400);
372 +   plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
373 +   hsig_MuIso_mht->SetMinimum(0.0);
374 +   hsig_MuIso_mht->SetMaximum(0.05);
375 +   hsig_MuIso_mht->SetContour(5);
376 +   hsig_MuIso_mht->Draw("colz");
377 +   if (gobsexcl)    gobsexcl->Draw("l");
378 +   c1->SaveAs("results_tb50/SigMuIso_MHT_m0_m12_tb50.pdf");
379 +  
380 +   // Signal MuIso Uncertainty  HT  
381 +   c1->SetLogz(0);
382 +   TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
383 +                     40,200,600,26,140,400);
384 +   plotTools->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
385 +   hsig_MuIso_ht->SetMinimum(0.0);
386 +   hsig_MuIso_ht->SetMaximum(0.05);
387 +   hsig_MuIso_ht->SetContour(5);
388 +   hsig_MuIso_ht->Draw("colz");
389 +   if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
390 +   c1->SaveAs("results_tb50/SigMuIso_HT_m0_m12_tb50.pdf");
391 +  
392 +
393 +   // Signal kFactor Uncertainty  MHT  
394 +   c1->SetLogz(0);
395 +   TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
396 +                     40,200,600,26,140,400);
397 +   plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalUncertKfactor);
398 +   hsig_kFactor_mht->SetMinimum(0.1);
399 +   hsig_kFactor_mht->SetMaximum(0.2);
400 +   hsig_kFactor_mht->SetContour(10);
401 +   hsig_kFactor_mht->Draw("colz");
402 +   if (gobsexcl)    gobsexcl->Draw("l");
403 +   c1->SaveAs("results_tb50/SigkFactor_MHT_m0_m12_tb50.pdf");
404 +  
405 +   // Signal kFactor Uncertainty  HT  
406 +   c1->SetLogz(0);
407 +   TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
408 +                     40,200,600,26,140,400);
409 +   plotTools->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalUncertKfactor);
410 +   hsig_kFactor_ht->SetMinimum(0.1);
411 +   hsig_kFactor_ht->SetMaximum(0.2);
412 +   hsig_kFactor_ht->SetContour(10);
413 +   hsig_kFactor_ht->Draw("colz");
414 +   if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
415 +   c1->SaveAs("results_tb50/SigkFactor_HT_m0_m12_tb50.pdf");
416 +  
417  
418     //-----------------------------------------------------------------------------------
419     c1->SetLogz(1);
# Line 463 | Line 538 | int plot(int argc, char** argv)
538     // Exclusion in squark - gluino mass ----------------------------------------
539     TH2F*hexcl_qg = new TH2F("hexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
540                       40,200,900,30,200,900);
541 <   TGraph * gexpexcl_qg         = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL,       3,0, 2,2);
542 <   TGraph * gobsexcl_qg         = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ObsExclCL,       3,0, 2,1);
541 >   TGraph * gexpexcl_qg         = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,  NLOExpExclCL,       3,0, 2,2);
542 >   TGraph * gexpexcl_qg_LO      = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,  ExpExclCL,       3,0, 2,4);
543 >   TGraph * gobsexcl_qg         = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,  NLOObsExclCL,       3,0, 2,1);
544     TGraph * gMCMCexpexcl_qg     = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCExpExclusion,3,0, 3,2);
545     TGraph * gMCMCobsexcl_qg     = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCObsExclusion,3,0, 3,1);
546 <   TGraph * gCLsExpExclHT_qg    = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL,       3,0, 1,2);
547 <   TGraph * gCLsObsExclHT_qg    = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ObsExclCL,       3,0, 1,1);
548 <   TGraph * gCLsExpExclHTm1_qg  = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCLm1sigma,3,0, 5,2);
549 <   TGraph * gCLsExpExclHTp1_qg  = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCLp1sigma,3,0, 5,2);
546 >   TGraph * gCLsExpExclHT_qg    = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL,       3,0, 1,2);
547 >   TGraph * gCLsExpExclHT_qg_LO = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL,       3,0, 1,4);
548 >   TGraph * gCLsObsExclHT_qg    = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL,       3,0, 1,1);
549 >   TGraph * gCLsExpExclHTm1_qg  = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLm1sigma,3,0, 5,2);
550 >   TGraph * gCLsExpExclHTp1_qg  = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLp1sigma,3,0, 5,2);
551     TGraph * gCLsExp1Sigma_qg    = MakeBand(gCLsExpExclHTp1_qg, gCLsExpExclHTm1_qg);
552     hexcl_qg->Draw("");
553     //
# Line 478 | Line 555 | int plot(int argc, char** argv)
555     if (gobsexcl_qg)      gobsexcl_qg->Draw("l");
556     if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->Draw("l");
557     if (gexpexcl_qg)      gexpexcl_qg->Draw("l");
558 +   if (gexpexcl_qg_LO)   gexpexcl_qg_LO->Draw("l");
559     if (gCLsExpExclHT_qg) gCLsExpExclHT_qg->Draw("l");
560 +   if (gCLsExpExclHT_qg_LO) gCLsExpExclHT_qg_LO->Draw("l");
561     //
562     TLegend * leg_qg = new TLegend(0.45,0.78,0.85,0.89);
563     leg_qg->SetBorderSize(0);leg_qg->SetFillColor(0);
564 <   if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"LO Observed (MHT, CLs)","l");
565 <   if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"LO Expected (MHT, CLs)","l");
566 <   if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"LO Observed (HT, CLs)","l");
567 <   if (gCLsExpExclHT_qg) leg_qg->AddEntry(gCLsExpExclHT_qg,"LO Expected (HT, CLs)","l");
564 >   TGraph * expLeg_qg = (TGraph*)gCLsExpExclHT_qg->Clone();expLeg_qg->SetFillStyle(gCLsExp1Sigma_qg->GetFillStyle());expLeg_qg->SetFillColor(gCLsExp1Sigma_qg->GetFillColor());
565 >   if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"NLO Observed (MHT, CLs)","l");
566 >   if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"NLO Expected (MHT, CLs)","l");
567 >   if (gexpexcl_qg_LO) leg_qg->AddEntry(gexpexcl_qg_LO,"LO Expected (MHT, CLs)","l");
568 >   if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"NLO Observed (HT, CLs)","l");
569 >   if (gCLsExpExclHT_qg) leg_qg->AddEntry(expLeg_qg,"NLO Expected (HT, CLs)","lf");
570 >   if (gCLsExpExclHT_qg_LO) leg_qg->AddEntry(gCLsExpExclHT_qg_LO,"LO Expected (HT, CLs)","l");
571     //if (gPLobsexcl_qg) leg_qg->AddEntry(gPLobsexcl_qg,"Observed (PL, RooStat)","l");
572     //if (gPLexpexcl_qg) leg_qg->AddEntry(gPLexpexcl_qg,"Expected (PL, RooStat)","l");
573     leg_qg->Draw();

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines