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.7 by auterman, Thu Mar 10 20:59:52 2011 UTC vs.
Revision 1.8 by auterman, Fri May 20 07:39:34 2011 UTC

# Line 86 | Line 86 | int plot(int argc, char** argv)
86     TheLimits * genpoints = new TheLimits();
87     //genpoints->Fill(argc, argv);
88     //genpoints->Fill("limits_MHT_tb50_20110303/filelist.txt");
89 <   genpoints->Fill("limits_MHT_20110310/filelist_tb50.txt");
89 >   genpoints->Fill("limits_moriond_MHT/filelist_tb50.txt");
90  
91     TheLimits * genpointsHT = new TheLimits();
92     //genpointsHT->Fill("limits_HT_tb50_20110303/filelist.txt");
93 <   genpointsHT->Fill("limits_HT_20110310/filelist_tb50.txt");
93 >   genpointsHT->Fill("limits_moriond_HT/filelist_tb50.txt");
94  
95    
96     //Replace read limits with specific numbers
# Line 101 | Line 101 | int plot(int argc, char** argv)
101     genpointsHT->FillGeneratorMasses("GenScan_tb50.dat");
102     genpointsHT->match();
103  
104 +   //genpoints->ExpandGrid(1);
105 +   //genpoints->ExpandGrid(1);
106 +   //genpoints->ExpandGrid(1);
107 +   //genpointsHT->ExpandGrid(1);
108 +   //genpointsHT->ExpandGrid(1);
109 +   //genpointsHT->ExpandGrid(1);
110 +
111     //the plotting ----------------------------------------------------------------------
112     //plotting helper functions
113     PlotTools<SusyScan> * plotTools = new PlotTools<SusyScan>(genpoints->GetScan());
# Line 382 | Line 389 | int plot(int argc, char** argv)
389     TF1* lngl[4];
390     TLatex sqt; sqt.SetTextSize(0.02); sqt.SetTextAngle(-22);sqt.SetTextColor(kGray+2);
391     sqt.DrawLatex(270,180,"#font[92]{#tilde{q}(500)GeV}");
392 <   sqt.DrawLatex(327,255,"#font[92]{#tilde{q}(650)GeV}");
392 >   sqt.DrawLatex(327,245,"#font[92]{#tilde{q}(650)GeV}");
393     sqt.DrawLatex(390,348,"#font[92]{#tilde{q}(800)GeV}");
394     TLatex glt; glt.SetTextSize(0.02); sqt.SetTextAngle(-4); glt.SetTextColor(kGray+2);
395     glt.DrawLatex(528,195,"#font[92]{#tilde{g}(500)GeV}");
# Line 421 | Line 428 | int plot(int argc, char** argv)
428  
429  
430  
431 +   hexcl->Draw("colz");
432 +   LEP_ch->Draw("fsame");
433 +   LEP_sl->Draw("fsame");
434 +   if (gCLsExp1Sigma)    gCLsExp1Sigma->Draw("lf");
435 +   sFirst->Draw("same");
436 +   Atlas->Draw("c,same");
437 +   ms.DrawLatex(490,458,"tan#beta=50, #mu>0, A_{0}=0");
438 +   CLsObsNLO->Draw("l");
439 +   //CLsObsLO->Draw("l");
440 +   CLsExpNLO->Draw("l");
441 +   //FCExpNLO->Draw("l");
442 +   //Jim's limits///////////////////////////////////////////////////////////////////////////////
443 +   TGraph* JimObsHT  = Jim_ht_tb50(0);
444 +   TGraph* JimObsMHT = Jim_mht_tb50(0);
445 +   TGraph* JimExpHT  = Jim_ht_tb50(1);
446 +   TGraph* JimExpMHT = Jim_mht_tb50(1);
447 +   TGraph* JimExpHTup  = Jim_ht_tb50(2);
448 +   TGraph* JimExpMHTup = Jim_mht_tb50(2);
449 +   TGraph* JimExpHTdn  = Jim_ht_tb50(3);
450 +   TGraph* JimExpMHTdn = Jim_mht_tb50(3);
451 +   Smooth(JimObsHT, 2);
452 +   Smooth(JimObsMHT,2);
453 +   Smooth(JimExpHT, 2);
454 +   Smooth(JimExpMHT,2);
455 +   Smooth(JimExpHTup, 2);
456 +   Smooth(JimExpMHTup,2);
457 +   Smooth(JimExpHTdn, 2);
458 +   Smooth(JimExpMHTdn,2);
459 +   TGraph * JimObs = plotToolsHT->ChooseBest(JimObsHT,JimObsMHT,JimObsHT,JimObsMHT);
460 +   TGraph * JimExp = plotToolsHT->ChooseBest(JimExpHT,JimExpMHT,JimExpHT,JimExpMHT);
461 +   TGraph * JimExpup = plotToolsHT->ChooseBest(JimExpHTup,JimExpMHTup,JimExpHTup,JimExpMHTup);
462 +   TGraph * JimExpdn = plotToolsHT->ChooseBest(JimExpHTdn,JimExpMHTdn,JimExpHTdn,JimExpMHTdn);
463 +   TGraph * JimExp1Sigma = MakeBand(JimExpup, JimExpdn);JimExp1Sigma->SetFillStyle(3001);
464 +   TGraph * JimLeg = (TGraph*)JimExp->Clone();JimLeg->SetFillStyle(JimExp1Sigma->GetFillStyle());JimLeg->SetFillColor(JimExp1Sigma->GetFillColor());
465 +   JimExp1Sigma->Draw("f");
466 +   JimObs->Draw("c");
467 +   JimExp->Draw("c");
468 +   stau->Draw("fsame");
469 +   b.DrawLatex( 220,420,"#tilde{#tau} LSP");
470 +   //constant ssqquark and gluino lines
471 +   sqt.DrawLatex(270,180,"#font[92]{#tilde{q}(500)GeV}");
472 +   sqt.DrawLatex(327,245,"#font[92]{#tilde{q}(650)GeV}");
473 +   sqt.DrawLatex(390,348,"#font[92]{#tilde{q}(800)GeV}");
474 +   glt.DrawLatex(528,195,"#font[92]{#tilde{g}(500)GeV}");
475 +   glt.DrawLatex(528,258,"#font[92]{#tilde{g}(650)GeV}");
476 +   glt.DrawLatex(528,322,"#font[92]{#tilde{g}(800)GeV}");
477 +   for(int i = 0; i < 4; i++){
478 +    lngl[i]->Draw("same");  
479 +    lnsq[i]->Draw("same");
480 +   }
481 +   legexp->Draw();
482 +   TLegend * legBayes = new TLegend(0.38,0.8,0.73,0.92);
483 +   legBayes->SetBorderSize(0);legBayes->SetFillColor(0);legBayes->SetFillStyle(4000);legBayes->SetTextFont(42);
484 +   legBayes->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
485 +   legBayes->AddEntry(JimObs,   "Observed, Bayes","l");
486 +   legBayes->AddEntry(JimLeg,   "Expected #pm 1#sigma, Bayes","lf");
487 +   legBayes->AddEntry(CLsObsNLO,"Observed, CLs","l");
488 +   legBayes->AddEntry(expLeg,   "Expected #pm 1#sigma, CLs","lf");
489 +   legBayes->Draw();
490 +   gPad->RedrawAxis();
491 +   c1->SaveAs("results/Exclusion_m0_m12_tb50_Bayes.pdf");
492 +
493 +
494 +   //The RA1-style Expected limits - no-signal hypothesis *only* for pseudo data ///////////////////////////////
495 +   TGraph * gCLsExpNoSExclMHT   = plotTools  ->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL,  3,0, 1,2);
496 +   TGraph * gCLsExpNoSExclHT    = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL,  3,0, 1,2);
497 +   TGraph * gCLsExpNoSExclHTm1  = plotToolsHT->ModifyExpSigma(gCLsExpExclHTm1, gCLsExpExclHT, gCLsExpNoSExclHT);
498 +   TGraph * gCLsExpNoSExclHTp1  = plotToolsHT->ModifyExpSigma(gCLsExpExclHTp1, gCLsExpExclHT, gCLsExpNoSExclHT);
499 +   TGraph * gCLsExpNoSExclMHTm1 = plotTools  ->ModifyExpSigma(gCLsExpExclMHTm1,gCLsExpExclHT, gCLsExpNoSExclHT);
500 +   TGraph * gCLsExpNoSExclMHTp1 = plotTools  ->ModifyExpSigma(gCLsExpExclMHTp1,gCLsExpExclHT, gCLsExpNoSExclHT);
501 +   Smooth( gCLsExpNoSExclMHT, 27 ); gCLsExpNoSExclMHT->SetLineWidth( 3 );
502 +   Smooth( gCLsExpNoSExclHT, 27 );  gCLsExpNoSExclHT->SetLineWidth( 3 );
503 +   Smooth( gCLsExpNoSExclHTm1, 27 );
504 +   Smooth( gCLsExpNoSExclHTp1, 27 );
505 +   Smooth( gCLsExpNoSExclMHTm1, 27 );
506 +   Smooth( gCLsExpNoSExclMHTp1, 27 );
507 +   TGraph * CLsExpNoSNLO = plotToolsHT->ChooseBest(gCLsExpNoSExclHT,gCLsExpNoSExclMHT, gCLsExpNoSExclHT,gCLsExpNoSExclMHT);
508 +   TGraph * gCLsExpNoSExclp1 = plotToolsHT->ChooseBest(gCLsExpNoSExclHTp1,gCLsExpNoSExclMHTp1, gCLsExpNoSExclHTp1,gCLsExpNoSExclMHTp1);
509 +   TGraph * gCLsExpNoSExclm1 = plotToolsHT->ChooseBest(gCLsExpNoSExclHTm1,gCLsExpNoSExclMHTm1, gCLsExpNoSExclHTm1,gCLsExpNoSExclMHTm1);
510 +   TGraph * gCLsExpNoS1Sigma = MakeBand(gCLsExpNoSExclp1, gCLsExpNoSExclm1);gCLsExpNoS1Sigma->SetFillStyle(4010);
511 +   hexcl->Draw("colz");
512 +   LEP_ch->Draw("fsame");
513 +   LEP_sl->Draw("fsame");
514 +   if (gCLsExp1Sigma)    gCLsExp1Sigma->Draw("lf");
515 +   sFirst->Draw("same");
516 +   Atlas->Draw("c,same");
517 +   ms.DrawLatex(490,458,"tan#beta=50, #mu>0, A_{0}=0");
518 +   gCLsExpNoS1Sigma->Draw("lf,same");
519 +   sFirst->Draw("same");
520 +   Atlas->Draw("c,same");
521 +   CLsObsNLO->Draw("l,same");
522 +   CLsObsLO->Draw("l,same");
523 +   CLsExpNoSNLO->Draw("l,same");
524 +   //gCLsExpNoSExclMHT->Draw("l,same");
525 +   //gCLsExpNoSExclHT->Draw("l,same");
526 +   //gCLsObsExclHT->Draw("l,same");
527 +   //gobsexcl->Draw("l,same");
528 +   //FCExpNLO->Draw("l");
529 +   stau->Draw("fsame");
530 +   b.DrawLatex( 220,420,"#tilde{#tau} LSP");
531 +   //constant ssqquark and gluino lines
532 +   sqt.DrawLatex(270,180,"#font[92]{#tilde{q}(500)GeV}");
533 +   sqt.DrawLatex(327,245,"#font[92]{#tilde{q}(650)GeV}");
534 +   sqt.DrawLatex(390,348,"#font[92]{#tilde{q}(800)GeV}");
535 +   glt.DrawLatex(528,195,"#font[92]{#tilde{g}(500)GeV}");
536 +   glt.DrawLatex(528,258,"#font[92]{#tilde{g}(650)GeV}");
537 +   glt.DrawLatex(528,322,"#font[92]{#tilde{g}(800)GeV}");
538 +   for(int i = 0; i < 4; i++){
539 +    lngl[i]->Draw("same");  
540 +    lnsq[i]->Draw("same");
541 +   }
542 +   legexp->Draw();
543 +   leg->Draw();
544 +   gPad->RedrawAxis();
545 +   c1->SaveAs("results/Exclusion_m0_m12_tb50_NoSigHypPseudoData.pdf");
546 +
547 +
548 +
549     // Signal Contamination in M0 - M1/2
550     c1->SetLogz(0);
551     TH2F*hsigcont = new TH2F("sigcont",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
# Line 678 | Line 803 | int plot(int argc, char** argv)
803     //c1->SaveAs("plot_tb10.pdf");
804  
805  
806 <  
806 > /*  
807     //c1->SaveAs("plot_tb50.pdf");
808     c1->SetLogy(1);
809     c1->SetLogx(1);
# Line 737 | Line 862 | int plot(int argc, char** argv)
862     leg_S->AddEntry(gMCMCExpUncert,"Expected (MCMC, RooStat)","l");
863     leg_S->Draw();
864     c1->SaveAs("results_tb50/UncertaintyScan.pdf");
865 <
865 > */
866     //theApp.Run();
867   }
868  

Diff Legend

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