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

Comparing UserCode/auterman/SusyScan/PlotScript/plot3tb.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 201 | Line 201 | int plot(int argc, char** argv)
201     TH2F*texcl = new TH2F("texcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
202                       50,0,509.9,35,100,450);
203     TH2F*texpexcl=(TH2F*)texcl->Clone();
204 <   plotTools->Area(texpexcl, Mzero, Mhalf, MCMCExpExclusion);
204 >   plotTools->Area(texpexcl, Mzero, Mhalf, NLOExpExclCL);
205     TH2F*tobsexcl=(TH2F*)texcl->Clone();
206     plotTools->Area(tobsexcl, Mzero, Mhalf, MCMCObsExclusion);
207     std::vector<TGraph*> contours = plotTools->GetContours(texpexcl,3);
# 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                       50,0,509.9,35,100,450);
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);
233 >   TGraph * gexpexcl         = plotTools     ->GetContour(hexcl,Mzero,Mhalf,NLOExpExclCL,       3,2, 2,2);
234 >   TGraph * gobsexcl         = plotTools     ->GetContour(hexcl,Mzero,Mhalf,NLOObsExclCL,       3,2, 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);
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 * gCLsExpExclHT_LO = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ExpExclCL,       3,0, 1,4);
240 >   TGraph * gCLsObsExclHT_LO = plotToolsHT   ->GetContour(hexcl,Mzero,Mhalf,ObsExclCL,       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 >   //TGraph * tt=(TGraph * )gCLsObsExclHT->Clone();
244 >   Smooth( gCLsObsExclHT_LO, 33); gCLsExpExclHT_LO->SetLineWidth( 3 );
245 >   Smooth( gCLsExpExclHT_LO, 33 );
246 >   Smooth( gCLsObsExclHT, 33); gCLsObsExclHT->SetLineWidth( 3 );
247 >   Smooth( gCLsExpExclHT, 33 ); gCLsExpExclHT->SetLineWidth( 3 );
248 >   Smooth( gexpexcl, 33 );
249 >   Smooth( gobsexcl, 33 );
250 >   Smooth( gMCMCexpexcl, 33 );
251 >   Smooth( gCLsExpExclHTm1, 33 );
252 >   Smooth( gCLsExpExclHTp1, 33 );
253     TGraph * gCLsExp1Sigma    = MakeBand(gCLsExpExclHTp1, gCLsExpExclHTm1);
254     hexcl->Draw("colz");  
243 –   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");
255     //set old exclusion Limits
256     TGraph* LEP_ch = set_lep_ch(3);
257     TGraph* LEP_sl = set_lep_sl(3);//slepton curve
# Line 271 | Line 271 | int plot(int argc, char** argv)
271     First->GetYaxis()->SetTitle("m_{1/2} (GeV)");
272     TSpline3 *sFirst = new TSpline3("sFirst",First);
273     sFirst->SetLineColor(kBlue);
274 <   sFirst->SetLineWidth(3);
274 >   //sFirst->SetLineWidth(3);
275     TSpline3 *sSecond = new TSpline3("sSecond",Second);
276     sSecond->SetLineColor(kBlue);
277     sSecond->SetLineStyle(2);
278 <   sSecond->SetLineWidth(3);
278 >   //sSecond->SetLineWidth(3);
279     TSpline3 *sThird = new TSpline3("sThird",Third);
280 <   sThird->SetLineColor(kGreen+2);
280 >   sThird->SetLineColor(kBlue);
281     sThird->SetLineStyle(4);
282 <   sThird->SetLineWidth(3);
282 >   //sThird->SetLineWidth(3);
283 >   if (gCLsExp1Sigma)    gCLsExp1Sigma->Draw("lf");
284     TEV_sn_d0_1->Draw("fsame");
285     TEV_sn_d0_2->Draw("fsame");
286     TEV_sg_d0->Draw("fsame");
# Line 290 | Line 291 | int plot(int argc, char** argv)
291     sFirst->Draw("same");
292     sSecond->Draw("same");
293     sThird->Draw("same");
294 +   if (gexpexcl)         gexpexcl->Draw("l");
295 +   if (gobsexcl)         gobsexcl->Draw("l");
296 +   if (gCLsExpExclHT)    gCLsExpExclHT->Draw("l");
297 +   if (gMCMCexpexcl)     gMCMCexpexcl->Draw("l");
298 +   if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
299 +   if (gCLsExpExclHT_LO)    gCLsExpExclHT_LO->Draw("l");
300 +
301     TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
302     b.DrawLatex( 10,420,"#tilde{#tau} LSP");
303     b.DrawLatex(  6,150,"LEP2");
# Line 300 | Line 308 | int plot(int argc, char** argv)
308     b.DrawLatex( 80,180,"D0 #tilde{#nu}");
309     TLegend * leg = new TLegend(0.45,0.7,0.85,0.89);
310     leg->SetBorderSize(0);leg->SetFillColor(0);
311 <   //if (gobsexcl) leg->AddEntry(gobsexcl,"Observed (MHT, CLs, TLimit)","l");
312 <   if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,"LO Observed (HT, CLs)","l");
313 <   if (gCLsExpExclHT) leg->AddEntry(gCLsExpExclHT,"LO Expected (HT, CLs)","l");
314 <   if (gexpexcl) leg->AddEntry(gexpexcl,"LO Expected (MHT, CLs, TLimit)","l");
311 >   TGraph * expLeg = (TGraph*)gCLsExpExclHT->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor());
312 >   if (gCLsExpExclHT) leg->AddEntry(expLeg       ,"NLO Expected (HT, CLs)","lf");
313 >   if (gCLsObsExclHT) leg->AddEntry(gCLsObsExclHT,"NLO Observed (HT, CLs)","l");
314 >   if (gCLsExpExclHT_LO) leg->AddEntry(gCLsExpExclHT_LO,"LO Expected (HT, CLs)","l");
315     //if (gPLobsexcl) leg->AddEntry(gPLobsexcl,"Observed (PL, RooStat)","l");
316     //if (gPLexpexcl) leg->AddEntry(gPLexpexcl,"Expected (PL, RooStat)","l");
317     //if (gFCobsexcl) leg->AddEntry(gFCobsexcl,"Observed (FC, RooStat)","l");
318     //if (gFCexpexcl) leg->AddEntry(gFCexpexcl,"Expected (FC, RooStat)","l");
319     //if (gMCMCobsexcl) leg->AddEntry(gMCMCobsexcl,"Observed (MHT, MCMC, RooStat)","l");
320     if (gMCMCexpexcl) leg->AddEntry(gMCMCexpexcl,"LO Expected (HT, MCMC, RooStat)","l");
321 <   if (sFirst)  leg->AddEntry(sFirst, "NLO RA1 Observed");
322 <   if (sSecond) leg->AddEntry(sSecond,"NLO RA1 Expected");
323 <   if (sThird) leg->AddEntry(sThird,"LO RA1 Observed");
321 >   if (gexpexcl) leg->AddEntry(gexpexcl,"NLO Expected (MHT, CLs, TLimit)","l");
322 >   if (gobsexcl) leg->AddEntry(gobsexcl,"NLO Observed (MHT, CLs, TLimit)","l");
323 >   if (sSecond) leg->AddEntry(sSecond,"NLO #alpha_{T} Expected");
324 >   if (sFirst)  leg->AddEntry(sFirst, "NLO #alpha_{T} Observed");
325 >   if (sThird) leg->AddEntry(sThird,"LO #alpha_{T} Expected");
326     leg->Draw();
327     gPad->RedrawAxis();
328     c1->SaveAs("results_tb3/Exclusion_m0_m12_tb3.pdf");
# Line 341 | Line 351 | int plot(int argc, char** argv)
351     if (gobsexcl)    gobsexcl->Draw("l");
352     c1->SaveAs("results_tb3/SignalContamination_MHT_m0_m12_tb3.pdf");
353    
354 +  
355 +   // Signal JEC Uncertainty  MHT  
356 +   c1->SetLogz(0);
357 +   TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
358 +                     50,0,509.9,25,100,350);
359 +   plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
360 +   hsig_jec_mht->SetMinimum(0.0);
361 +   hsig_jec_mht->SetMaximum(0.25);
362 +   hsig_jec_mht->SetContour(10);
363 +   hsig_jec_mht->Draw("colz");
364 +   if (gobsexcl)    gobsexcl->Draw("l");
365 +   c1->SaveAs("results_tb3/SigJEC_MHT_m0_m12_tb3.pdf");
366 +  
367 +   // Signal JEC Uncertainty  HT  
368 +   c1->SetLogz(0);
369 +   TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
370 +                    50,0,509.9,25,100,350);
371 +   plotTools->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
372 +   hsig_jec_ht->SetMinimum(0.0);
373 +   hsig_jec_ht->SetMaximum(0.25);
374 +   hsig_jec_ht->SetContour(10);
375 +   hsig_jec_ht->Draw("colz");
376 +   if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
377 +   c1->SaveAs("results_tb3/SigJEC_HT_m0_m12_tb3.pdf");
378 +  
379  
380 +   // Signal MuIso Uncertainty  MHT  
381 +   c1->SetLogz(0);
382 +   TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
383 +                    50,0,509.9,25,100,350);
384 +   plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
385 +   hsig_MuIso_mht->SetMinimum(0.0);
386 +   hsig_MuIso_mht->SetMaximum(0.05);
387 +   hsig_MuIso_mht->SetContour(5);
388 +   hsig_MuIso_mht->Draw("colz");
389 +   if (gobsexcl)    gobsexcl->Draw("l");
390 +   c1->SaveAs("results_tb3/SigMuIso_MHT_m0_m12_tb3.pdf");
391 +  
392 +   // Signal MuIso Uncertainty  HT  
393 +   c1->SetLogz(0);
394 +   TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
395 +                    50,0,509.9,25,100,350);
396 +   plotTools->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
397 +   hsig_MuIso_ht->SetMinimum(0.0);
398 +   hsig_MuIso_ht->SetMaximum(0.05);
399 +   hsig_MuIso_ht->SetContour(5);
400 +   hsig_MuIso_ht->Draw("colz");
401 +   if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
402 +   c1->SaveAs("results_tb3/SigMuIso_HT_m0_m12_tb3.pdf");
403 +  
404 +
405 +   // Signal kFactor Uncertainty  MHT  
406 +   c1->SetLogz(0);
407 +   TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
408 +                    50,0,509.9,25,100,350);
409 +   plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalUncertKfactor);
410 +   hsig_kFactor_mht->SetMinimum(0.1);
411 +   hsig_kFactor_mht->SetMaximum(0.2);
412 +   hsig_kFactor_mht->SetContour(10);
413 +   hsig_kFactor_mht->Draw("colz");
414 +   if (gobsexcl)    gobsexcl->Draw("l");
415 +   c1->SaveAs("results_tb3/SigkFactor_MHT_m0_m12_tb3.pdf");
416 +  
417 +   // Signal kFactor Uncertainty  HT  
418 +   c1->SetLogz(0);
419 +   TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
420 +                    50,0,509.9,25,100,350);
421 +   plotTools->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalUncertKfactor);
422 +   hsig_kFactor_ht->SetMinimum(0.1);
423 +   hsig_kFactor_ht->SetMaximum(0.2);
424 +   hsig_kFactor_ht->SetContour(10);
425 +   hsig_kFactor_ht->Draw("colz");
426 +   if (gCLsObsExclHT)    gCLsObsExclHT->Draw("l");
427 +   c1->SaveAs("results_tb3/SigkFactor_HT_m0_m12_tb3.pdf");
428 +  
429     //-----------------------------------------------------------------------------------
430     c1->SetLogz(1);
431  
# Line 464 | Line 548 | int plot(int argc, char** argv)
548  
549     TH2F*hexcl_qg = new TH2F("hexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
550                       50,200,800,30,200,900);
551 <   TGraph * gexpexcl_qg         = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL,       3,0, 2,2);
552 <   TGraph * gobsexcl_qg         = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ObsExclCL,       3,0, 2,1);
551 >   TGraph * gexpexcl_qg         = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL,       3,0, 2,2);
552 >   TGraph * gexpexcl_qg_LO      = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL,       3,0, 2,4);
553 >   TGraph * gobsexcl_qg         = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL,       3,0, 2,1);
554     TGraph * gMCMCexpexcl_qg     = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCExpExclusion,3,0, 3,2);
555     TGraph * gMCMCobsexcl_qg     = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCObsExclusion,3,0, 3,1);
556 <   TGraph * gCLsExpExclHT_qg    = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL,       3,0, 1,2);
557 <   TGraph * gCLsObsExclHT_qg    = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ObsExclCL,       3,0, 1,1);
558 <   TGraph * gCLsExpExclHTm1_qg  = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCLm1sigma,3,0, 5,2);
559 <   TGraph * gCLsExpExclHTp1_qg  = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCLp1sigma,3,0, 5,2);
556 >   TGraph * gCLsExpExclHT_qg    = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL,       3,0, 1,2);
557 >   TGraph * gCLsExpExclHT_qg_LO = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL,       3,0, 1,4);
558 >   TGraph * gCLsObsExclHT_qg    = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL,       3,0, 1,1);
559 >   TGraph * gCLsExpExclHTm1_qg  = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLm1sigma,3,0, 5,2);
560 >   TGraph * gCLsExpExclHTp1_qg  = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLp1sigma,3,0, 5,2);
561 >   Smooth(gexpexcl_qg, 11);
562 >   Smooth(gobsexcl_qg, 11);
563 >   Smooth(gCLsExpExclHT_qg, 11);
564 >   Smooth(gCLsObsExclHT_qg, 11);
565 >   Smooth(gCLsExpExclHTm1_qg, 11);
566 >   Smooth(gCLsExpExclHTp1_qg, 11);
567     TGraph * gCLsExp1Sigma_qg    = MakeBand(gCLsExpExclHTp1_qg, gCLsExpExclHTm1_qg);
568     hexcl_qg->Draw("");
569     //
# Line 480 | Line 572 | int plot(int argc, char** argv)
572     if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->Draw("l");
573     if (gexpexcl_qg)      gexpexcl_qg->Draw("l");
574     if (gCLsExpExclHT_qg) gCLsExpExclHT_qg->Draw("l");
575 +   if (gexpexcl_qg_LO)      gexpexcl_qg_LO->Draw("l");
576 +   if (gCLsExpExclHT_qg_LO) gCLsExpExclHT_qg_LO->Draw("l");
577     //
578     TLegend * leg_qg = new TLegend(0.45,0.78,0.85,0.89);
579     leg_qg->SetBorderSize(0);leg_qg->SetFillColor(0);
580 <   if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"LO Observed (MHT, CLs)","l");
581 <   if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"LO Expected (MHT, CLs)","l");
582 <   if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"LO Observed (HT, CLs)","l");
583 <   if (gCLsExpExclHT_qg) leg_qg->AddEntry(gCLsExpExclHT_qg,"LO Expected (HT, CLs)","l");
580 >   TGraph * expLeg_qg = (TGraph*)gCLsExpExclHT_qg->Clone();expLeg_qg->SetFillStyle(gCLsExp1Sigma_qg->GetFillStyle());expLeg_qg->SetFillColor(gCLsExp1Sigma_qg->GetFillColor());
581 >   if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"NLO Observed (MHT, CLs)","l");
582 >   if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"NLO Expected (MHT, CLs)","l");
583 >   if (gexpexcl_qg_LO) leg_qg->AddEntry(gexpexcl_qg_LO,"LO Expected (MHT, CLs)","l");
584 >   if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"NLO Observed (HT, CLs)","l");
585 >   if (gCLsExpExclHT_qg) leg_qg->AddEntry(expLeg_qg,"NLO Expected (HT, CLs)","lf");
586 >   if (gCLsExpExclHT_qg_LO) leg_qg->AddEntry(gCLsExpExclHT_qg_LO,"LO Expected (HT, CLs)","l");
587     //if (gPLobsexcl_qg) leg_qg->AddEntry(gPLobsexcl_qg,"Observed (PL, RooStat)","l");
588     //if (gPLexpexcl_qg) leg_qg->AddEntry(gPLexpexcl_qg,"Expected (PL, RooStat)","l");
589     leg_qg->Draw();

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines