| 1 |
#include "plot.h"
|
| 2 |
#include "SusyScan.h"
|
| 3 |
#include "GeneratorMasses.h"
|
| 4 |
#include "PlotTools.h"
|
| 5 |
#include "PrintTools.h"
|
| 6 |
#include "TheLimits.h"
|
| 7 |
#include "GlobalFunctions.h"
|
| 8 |
#include "StyleSettings.h"
|
| 9 |
|
| 10 |
#include "TRint.h"
|
| 11 |
#include "TROOT.h"
|
| 12 |
#include "TObjArray.h"
|
| 13 |
#include "TStyle.h"
|
| 14 |
|
| 15 |
#include "TChain.h"
|
| 16 |
#include "TFile.h"
|
| 17 |
#include "TGraph.h"
|
| 18 |
#include "TH1.h"
|
| 19 |
#include "TH2.h"
|
| 20 |
#include "TH2F.h"
|
| 21 |
#include "TTree.h"
|
| 22 |
#include "TKey.h"
|
| 23 |
#include "Riostream.h"
|
| 24 |
#include "TCanvas.h"
|
| 25 |
#include "TLegend.h"
|
| 26 |
#include "TLatex.h"
|
| 27 |
#include "TMarker.h"
|
| 28 |
#include "TPaveText.h"
|
| 29 |
|
| 30 |
#include <string>
|
| 31 |
#include <cmath>
|
| 32 |
#include <stdio.h>
|
| 33 |
|
| 34 |
int plot(int argc, char** argv)
|
| 35 |
{
|
| 36 |
//interactive root session
|
| 37 |
//TApplication theApp("App", 0, 0);
|
| 38 |
if (gROOT->IsBatch()) {
|
| 39 |
fprintf(stderr, "%s: cannot run in batch mode\n", argv[0]);
|
| 40 |
return 1;
|
| 41 |
}
|
| 42 |
|
| 43 |
//Style stuff
|
| 44 |
/*
|
| 45 |
gStyle->SetHistFillColor(0);
|
| 46 |
gStyle->SetPalette(1);
|
| 47 |
gStyle->SetCanvasColor(0);
|
| 48 |
gStyle->SetCanvasBorderMode(0);
|
| 49 |
gStyle->SetPadColor(0);
|
| 50 |
gStyle->SetPadBorderMode(0);
|
| 51 |
gStyle->SetFrameBorderMode(0);
|
| 52 |
|
| 53 |
gStyle->SetTitleFillColor(0);
|
| 54 |
gStyle->SetTitleBorderSize(0);
|
| 55 |
gStyle->SetTitleX(0.10);
|
| 56 |
gStyle->SetTitleY(0.98);
|
| 57 |
gStyle->SetTitleW(0.8);
|
| 58 |
gStyle->SetTitleH(0.06);
|
| 59 |
|
| 60 |
gStyle->SetErrorX(0);
|
| 61 |
gStyle->SetStatColor(0);
|
| 62 |
gStyle->SetStatBorderSize(0);
|
| 63 |
gStyle->SetStatX(0);
|
| 64 |
gStyle->SetStatY(0);
|
| 65 |
gStyle->SetStatW(0);
|
| 66 |
gStyle->SetStatH(0);
|
| 67 |
|
| 68 |
gStyle->SetTitleFont(22);
|
| 69 |
gStyle->SetLabelFont(22,"X");
|
| 70 |
gStyle->SetLabelFont(22,"Y");
|
| 71 |
gStyle->SetLabelFont(22,"Z");
|
| 72 |
gStyle->SetLabelSize(0.03,"X");
|
| 73 |
gStyle->SetLabelSize(0.03,"Y");
|
| 74 |
gStyle->SetLabelSize(0.03,"Z");
|
| 75 |
*/
|
| 76 |
util::StyleSettings::paperNoTitle();
|
| 77 |
gStyle->SetPadBottomMargin(0.18);
|
| 78 |
|
| 79 |
//gROOT->SetStyle("MyStyle");
|
| 80 |
TCanvas * c1 = new TCanvas("c1","c1",600,600);
|
| 81 |
//c1->SetFillStyle ( 4000 );
|
| 82 |
//c1->SetLeftMargin ( 0.15 );
|
| 83 |
//c1->SetRightMargin ( 0.15 );
|
| 84 |
//c1->SetBottomMargin( 0.10 );
|
| 85 |
c1->cd();
|
| 86 |
|
| 87 |
//Get limits and generator masses ---------------------------------------------------
|
| 88 |
TheLimits * genpoints = new TheLimits();
|
| 89 |
//genpoints->Fill(argc, argv);
|
| 90 |
//genpoints->Fill("limits_MHT_tb10_20110303/filelist.txt");
|
| 91 |
genpoints->Fill("limits_MHT_20110609/filelist_tb10.txt");
|
| 92 |
|
| 93 |
TheLimits * genpointsHT = new TheLimits();
|
| 94 |
//genpointsHT->Fill("limits_HT_tb10_20110303/filelist.txt");
|
| 95 |
genpointsHT->Fill("limits_HT_20110609/filelist_tb10.txt");
|
| 96 |
|
| 97 |
//TheLimits * genpointsHT = new TheLimits();
|
| 98 |
//genpointsHT->Fill("limits_HT_SigCont/filelist.txt");
|
| 99 |
|
| 100 |
//TheLimits * genpointsMHT = new TheLimits();
|
| 101 |
//genpointsMHT->Fill("limits_MHT_SigCont/filelist.txt");
|
| 102 |
|
| 103 |
|
| 104 |
//Replace read limits with specific numbers
|
| 105 |
//genpoints->OverwriteLimits("ABCD_MHT");
|
| 106 |
|
| 107 |
genpoints->FillGeneratorMasses("GenScan_tb10.dat");
|
| 108 |
genpoints->match();
|
| 109 |
genpointsHT->FillGeneratorMasses("GenScan_tb10.dat");
|
| 110 |
genpointsHT->match();
|
| 111 |
//genpointsHT->FillGeneratorMasses("GenScan_tb10.dat");
|
| 112 |
//genpointsHT->match();
|
| 113 |
//genpointsMHT->FillGeneratorMasses("GenScan_tb10.dat");
|
| 114 |
//genpointsMHT->match();
|
| 115 |
|
| 116 |
//genpoints->ExpandGrid(1);
|
| 117 |
//genpoints->ExpandGrid(1);
|
| 118 |
//genpoints->ExpandGrid(1);
|
| 119 |
//genpointsHT->ExpandGrid(1);
|
| 120 |
//genpointsHT->ExpandGrid(1);
|
| 121 |
//genpointsHT->ExpandGrid(1);
|
| 122 |
|
| 123 |
//the plotting ----------------------------------------------------------------------
|
| 124 |
//plotting helper functions
|
| 125 |
PlotTools<SusyScan> * plotTools = new PlotTools<SusyScan>(genpoints->GetScan());
|
| 126 |
PlotTools<SusyScan> * plotToolsHT = new PlotTools<SusyScan>(genpointsHT->GetScan());
|
| 127 |
//PlotTools<SusyScan> * plotToolsHT = new PlotTools<SusyScan>(genpointsHT->GetScan());
|
| 128 |
//PlotTools<SusyScan> * plotToolsMHT = new PlotTools<SusyScan>(genpointsMHT->GetScan());
|
| 129 |
PlotTools<GeneratorMasses> * plotMasses = new PlotTools<GeneratorMasses>(genpoints->GetGeneratorMasses());
|
| 130 |
|
| 131 |
//iso mass lines
|
| 132 |
TGraph * gl500 = plotMasses->Line(Mzero, Mhalf, MGluino, 500, 0.8 );
|
| 133 |
TGraph * gl600 = plotMasses->Line(Mzero, Mhalf, MGluino, 600, 0.8 );
|
| 134 |
TGraph * gl700 = plotMasses->Line(Mzero, Mhalf, MGluino, 700, 0.8 );
|
| 135 |
TGraph * gl800 = plotMasses->Line(Mzero, Mhalf, MGluino, 800, 0.8 );
|
| 136 |
TGraph * gl900 = plotMasses->Line(Mzero, Mhalf, MGluino, 900, 0.8 );
|
| 137 |
TGraph * gl1000 = plotMasses->Line(Mzero, Mhalf, MGluino, 1000, 0.8);
|
| 138 |
TGraph * sq500 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 500, 1);
|
| 139 |
TGraph * sq600 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 600, 1);
|
| 140 |
TGraph * sq700 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 700, 1);
|
| 141 |
TGraph * sq800 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 800, 1);
|
| 142 |
TGraph * sq900 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 900, 1);
|
| 143 |
TGraph * sq1000 = plotMasses->Line(Mzero, Mhalf, MSquarkL, 1000, 1);
|
| 144 |
|
| 145 |
TGraph * chi100 = plotMasses->Line(Mzero, Mhalf, MChi1, 50, 20);
|
| 146 |
TGraph * cha200 = plotMasses->Line(Mzero, Mhalf, MCha1, 200, 20);
|
| 147 |
cha200->SetLineColor(2);
|
| 148 |
|
| 149 |
//the histograms
|
| 150 |
c1->SetLogz(1);
|
| 151 |
//h->SetMaximum(27);
|
| 152 |
//h->SetMinimum(0.01);
|
| 153 |
|
| 154 |
/**/
|
| 155 |
// cross-section in M0 - M1/2
|
| 156 |
c1->SetRightMargin ( 0.2 );
|
| 157 |
TH2F*hxsec = new TH2F("xsec",";m_{0} [GeV]; m_{1/2} [GeV]; cross section [pb]",
|
| 158 |
200,0,2009.9,35,100,450);
|
| 159 |
plotTools->Area(hxsec, Mzero, Mhalf, Xsection);
|
| 160 |
hxsec->SetMinimum(0.01);
|
| 161 |
//sq500->Draw();
|
| 162 |
//gl500->Draw();
|
| 163 |
hxsec->Draw("colz");
|
| 164 |
c1->SaveAs("results_tb10/Xsection_m0_m12_tb10.pdf");
|
| 165 |
std::string wait;
|
| 166 |
//std::cin>>wait;
|
| 167 |
|
| 168 |
// Observed Limit in M0 - M1/2
|
| 169 |
TH2F*hobslimit = new TH2F("obslimit",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Limit [pb]",
|
| 170 |
200,0,2009.9,35,100,450);
|
| 171 |
plotTools->Area(hobslimit, Mzero, Mhalf, ObsXsecLimit);
|
| 172 |
hobslimit->SetMinimum(0.01);
|
| 173 |
hobslimit->Draw("colz");
|
| 174 |
c1->SaveAs("results_tb10/ObsLimit_m0_m12_tb10.pdf");
|
| 175 |
|
| 176 |
|
| 177 |
// Expected Limit in M0 - M1/2
|
| 178 |
TH2F*hexplimit = new TH2F("explimit",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Limit [pb]",
|
| 179 |
200,0,2009.9,35,100,450);
|
| 180 |
plotTools->Area(hexplimit, Mzero, Mhalf, ExpXsecLimit);
|
| 181 |
hexplimit->SetMinimum(0.01);
|
| 182 |
hexplimit->Draw("colz");
|
| 183 |
c1->SaveAs("results_tb10/ExpLimit_m0_m12_tb10.pdf");
|
| 184 |
|
| 185 |
// Signal Acceptance in M0 - M1/2
|
| 186 |
c1->SetRightMargin ( 0.2 );
|
| 187 |
c1->SetTopMargin ( 0.11 );
|
| 188 |
c1->SetLogz(0);
|
| 189 |
TH2F*hsigacc = new TH2F("sigacc",";m_{0} [GeV]; m_{1/2} [GeV]; Signal Acceptance",
|
| 190 |
100,0,2009.9,35,100,450);
|
| 191 |
hsigacc->SetNdivisions(505);
|
| 192 |
plotTools->Area(hsigacc, Mzero, Mhalf, NLOSignalAcceptance);
|
| 193 |
hsigacc->SetMinimum(0.0);
|
| 194 |
hsigacc->SetMaximum(0.35);
|
| 195 |
hsigacc->SetContour(14);
|
| 196 |
hsigacc->GetZaxis()->SetTitleOffset(1.5);
|
| 197 |
hsigacc->Draw("colz");
|
| 198 |
TLatex as; as.SetTextSize(0.025); as.SetTextFont(42);//ms.SetTextColor(12);
|
| 199 |
as.DrawLatex(520,461,"MHT selection, tan#beta=10, #mu>0, A_{0}=0");
|
| 200 |
as.SetTextSize(0.04);
|
| 201 |
as.DrawLatex( 50,461,"CMS");
|
| 202 |
c1->SaveAs("results_tb10/SigAccMHT_m0_m12_tb10.pdf");
|
| 203 |
|
| 204 |
|
| 205 |
// Signal Acceptance in M0 - M1/2
|
| 206 |
TH2F*hsigaccHT = new TH2F("sigaccHT",";m_{0} [GeV]; m_{1/2} [GeV]; Signal Acceptance",
|
| 207 |
100,0,2009.9,35,100,450);
|
| 208 |
hsigaccHT->SetNdivisions(505);
|
| 209 |
plotToolsHT->Area(hsigaccHT, Mzero, Mhalf, NLOSignalAcceptance);
|
| 210 |
hsigaccHT->SetMinimum(0.0);
|
| 211 |
hsigaccHT->SetMaximum(0.35);
|
| 212 |
hsigaccHT->SetContour(14);
|
| 213 |
hsigaccHT->GetZaxis()->SetTitleOffset(1.5);
|
| 214 |
hsigaccHT->Draw("colz");
|
| 215 |
as.SetTextSize(0.025); as.SetTextFont(42);//ms.SetTextColor(12);
|
| 216 |
as.DrawLatex(525,461,"HT selection, tan#beta=10, #mu>0, A_{0}=0");
|
| 217 |
as.SetTextSize(0.04);
|
| 218 |
as.DrawLatex( 50,461,"CMS");
|
| 219 |
c1->SaveAs("results_tb10/SigAccHT_m0_m12_tb10.pdf");
|
| 220 |
c1->SetTopMargin ( 0.1 );
|
| 221 |
|
| 222 |
|
| 223 |
// Exp. Limit on Number of Signal Events in M0 - M1/2
|
| 224 |
c1->SetLogz(0);
|
| 225 |
TH2F*hexplimitnsig = new TH2F("explimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL exp. limit signal events [# ]",
|
| 226 |
200,0,2009.9,35,100,450);
|
| 227 |
plotTools->Area(hexplimitnsig, Mzero, Mhalf, ExpNSignLimit);
|
| 228 |
hexplimitnsig->SetMinimum(0.0);
|
| 229 |
hexplimitnsig->SetMaximum(20);
|
| 230 |
hexplimitnsig->Draw("colz");
|
| 231 |
c1->SaveAs("results_tb10/ExpLimitOnNSig_m0_m12_tb10.pdf");
|
| 232 |
|
| 233 |
// Obs. Limit on Number of Signal Events in M0 - M1/2
|
| 234 |
TH2F*hobslimitnsig = new TH2F("obslimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL obs. limit signal events [# ]",
|
| 235 |
200,0,2009.9,35,100,450);
|
| 236 |
plotTools->Area(hobslimitnsig, Mzero, Mhalf, ObsNSignLimit);
|
| 237 |
hobslimitnsig->SetMinimum(0.0);
|
| 238 |
hobslimitnsig->SetMaximum(20);
|
| 239 |
hobslimitnsig->Draw("colz");
|
| 240 |
c1->SaveAs("results_tb10/ObsLimitOnNSig_m0_m12_tb10.pdf");
|
| 241 |
|
| 242 |
c1->SetLogz(0);
|
| 243 |
// Expected Exclusion in M0 - M1/2
|
| 244 |
TH2F*hexpexcl = new TH2F("expexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
|
| 245 |
200,0,2009.9,35,100,450);
|
| 246 |
plotTools->Area(hexpexcl, Mzero, Mhalf, NLOExpExclCL);
|
| 247 |
hexpexcl->Draw("colz");
|
| 248 |
c1->SaveAs("results_tb10/ExpExclusionMHT_m0_m12_tb10.pdf");
|
| 249 |
|
| 250 |
// Expected Exclusion in M0 - M1/2
|
| 251 |
TH2F*hexpexclHT = new TH2F("expexclHT",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
|
| 252 |
200,0,2009.9,35,100,450);
|
| 253 |
plotToolsHT->Area(hexpexclHT, Mzero, Mhalf, NLOExpExclCL);
|
| 254 |
hexpexclHT->Draw("colz");
|
| 255 |
c1->SaveAs("results_tb10/ExpExclusionHT_m0_m12_tb10.pdf");
|
| 256 |
|
| 257 |
// Observed Exclusion in M0 - M1/2
|
| 258 |
TH2F*hobsexcl = new TH2F("obsexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Exclusion",
|
| 259 |
200,0,2009.9,35,100,450);
|
| 260 |
plotTools->Area(hobsexcl, Mzero, Mhalf, NLOObsExclCL);
|
| 261 |
hobsexcl->Draw("colz");
|
| 262 |
c1->SaveAs("results_tb10/ObsExclusionMHT_m0_m12_tb10.pdf");
|
| 263 |
|
| 264 |
// Observed Exclusion in M0 - M1/2
|
| 265 |
TH2F*hobsexclHT = new TH2F("obsexclHT",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Exclusion",
|
| 266 |
200,0,2009.9,35,100,450);
|
| 267 |
plotToolsHT->Area(hobsexclHT, Mzero, Mhalf, NLOObsExclCL);
|
| 268 |
hobsexclHT->Draw("colz");
|
| 269 |
c1->SaveAs("results_tb10/ObsExclusionHT_m0_m12_tb10.pdf");
|
| 270 |
|
| 271 |
|
| 272 |
|
| 273 |
// TestContours in M0 - M1/2 /////////////////////////////////////////////////////////////////////
|
| 274 |
c1->SetRightMargin ( 0.1 );
|
| 275 |
TH2F*texcl = new TH2F("texcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
|
| 276 |
200,0,2009.9,35,100,450);
|
| 277 |
TH2F*testExclold=(TH2F*)texcl->Clone();
|
| 278 |
TH2F*testExclnew=(TH2F*)texcl->Clone();
|
| 279 |
plotTools->Area(testExclold, Mzero, Mhalf, NLOHybridObsExclCL);
|
| 280 |
plotTools->Area(testExclnew, Mzero, Mhalf, NLOHybridObsCL );
|
| 281 |
|
| 282 |
//std::vector<TGraph*> contours = plotTools->GetContours005(testExclnew,3);
|
| 283 |
std::vector<TGraph*> contours = plotToolsHT->GetContours(testExclold,3);
|
| 284 |
testExclold->SetNdivisions(505);
|
| 285 |
testExclold->Draw("colz");
|
| 286 |
//testExclnew->Draw("colz");
|
| 287 |
int col=kBlue-10;
|
| 288 |
for (std::vector<TGraph*>::iterator cont=contours.begin(); cont!=contours.end(); ++cont){
|
| 289 |
if (! *cont) continue;
|
| 290 |
double x, y;
|
| 291 |
(*cont)->GetPoint(0, x, y);
|
| 292 |
(*cont)->SetLineColor(col);
|
| 293 |
(*cont)->Draw("l");
|
| 294 |
TLatex l; l.SetTextSize(0.04); l.SetTextColor(col++);
|
| 295 |
char val[20];
|
| 296 |
sprintf(val,"%d",(int)(cont-contours.begin()));
|
| 297 |
l.DrawLatex(x,y,val);
|
| 298 |
//if (cont-contours.begin()>3) break;
|
| 299 |
}
|
| 300 |
c1->SaveAs("results_tb10/ExclusionTestContours_m0_m12_tb10.pdf");
|
| 301 |
|
| 302 |
|
| 303 |
// Exclusion in M0 - M1/2
|
| 304 |
TH2F*hexcl = new TH2F("hexcl",";m_{0} (GeV); m_{1/2} (GeV); 95% CL Exclusion",
|
| 305 |
200,0,2009.9,35,110,450);
|
| 306 |
hexcl->SetNdivisions(505);
|
| 307 |
TH2F*hs = new TH2F("hs","",200,0,2009.9,35,100,450);
|
| 308 |
TGraph * gexpexcl = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCL, 3,0, 2,2);
|
| 309 |
TGraph * gexpexcl_LO = plotTools ->GetContour(hs,Mzero,Mhalf,ExpExclCL, 3,0, 2,4);
|
| 310 |
TGraph * gobsexcl = plotTools ->GetContour(hs,Mzero,Mhalf,NLOObsExclCL, 3,0, 2,1);
|
| 311 |
TGraph * gobsexcl_LO = plotTools ->GetContour(hs,Mzero,Mhalf,ObsExclCL, 3,0, 2,1);
|
| 312 |
TGraph * gMCMCexpexcl = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCExpExclusion, 3,0, 3,2);
|
| 313 |
TGraph * gMCMCobsexcl = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCObsExclusion, 3,0, 3,1);
|
| 314 |
TGraph * gFCexpexclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion, 3,0, 3,2);
|
| 315 |
TGraph * gFCexpexclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,FCExpExclusion, 3,0, 3,4);
|
| 316 |
TGraph * gFCexpexclMHT = plotTools->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion, 3,0, 3,2);
|
| 317 |
TGraph * gFCexpexclMHT_LO = plotTools->GetContour(hs,Mzero,Mhalf,FCExpExclusion, 3,0, 3,4);
|
| 318 |
TGraph * gFCobsexclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion, 3,0, 3,1);
|
| 319 |
TGraph * gFCobsexclMHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion, 3,0, 3,1);
|
| 320 |
TGraph * gCLsExpExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCL, 3,0, 1,2);
|
| 321 |
//TGraph * gCLsExpNoSExclMHT= plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL, 3,0, 3,2);
|
| 322 |
//TGraph * gCLsExpNoSExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL, 3,0, 3,2);
|
| 323 |
TGraph * gCLsObsExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOObsExclCL, 3,0, 1,1);
|
| 324 |
TGraph * gCLsExpExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ExpExclCL, 3,0, 1,4);
|
| 325 |
TGraph * gCLsObsExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ObsExclCL, 3,0, 2,1);
|
| 326 |
TGraph * gCLsExpExclHTm1 = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma,3,0, 5,2);
|
| 327 |
TGraph * gCLsExpExclHTp1 = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma,3,0, 5,2);
|
| 328 |
TGraph * gCLsExpExclMHTm1 = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma, 3,0, 5,2);
|
| 329 |
TGraph * gCLsExpExclMHTp1 = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma, 3,0, 5,2);
|
| 330 |
TGraph * gHybridCLsObsExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOHybridObsExclCL, 3,1, 6,1);
|
| 331 |
TGraph * gHybridCLsObsExclMHT= plotTools ->GetContour(hs,Mzero,Mhalf,NLOHybridObsExclCL, 3,1, 6,1);
|
| 332 |
TGraph * gCheck1 = (TGraph*)gCLsObsExclHT->Clone();
|
| 333 |
TGraph * gCheck2 = (TGraph*)gCLsExpExclHT->Clone();
|
| 334 |
Smooth( gCLsObsExclHT, 27); gCLsObsExclHT->SetLineWidth( 3 );
|
| 335 |
//Smooth( gCLsExpNoSExclMHT, 27 ); gCLsExpNoSExclMHT->SetLineWidth( 3 );
|
| 336 |
//Smooth( gCLsExpNoSExclHT, 27 ); gCLsExpNoSExclHT->SetLineWidth( 3 );
|
| 337 |
Smooth( gCLsExpExclHT, 27 ); gCLsExpExclHT->SetLineWidth( 3 );
|
| 338 |
Smooth( gCLsObsExclHT_LO, 27); gCLsObsExclHT_LO->SetLineWidth( 2 );
|
| 339 |
Smooth( gCLsExpExclHT_LO, 27 ); gCLsExpExclHT_LO->SetLineWidth( 3 );
|
| 340 |
Smooth( gexpexcl, 27 );
|
| 341 |
Smooth( gobsexcl_LO, 27 );
|
| 342 |
Smooth( gobsexcl, 27 );
|
| 343 |
Smooth( gexpexcl_LO, 27 );
|
| 344 |
Smooth( gFCexpexclHT, 27 );
|
| 345 |
Smooth( gFCexpexclHT_LO, 27 );gFCexpexclHT_LO->SetLineWidth( 3 );
|
| 346 |
Smooth( gFCexpexclMHT, 27 );
|
| 347 |
Smooth( gFCexpexclMHT_LO, 27 );
|
| 348 |
Smooth( gFCobsexclHT, 27 );gFCobsexclHT->SetLineWidth( 3 );
|
| 349 |
Smooth( gFCobsexclMHT, 27 );
|
| 350 |
Smooth( gCLsExpExclHTm1, 27 );
|
| 351 |
Smooth( gCLsExpExclHTp1, 27 );
|
| 352 |
Smooth( gCLsExpExclMHTm1, 27 );
|
| 353 |
Smooth( gCLsExpExclMHTp1, 27 );
|
| 354 |
// Smooth( gHybridCLsObsExclHT, 27); gHybridCLsObsExclHT->SetLineWidth( 3 );
|
| 355 |
// Smooth( gHybridCLsObsExclMHT, 27);gHybridCLsObsExclMHT->SetLineWidth( 3 );
|
| 356 |
TGraph * CLsObsNLO = plotToolsHT->ChooseBest(gCLsObsExclHT,gobsexcl, gCLsObsExclHT,gobsexcl);
|
| 357 |
TGraph * CLsObsLO = plotToolsHT->ChooseBest(gCLsObsExclHT_LO,gobsexcl_LO, gCLsObsExclHT_LO,gobsexcl_LO);
|
| 358 |
TGraph * CLsExpNLO = plotToolsHT->ChooseBest(gCLsExpExclHT,gexpexcl, gCLsExpExclHT,gexpexcl);
|
| 359 |
//TGraph * CLsExpNoSNLO = plotToolsHT->ChooseBest(gCLsExpNoSExclHT,gCLsExpNoSExclMHT, gCLsExpNoSExclHT,gCLsExpNoSExclMHT);
|
| 360 |
TGraph * CLsExpLO = plotToolsHT->ChooseBest(gCLsExpExclHT_LO,gexpexcl_LO, gCLsExpExclHT_LO,gexpexcl_LO);
|
| 361 |
TGraph * FCExpLO = plotToolsHT->ChooseBest(gFCexpexclHT_LO,gFCexpexclMHT_LO,gFCexpexclHT_LO,gFCexpexclMHT_LO);
|
| 362 |
TGraph * FCExpNLO = plotToolsHT->ChooseBest(gFCexpexclHT,gFCexpexclMHT,gFCexpexclHT,gFCexpexclMHT);
|
| 363 |
TGraph * FCObsNLO = plotToolsHT->ChooseBest(gFCobsexclHT,gFCobsexclMHT,gFCobsexclHT,gFCobsexclMHT);
|
| 364 |
TGraph * HybridCLsObsNLO = plotToolsHT->ChooseBest(gHybridCLsObsExclHT,gHybridCLsObsExclMHT, gHybridCLsObsExclHT,gHybridCLsObsExclMHT);
|
| 365 |
TGraph * gCLsExpExclp1 = plotToolsHT->ChooseBest(gCLsExpExclHTp1,gCLsExpExclMHTp1, gCLsExpExclHTp1,gCLsExpExclMHTp1);
|
| 366 |
TGraph * gCLsExpExclm1 = plotToolsHT->ChooseBest(gCLsExpExclHTm1,gCLsExpExclMHTm1, gCLsExpExclHTm1,gCLsExpExclMHTm1);
|
| 367 |
TGraph * gCLsExp1Sigma = MakeBand(gCLsExpExclp1, gCLsExpExclm1);gCLsExp1Sigma->SetFillStyle(4010);
|
| 368 |
hexcl->GetYaxis()->SetTitleOffset(1.3);
|
| 369 |
hexcl->GetXaxis()->SetTitleOffset(0.92);
|
| 370 |
hexcl->Draw("colz");
|
| 371 |
//set old exclusion Limits
|
| 372 |
TGraph* LEP_ch = set_lep_ch(10);
|
| 373 |
TGraph* LEP_sl = set_lep_sl(10);//slepton curve
|
| 374 |
TGraph* TEV_sg_cdf = set_tev_sg_cdf(10);//squark gluino cdf
|
| 375 |
TGraph* TEV_sg_d0 = set_tev_sg_d0(10);//squark gluino d0
|
| 376 |
//TGraph* TEV_tlp_cdf = set_tev_tlp_cdf(10);//trilepton cdf
|
| 377 |
//TGraph* TEV_tlp_d0 = set_tev_tlp_d0(10);//trilepton d0
|
| 378 |
TGraph* Atlas = Atlas_m0_m12_tb3_obs();
|
| 379 |
Atlas->SetLineColor( 28 );
|
| 380 |
TGraph* stau = set_tev_stau(10);//stau
|
| 381 |
TGraph* TEV_sn_d0_1 = set_sneutrino_d0_1(10);
|
| 382 |
TGraph* TEV_sn_d0_2 = set_sneutrino_d0_2(10);
|
| 383 |
TGraphErrors* First = getObserved_NLO_tanBeta10();
|
| 384 |
TGraphErrors* Second = getExpected_NLO_tanBeta10();//getLO_jetMultis();
|
| 385 |
TGraphErrors* Third = getLO_tanBeta10();
|
| 386 |
First->GetXaxis()->SetRangeUser(0,505);
|
| 387 |
First->GetYaxis()->SetRangeUser(80,500);
|
| 388 |
First->GetXaxis()->SetTitle("m_{0} (GeV)");
|
| 389 |
First->GetYaxis()->SetTitle("m_{1/2} (GeV)");
|
| 390 |
TSpline3 *sFirst = new TSpline3("sFirst",First); sFirst->SetLineColor(kBlue);sFirst->SetLineWidth(1);
|
| 391 |
TSpline3 *sSecond = new TSpline3("sSecond",Second);
|
| 392 |
sSecond->SetLineColor(kBlue); sSecond->SetLineStyle(2); sSecond->SetLineWidth(1);
|
| 393 |
TSpline3 *sThird = new TSpline3("sThird",Third);
|
| 394 |
sThird->SetLineColor(kBlue); sThird->SetLineStyle(4); sThird->SetLineWidth(1);
|
| 395 |
//TEV_sn_d0_1->Draw("fsame");
|
| 396 |
//TEV_sn_d0_2->Draw("fsame"); //only for tb=3
|
| 397 |
LEP_ch->Draw("fsame");
|
| 398 |
LEP_sl->Draw("fsame");
|
| 399 |
TEV_sg_cdf->Draw("fsame");
|
| 400 |
TEV_sg_cdf->Draw("lsame");
|
| 401 |
TEV_sg_d0->Draw("fsame");
|
| 402 |
TEV_sg_d0->Draw("lsame");
|
| 403 |
TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
|
| 404 |
//b.DrawLatex( 6,150,"LEP2");
|
| 405 |
//b.DrawLatex( 22,138,"#tilde{l}^{#pm}");
|
| 406 |
//b.DrawLatex(360,110,"LEP2 #tilde{#chi}^{#pm}");
|
| 407 |
//b.DrawLatex(90,160,"CDF");
|
| 408 |
//b.DrawLatex(248,129,"D0 #tilde{g}, #tilde{q}");
|
| 409 |
//b.DrawLatex( 80,180,"D0 #tilde{#nu}");
|
| 410 |
// if (gCLsExp1Sigma) gCLsExp1Sigma->Draw("lf");
|
| 411 |
sFirst->Draw("same");
|
| 412 |
//sSecond->Draw("same");
|
| 413 |
//sThird->Draw("same");
|
| 414 |
Atlas->Draw("c,same");
|
| 415 |
TLatex ms; ms.SetTextSize(0.025); ms.SetTextFont(42);//ms.SetTextColor(12);
|
| 416 |
ms.DrawLatex(750,458,"tan#beta=10, #mu>0, A_{0}=0");
|
| 417 |
|
| 418 |
//gCheck1->Draw("same");
|
| 419 |
//gCheck2->Draw("same");
|
| 420 |
|
| 421 |
//LM points
|
| 422 |
TMarker* LM0 = new TMarker(200.,160.,20);
|
| 423 |
TMarker* LM1 = new TMarker(60.,250.,20);
|
| 424 |
LM0->SetMarkerSize(0.7);
|
| 425 |
LM1->SetMarkerSize(0.7);
|
| 426 |
TLatex* tLM0 = new TLatex(205.,160.,"LM0");
|
| 427 |
tLM0->SetTextSize(0.02);
|
| 428 |
TLatex* tLM1 = new TLatex(65.,243.,"LM1");
|
| 429 |
tLM1->SetTextSize(0.02);
|
| 430 |
LM0->Draw("same");
|
| 431 |
tLM0->Draw("same");
|
| 432 |
LM1->Draw("same");
|
| 433 |
tLM1->Draw("same");
|
| 434 |
TLegend* legexp = new TLegend(0.61,0.61,0.93,0.88,NULL,"brNDC");
|
| 435 |
legexp->SetFillColor(0);legexp->SetShadowColor(0);legexp->SetFillStyle(4000);legexp->SetTextFont(42);legexp->SetTextSize(0.025);legexp->SetBorderSize(0);
|
| 436 |
//TEV_sg_cdf.SetLineColor(1);
|
| 437 |
legexp->SetHeader("CMS");
|
| 438 |
legexp->AddEntry(TEV_sg_cdf,"CDF #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan#beta=5, #mu<0}","f");
|
| 439 |
legexp->AddEntry(TEV_sg_d0,"D0 #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan#beta=3, #mu<0}","f");
|
| 440 |
//ch_gr.SetLineColor(1);
|
| 441 |
legexp->AddEntry(LEP_ch,"LEP2 #tilde{#chi}_{1}^{#pm}","f");
|
| 442 |
//sl_gr.SetLineColor(1);
|
| 443 |
legexp->AddEntry(LEP_sl,"LEP2 #tilde{#font[12]{l}}^{#pm}","f"); //NOT FOR tb=50!
|
| 444 |
//if(tanbeta == 3)
|
| 445 |
//legexp->AddEntry(TEV_sn_d0_1,"D0 #chi^{#pm}_{1}, #chi^{0}_{2}","f");
|
| 446 |
legexp->AddEntry(sFirst, "CMS #alpha_{T}");
|
| 447 |
//legexp->AddEntry(Atlas, "Atlas, #scale[0.8]{tan#beta=3}","l");
|
| 448 |
legexp->Draw();
|
| 449 |
|
| 450 |
CLsObsNLO->Draw("l");
|
| 451 |
CLsObsLO->Draw("l");
|
| 452 |
CLsExpNLO->Draw("l");
|
| 453 |
//FCExpNLO->Draw("l");
|
| 454 |
//CLsExpNoSNLO->Draw("l");
|
| 455 |
// HybridCLsObsNLO->Draw("l");
|
| 456 |
//gobsexcl->Draw("l");
|
| 457 |
//gHybridCLsObsExclMHT->Draw("l");
|
| 458 |
|
| 459 |
|
| 460 |
stau->Draw("fsame");
|
| 461 |
b.DrawLatex( 10,380,"#tilde{#tau} LSP");
|
| 462 |
//constant ssqquark and gluino lines
|
| 463 |
TF1* lnsq[4];
|
| 464 |
TF1* lngl[4];
|
| 465 |
TLatex sqt; sqt.SetTextSize(0.02); sqt.SetTextAngle(-22);sqt.SetTextColor(kGray+2);
|
| 466 |
sqt.DrawLatex(148,220,"#font[92]{#tilde{q}(500)GeV}");
|
| 467 |
sqt.DrawLatex(390,350,"#font[92]{#tilde{q}(800)GeV}");
|
| 468 |
TLatex glt; glt.SetTextSize(0.02); sqt.SetTextAngle(-4); glt.SetTextColor(kGray+2);
|
| 469 |
glt.DrawLatex(840,172,"#font[92]{#tilde{g}(500)GeV}");
|
| 470 |
glt.DrawLatex(840,235,"#font[92]{#tilde{g}(650)GeV}");
|
| 471 |
glt.DrawLatex(838,297,"#font[92]{#tilde{g}(800)GeV}");
|
| 472 |
int tanBeta_=10;
|
| 473 |
for(int i = 0; i < 4; i++){
|
| 474 |
lnsq[i] = constant_squark(tanBeta_,i);
|
| 475 |
lngl[i] = constant_gluino(tanBeta_,i);
|
| 476 |
lngl[i]->Draw("same");
|
| 477 |
lnsq[i]->Draw("same");
|
| 478 |
}
|
| 479 |
TLegend * leg = new TLegend(0.25,0.76,0.6,0.88);
|
| 480 |
leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42);legexp->SetTextSize(0.025);
|
| 481 |
TGraph * expLeg = (TGraph*)CLsExpNLO->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor());
|
| 482 |
leg->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
|
| 483 |
leg->AddEntry(CLsObsNLO,"Observed, NLO","l");
|
| 484 |
leg->AddEntry(CLsObsLO, "Observed, LO","l");
|
| 485 |
leg->AddEntry(expLeg, "Expected #pm 1#sigma, NLO","lf");
|
| 486 |
//leg->AddEntry(CLsExpNoSNLO, "Expected, no-signal hyp., NLO","l");
|
| 487 |
leg->Draw();
|
| 488 |
gPad->RedrawAxis();
|
| 489 |
c1->SaveAs("results_tb10/Exclusion_m0_m12_tb10.pdf");
|
| 490 |
|
| 491 |
|
| 492 |
//plotToolsHT->Print(NLOSignal, Mzero,Mhalf, gCLsObsExclHT, 5);
|
| 493 |
// plotToolsHT->Print(NLOXsection, Mzero,Mhalf, MGluino, MSquarkL, gCLsObsExclHT, 5);
|
| 494 |
|
| 495 |
hexcl->Draw("");
|
| 496 |
gCLsExp1Sigma->SetName("CLsNLO_Exp1SigmaUncertaintyTb10");
|
| 497 |
CLsObsNLO->SetName("CLsNLO_ObservedTb10");
|
| 498 |
CLsExpNLO->SetName("CLsNLO_ExpectedTb10");
|
| 499 |
CLsExpLO->SetName("CLsLO_ExpectedTb10");
|
| 500 |
gCLsExp1Sigma->Draw("l");
|
| 501 |
CLsObsNLO->Draw("l");
|
| 502 |
CLsExpNLO->Draw("l");
|
| 503 |
CLsObsLO->Draw("l");
|
| 504 |
c1->SaveAs("results_tb10/LimitContours_tb10.C");
|
| 505 |
// ===============================================================================================================
|
| 506 |
|
| 507 |
|
| 508 |
//The RA1-style Expected limits - no-signal hypothesis *only* for pseudo data ///////////////////////////////
|
| 509 |
TGraph * gCLsExpNoSExclMHT = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL, 3,0, 1,2);
|
| 510 |
TGraph * gCLsExpNoSExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL, 3,0, 1,2);
|
| 511 |
TGraph * gCLsExpNoSExclHTm1 = plotToolsHT->ModifyExpSigma(gCLsExpExclHTm1, gCLsExpExclHT, gCLsExpNoSExclHT);
|
| 512 |
TGraph * gCLsExpNoSExclHTp1 = plotToolsHT->ModifyExpSigma(gCLsExpExclHTp1, gCLsExpExclHT, gCLsExpNoSExclHT);
|
| 513 |
TGraph * gCLsExpNoSExclMHTm1 = plotTools ->ModifyExpSigma(gCLsExpExclMHTm1,gCLsExpExclHT, gCLsExpNoSExclHT);
|
| 514 |
TGraph * gCLsExpNoSExclMHTp1 = plotTools ->ModifyExpSigma(gCLsExpExclMHTp1,gCLsExpExclHT, gCLsExpNoSExclHT);
|
| 515 |
Smooth( gCLsExpNoSExclMHT, 27 ); gCLsExpNoSExclMHT->SetLineWidth( 3 );
|
| 516 |
Smooth( gCLsExpNoSExclHT, 27 ); gCLsExpNoSExclHT->SetLineWidth( 3 );
|
| 517 |
Smooth( gCLsExpNoSExclHTm1, 27 );
|
| 518 |
Smooth( gCLsExpNoSExclHTp1, 27 );
|
| 519 |
Smooth( gCLsExpNoSExclMHTm1, 27 );
|
| 520 |
Smooth( gCLsExpNoSExclMHTp1, 27 );
|
| 521 |
TGraph * CLsExpNoSNLO = plotToolsHT->ChooseBest(gCLsExpNoSExclHT,gCLsExpNoSExclMHT, gCLsExpNoSExclHT,gCLsExpNoSExclMHT);
|
| 522 |
TGraph * gCLsExpNoSExclp1 = plotToolsHT->ChooseBest(gCLsExpNoSExclHTp1,gCLsExpNoSExclMHTp1, gCLsExpNoSExclHTp1,gCLsExpNoSExclMHTp1);
|
| 523 |
TGraph * gCLsExpNoSExclm1 = plotToolsHT->ChooseBest(gCLsExpNoSExclHTm1,gCLsExpNoSExclMHTm1, gCLsExpNoSExclHTm1,gCLsExpNoSExclMHTm1);
|
| 524 |
TGraph * gCLsExpNoS1Sigma = MakeBand(gCLsExpNoSExclp1, gCLsExpNoSExclm1);gCLsExpNoS1Sigma->SetFillStyle(4010);
|
| 525 |
hexcl->Draw("colz");
|
| 526 |
LEP_ch->Draw("fsame");
|
| 527 |
LEP_sl->Draw("fsame");
|
| 528 |
TEV_sg_cdf->Draw("fsame");
|
| 529 |
TEV_sg_cdf->Draw("lsame");
|
| 530 |
TEV_sg_d0->Draw("fsame");
|
| 531 |
TEV_sg_d0->Draw("lsame");
|
| 532 |
ms.DrawLatex(1500,458,"tan#beta=10, #mu>0, A_{0}=0");
|
| 533 |
LM0->SetMarkerSize(0.7);
|
| 534 |
LM1->SetMarkerSize(0.7);
|
| 535 |
tLM0->SetTextSize(0.02);
|
| 536 |
tLM1->SetTextSize(0.02);
|
| 537 |
LM0->Draw("same");
|
| 538 |
tLM0->Draw("same");
|
| 539 |
LM1->Draw("same");
|
| 540 |
tLM1->Draw("same");
|
| 541 |
gCLsExpNoS1Sigma->Draw("lf,same");
|
| 542 |
sFirst->Draw("same");
|
| 543 |
//Atlas->Draw("c,same");
|
| 544 |
CLsObsNLO->Draw("l,same");
|
| 545 |
CLsObsLO->Draw("l,same");
|
| 546 |
CLsExpNoSNLO->Draw("l,same");
|
| 547 |
//FCExpNLO->Draw("l");
|
| 548 |
stau->Draw("fsame");
|
| 549 |
b.DrawLatex( 10,380,"#tilde{#tau} LSP");
|
| 550 |
//constant ssqquark and gluino lines
|
| 551 |
sqt.DrawLatex(148,220,"#font[92]{#tilde{q}(500)GeV}");
|
| 552 |
sqt.DrawLatex(390,350,"#font[92]{#tilde{q}(800)GeV}");
|
| 553 |
glt.DrawLatex(840,172,"#font[92]{#tilde{g}(500)GeV}");
|
| 554 |
glt.DrawLatex(840,235,"#font[92]{#tilde{g}(650)GeV}");
|
| 555 |
glt.DrawLatex(838,297,"#font[92]{#tilde{g}(800)GeV}");
|
| 556 |
for(int i = 0; i < 4; i++){
|
| 557 |
lngl[i]->Draw("same");
|
| 558 |
lnsq[i]->Draw("same");
|
| 559 |
}
|
| 560 |
legexp->Draw();
|
| 561 |
leg->Draw();
|
| 562 |
gPad->RedrawAxis();
|
| 563 |
c1->SaveAs("results_tb10/Exclusion_m0_m12_tb10_NoSigHypPseudoData.pdf");
|
| 564 |
|
| 565 |
|
| 566 |
hexcl->Draw("");
|
| 567 |
gCLsExp1Sigma->SetName("CLsNLO_Exp1SigmaUncertaintyTb10");
|
| 568 |
CLsObsNLO->SetName("CLsNLO_ObservedTb10");
|
| 569 |
CLsExpNoSNLO->SetName("CLsNLO_ExpectedTb10");
|
| 570 |
CLsObsLO->SetName("CLsLO_ObservedTb10");
|
| 571 |
//gCLsExpNoS1Sigma->Draw("l");
|
| 572 |
CLsObsNLO->Draw("l,same");
|
| 573 |
//CLsObsLO->Draw("l,same");
|
| 574 |
CLsExpNoSNLO->Draw("l,same");
|
| 575 |
c1->SaveAs("results_tb10/LimitContours_tb10_NoSigHypPseudoData.C");
|
| 576 |
|
| 577 |
hexcl->Draw("colz");
|
| 578 |
LEP_ch->Draw("fsame");
|
| 579 |
LEP_sl->Draw("fsame");
|
| 580 |
TEV_sg_cdf->Draw("fsame");
|
| 581 |
TEV_sg_cdf->Draw("lsame");
|
| 582 |
TEV_sg_d0->Draw("fsame");
|
| 583 |
TEV_sg_d0->Draw("lsame");
|
| 584 |
gCLsExpNoS1Sigma->Draw("lf,same");
|
| 585 |
sFirst->Draw("same");
|
| 586 |
Atlas->Draw("c,same");
|
| 587 |
CLsObsNLO->Draw("l,same");
|
| 588 |
CLsExpNoSNLO->Draw("l,same");
|
| 589 |
ms.DrawLatex(750,458,"tan#beta=10, #mu>0, A_{0}=0");
|
| 590 |
LM0->SetMarkerSize(0.7);
|
| 591 |
LM1->SetMarkerSize(0.7);
|
| 592 |
tLM0->SetTextSize(0.02);
|
| 593 |
tLM1->SetTextSize(0.02);
|
| 594 |
LM0->Draw("same");
|
| 595 |
tLM0->Draw("same");
|
| 596 |
LM1->Draw("same");
|
| 597 |
tLM1->Draw("same");
|
| 598 |
stau->Draw("fsame");
|
| 599 |
b.DrawLatex( 10,380,"#tilde{#tau} LSP");
|
| 600 |
//Jim's limits///////////////////////////////////////////////////////////////////////////////
|
| 601 |
TGraph* JimObsHT = Jim_ht_tb10(0);
|
| 602 |
TGraph* JimObsMHT = Jim_mht_tb10(0);
|
| 603 |
TGraph* JimExpHT = Jim_ht_tb10(1);
|
| 604 |
TGraph* JimExpMHT = Jim_mht_tb10(1);
|
| 605 |
TGraph* JimExpHTup = Jim_ht_tb10(2);
|
| 606 |
TGraph* JimExpMHTup = Jim_mht_tb10(2);
|
| 607 |
TGraph* JimExpHTdn = Jim_ht_tb10(3);
|
| 608 |
TGraph* JimExpMHTdn = Jim_mht_tb10(3);
|
| 609 |
Smooth(JimObsHT, 2);
|
| 610 |
Smooth(JimObsMHT,2);
|
| 611 |
Smooth(JimExpHT, 2);
|
| 612 |
Smooth(JimExpMHT,2);
|
| 613 |
Smooth(JimExpHTup, 2);
|
| 614 |
Smooth(JimExpMHTup,2);
|
| 615 |
Smooth(JimExpHTdn, 2);
|
| 616 |
Smooth(JimExpMHTdn,2);
|
| 617 |
TGraph * JimObs = plotToolsHT->ChooseBest(JimObsHT,JimObsMHT,JimObsHT,JimObsMHT);
|
| 618 |
TGraph * JimExp = plotToolsHT->ChooseBest(JimExpHT,JimExpMHT,JimExpHT,JimExpMHT);
|
| 619 |
TGraph * JimExpup = plotToolsHT->ChooseBest(JimExpHTup,JimExpMHTup,JimExpHTup,JimExpMHTup);
|
| 620 |
TGraph * JimExpdn = plotToolsHT->ChooseBest(JimExpHTdn,JimExpMHTdn,JimExpHTdn,JimExpMHTdn);
|
| 621 |
TGraph * JimExp1Sigma = MakeBand(JimExpup, JimExpdn);JimExp1Sigma->SetFillStyle(3001);
|
| 622 |
TGraph * JimLeg = (TGraph*)JimExp->Clone();JimLeg->SetFillStyle(JimExp1Sigma->GetFillStyle());JimLeg->SetFillColor(JimExp1Sigma->GetFillColor());
|
| 623 |
JimExp1Sigma->Draw("f");
|
| 624 |
JimObs->Draw("c");
|
| 625 |
JimExp->Draw("c");
|
| 626 |
stau->Draw("fsame");
|
| 627 |
b.DrawLatex( 10,380,"#tilde{#tau} LSP");
|
| 628 |
//constant ssqquark and gluino lines
|
| 629 |
sqt.DrawLatex(148,220,"#font[92]{#tilde{q}(500)GeV}");
|
| 630 |
sqt.DrawLatex(390,350,"#font[92]{#tilde{q}(800)GeV}");
|
| 631 |
glt.DrawLatex(840,172,"#font[92]{#tilde{g}(500)GeV}");
|
| 632 |
glt.DrawLatex(840,235,"#font[92]{#tilde{g}(650)GeV}");
|
| 633 |
glt.DrawLatex(838,297,"#font[92]{#tilde{g}(800)GeV}");
|
| 634 |
for(int i = 0; i < 4; i++){
|
| 635 |
lngl[i]->Draw("same");
|
| 636 |
lnsq[i]->Draw("same");
|
| 637 |
}
|
| 638 |
legexp->Draw();
|
| 639 |
TLegend * legBayes = new TLegend(0.25,0.730,0.6,0.88);
|
| 640 |
legBayes->SetBorderSize(0);legBayes->SetFillColor(0);legBayes->SetFillStyle(4000);legBayes->SetTextFont(42);
|
| 641 |
legBayes->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
|
| 642 |
legBayes->AddEntry(JimObs, "Observed, Bayes","l");
|
| 643 |
legBayes->AddEntry(JimLeg, "Expected #pm 1#sigma, Bayes","lf");
|
| 644 |
legBayes->AddEntry(CLsObsNLO,"Observed, CLs","l");
|
| 645 |
legBayes->AddEntry(expLeg, "Expected #pm 1#sigma, CLs","lf");
|
| 646 |
legBayes->Draw();
|
| 647 |
gPad->RedrawAxis();
|
| 648 |
c1->SaveAs("results_tb10/Exclusion_m0_m12_tb10_Bayes.pdf");
|
| 649 |
|
| 650 |
|
| 651 |
|
| 652 |
// Signal Contamination in M0 - M1/2
|
| 653 |
c1->SetLogz(0);
|
| 654 |
c1->SetRightMargin ( 0.2 );
|
| 655 |
TH2F*hsigcont = new TH2F("sigcont",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
|
| 656 |
200,0,2009.9,35,100,450);
|
| 657 |
plotToolsHT->Area(hsigcont, Mzero, Mhalf, SignalContamination);
|
| 658 |
hsigcont->SetMinimum(0.01);
|
| 659 |
hsigcont->SetMaximum(20);
|
| 660 |
hsigcont->Draw("colz");
|
| 661 |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 662 |
c1->SaveAs("results_tb10/SignalContamination_HT_m0_m12_tb10.pdf");
|
| 663 |
|
| 664 |
// Signal Contamination in M0 - M1/2
|
| 665 |
c1->SetLogz(0);
|
| 666 |
TH2F*hsigcontMHT = new TH2F("sigcontMHT",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
|
| 667 |
200,0,2009.9,35,100,450);
|
| 668 |
plotTools->Area(hsigcontMHT, Mzero, Mhalf, SignalContamination);
|
| 669 |
hsigcontMHT->SetMinimum(0.01);
|
| 670 |
hsigcontMHT->SetMaximum(20);
|
| 671 |
hsigcontMHT->Draw("colz");
|
| 672 |
if (gobsexcl) gobsexcl->Draw("l");
|
| 673 |
c1->SaveAs("results_tb10/SignalContamination_MHT_m0_m12_tb10.pdf");
|
| 674 |
|
| 675 |
// Signal JEC Uncertainty MHT
|
| 676 |
c1->SetLogz(0);
|
| 677 |
TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
|
| 678 |
200,0,2009.9,35,100,450);
|
| 679 |
plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
|
| 680 |
hsig_jec_mht->SetMinimum(0.0);
|
| 681 |
hsig_jec_mht->SetMaximum(0.25);
|
| 682 |
hsig_jec_mht->SetContour(10);
|
| 683 |
hsig_jec_mht->Draw("colz");
|
| 684 |
if (gobsexcl) gobsexcl->Draw("l");
|
| 685 |
c1->SaveAs("results_tb10/SigJEC_MHT_m0_m12_tb10.pdf");
|
| 686 |
|
| 687 |
// Signal JEC Uncertainty HT
|
| 688 |
c1->SetLogz(0);
|
| 689 |
TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
|
| 690 |
200,0,2009.9,35,100,450);
|
| 691 |
plotToolsHT->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
|
| 692 |
hsig_jec_ht->SetMinimum(0.0);
|
| 693 |
hsig_jec_ht->SetMaximum(0.25);
|
| 694 |
hsig_jec_ht->SetContour(10);
|
| 695 |
hsig_jec_ht->Draw("colz");
|
| 696 |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 697 |
c1->SaveAs("results_tb10/SigJEC_HT_m0_m12_tb10.pdf");
|
| 698 |
|
| 699 |
|
| 700 |
// Signal MuIso Uncertainty MHT
|
| 701 |
c1->SetLogz(0);
|
| 702 |
TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
|
| 703 |
200,0,2009.9,35,100,450);
|
| 704 |
plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
|
| 705 |
hsig_MuIso_mht->SetMinimum(0.0);hsig_MuIso_mht->GetZaxis()->SetTitleOffset(1.3);
|
| 706 |
hsig_MuIso_mht->SetMaximum(0.05);
|
| 707 |
hsig_MuIso_mht->SetContour(5);
|
| 708 |
hsig_MuIso_mht->Draw("colz");
|
| 709 |
if (gobsexcl) gobsexcl->Draw("l");
|
| 710 |
c1->SaveAs("results_tb10/SigMuIso_MHT_m0_m12_tb10.pdf");
|
| 711 |
|
| 712 |
// Signal MuIso Uncertainty HT
|
| 713 |
c1->SetLogz(0);
|
| 714 |
TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
|
| 715 |
200,0,2009.9,35,100,450);
|
| 716 |
plotToolsHT->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
|
| 717 |
hsig_MuIso_ht->SetMinimum(0.0);hsig_MuIso_ht->GetZaxis()->SetTitleOffset(1.3);
|
| 718 |
hsig_MuIso_ht->SetMaximum(0.05);
|
| 719 |
hsig_MuIso_ht->SetContour(5);
|
| 720 |
hsig_MuIso_ht->Draw("colz");
|
| 721 |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 722 |
c1->SaveAs("results_tb10/SigMuIso_HT_m0_m12_tb10.pdf");
|
| 723 |
|
| 724 |
|
| 725 |
|
| 726 |
// Signal PDF Uncertainty MHT
|
| 727 |
c1->SetLogz(0);
|
| 728 |
TH2F*hsig_PDF_mht = new TH2F("sigPDFmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal PDF uncertainty",
|
| 729 |
200,0,2009.9,35,100,450);
|
| 730 |
plotTools->Area(hsig_PDF_mht, Mzero, Mhalf, SignalUncertPDF);
|
| 731 |
hsig_PDF_mht->SetMinimum(0.0);hsig_PDF_mht->GetZaxis()->SetTitleOffset(1.3);
|
| 732 |
hsig_PDF_mht->SetMaximum(0.20);
|
| 733 |
hsig_PDF_mht->SetContour(10);
|
| 734 |
hsig_PDF_mht->Draw("colz");
|
| 735 |
if (gobsexcl) gobsexcl->Draw("l");
|
| 736 |
c1->SaveAs("results_tb10/SigPDF_MHT_m0_m12_tb10.pdf");
|
| 737 |
|
| 738 |
// Signal PDF Uncertainty HT
|
| 739 |
c1->SetLogz(0);
|
| 740 |
TH2F*hsig_PDF_ht = new TH2F("sigPDFht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal PDF uncertainty",
|
| 741 |
200,0,2009.9,35,100,450);
|
| 742 |
plotToolsHT->Area(hsig_PDF_ht, Mzero, Mhalf, SignalUncertPDF);
|
| 743 |
hsig_PDF_ht->SetMinimum(0.0);hsig_PDF_ht->GetZaxis()->SetTitleOffset(1.3);
|
| 744 |
hsig_PDF_ht->SetMaximum(0.20);
|
| 745 |
hsig_PDF_ht->SetContour(10);
|
| 746 |
hsig_PDF_ht->Draw("colz");
|
| 747 |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 748 |
c1->SaveAs("results_tb10/SigPDF_HT_m0_m12_tb10.pdf");
|
| 749 |
|
| 750 |
|
| 751 |
// Signal kFactor Uncertainty MHT
|
| 752 |
c1->SetLogz(0);
|
| 753 |
TH2F*hsig_kFactorUncert_mht = new TH2F("sigkFactorUncertmht",";m_{0} [GeV]; m_{1/2} [GeV]; signal NLO scale uncertainty",
|
| 754 |
200,0,2009.9,35,100,450);
|
| 755 |
plotTools->Area(hsig_kFactorUncert_mht, Mzero, Mhalf, SignalUncertKfactor);
|
| 756 |
hsig_kFactorUncert_mht->SetMinimum(0.1);hsig_kFactorUncert_mht->GetZaxis()->SetTitleOffset(1.3);
|
| 757 |
hsig_kFactorUncert_mht->SetMaximum(0.2);
|
| 758 |
hsig_kFactorUncert_mht->SetContour(10);
|
| 759 |
hsig_kFactorUncert_mht->Draw("colz");
|
| 760 |
if (gobsexcl) gobsexcl->Draw("l");
|
| 761 |
c1->SaveAs("results_tb10/SigkFactorUncert_MHT_m0_m12_tb10.pdf");
|
| 762 |
|
| 763 |
// Signal kFactor Uncertainty HT
|
| 764 |
c1->SetLogz(0);
|
| 765 |
TH2F*hsig_kFactorUncert_ht = new TH2F("sigkFactorUncertht",";m_{0} [GeV]; m_{1/2} [GeV]; signal NLO scale uncertainty",
|
| 766 |
200,0,2009.9,35,100,450);
|
| 767 |
plotToolsHT->Area(hsig_kFactorUncert_ht, Mzero, Mhalf, SignalUncertKfactor);
|
| 768 |
hsig_kFactorUncert_ht->SetMinimum(0.1);hsig_kFactorUncert_ht->GetZaxis()->SetTitleOffset(1.3);
|
| 769 |
hsig_kFactorUncert_ht->SetMaximum(0.2);
|
| 770 |
hsig_kFactorUncert_ht->SetContour(10);
|
| 771 |
hsig_kFactorUncert_ht->Draw("colz");
|
| 772 |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 773 |
c1->SaveAs("results_tb10/SigkFactorUncert_HT_m0_m12_tb10.pdf");
|
| 774 |
|
| 775 |
// Signal kFactor HT
|
| 776 |
c1->SetLogz(0);
|
| 777 |
TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{0} [GeV]; m_{1/2} [GeV]; NLO signal k-factor",
|
| 778 |
200,0,2009.9,35,100,450);
|
| 779 |
plotToolsHT->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalKfactor);
|
| 780 |
hsig_kFactor_ht->SetMinimum(1.0);hsig_kFactor_ht->GetZaxis()->SetTitleOffset(1.3);
|
| 781 |
hsig_kFactor_ht->SetMaximum(2.4);
|
| 782 |
hsig_kFactor_ht->SetContour(7);
|
| 783 |
hsig_kFactor_ht->Draw("colz");
|
| 784 |
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 785 |
c1->SaveAs("results_tb10/SigkFactor_HT_m0_m12_tb10.pdf");
|
| 786 |
|
| 787 |
// Signal kFactor MHT
|
| 788 |
c1->SetLogz(0);
|
| 789 |
TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{0} [GeV]; m_{1/2} [GeV]; NLO signal k-factor",
|
| 790 |
200,0,2009.9,35,100,450);
|
| 791 |
plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalKfactor);
|
| 792 |
hsig_kFactor_mht->SetMinimum(1.0);hsig_kFactor_ht->GetZaxis()->SetTitleOffset(1.3);
|
| 793 |
hsig_kFactor_mht->SetMaximum(2.4);
|
| 794 |
hsig_kFactor_mht->SetContour(7);
|
| 795 |
hsig_kFactor_mht->Draw("colz");
|
| 796 |
if (gobsexcl) gobsexcl->Draw("l");
|
| 797 |
c1->SaveAs("results_tb10/SigkFactor_MHT_m0_m12_tb10.pdf");
|
| 798 |
|
| 799 |
|
| 800 |
|
| 801 |
|
| 802 |
|
| 803 |
|
| 804 |
//theApp.Run();
|
| 805 |
}
|
| 806 |
|
| 807 |
|
| 808 |
int main(int argc, char** argv)
|
| 809 |
{
|
| 810 |
return plot(argc, argv);
|
| 811 |
}
|
