| 1 |
auterman |
1.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_moriond_MHT.old/filelist_tb10.txt");
|
| 92 |
|
|
|
| 93 |
|
|
TheLimits * genpointsHT = new TheLimits();
|
| 94 |
|
|
//genpointsHT->Fill("limits_HT_tb10_20110303/filelist.txt");
|
| 95 |
|
|
genpointsHT->Fill("limits_moriond_HT.old/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 |
|
|
100,0,1009.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/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 |
|
|
100,0,1009.9,35,100,450);
|
| 171 |
|
|
plotTools->Area(hobslimit, Mzero, Mhalf, ObsXsecLimit);
|
| 172 |
|
|
hobslimit->SetMinimum(0.01);
|
| 173 |
|
|
hobslimit->Draw("colz");
|
| 174 |
|
|
c1->SaveAs("results/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 |
|
|
100,0,1009.9,35,100,450);
|
| 180 |
|
|
plotTools->Area(hexplimit, Mzero, Mhalf, ExpXsecLimit);
|
| 181 |
|
|
hexplimit->SetMinimum(0.01);
|
| 182 |
|
|
hexplimit->Draw("colz");
|
| 183 |
|
|
c1->SaveAs("results/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 #times Efficiency",
|
| 190 |
|
|
50,0,1009.9,35,100,450);
|
| 191 |
|
|
plotTools->Area(hsigacc, Mzero, Mhalf, NLOSignalAcceptance);
|
| 192 |
|
|
hsigacc->SetMinimum(0.0);
|
| 193 |
|
|
hsigacc->SetMaximum(0.35);
|
| 194 |
|
|
hsigacc->SetContour(14);
|
| 195 |
|
|
hsigacc->GetZaxis()->SetTitleOffset(1.5);
|
| 196 |
|
|
hsigacc->Draw("colz");
|
| 197 |
|
|
TLatex as; as.SetTextSize(0.035); as.SetTextFont(42);//ms.SetTextColor(12);
|
| 198 |
|
|
as.DrawLatex(330,461,"#slash H_{T} selection, tan #beta=10, #mu>0, A_{0}=0");
|
| 199 |
|
|
as.SetTextSize(0.04);
|
| 200 |
|
|
as.DrawLatex( 50,461,"CMS");
|
| 201 |
|
|
c1->SaveAs("results/SigAccMHT_m0_m12_tb10.pdf");
|
| 202 |
|
|
|
| 203 |
|
|
|
| 204 |
|
|
// Signal Acceptance in M0 - M1/2
|
| 205 |
|
|
TH2F*hsigaccHT = new TH2F("sigaccHT",";m_{0} (GeV); m_{1/2} (GeV); Signal Acceptance #times Efficiency",
|
| 206 |
|
|
50,0,1009.9,35,100,450);
|
| 207 |
|
|
plotToolsHT->Area(hsigaccHT, Mzero, Mhalf, NLOSignalAcceptance);
|
| 208 |
|
|
hsigaccHT->SetMinimum(0.0);
|
| 209 |
|
|
hsigaccHT->SetMaximum(0.35);
|
| 210 |
|
|
hsigaccHT->SetContour(14);
|
| 211 |
|
|
hsigaccHT->GetZaxis()->SetTitleOffset(1.5);
|
| 212 |
|
|
hsigaccHT->Draw("colz");
|
| 213 |
|
|
as.SetTextSize(0.035); as.SetTextFont(42);//ms.SetTextColor(12);
|
| 214 |
|
|
as.DrawLatex(335,461,"H_{T} selection, tan #beta=10, #mu>0, A_{0}=0");
|
| 215 |
|
|
as.SetTextSize(0.04);
|
| 216 |
|
|
as.DrawLatex( 50,461,"CMS");
|
| 217 |
|
|
c1->SaveAs("results/SigAccHT_m0_m12_tb10.pdf");
|
| 218 |
|
|
c1->SetTopMargin ( 0.1 );
|
| 219 |
|
|
|
| 220 |
|
|
|
| 221 |
|
|
// Exp. Limit on Number of Signal Events in M0 - M1/2
|
| 222 |
|
|
c1->SetLogz(0);
|
| 223 |
|
|
TH2F*hexplimitnsig = new TH2F("explimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL exp. limit signal events [# ]",
|
| 224 |
|
|
100,0,1009.9,35,100,450);
|
| 225 |
|
|
plotTools->Area(hexplimitnsig, Mzero, Mhalf, ExpNSignLimit);
|
| 226 |
|
|
hexplimitnsig->SetMinimum(0.0);
|
| 227 |
|
|
hexplimitnsig->SetMaximum(20);
|
| 228 |
|
|
hexplimitnsig->Draw("colz");
|
| 229 |
|
|
c1->SaveAs("results/ExpLimitOnNSig_m0_m12_tb10.pdf");
|
| 230 |
|
|
|
| 231 |
|
|
// Obs. Limit on Number of Signal Events in M0 - M1/2
|
| 232 |
|
|
TH2F*hobslimitnsig = new TH2F("obslimitnsig",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL obs. limit signal events [# ]",
|
| 233 |
|
|
100,0,1009.9,35,100,450);
|
| 234 |
|
|
plotTools->Area(hobslimitnsig, Mzero, Mhalf, ObsNSignLimit);
|
| 235 |
|
|
hobslimitnsig->SetMinimum(0.0);
|
| 236 |
|
|
hobslimitnsig->SetMaximum(20);
|
| 237 |
|
|
hobslimitnsig->Draw("colz");
|
| 238 |
|
|
c1->SaveAs("results/ObsLimitOnNSig_m0_m12_tb10.pdf");
|
| 239 |
|
|
|
| 240 |
|
|
c1->SetLogz(0);
|
| 241 |
|
|
// Expected Exclusion in M0 - M1/2
|
| 242 |
|
|
TH2F*hexpexcl = new TH2F("expexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
|
| 243 |
|
|
100,0,1009.9,35,100,450);
|
| 244 |
|
|
plotTools->Area(hexpexcl, Mzero, Mhalf, ExpExclCL);
|
| 245 |
|
|
hexpexcl->Draw("colz");
|
| 246 |
|
|
c1->SaveAs("results/ExpExclusion_m0_m12_tb10.pdf");
|
| 247 |
|
|
|
| 248 |
|
|
// Observed Exclusion in M0 - M1/2
|
| 249 |
|
|
TH2F*hobsexcl = new TH2F("obsexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Exclusion",
|
| 250 |
|
|
100,0,1009.9,35,100,450);
|
| 251 |
|
|
plotTools->Area(hobsexcl, Mzero, Mhalf, ObsExclCL);
|
| 252 |
|
|
hobsexcl->Draw("colz");
|
| 253 |
|
|
c1->SaveAs("results/ObsExclusion_m0_m12_tb10.pdf");
|
| 254 |
|
|
|
| 255 |
|
|
// Observed Exclusion in M0 - M1/2
|
| 256 |
|
|
TH2F*hPLobsexcl = new TH2F("plobsexcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Observed Exclusion",
|
| 257 |
|
|
100,0,1009.9,35,100,450);
|
| 258 |
|
|
plotTools->Area(hPLobsexcl, Mzero, Mhalf, PLObsExclusion);
|
| 259 |
|
|
hPLobsexcl->Draw("colz");
|
| 260 |
|
|
c1->SaveAs("results/PL_ObsExclusion_m0_m12_tb10.pdf");
|
| 261 |
|
|
|
| 262 |
|
|
// TestContours in M0 - M1/2 /////////////////////////////////////////////////////////////////////
|
| 263 |
|
|
c1->SetRightMargin ( 0.1 );
|
| 264 |
|
|
TH2F*texcl = new TH2F("texcl",";m_{0} [GeV]; m_{1/2} [GeV]; 95% CL Expected Exclusion",
|
| 265 |
|
|
100,0,1009.9,35,100,450);
|
| 266 |
|
|
TH2F*testExclold=(TH2F*)texcl->Clone();
|
| 267 |
|
|
TH2F*testExclnew=(TH2F*)texcl->Clone();
|
| 268 |
|
|
plotTools->Area(testExclold, Mzero, Mhalf, NLOHybridObsExclCL);
|
| 269 |
|
|
plotTools->Area(testExclnew, Mzero, Mhalf, NLOHybridObsCL );
|
| 270 |
|
|
|
| 271 |
|
|
//std::vector<TGraph*> contours = plotTools->GetContours005(testExclnew,3);
|
| 272 |
|
|
std::vector<TGraph*> contours = plotToolsHT->GetContours(testExclold,3);
|
| 273 |
|
|
testExclold->Draw("colz");
|
| 274 |
|
|
//testExclnew->Draw("colz");
|
| 275 |
|
|
int col=kBlue-10;
|
| 276 |
|
|
for (std::vector<TGraph*>::iterator cont=contours.begin(); cont!=contours.end(); ++cont){
|
| 277 |
|
|
if (! *cont) continue;
|
| 278 |
|
|
double x, y;
|
| 279 |
|
|
(*cont)->GetPoint(0, x, y);
|
| 280 |
|
|
(*cont)->SetLineColor(col);
|
| 281 |
|
|
(*cont)->Draw("l");
|
| 282 |
|
|
TLatex l; l.SetTextSize(0.04); l.SetTextColor(col++);
|
| 283 |
|
|
char val[20];
|
| 284 |
|
|
sprintf(val,"%d",(int)(cont-contours.begin()));
|
| 285 |
|
|
l.DrawLatex(x,y,val);
|
| 286 |
|
|
//if (cont-contours.begin()>3) break;
|
| 287 |
|
|
}
|
| 288 |
|
|
c1->SaveAs("results/ExclusionTestContours_m0_m12_tb10.pdf");
|
| 289 |
|
|
|
| 290 |
|
|
// Exclusion in M0 - M1/2
|
| 291 |
|
|
TH2F*hexcl = new TH2F("hexcl",";m_{0} (GeV); m_{1/2} (GeV); 95% CL Exclusion",
|
| 292 |
|
|
100,0,1009.9,35,110,450);
|
| 293 |
|
|
TH2F*hs = new TH2F("hs","",100,0,1009.9,35,100,450);
|
| 294 |
|
|
TGraph * gexpexcl = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCL, 3,0, 2,2);
|
| 295 |
|
|
TGraph * gexpexcl_LO = plotTools ->GetContour(hs,Mzero,Mhalf,ExpExclCL, 3,0, 2,4);
|
| 296 |
|
|
TGraph * gobsexcl = plotTools ->GetContour(hs,Mzero,Mhalf,NLOObsExclCL, 3,0, 2,1);
|
| 297 |
|
|
TGraph * gobsexcl_LO = plotTools ->GetContour(hs,Mzero,Mhalf,ObsExclCL, 3,0, 2,1);
|
| 298 |
|
|
TGraph * gMCMCexpexcl = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCExpExclusion, 3,0, 3,2);
|
| 299 |
|
|
TGraph * gMCMCobsexcl = plotToolsHT->GetContour(hs,Mzero,Mhalf,MCMCObsExclusion, 3,0, 3,1);
|
| 300 |
|
|
TGraph * gFCexpexclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion, 3,0, 3,2);
|
| 301 |
|
|
TGraph * gFCexpexclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,FCExpExclusion, 3,0, 3,4);
|
| 302 |
|
|
TGraph * gFCexpexclMHT = plotTools->GetContour(hs,Mzero,Mhalf,NLOFCExpExclusion, 3,0, 3,2);
|
| 303 |
|
|
TGraph * gFCexpexclMHT_LO = plotTools->GetContour(hs,Mzero,Mhalf,FCExpExclusion, 3,0, 3,4);
|
| 304 |
|
|
TGraph * gFCobsexclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion, 3,0, 3,1);
|
| 305 |
|
|
TGraph * gFCobsexclMHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOFCObsExclusion, 3,0, 3,1);
|
| 306 |
|
|
TGraph * gCLsExpExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCL, 3,0, 1,2);
|
| 307 |
|
|
//TGraph * gCLsExpNoSExclMHT= plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL, 3,0, 3,2);
|
| 308 |
|
|
//TGraph * gCLsExpNoSExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpNoSigExclCL, 3,0, 3,2);
|
| 309 |
|
|
TGraph * gCLsObsExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOObsExclCL, 3,0, 1,1);
|
| 310 |
|
|
TGraph * gCLsExpExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ExpExclCL, 3,0, 1,4);
|
| 311 |
|
|
TGraph * gCLsObsExclHT_LO = plotToolsHT->GetContour(hs,Mzero,Mhalf,ObsExclCL, 3,0, 2,1);
|
| 312 |
|
|
TGraph * gCLsExpExclHTm1 = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma,3,0, 5,2);
|
| 313 |
|
|
TGraph * gCLsExpExclHTp1 = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma,3,0, 5,2);
|
| 314 |
|
|
TGraph * gCLsExpExclMHTm1 = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLm1sigma, 3,0, 5,2);
|
| 315 |
|
|
TGraph * gCLsExpExclMHTp1 = plotTools ->GetContour(hs,Mzero,Mhalf,NLOExpExclCLp1sigma, 3,0, 5,2);
|
| 316 |
|
|
TGraph * gHybridCLsObsExclHT = plotToolsHT->GetContour(hs,Mzero,Mhalf,NLOHybridObsExclCL, 3,1, 6,1);
|
| 317 |
|
|
TGraph * gHybridCLsObsExclMHT= plotTools ->GetContour(hs,Mzero,Mhalf,NLOHybridObsExclCL, 3,1, 6,1);
|
| 318 |
|
|
TGraph * gCheck1 = (TGraph*)gCLsObsExclHT->Clone();
|
| 319 |
|
|
TGraph * gCheck2 = (TGraph*)gCLsExpExclHT->Clone();
|
| 320 |
|
|
Smooth( gCLsObsExclHT, 27); gCLsObsExclHT->SetLineWidth( 3 );
|
| 321 |
|
|
//Smooth( gCLsExpNoSExclMHT, 27 ); gCLsExpNoSExclMHT->SetLineWidth( 3 );
|
| 322 |
|
|
//Smooth( gCLsExpNoSExclHT, 27 ); gCLsExpNoSExclHT->SetLineWidth( 3 );
|
| 323 |
|
|
Smooth( gCLsExpExclHT, 27 ); gCLsExpExclHT->SetLineWidth( 3 );
|
| 324 |
|
|
Smooth( gCLsObsExclHT_LO, 27); gCLsObsExclHT_LO->SetLineWidth( 2 );
|
| 325 |
|
|
Smooth( gCLsExpExclHT_LO, 27 ); gCLsExpExclHT_LO->SetLineWidth( 3 );
|
| 326 |
|
|
Smooth( gexpexcl, 27 );
|
| 327 |
|
|
Smooth( gobsexcl_LO, 27 );
|
| 328 |
|
|
Smooth( gobsexcl, 27 );
|
| 329 |
|
|
Smooth( gexpexcl_LO, 27 );
|
| 330 |
|
|
Smooth( gFCexpexclHT, 27 );
|
| 331 |
|
|
Smooth( gFCexpexclHT_LO, 27 );gFCexpexclHT_LO->SetLineWidth( 3 );
|
| 332 |
|
|
Smooth( gFCexpexclMHT, 27 );
|
| 333 |
|
|
Smooth( gFCexpexclMHT_LO, 27 );
|
| 334 |
|
|
Smooth( gFCobsexclHT, 27 );gFCobsexclHT->SetLineWidth( 3 );
|
| 335 |
|
|
Smooth( gFCobsexclMHT, 27 );
|
| 336 |
|
|
Smooth( gCLsExpExclHTm1, 27 );
|
| 337 |
|
|
Smooth( gCLsExpExclHTp1, 27 );
|
| 338 |
|
|
Smooth( gCLsExpExclMHTm1, 27 );
|
| 339 |
|
|
Smooth( gCLsExpExclMHTp1, 27 );
|
| 340 |
|
|
// Smooth( gHybridCLsObsExclHT, 27); gHybridCLsObsExclHT->SetLineWidth( 3 );
|
| 341 |
|
|
// Smooth( gHybridCLsObsExclMHT, 27);gHybridCLsObsExclMHT->SetLineWidth( 3 );
|
| 342 |
|
|
TGraph * CLsObsNLO = plotToolsHT->ChooseBest(gCLsObsExclHT,gobsexcl, gCLsObsExclHT,gobsexcl);
|
| 343 |
|
|
TGraph * CLsObsLO = plotToolsHT->ChooseBest(gCLsObsExclHT_LO,gobsexcl_LO, gCLsObsExclHT_LO,gobsexcl_LO);
|
| 344 |
|
|
TGraph * CLsExpNLO = plotToolsHT->ChooseBest(gCLsExpExclHT,gexpexcl, gCLsExpExclHT,gexpexcl);
|
| 345 |
|
|
//TGraph * CLsExpNoSNLO = plotToolsHT->ChooseBest(gCLsExpNoSExclHT,gCLsExpNoSExclMHT, gCLsExpNoSExclHT,gCLsExpNoSExclMHT);
|
| 346 |
|
|
TGraph * CLsExpLO = plotToolsHT->ChooseBest(gCLsExpExclHT_LO,gexpexcl_LO, gCLsExpExclHT_LO,gexpexcl_LO);
|
| 347 |
|
|
TGraph * FCExpLO = plotToolsHT->ChooseBest(gFCexpexclHT_LO,gFCexpexclMHT_LO,gFCexpexclHT_LO,gFCexpexclMHT_LO);
|
| 348 |
|
|
TGraph * FCExpNLO = plotToolsHT->ChooseBest(gFCexpexclHT,gFCexpexclMHT,gFCexpexclHT,gFCexpexclMHT);
|
| 349 |
|
|
TGraph * FCObsNLO = plotToolsHT->ChooseBest(gFCobsexclHT,gFCobsexclMHT,gFCobsexclHT,gFCobsexclMHT);
|
| 350 |
|
|
TGraph * HybridCLsObsNLO = plotToolsHT->ChooseBest(gHybridCLsObsExclHT,gHybridCLsObsExclMHT, gHybridCLsObsExclHT,gHybridCLsObsExclMHT);
|
| 351 |
|
|
TGraph * gCLsExpExclp1 = plotToolsHT->ChooseBest(gCLsExpExclHTp1,gCLsExpExclMHTp1, gCLsExpExclHTp1,gCLsExpExclMHTp1);
|
| 352 |
|
|
TGraph * gCLsExpExclm1 = plotToolsHT->ChooseBest(gCLsExpExclHTm1,gCLsExpExclMHTm1, gCLsExpExclHTm1,gCLsExpExclMHTm1);
|
| 353 |
|
|
TGraph * gCLsExp1Sigma = MakeBand(gCLsExpExclp1, gCLsExpExclm1);gCLsExp1Sigma->SetFillStyle(4010);
|
| 354 |
|
|
hexcl->GetYaxis()->SetTitleOffset(1.3);
|
| 355 |
|
|
hexcl->GetXaxis()->SetTitleOffset(0.92);
|
| 356 |
|
|
hexcl->Draw("colz");
|
| 357 |
|
|
//set old exclusion Limits
|
| 358 |
|
|
TGraph* LEP_ch = set_lep_ch(10);
|
| 359 |
|
|
TGraph* LEP_sl = set_lep_sl(10);//slepton curve
|
| 360 |
|
|
TGraph* TEV_sg_cdf = set_tev_sg_cdf(10);//squark gluino cdf
|
| 361 |
|
|
TGraph* TEV_sg_d0 = set_tev_sg_d0(10);//squark gluino d0
|
| 362 |
|
|
//TGraph* TEV_tlp_cdf = set_tev_tlp_cdf(10);//trilepton cdf
|
| 363 |
|
|
//TGraph* TEV_tlp_d0 = set_tev_tlp_d0(10);//trilepton d0
|
| 364 |
|
|
TGraph* Atlas = Atlas_m0_m12_tb3_obs();
|
| 365 |
|
|
Atlas->SetLineColor( 28 );
|
| 366 |
|
|
TGraph* stau = set_tev_stau(10);//stau
|
| 367 |
|
|
TGraph* TEV_sn_d0_1 = set_sneutrino_d0_1(10);
|
| 368 |
|
|
TGraph* TEV_sn_d0_2 = set_sneutrino_d0_2(10);
|
| 369 |
|
|
TGraphErrors* First = getObserved_NLO_tanBeta10();
|
| 370 |
|
|
TGraphErrors* Second = getExpected_NLO_tanBeta10();//getLO_jetMultis();
|
| 371 |
|
|
TGraphErrors* Third = getLO_tanBeta10();
|
| 372 |
|
|
First->GetXaxis()->SetRangeUser(0,505);
|
| 373 |
|
|
First->GetYaxis()->SetRangeUser(80,500);
|
| 374 |
|
|
First->GetXaxis()->SetTitle("m_{0} (GeV)");
|
| 375 |
|
|
First->GetYaxis()->SetTitle("m_{1/2} (GeV)");
|
| 376 |
|
|
TSpline3 *sFirst = new TSpline3("sFirst",First); sFirst->SetLineColor(kBlue);sFirst->SetLineWidth(1);
|
| 377 |
|
|
TSpline3 *sSecond = new TSpline3("sSecond",Second);
|
| 378 |
|
|
sSecond->SetLineColor(kBlue); sSecond->SetLineStyle(2); sSecond->SetLineWidth(1);
|
| 379 |
|
|
TSpline3 *sThird = new TSpline3("sThird",Third);
|
| 380 |
|
|
sThird->SetLineColor(kBlue); sThird->SetLineStyle(4); sThird->SetLineWidth(1);
|
| 381 |
|
|
//TEV_sn_d0_1->Draw("fsame");
|
| 382 |
|
|
//TEV_sn_d0_2->Draw("fsame"); //only for tb=3
|
| 383 |
|
|
LEP_ch->Draw("fsame");
|
| 384 |
|
|
LEP_sl->Draw("fsame");
|
| 385 |
|
|
TEV_sg_cdf->Draw("fsame");
|
| 386 |
|
|
TEV_sg_cdf->Draw("lsame");
|
| 387 |
|
|
TEV_sg_d0->Draw("fsame");
|
| 388 |
|
|
TEV_sg_d0->Draw("lsame");
|
| 389 |
|
|
TLatex b; b.SetTextSize(0.02); b.SetTextColor(1);
|
| 390 |
|
|
//b.DrawLatex( 6,150,"LEP2");
|
| 391 |
|
|
//b.DrawLatex( 22,138,"#tilde{l}^{#pm}");
|
| 392 |
|
|
//b.DrawLatex(360,110,"LEP2 #tilde{#chi}^{#pm}");
|
| 393 |
|
|
//b.DrawLatex(90,160,"CDF");
|
| 394 |
|
|
//b.DrawLatex(248,129,"D0 #tilde{g}, #tilde{q}");
|
| 395 |
|
|
//b.DrawLatex( 80,180,"D0 #tilde{#nu}");
|
| 396 |
|
|
// if (gCLsExp1Sigma) gCLsExp1Sigma->Draw("lf");
|
| 397 |
|
|
sFirst->Draw("same");
|
| 398 |
|
|
//sSecond->Draw("same");
|
| 399 |
|
|
//sThird->Draw("same");
|
| 400 |
|
|
Atlas->Draw("c,same");
|
| 401 |
|
|
TLatex ms; ms.SetTextSize(0.025); ms.SetTextFont(42);//ms.SetTextColor(12);
|
| 402 |
|
|
ms.DrawLatex(750,458,"tan #beta=10, #mu>0, A_{0}=0");
|
| 403 |
|
|
|
| 404 |
|
|
//gCheck1->Draw("same");
|
| 405 |
|
|
//gCheck2->Draw("same");
|
| 406 |
|
|
|
| 407 |
|
|
//LM points
|
| 408 |
|
|
TMarker* LM0 = new TMarker(200.,160.,20);
|
| 409 |
|
|
TMarker* LM1 = new TMarker(60.,250.,20);
|
| 410 |
|
|
LM0->SetMarkerSize(0.7);
|
| 411 |
|
|
LM1->SetMarkerSize(0.7);
|
| 412 |
|
|
TLatex* tLM0 = new TLatex(205.,160.,"LM0");
|
| 413 |
|
|
tLM0->SetTextSize(0.02);
|
| 414 |
|
|
TLatex* tLM1 = new TLatex(65.,243.,"LM1");
|
| 415 |
|
|
tLM1->SetTextSize(0.02);
|
| 416 |
|
|
LM0->Draw("same");
|
| 417 |
|
|
tLM0->Draw("same");
|
| 418 |
|
|
LM1->Draw("same");
|
| 419 |
|
|
tLM1->Draw("same");
|
| 420 |
|
|
TLegend* legexp = new TLegend(0.61,0.61,0.93,0.88,NULL,"brNDC");
|
| 421 |
|
|
legexp->SetFillColor(0);legexp->SetShadowColor(0);legexp->SetFillStyle(4000);legexp->SetTextFont(42);legexp->SetTextSize(0.025);legexp->SetBorderSize(0);
|
| 422 |
|
|
//TEV_sg_cdf.SetLineColor(1);
|
| 423 |
|
|
legexp->SetHeader("CMS tan #beta=10, #mu>0, A_{0}=0");
|
| 424 |
|
|
legexp->AddEntry(TEV_sg_cdf,"CDF #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan #beta=5, #mu<0}","f");
|
| 425 |
|
|
legexp->AddEntry(TEV_sg_d0,"D0 #tilde{#font[12]{g}}, #tilde{#font[12]{q}}, #scale[0.8]{tan #beta=3, #mu<0}","f");
|
| 426 |
|
|
//ch_gr.SetLineColor(1);
|
| 427 |
|
|
legexp->AddEntry(LEP_ch,"LEP2 #tilde{#chi}_{1}^{#pm}","f");
|
| 428 |
|
|
//sl_gr.SetLineColor(1);
|
| 429 |
|
|
legexp->AddEntry(LEP_sl,"LEP2 #tilde{#font[12]{l}}^{#pm}","f"); //NOT FOR tb=50!
|
| 430 |
|
|
//if(tanbeta == 3)
|
| 431 |
|
|
//legexp->AddEntry(TEV_sn_d0_1,"D0 #chi^{#pm}_{1}, #chi^{0}_{2}","f");
|
| 432 |
|
|
legexp->AddEntry(sFirst, "CMS #alpha_{T}");
|
| 433 |
|
|
//legexp->AddEntry(Atlas, "Atlas, #scale[0.8]{tan#beta=3}","l");
|
| 434 |
|
|
legexp->Draw();
|
| 435 |
|
|
|
| 436 |
|
|
CLsObsNLO->Draw("l");
|
| 437 |
|
|
CLsObsLO->Draw("l");
|
| 438 |
|
|
CLsExpNLO->Draw("l");
|
| 439 |
|
|
//FCExpNLO->Draw("l");
|
| 440 |
|
|
//CLsExpNoSNLO->Draw("l");
|
| 441 |
|
|
// HybridCLsObsNLO->Draw("l");
|
| 442 |
|
|
//gobsexcl->Draw("l");
|
| 443 |
|
|
//gHybridCLsObsExclMHT->Draw("l");
|
| 444 |
|
|
|
| 445 |
|
|
|
| 446 |
|
|
stau->Draw("fsame");
|
| 447 |
|
|
b.DrawLatex( 10,380,"#tilde{#tau} LSP");
|
| 448 |
|
|
//constant ssqquark and gluino lines
|
| 449 |
|
|
TF1* lnsq[4];
|
| 450 |
|
|
TF1* lngl[4];
|
| 451 |
|
|
TLatex sqt; sqt.SetTextSize(0.02); sqt.SetTextAngle(-22);sqt.SetTextColor(kGray+2);
|
| 452 |
|
|
sqt.DrawLatex(148,220,"#font[92]{#tilde{q}(500)GeV}");
|
| 453 |
|
|
sqt.DrawLatex(390,350,"#font[92]{#tilde{q}(800)GeV}");
|
| 454 |
|
|
TLatex glt; glt.SetTextSize(0.02); sqt.SetTextAngle(-4); glt.SetTextColor(kGray+2);
|
| 455 |
|
|
glt.DrawLatex(840,172,"#font[92]{#tilde{g}(500)GeV}");
|
| 456 |
|
|
glt.DrawLatex(840,235,"#font[92]{#tilde{g}(650)GeV}");
|
| 457 |
|
|
glt.DrawLatex(838,297,"#font[92]{#tilde{g}(800)GeV}");
|
| 458 |
|
|
int tanBeta_=10;
|
| 459 |
|
|
for(int i = 0; i < 4; i++){
|
| 460 |
|
|
lnsq[i] = constant_squark(tanBeta_,i);
|
| 461 |
|
|
lngl[i] = constant_gluino(tanBeta_,i);
|
| 462 |
|
|
lngl[i]->Draw("same");
|
| 463 |
|
|
lnsq[i]->Draw("same");
|
| 464 |
|
|
}
|
| 465 |
|
|
TLegend * leg = new TLegend(0.25,0.72,0.68,0.88);
|
| 466 |
|
|
leg->SetBorderSize(0);leg->SetFillColor(0);leg->SetFillStyle(4000);leg->SetTextFont(42);leg->SetTextSize(0.025);
|
| 467 |
|
|
TGraph * expLeg = (TGraph*)CLsExpNLO->Clone();expLeg->SetFillStyle(gCLsExp1Sigma->GetFillStyle());expLeg->SetFillColor(gCLsExp1Sigma->GetFillColor());
|
| 468 |
|
|
leg->SetHeader("L_{int} = 36 pb^{-1}, #sqrt{s} = 7 TeV");
|
| 469 |
|
|
leg->AddEntry(CLsObsNLO,"Observed, NLO","l");
|
| 470 |
|
|
leg->AddEntry(CLsObsLO, "Observed, LO","l");
|
| 471 |
|
|
leg->AddEntry(expLeg, "Expected #pm 1#sigma, NLO","lf");
|
| 472 |
|
|
//leg->AddEntry(CLsExpNoSNLO, "Expected, no-signal hyp., NLO","l");
|
| 473 |
|
|
leg->Draw();
|
| 474 |
|
|
gPad->RedrawAxis();
|
| 475 |
|
|
c1->SaveAs("results/Exclusion_m0_m12_tb10.pdf");
|
| 476 |
|
|
|
| 477 |
|
|
|
| 478 |
|
|
//plotToolsHT->Print(NLOSignal, Mzero,Mhalf, gCLsObsExclHT, 5);
|
| 479 |
|
|
// plotToolsHT->Print(NLOXsection, Mzero,Mhalf, MGluino, MSquarkL, gCLsObsExclHT, 5);
|
| 480 |
|
|
|
| 481 |
|
|
hexcl->Draw("");
|
| 482 |
|
|
gCLsExp1Sigma->SetName("CLsNLO_Exp1SigmaUncertaintyTb10");
|
| 483 |
|
|
CLsObsNLO->SetName("CLsNLO_ObservedTb10");
|
| 484 |
|
|
CLsExpNLO->SetName("CLsNLO_ExpectedTb10");
|
| 485 |
|
|
CLsExpLO->SetName("CLsLO_ExpectedTb10");
|
| 486 |
|
|
gCLsExp1Sigma->Draw("l");
|
| 487 |
|
|
CLsObsNLO->Draw("l");
|
| 488 |
|
|
CLsExpNLO->Draw("l");
|
| 489 |
|
|
CLsObsLO->Draw("l");
|
| 490 |
|
|
c1->SaveAs("results/LimitContours_tb10.C");
|
| 491 |
|
|
// ===============================================================================================================
|
| 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 |
|
|
TEV_sg_cdf->Draw("fsame");
|
| 515 |
|
|
TEV_sg_cdf->Draw("lsame");
|
| 516 |
|
|
TEV_sg_d0->Draw("fsame");
|
| 517 |
|
|
TEV_sg_d0->Draw("lsame");
|
| 518 |
|
|
//ms.DrawLatex(750,458,"tan#beta=10, #mu>0, A_{0}=0");
|
| 519 |
|
|
//ms.DrawLatex(10,458,"L_{int} = 36/pb, #sqrt{s} = 7 TeV");
|
| 520 |
|
|
LM0->SetMarkerSize(0.7);
|
| 521 |
|
|
LM1->SetMarkerSize(0.7);
|
| 522 |
|
|
tLM0->SetTextSize(0.02);
|
| 523 |
|
|
tLM1->SetTextSize(0.02);
|
| 524 |
|
|
//LM0->Draw("same");
|
| 525 |
|
|
//tLM0->Draw("same");
|
| 526 |
|
|
LM1->Draw("same");
|
| 527 |
|
|
tLM1->Draw("same");
|
| 528 |
|
|
gCLsExpNoS1Sigma->Draw("lf,same");
|
| 529 |
|
|
sFirst->Draw("same");
|
| 530 |
|
|
//Atlas->Draw("c,same");
|
| 531 |
|
|
CLsObsNLO->Draw("l,same");
|
| 532 |
|
|
CLsObsLO->Draw("l,same");
|
| 533 |
|
|
CLsExpNoSNLO->Draw("l,same");
|
| 534 |
|
|
//FCExpNLO->Draw("l");
|
| 535 |
|
|
stau->Draw("fsame");
|
| 536 |
|
|
b.DrawLatex( 10,380,"#tilde{#tau} LSP");
|
| 537 |
|
|
//constant ssqquark and gluino lines
|
| 538 |
|
|
sqt.DrawLatex(148,220,"#font[92]{#tilde{q}(500)GeV}");
|
| 539 |
|
|
sqt.DrawLatex(390,350,"#font[92]{#tilde{q}(800)GeV}");
|
| 540 |
|
|
glt.DrawLatex(840,172,"#font[92]{#tilde{g}(500)GeV}");
|
| 541 |
|
|
glt.DrawLatex(840,235,"#font[92]{#tilde{g}(650)GeV}");
|
| 542 |
|
|
glt.DrawLatex(838,297,"#font[92]{#tilde{g}(800)GeV}");
|
| 543 |
|
|
for(int i = 0; i < 4; i++){
|
| 544 |
|
|
lngl[i]->Draw("same");
|
| 545 |
|
|
lnsq[i]->Draw("same");
|
| 546 |
|
|
}
|
| 547 |
|
|
legexp->Draw();
|
| 548 |
|
|
leg->Draw();
|
| 549 |
|
|
gPad->RedrawAxis();
|
| 550 |
|
|
c1->SaveAs("results/Exclusion_m0_m12_tb10_NoSigHypPseudoData.pdf");
|
| 551 |
|
|
|
| 552 |
|
|
|
| 553 |
|
|
hexcl->Draw("");
|
| 554 |
|
|
gCLsExp1Sigma->SetName("CLsNLO_Exp1SigmaUncertaintyTb10");
|
| 555 |
|
|
CLsObsNLO->SetName("CLsNLO_ObservedTb10");
|
| 556 |
|
|
CLsExpNoSNLO->SetName("CLsNLO_ExpectedTb10");
|
| 557 |
|
|
CLsObsLO->SetName("CLsLO_ObservedTb10");
|
| 558 |
|
|
//gCLsExpNoS1Sigma->Draw("l");
|
| 559 |
|
|
CLsObsNLO->Draw("l,same");
|
| 560 |
|
|
//CLsObsLO->Draw("l,same");
|
| 561 |
|
|
CLsExpNoSNLO->Draw("l,same");
|
| 562 |
|
|
c1->SaveAs("results/LimitContours_tb10_NoSigHypPseudoData.C");
|
| 563 |
|
|
|
| 564 |
|
|
hexcl->Draw("colz");
|
| 565 |
|
|
LEP_ch->Draw("fsame");
|
| 566 |
|
|
LEP_sl->Draw("fsame");
|
| 567 |
|
|
TEV_sg_cdf->Draw("fsame");
|
| 568 |
|
|
TEV_sg_cdf->Draw("lsame");
|
| 569 |
|
|
TEV_sg_d0->Draw("fsame");
|
| 570 |
|
|
TEV_sg_d0->Draw("lsame");
|
| 571 |
|
|
gCLsExpNoS1Sigma->Draw("lf,same");
|
| 572 |
|
|
sFirst->Draw("same");
|
| 573 |
|
|
Atlas->Draw("c,same");
|
| 574 |
|
|
CLsObsNLO->Draw("l,same");
|
| 575 |
|
|
CLsExpNoSNLO->Draw("l,same");
|
| 576 |
|
|
ms.DrawLatex(750,458,"tan #beta=10, #mu>0, A_{0}=0");
|
| 577 |
|
|
LM0->SetMarkerSize(0.7);
|
| 578 |
|
|
LM1->SetMarkerSize(0.7);
|
| 579 |
|
|
tLM0->SetTextSize(0.02);
|
| 580 |
|
|
tLM1->SetTextSize(0.02);
|
| 581 |
|
|
LM0->Draw("same");
|
| 582 |
|
|
tLM0->Draw("same");
|
| 583 |
|
|
LM1->Draw("same");
|
| 584 |
|
|
tLM1->Draw("same");
|
| 585 |
|
|
stau->Draw("fsame");
|
| 586 |
|
|
b.DrawLatex( 10,380,"#tilde{#tau} LSP");
|
| 587 |
|
|
//Jim's limits///////////////////////////////////////////////////////////////////////////////
|
| 588 |
|
|
TGraph* JimObsHT = Jim_ht_tb10(0);
|
| 589 |
|
|
TGraph* JimObsMHT = Jim_mht_tb10(0);
|
| 590 |
|
|
TGraph* JimExpHT = Jim_ht_tb10(1);
|
| 591 |
|
|
TGraph* JimExpMHT = Jim_mht_tb10(1);
|
| 592 |
|
|
TGraph* JimExpHTup = Jim_ht_tb10(2);
|
| 593 |
|
|
TGraph* JimExpMHTup = Jim_mht_tb10(2);
|
| 594 |
|
|
TGraph* JimExpHTdn = Jim_ht_tb10(3);
|
| 595 |
|
|
TGraph* JimExpMHTdn = Jim_mht_tb10(3);
|
| 596 |
|
|
Smooth(JimObsHT, 2);
|
| 597 |
|
|
Smooth(JimObsMHT,2);
|
| 598 |
|
|
Smooth(JimExpHT, 2);
|
| 599 |
|
|
Smooth(JimExpMHT,2);
|
| 600 |
|
|
Smooth(JimExpHTup, 2);
|
| 601 |
|
|
Smooth(JimExpMHTup,2);
|
| 602 |
|
|
Smooth(JimExpHTdn, 2);
|
| 603 |
|
|
Smooth(JimExpMHTdn,2);
|
| 604 |
|
|
TGraph * JimObs = plotToolsHT->ChooseBest(JimObsHT,JimObsMHT,JimObsHT,JimObsMHT);
|
| 605 |
|
|
TGraph * JimExp = plotToolsHT->ChooseBest(JimExpHT,JimExpMHT,JimExpHT,JimExpMHT);
|
| 606 |
|
|
TGraph * JimExpup = plotToolsHT->ChooseBest(JimExpHTup,JimExpMHTup,JimExpHTup,JimExpMHTup);
|
| 607 |
|
|
TGraph * JimExpdn = plotToolsHT->ChooseBest(JimExpHTdn,JimExpMHTdn,JimExpHTdn,JimExpMHTdn);
|
| 608 |
|
|
TGraph * JimExp1Sigma = MakeBand(JimExpup, JimExpdn);JimExp1Sigma->SetFillStyle(3001);
|
| 609 |
|
|
TGraph * JimLeg = (TGraph*)JimExp->Clone();JimLeg->SetFillStyle(JimExp1Sigma->GetFillStyle());JimLeg->SetFillColor(JimExp1Sigma->GetFillColor());
|
| 610 |
|
|
JimExp1Sigma->Draw("f");
|
| 611 |
|
|
JimObs->Draw("c");
|
| 612 |
|
|
JimExp->Draw("c");
|
| 613 |
|
|
stau->Draw("fsame");
|
| 614 |
|
|
b.DrawLatex( 10,380,"#tilde{#tau} LSP");
|
| 615 |
|
|
//constant ssqquark and gluino lines
|
| 616 |
|
|
sqt.DrawLatex(148,220,"#font[92]{#tilde{q}(500)GeV}");
|
| 617 |
|
|
sqt.DrawLatex(390,350,"#font[92]{#tilde{q}(800)GeV}");
|
| 618 |
|
|
glt.DrawLatex(840,172,"#font[92]{#tilde{g}(500)GeV}");
|
| 619 |
|
|
glt.DrawLatex(840,235,"#font[92]{#tilde{g}(650)GeV}");
|
| 620 |
|
|
glt.DrawLatex(838,297,"#font[92]{#tilde{g}(800)GeV}");
|
| 621 |
|
|
for(int i = 0; i < 4; i++){
|
| 622 |
|
|
lngl[i]->Draw("same");
|
| 623 |
|
|
lnsq[i]->Draw("same");
|
| 624 |
|
|
}
|
| 625 |
|
|
legexp->Draw();
|
| 626 |
|
|
TLegend * legBayes = new TLegend(0.25,0.730,0.6,0.88);
|
| 627 |
|
|
legBayes->SetBorderSize(0);legBayes->SetFillColor(0);legBayes->SetFillStyle(4000);legBayes->SetTextFont(42);
|
| 628 |
|
|
legBayes->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
|
| 629 |
|
|
legBayes->AddEntry(JimObs, "Observed, Bayes","l");
|
| 630 |
|
|
legBayes->AddEntry(JimLeg, "Expected #pm 1#sigma, Bayes","lf");
|
| 631 |
|
|
legBayes->AddEntry(CLsObsNLO,"Observed, CLs","l");
|
| 632 |
|
|
legBayes->AddEntry(expLeg, "Expected #pm 1#sigma, CLs","lf");
|
| 633 |
|
|
legBayes->Draw();
|
| 634 |
|
|
gPad->RedrawAxis();
|
| 635 |
|
|
c1->SaveAs("results/Exclusion_m0_m12_tb10_Bayes.pdf");
|
| 636 |
|
|
|
| 637 |
|
|
|
| 638 |
|
|
|
| 639 |
|
|
// Signal Contamination in M0 - M1/2
|
| 640 |
|
|
c1->SetLogz(0);
|
| 641 |
|
|
c1->SetRightMargin ( 0.2 );
|
| 642 |
|
|
TH2F*hsigcont = new TH2F("sigcont",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
|
| 643 |
|
|
100,0,1009.9,35,100,450);
|
| 644 |
|
|
plotToolsHT->Area(hsigcont, Mzero, Mhalf, SignalContamination);
|
| 645 |
|
|
hsigcont->SetMinimum(0.01);
|
| 646 |
|
|
hsigcont->SetMaximum(20);
|
| 647 |
|
|
hsigcont->Draw("colz");
|
| 648 |
|
|
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 649 |
|
|
c1->SaveAs("results/SignalContamination_HT_m0_m12_tb10.pdf");
|
| 650 |
|
|
|
| 651 |
|
|
// Signal Contamination in M0 - M1/2
|
| 652 |
|
|
c1->SetLogz(0);
|
| 653 |
|
|
TH2F*hsigcontMHT = new TH2F("sigcontMHT",";m_{0} [GeV]; m_{1/2}; number of signal in bkg yield [events]",
|
| 654 |
|
|
100,0,1009.9,35,100,450);
|
| 655 |
|
|
plotTools->Area(hsigcontMHT, Mzero, Mhalf, SignalContamination);
|
| 656 |
|
|
hsigcontMHT->SetMinimum(0.01);
|
| 657 |
|
|
hsigcontMHT->SetMaximum(20);
|
| 658 |
|
|
hsigcontMHT->Draw("colz");
|
| 659 |
|
|
if (gobsexcl) gobsexcl->Draw("l");
|
| 660 |
|
|
c1->SaveAs("results/SignalContamination_MHT_m0_m12_tb10.pdf");
|
| 661 |
|
|
|
| 662 |
|
|
// Signal JEC Uncertainty MHT
|
| 663 |
|
|
c1->SetLogz(0);
|
| 664 |
|
|
TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
|
| 665 |
|
|
100,0,1009.9,35,100,450);
|
| 666 |
|
|
plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
|
| 667 |
|
|
hsig_jec_mht->SetMinimum(0.0);
|
| 668 |
|
|
hsig_jec_mht->SetMaximum(0.25);
|
| 669 |
|
|
hsig_jec_mht->SetContour(10);
|
| 670 |
|
|
hsig_jec_mht->Draw("colz");
|
| 671 |
|
|
if (gobsexcl) gobsexcl->Draw("l");
|
| 672 |
|
|
c1->SaveAs("results/SigJEC_MHT_m0_m12_tb10.pdf");
|
| 673 |
|
|
|
| 674 |
|
|
// Signal JEC Uncertainty HT
|
| 675 |
|
|
c1->SetLogz(0);
|
| 676 |
|
|
TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal JEC uncertainty",
|
| 677 |
|
|
100,0,1009.9,35,100,450);
|
| 678 |
|
|
plotToolsHT->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
|
| 679 |
|
|
hsig_jec_ht->SetMinimum(0.0);
|
| 680 |
|
|
hsig_jec_ht->SetMaximum(0.25);
|
| 681 |
|
|
hsig_jec_ht->SetContour(10);
|
| 682 |
|
|
hsig_jec_ht->Draw("colz");
|
| 683 |
|
|
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 684 |
|
|
c1->SaveAs("results/SigJEC_HT_m0_m12_tb10.pdf");
|
| 685 |
|
|
|
| 686 |
|
|
|
| 687 |
|
|
// Signal MuIso Uncertainty MHT
|
| 688 |
|
|
c1->SetLogz(0);
|
| 689 |
|
|
TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
|
| 690 |
|
|
100,0,1009.9,35,100,450);
|
| 691 |
|
|
plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
|
| 692 |
|
|
hsig_MuIso_mht->SetMinimum(0.0);hsig_MuIso_mht->GetZaxis()->SetTitleOffset(1.3);
|
| 693 |
|
|
hsig_MuIso_mht->SetMaximum(0.05);
|
| 694 |
|
|
hsig_MuIso_mht->SetContour(5);
|
| 695 |
|
|
hsig_MuIso_mht->Draw("colz");
|
| 696 |
|
|
if (gobsexcl) gobsexcl->Draw("l");
|
| 697 |
|
|
c1->SaveAs("results/SigMuIso_MHT_m0_m12_tb10.pdf");
|
| 698 |
|
|
|
| 699 |
|
|
// Signal MuIso Uncertainty HT
|
| 700 |
|
|
c1->SetLogz(0);
|
| 701 |
|
|
TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal #mu_{iso} uncertainty",
|
| 702 |
|
|
100,0,1009.9,35,100,450);
|
| 703 |
|
|
plotToolsHT->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
|
| 704 |
|
|
hsig_MuIso_ht->SetMinimum(0.0);hsig_MuIso_ht->GetZaxis()->SetTitleOffset(1.3);
|
| 705 |
|
|
hsig_MuIso_ht->SetMaximum(0.05);
|
| 706 |
|
|
hsig_MuIso_ht->SetContour(5);
|
| 707 |
|
|
hsig_MuIso_ht->Draw("colz");
|
| 708 |
|
|
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 709 |
|
|
c1->SaveAs("results/SigMuIso_HT_m0_m12_tb10.pdf");
|
| 710 |
|
|
|
| 711 |
|
|
|
| 712 |
|
|
|
| 713 |
|
|
// Signal PDF Uncertainty MHT
|
| 714 |
|
|
c1->SetLogz(0);
|
| 715 |
|
|
TH2F*hsig_PDF_mht = new TH2F("sigPDFmht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal PDF uncertainty",
|
| 716 |
|
|
100,0,1009.9,35,100,450);
|
| 717 |
|
|
plotTools->Area(hsig_PDF_mht, Mzero, Mhalf, SignalUncertPDF);
|
| 718 |
|
|
hsig_PDF_mht->SetMinimum(0.0);hsig_PDF_mht->GetZaxis()->SetTitleOffset(1.3);
|
| 719 |
|
|
hsig_PDF_mht->SetMaximum(0.20);
|
| 720 |
|
|
hsig_PDF_mht->SetContour(10);
|
| 721 |
|
|
hsig_PDF_mht->Draw("colz");
|
| 722 |
|
|
if (gobsexcl) gobsexcl->Draw("l");
|
| 723 |
|
|
c1->SaveAs("results/SigPDF_MHT_m0_m12_tb10.pdf");
|
| 724 |
|
|
|
| 725 |
|
|
// Signal PDF Uncertainty HT
|
| 726 |
|
|
c1->SetLogz(0);
|
| 727 |
|
|
TH2F*hsig_PDF_ht = new TH2F("sigPDFht",";m_{0} [GeV]; m_{1/2} [GeV]; Signal PDF uncertainty",
|
| 728 |
|
|
100,0,1009.9,35,100,450);
|
| 729 |
|
|
plotToolsHT->Area(hsig_PDF_ht, Mzero, Mhalf, SignalUncertPDF);
|
| 730 |
|
|
hsig_PDF_ht->SetMinimum(0.0);hsig_PDF_ht->GetZaxis()->SetTitleOffset(1.3);
|
| 731 |
|
|
hsig_PDF_ht->SetMaximum(0.20);
|
| 732 |
|
|
hsig_PDF_ht->SetContour(10);
|
| 733 |
|
|
hsig_PDF_ht->Draw("colz");
|
| 734 |
|
|
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 735 |
|
|
c1->SaveAs("results/SigPDF_HT_m0_m12_tb10.pdf");
|
| 736 |
|
|
|
| 737 |
|
|
|
| 738 |
|
|
// Signal kFactor Uncertainty MHT
|
| 739 |
|
|
c1->SetLogz(0);
|
| 740 |
|
|
TH2F*hsig_kFactorUncert_mht = new TH2F("sigkFactorUncertmht",";m_{0} [GeV]; m_{1/2} [GeV]; signal scale uncertainty",
|
| 741 |
|
|
100,0,1009.9,35,100,450);
|
| 742 |
|
|
plotTools->Area(hsig_kFactorUncert_mht, Mzero, Mhalf, SignalUncertKfactor);
|
| 743 |
|
|
hsig_kFactorUncert_mht->SetMinimum(0.1);hsig_kFactorUncert_mht->GetZaxis()->SetTitleOffset(1.3);
|
| 744 |
|
|
hsig_kFactorUncert_mht->SetMaximum(0.2);
|
| 745 |
|
|
hsig_kFactorUncert_mht->SetContour(10);
|
| 746 |
|
|
hsig_kFactorUncert_mht->Draw("colz");
|
| 747 |
|
|
if (gobsexcl) gobsexcl->Draw("l");
|
| 748 |
|
|
c1->SaveAs("results/SigkFactorUncert_MHT_m0_m12_tb10.pdf");
|
| 749 |
|
|
|
| 750 |
|
|
// Signal kFactor Uncertainty HT
|
| 751 |
|
|
c1->SetLogz(0);
|
| 752 |
|
|
TH2F*hsig_kFactorUncert_ht = new TH2F("sigkFactorUncertht",";m_{0} [GeV]; m_{1/2} [GeV]; signal scale uncertainty",
|
| 753 |
|
|
100,0,1009.9,35,100,450);
|
| 754 |
|
|
plotToolsHT->Area(hsig_kFactorUncert_ht, Mzero, Mhalf, SignalUncertKfactor);
|
| 755 |
|
|
hsig_kFactorUncert_ht->SetMinimum(0.1);hsig_kFactorUncert_ht->GetZaxis()->SetTitleOffset(1.3);
|
| 756 |
|
|
hsig_kFactorUncert_ht->SetMaximum(0.2);
|
| 757 |
|
|
hsig_kFactorUncert_ht->SetContour(10);
|
| 758 |
|
|
hsig_kFactorUncert_ht->Draw("colz");
|
| 759 |
|
|
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 760 |
|
|
c1->SaveAs("results/SigkFactorUncert_HT_m0_m12_tb10.pdf");
|
| 761 |
|
|
|
| 762 |
|
|
// Signal kFactor HT
|
| 763 |
|
|
c1->SetLogz(0);
|
| 764 |
|
|
TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{0} [GeV]; m_{1/2} [GeV]; NLO signal k-factor",
|
| 765 |
|
|
100,0,1009.9,35,100,450);
|
| 766 |
|
|
plotToolsHT->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalKfactor);
|
| 767 |
|
|
hsig_kFactor_ht->SetMinimum(1.0);hsig_kFactor_ht->GetZaxis()->SetTitleOffset(1.3);
|
| 768 |
|
|
hsig_kFactor_ht->SetMaximum(2.4);
|
| 769 |
|
|
hsig_kFactor_ht->SetContour(7);
|
| 770 |
|
|
hsig_kFactor_ht->Draw("colz");
|
| 771 |
|
|
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 772 |
|
|
c1->SaveAs("results/SigkFactor_HT_m0_m12_tb10.pdf");
|
| 773 |
|
|
|
| 774 |
|
|
// Signal kFactor MHT
|
| 775 |
|
|
c1->SetLogz(0);
|
| 776 |
|
|
TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{0} [GeV]; m_{1/2} [GeV]; NLO signal k-factor",
|
| 777 |
|
|
100,0,1009.9,35,100,450);
|
| 778 |
|
|
plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalKfactor);
|
| 779 |
|
|
hsig_kFactor_mht->SetMinimum(1.0);hsig_kFactor_ht->GetZaxis()->SetTitleOffset(1.3);
|
| 780 |
|
|
hsig_kFactor_mht->SetMaximum(2.4);
|
| 781 |
|
|
hsig_kFactor_mht->SetContour(7);
|
| 782 |
|
|
hsig_kFactor_mht->Draw("colz");
|
| 783 |
|
|
if (gobsexcl) gobsexcl->Draw("l");
|
| 784 |
|
|
c1->SaveAs("results/SigkFactor_MHT_m0_m12_tb10.pdf");
|
| 785 |
|
|
|
| 786 |
|
|
|
| 787 |
|
|
|
| 788 |
|
|
//-----------------------------------------------------------------------------------
|
| 789 |
|
|
c1->SetLogz(1);
|
| 790 |
|
|
|
| 791 |
|
|
TGraph * mz500 = plotMasses->Line( MSquarkL, MGluino, Mzero, 500, 1);
|
| 792 |
|
|
TGraph * mh250 = plotMasses->Line( MSquarkL, MGluino, Mhalf, 250, 1);
|
| 793 |
|
|
|
| 794 |
|
|
// cross-section in squark - gluino mass
|
| 795 |
|
|
TH2F*hxsec_qg = new TH2F("xsec_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; cross section [pb]",
|
| 796 |
|
|
60,200,1400,50,200,1200);
|
| 797 |
|
|
plotTools->Area(hxsec_qg, MSquarkL, MGluino, Xsection);
|
| 798 |
|
|
hxsec_qg->SetMinimum(0.01);
|
| 799 |
|
|
hxsec_qg->Draw("colz");
|
| 800 |
|
|
mz500->Draw();
|
| 801 |
|
|
mh250->Draw();
|
| 802 |
|
|
c1->SaveAs("results/Xsection_mSql_mGl_tb10.pdf");
|
| 803 |
|
|
|
| 804 |
|
|
// Observed Limit in squark - gluino mass
|
| 805 |
|
|
TH2F*hobslimit_qg = new TH2F("obslimit_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Limit [pb]",
|
| 806 |
|
|
60,200,1400,50,200,1200);
|
| 807 |
|
|
plotTools->Area(hobslimit_qg, MSquarkL, MGluino, ObsXsecLimit);
|
| 808 |
|
|
hobslimit_qg->SetMinimum(0.01);
|
| 809 |
|
|
hobslimit_qg->Draw("colz");
|
| 810 |
|
|
c1->SaveAs("results/ObsLimit_mSql_mGl_tb10.pdf");
|
| 811 |
|
|
|
| 812 |
|
|
// Expected Limit in squark - gluino mass
|
| 813 |
|
|
TH2F*hexplimit_qg = new TH2F("explimit_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Expected Limit [pb]",
|
| 814 |
|
|
60,200,1400,50,200,1200);
|
| 815 |
|
|
plotTools->Area(hexplimit_qg, MSquarkL, MGluino, ExpXsecLimit);
|
| 816 |
|
|
hexplimit_qg->SetMinimum(0.01);
|
| 817 |
|
|
hexplimit_qg->Draw("colz");
|
| 818 |
|
|
c1->SaveAs("results/ExpLimit_mSql_mGl_tb10.pdf");
|
| 819 |
|
|
|
| 820 |
|
|
// Signal Acceptance in squark - gluino mass
|
| 821 |
|
|
TH2F*hsigacc_qg = new TH2F("sigacc_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Acceptance",
|
| 822 |
|
|
60,200,1400,50,200,1200);
|
| 823 |
|
|
plotTools->Area(hsigacc_qg, MSquarkL, MGluino, SignalAcceptance);
|
| 824 |
|
|
hsigacc_qg->SetMinimum(0.01);
|
| 825 |
|
|
hsigacc_qg->SetMaximum(1.0);
|
| 826 |
|
|
hsigacc_qg->Draw("colz");
|
| 827 |
|
|
c1->SaveAs("results/SigAccMHT_mSql_mGl_tb10.pdf");
|
| 828 |
|
|
|
| 829 |
|
|
// Signal Acceptance in squark - gluino mass
|
| 830 |
|
|
TH2F*hsigaccHT_qg = new TH2F("sigaccHT_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Acceptance",
|
| 831 |
|
|
60,200,1400,50,200,1200);
|
| 832 |
|
|
plotToolsHT->Area(hsigaccHT_qg, MSquarkL, MGluino, SignalAcceptance);
|
| 833 |
|
|
hsigaccHT_qg->SetMinimum(0.01);
|
| 834 |
|
|
hsigaccHT_qg->SetMaximum(1.0);
|
| 835 |
|
|
hsigaccHT_qg->Draw("colz");
|
| 836 |
|
|
c1->SaveAs("results/SigAccHT_mSql_mGl_tb10.pdf");
|
| 837 |
|
|
|
| 838 |
|
|
|
| 839 |
|
|
// Signal Contamination in squark - gluino mass
|
| 840 |
|
|
c1->SetLogz(0);
|
| 841 |
|
|
TH2F*hsigcont_qg = new TH2F("sigacc_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Contamination in bkg yield",
|
| 842 |
|
|
60,200,1400,50,200,1200);
|
| 843 |
|
|
plotTools->Area(hsigcont_qg, MSquarkL, MGluino, SignalContamination);
|
| 844 |
|
|
hsigcont_qg->SetMinimum(0.0);
|
| 845 |
|
|
hsigcont_qg->SetMaximum(20.0);
|
| 846 |
|
|
hsigcont_qg->Draw("colz");
|
| 847 |
|
|
c1->SaveAs("results/SignalContamination_MHT_mSql_mGl_tb10.pdf");
|
| 848 |
|
|
|
| 849 |
|
|
// Signal Contamination in squark - gluino mass
|
| 850 |
|
|
TH2F*hsigcontHT_qg = new TH2F("sigaccHT_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal Acceptance in bkg yield",
|
| 851 |
|
|
60,200,1400,50,200,1200);
|
| 852 |
|
|
plotToolsHT->Area(hsigcontHT_qg, MSquarkL, MGluino, SignalContamination);
|
| 853 |
|
|
hsigcontHT_qg->SetMinimum(0.0);
|
| 854 |
|
|
hsigcontHT_qg->SetMaximum(20.0);
|
| 855 |
|
|
hsigcontHT_qg->Draw("colz");
|
| 856 |
|
|
c1->SaveAs("results/SignalContamination_HT_mSql_mGl_tb10.pdf");
|
| 857 |
|
|
/*
|
| 858 |
|
|
// Signal JEC Uncertainty MHT
|
| 859 |
|
|
c1->SetLogz(0);
|
| 860 |
|
|
TH2F*hsig_jec_mht = new TH2F("sigjecmht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal JEC uncertainty",
|
| 861 |
|
|
100,0,1009.9,35,100,450);
|
| 862 |
|
|
plotTools->Area(hsig_jec_mht, Mzero, Mhalf, SignalUncertJEC);
|
| 863 |
|
|
hsig_jec_mht->SetMinimum(0.0);
|
| 864 |
|
|
hsig_jec_mht->SetMaximum(0.25);
|
| 865 |
|
|
hsig_jec_mht->SetContour(10);
|
| 866 |
|
|
hsig_jec_mht->Draw("colz");
|
| 867 |
|
|
if (gobsexcl) gobsexcl->Draw("l");
|
| 868 |
|
|
c1->SaveAs("results/SigJEC_MHT_m0_m12_tb10.pdf");
|
| 869 |
|
|
|
| 870 |
|
|
// Signal JEC Uncertainty HT
|
| 871 |
|
|
c1->SetLogz(0);
|
| 872 |
|
|
TH2F*hsig_jec_ht = new TH2F("sigjecht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal JEC uncertainty",
|
| 873 |
|
|
100,0,1009.9,35,100,450);
|
| 874 |
|
|
plotTools->Area(hsig_jec_ht, Mzero, Mhalf, SignalUncertJEC);
|
| 875 |
|
|
hsig_jec_ht->SetMinimum(0.0);
|
| 876 |
|
|
hsig_jec_ht->SetMaximum(0.25);
|
| 877 |
|
|
hsig_jec_ht->SetContour(10);
|
| 878 |
|
|
hsig_jec_ht->Draw("colz");
|
| 879 |
|
|
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 880 |
|
|
c1->SaveAs("results/SigJEC_HT_m0_m12_tb10.pdf");
|
| 881 |
|
|
|
| 882 |
|
|
|
| 883 |
|
|
// Signal MuIso Uncertainty MHT
|
| 884 |
|
|
c1->SetLogz(0);
|
| 885 |
|
|
TH2F*hsig_MuIso_mht = new TH2F("sigMuIsomht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
|
| 886 |
|
|
100,0,1009.9,35,100,450);
|
| 887 |
|
|
plotTools->Area(hsig_MuIso_mht, Mzero, Mhalf, SignalUncertMuIso);
|
| 888 |
|
|
hsig_MuIso_mht->SetMinimum(0.0);
|
| 889 |
|
|
hsig_MuIso_mht->SetMaximum(0.05);
|
| 890 |
|
|
hsig_MuIso_mht->SetContour(5);
|
| 891 |
|
|
hsig_MuIso_mht->Draw("colz");
|
| 892 |
|
|
if (gobsexcl) gobsexcl->Draw("l");
|
| 893 |
|
|
c1->SaveAs("results/SigMuIso_MHT_m0_m12_tb10.pdf");
|
| 894 |
|
|
|
| 895 |
|
|
// Signal MuIso Uncertainty HT
|
| 896 |
|
|
c1->SetLogz(0);
|
| 897 |
|
|
TH2F*hsig_MuIso_ht = new TH2F("sigMuIsoht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
|
| 898 |
|
|
100,0,1009.9,35,100,450);
|
| 899 |
|
|
plotTools->Area(hsig_MuIso_ht, Mzero, Mhalf, SignalUncertMuIso);
|
| 900 |
|
|
hsig_MuIso_ht->SetMinimum(0.0);
|
| 901 |
|
|
hsig_MuIso_ht->SetMaximum(0.05);
|
| 902 |
|
|
hsig_MuIso_ht->SetContour(5);
|
| 903 |
|
|
hsig_MuIso_ht->Draw("colz");
|
| 904 |
|
|
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 905 |
|
|
c1->SaveAs("results/SigMuIso_HT_m0_m12_tb10.pdf");
|
| 906 |
|
|
|
| 907 |
|
|
|
| 908 |
|
|
// Signal kFactor Uncertainty MHT
|
| 909 |
|
|
c1->SetLogz(0);
|
| 910 |
|
|
TH2F*hsig_kFactor_mht = new TH2F("sigkFactormht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
|
| 911 |
|
|
100,0,1009.9,35,100,450);
|
| 912 |
|
|
plotTools->Area(hsig_kFactor_mht, Mzero, Mhalf, SignalUncertKfactor);
|
| 913 |
|
|
hsig_kFactor_mht->SetMinimum(0.1);
|
| 914 |
|
|
hsig_kFactor_mht->SetMaximum(0.2);
|
| 915 |
|
|
hsig_kFactor_mht->SetContour(10);
|
| 916 |
|
|
hsig_kFactor_mht->Draw("colz");
|
| 917 |
|
|
if (gobsexcl) gobsexcl->Draw("l");
|
| 918 |
|
|
c1->SaveAs("results/SigkFactor_MHT_m0_m12_tb10.pdf");
|
| 919 |
|
|
|
| 920 |
|
|
// Signal kFactor Uncertainty HT
|
| 921 |
|
|
c1->SetLogz(0);
|
| 922 |
|
|
TH2F*hsig_kFactor_ht = new TH2F("sigkFactorht",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; Signal #mu_{iso} uncertainty",
|
| 923 |
|
|
100,0,1009.9,35,100,450);
|
| 924 |
|
|
plotTools->Area(hsig_kFactor_ht, Mzero, Mhalf, SignalUncertKfactor);
|
| 925 |
|
|
hsig_kFactor_ht->SetMinimum(0.1);
|
| 926 |
|
|
hsig_kFactor_ht->SetMaximum(0.2);
|
| 927 |
|
|
hsig_kFactor_ht->SetContour(10);
|
| 928 |
|
|
hsig_kFactor_ht->Draw("colz");
|
| 929 |
|
|
if (gCLsObsExclHT) gCLsObsExclHT->Draw("l");
|
| 930 |
|
|
c1->SaveAs("results/SigkFactor_HT_m0_m12_tb10.pdf");
|
| 931 |
|
|
|
| 932 |
|
|
*/
|
| 933 |
|
|
// Exp. Limit on Number of Signal Events in squark - gluino mass
|
| 934 |
|
|
c1->SetLogz(0);
|
| 935 |
|
|
TH2F*hexplimitnsig_qg = new TH2F("explimitnsig_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL exp. limit signal events [# ]",
|
| 936 |
|
|
60,200,1400,50,200,1200);
|
| 937 |
|
|
plotTools->Area(hexplimitnsig_qg, MSquarkL, MGluino, ExpNSignLimit);
|
| 938 |
|
|
hexplimitnsig_qg->SetMinimum(0.0);
|
| 939 |
|
|
hexplimitnsig_qg->SetMaximum(20);
|
| 940 |
|
|
hexplimitnsig_qg->Draw("colz");
|
| 941 |
|
|
c1->SaveAs("results/ExpLimitOnNSig_mSql_mGl_tb10.pdf");
|
| 942 |
|
|
|
| 943 |
|
|
// Obs. Limit on Number of Signal Events in squark - gluino mass
|
| 944 |
|
|
TH2F*hobslimitnsig_qg = new TH2F("obslimitnsig_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL obs. limit signal events [# ]",
|
| 945 |
|
|
60,200,1400,50,200,1200);
|
| 946 |
|
|
plotTools->Area(hobslimitnsig_qg, MSquarkL, MGluino, ObsNSignLimit);
|
| 947 |
|
|
hobslimitnsig_qg->SetMinimum(0.0);
|
| 948 |
|
|
hobslimitnsig_qg->SetMaximum(20);
|
| 949 |
|
|
hobslimitnsig_qg->Draw("colz");
|
| 950 |
|
|
c1->SaveAs("results/ObsLimitOnNSig_mSql_mGl_tb10.pdf");
|
| 951 |
|
|
|
| 952 |
|
|
c1->SetLogz(0);
|
| 953 |
|
|
// Expected Exclusion in squark - gluino mass
|
| 954 |
|
|
TH2F*hexpexcl_qg = new TH2F("expexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Expected Exclusion",
|
| 955 |
|
|
60,200,1400,50,200,1200);
|
| 956 |
|
|
plotTools->Area(hexpexcl_qg, MSquarkL, MGluino, ExpExclCL);
|
| 957 |
|
|
hexpexcl_qg->Draw("colz");
|
| 958 |
|
|
c1->SaveAs("results/ExpExclusion_mSql_mGl_tb10.pdf");
|
| 959 |
|
|
|
| 960 |
|
|
// Observed Exclusion in squark - gluino mass
|
| 961 |
|
|
TH2F*hobsexcl_qg = new TH2F("obsexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
|
| 962 |
|
|
60,200,1400,50,200,1200);
|
| 963 |
|
|
plotTools->Area(hobsexcl_qg, MSquarkL, MGluino, ObsExclCL);
|
| 964 |
|
|
//std::vector<TGraph *> vobsexcl_qg = plotTools->GetContours(hobsexcl_qg);
|
| 965 |
|
|
hobsexcl_qg->Draw("colz");
|
| 966 |
|
|
//for (std::vector<TGraph *>::iterator g=vobsexcl_qg.begin();g!=vobsexcl_qg.end();++g){
|
| 967 |
|
|
// if (! *g) continue;
|
| 968 |
|
|
// // (*g)->Draw("l");
|
| 969 |
|
|
//}
|
| 970 |
|
|
c1->SaveAs("results/ObsExclusion_mSql_mGl_tb10.pdf");
|
| 971 |
|
|
|
| 972 |
|
|
|
| 973 |
|
|
// Observed Exclusion in squark - gluino mass
|
| 974 |
|
|
TH2F*hPLobsexcl_qg = new TH2F("plobsexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
|
| 975 |
|
|
60,200,1400,50,200,1200);
|
| 976 |
|
|
plotTools->Area(hPLobsexcl_qg, MSquarkL, MGluino, PLObsExclusion);
|
| 977 |
|
|
//hPLobsexcl_qg->Draw("colz");
|
| 978 |
|
|
//c1->SaveAs("results/PL_ObsExclusion_mSql_mGl_tb10.pdf");
|
| 979 |
|
|
|
| 980 |
|
|
|
| 981 |
|
|
|
| 982 |
|
|
// Exclusion in squark - gluino mass ----------------------------------------
|
| 983 |
|
|
c1->SetRightMargin ( 0.05 );
|
| 984 |
|
|
TH2F*h_qg = new TH2F("hExcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
|
| 985 |
|
|
60,0,1200,50,0,1200);
|
| 986 |
|
|
TH2F*hExcl_qg = new TH2F("h_qg","",60,200,1400,50,200,1200);
|
| 987 |
|
|
TH2F*hexpFC_HT_LO_qg=(TH2F*)hExcl_qg->Clone();
|
| 988 |
|
|
plotToolsHT->Area(hexpFC_HT_LO_qg, MSquarkL, MGluino, FCExpExclusion);
|
| 989 |
|
|
TH2F*hexpFC_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
|
| 990 |
|
|
plotToolsHT->Area(hexpFC_HT_NLO_qg, MSquarkL, MGluino, NLOFCExpExclusion);
|
| 991 |
|
|
TH2F*hexpFC_MHT_LO_qg=(TH2F*)hExcl_qg->Clone();
|
| 992 |
|
|
plotTools->Area(hexpFC_MHT_LO_qg, MSquarkL, MGluino, FCExpExclusion);
|
| 993 |
|
|
TH2F*hexpFC_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
|
| 994 |
|
|
plotTools->Area(hexpFC_MHT_NLO_qg, MSquarkL, MGluino, NLOFCExpExclusion);
|
| 995 |
|
|
//
|
| 996 |
|
|
TH2F*hexpCLs_HT_LO_qg=(TH2F*)hExcl_qg->Clone();
|
| 997 |
|
|
plotToolsHT->Area(hexpCLs_HT_LO_qg, MSquarkL, MGluino, ExpExclCL);
|
| 998 |
|
|
TH2F*hexpCLs_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
|
| 999 |
|
|
plotToolsHT->Area(hexpCLs_HT_NLO_qg, MSquarkL, MGluino, NLOExpExclCL);
|
| 1000 |
|
|
TH2F*hexpCLs_MHT_LO_qg=(TH2F*)hExcl_qg->Clone();
|
| 1001 |
|
|
plotTools->Area(hexpCLs_MHT_LO_qg, MSquarkL, MGluino, ExpExclCL);
|
| 1002 |
|
|
TH2F*hexpCLs_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
|
| 1003 |
|
|
plotTools->Area(hexpCLs_MHT_NLO_qg, MSquarkL, MGluino, NLOExpExclCL);
|
| 1004 |
|
|
//
|
| 1005 |
|
|
TH2F*hexp1p_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
|
| 1006 |
|
|
plotToolsHT->Area(hexp1p_HT_NLO_qg, MSquarkL, MGluino, NLOExpExclCLp1sigma);
|
| 1007 |
|
|
TH2F*hexp1p_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
|
| 1008 |
|
|
plotTools->Area(hexp1p_MHT_NLO_qg, MSquarkL, MGluino, NLOExpExclCLp1sigma);
|
| 1009 |
|
|
TH2F*hexp1m_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
|
| 1010 |
|
|
plotToolsHT->Area(hexp1m_HT_NLO_qg, MSquarkL, MGluino, NLOExpExclCLm1sigma);
|
| 1011 |
|
|
TH2F*hexp1m_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
|
| 1012 |
|
|
plotTools->Area(hexp1m_MHT_NLO_qg, MSquarkL, MGluino, NLOExpExclCLm1sigma);
|
| 1013 |
|
|
|
| 1014 |
|
|
TH2F*hobsCLs_HT_NLO_qg=(TH2F*)hExcl_qg->Clone();
|
| 1015 |
|
|
plotToolsHT->Area(hobsCLs_HT_NLO_qg, MSquarkL, MGluino, NLOObsExclCL);
|
| 1016 |
|
|
TH2F*hobsCLs_MHT_NLO_qg=(TH2F*)hExcl_qg->Clone();
|
| 1017 |
|
|
plotTools->Area(hobsCLs_MHT_NLO_qg, MSquarkL, MGluino, NLOObsExclCL);
|
| 1018 |
|
|
//
|
| 1019 |
|
|
TH2 * hexpFC_LO_qg = plotToolsHT->BinWiseOr(hexpFC_HT_LO_qg, hexpFC_MHT_LO_qg);
|
| 1020 |
|
|
TH2 * hexpCLs_LO_qg = plotToolsHT->BinWiseOr(hexpCLs_HT_LO_qg, hexpCLs_MHT_LO_qg);
|
| 1021 |
|
|
TH2 * hexpCLs_NLO_qg = plotToolsHT->BinWiseOr(hexpCLs_HT_NLO_qg, hexpCLs_MHT_NLO_qg);
|
| 1022 |
|
|
TH2 * hobsCLs_NLO_qg = plotToolsHT->BinWiseOr(hobsCLs_HT_NLO_qg, hobsCLs_MHT_NLO_qg);
|
| 1023 |
|
|
TH2 * hexp1p_NLO_qg = plotToolsHT->BinWiseOr(hexp1p_HT_NLO_qg, hexp1p_MHT_NLO_qg);
|
| 1024 |
|
|
TH2 * hexp1m_NLO_qg = plotToolsHT->BinWiseOr(hexp1m_HT_NLO_qg, hexp1m_MHT_NLO_qg);
|
| 1025 |
|
|
//
|
| 1026 |
|
|
TGraph * gexpFC_LO_qg = plotToolsHT->GetContour(hexpFC_LO_qg,3,0);
|
| 1027 |
|
|
gexpFC_LO_qg->SetLineStyle(4); gexpFC_LO_qg->SetLineColor(3);gexpFC_LO_qg->SetLineWidth(2);
|
| 1028 |
|
|
TGraph * gexpCLs_LO_qg = plotToolsHT->GetContour(hexpCLs_LO_qg,3,0);
|
| 1029 |
|
|
gexpCLs_LO_qg->SetLineStyle(4); gexpCLs_LO_qg->SetLineColor(1);gexpCLs_LO_qg->SetLineWidth(2);
|
| 1030 |
|
|
TGraph * gexpCLs_NLO_qg = plotToolsHT->GetContour(hexpCLs_NLO_qg,3,0);
|
| 1031 |
|
|
gexpCLs_NLO_qg->SetLineStyle(2); gexpCLs_NLO_qg->SetLineColor(1);gexpCLs_NLO_qg->SetLineWidth(2);
|
| 1032 |
|
|
TGraph * gobsCLs_NLO_qg = plotToolsHT->GetContour(hobsCLs_NLO_qg,3,0);
|
| 1033 |
|
|
gobsCLs_NLO_qg->SetLineStyle(1); gobsCLs_NLO_qg->SetLineColor(1);gobsCLs_NLO_qg->SetLineWidth(2);
|
| 1034 |
|
|
TGraph * gexp1p_NLO_qg = plotToolsHT->GetContour(hexp1p_NLO_qg,3,0);
|
| 1035 |
|
|
TGraph * gexp1m_NLO_qg = plotToolsHT->GetContour(hexp1m_NLO_qg,3,0);
|
| 1036 |
|
|
Smooth(gexpFC_LO_qg, 11);
|
| 1037 |
|
|
Smooth(gexpCLs_LO_qg, 11);
|
| 1038 |
|
|
Smooth(gexpCLs_NLO_qg, 11);
|
| 1039 |
|
|
Smooth(gobsCLs_NLO_qg, 11);
|
| 1040 |
|
|
Smooth(gexp1p_NLO_qg, 11);
|
| 1041 |
|
|
Smooth(gexp1m_NLO_qg, 11);
|
| 1042 |
|
|
|
| 1043 |
|
|
TGraph * g1Sigma_qg = MakeBand(gexp1p_NLO_qg, gexp1m_NLO_qg);
|
| 1044 |
|
|
g1Sigma_qg->SetFillColor(5);
|
| 1045 |
|
|
h_qg->GetYaxis()->SetTitleOffset(1.5);
|
| 1046 |
|
|
h_qg->Draw("");
|
| 1047 |
|
|
TGraph * sqLEP = sq_LEP(); sqLEP->Draw("f");
|
| 1048 |
|
|
TGraph * glTEV = gl_TEV(); glTEV->Draw("f");
|
| 1049 |
|
|
TGraph * glCDF = gl_CDF(); glCDF->Draw("f");
|
| 1050 |
|
|
TGraph * glDEZ = gl_DEZ(); glDEZ->Draw("f");
|
| 1051 |
|
|
|
| 1052 |
|
|
g1Sigma_qg->Draw("f");
|
| 1053 |
|
|
gexpFC_LO_qg->Draw("l");
|
| 1054 |
|
|
gexpCLs_LO_qg->Draw("l");
|
| 1055 |
|
|
gexpCLs_NLO_qg->Draw("l");
|
| 1056 |
|
|
gobsCLs_NLO_qg->Draw("l");
|
| 1057 |
|
|
//TGraph * glWHT = gl_WHT(); glWHT->Draw("f");
|
| 1058 |
|
|
b.DrawLatex( 800, 50 ,"Tevatron RunI");
|
| 1059 |
|
|
b.DrawLatex( 800, 200,"CDF RunII");
|
| 1060 |
|
|
b.DrawLatex( 800, 240,"Dzero RunII");
|
| 1061 |
|
|
b.DrawLatex( 100, 600,"LEP #tilde{q}#tilde{q}");
|
| 1062 |
|
|
//
|
| 1063 |
|
|
TLegend * l_qg = new TLegend(0.45,0.7,0.85,0.89);
|
| 1064 |
|
|
l_qg->SetBorderSize(0);l_qg->SetFillColor(0);l_qg->SetTextFont(42);
|
| 1065 |
|
|
l_qg->SetHeader("L_{int} = 36/pb, #sqrt{s} = 7 TeV");
|
| 1066 |
|
|
TGraph * expL_qg = (TGraph*)gexpCLs_NLO_qg->Clone();expL_qg->SetFillStyle(g1Sigma_qg->GetFillStyle());expL_qg->SetFillColor(g1Sigma_qg->GetFillColor());
|
| 1067 |
|
|
if (expL_qg) l_qg->AddEntry(expL_qg, "Expected (CLs, NLO)","lf");
|
| 1068 |
|
|
if (gobsCLs_NLO_qg) l_qg->AddEntry(gobsCLs_NLO_qg,"Observed (CLs, NLO)","l");
|
| 1069 |
|
|
if (gexpCLs_LO_qg) l_qg->AddEntry(gexpCLs_LO_qg, "Expected (CLs, LO)","l");
|
| 1070 |
|
|
if (gexpFC_LO_qg) l_qg->AddEntry(gexpFC_LO_qg, "Expected (FC, LO)","l");
|
| 1071 |
|
|
l_qg->Draw();
|
| 1072 |
|
|
gPad->RedrawAxis();
|
| 1073 |
|
|
c1->SaveAs("results/Exclusion_mSql_mGl_tb10.pdf");
|
| 1074 |
|
|
|
| 1075 |
|
|
//====================================================================================================================
|
| 1076 |
|
|
// Exclusion in gluino - squark mass ----------------------------------------
|
| 1077 |
|
|
c1->SetTopMargin ( 0.05 );
|
| 1078 |
|
|
TH2F*h_gq = new TH2F("hExcl_gq",";m_{#tilde{g}} (GeV); m_{#tilde{q}} (GeV); 95% CL Observed Exclusion",
|
| 1079 |
|
|
50,0,900,60,0,1000);
|
| 1080 |
|
|
TH2F*hExcl_gq = new TH2F("h_gq","",50,200,1200,60,200,1400);
|
| 1081 |
|
|
/* TH2F*hexpFC_HT_LO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1082 |
|
|
plotToolsHT->Area(hexpFC_HT_LO_gq, MGluino, MSquarkL, FCExpExclusion);
|
| 1083 |
|
|
TH2F*hexpFC_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1084 |
|
|
plotToolsHT->Area(hexpFC_HT_NLO_gq, MGluino, MSquarkL, NLOFCExpExclusion);
|
| 1085 |
|
|
TH2F*hexpFC_MHT_LO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1086 |
|
|
plotTools->Area(hexpFC_MHT_LO_gq, MGluino, MSquarkL, FCExpExclusion);
|
| 1087 |
|
|
TH2F*hexpFC_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1088 |
|
|
plotTools->Area(hexpFC_MHT_NLO_gq, MGluino, MSquarkL, NLOFCExpExclusion);
|
| 1089 |
|
|
*/ //
|
| 1090 |
|
|
TH2F*hexpCLs_HT_LO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1091 |
|
|
plotToolsHT->Area(hexpCLs_HT_LO_gq, MGluino, MSquarkL, ExpExclCL);
|
| 1092 |
|
|
TH2F*hexpCLs_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1093 |
|
|
plotToolsHT->Area(hexpCLs_HT_NLO_gq, MGluino, MSquarkL, NLOExpExclCL);
|
| 1094 |
|
|
TH2F*hexpCLs_MHT_LO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1095 |
|
|
plotTools->Area(hexpCLs_MHT_LO_gq, MGluino, MSquarkL, ExpExclCL);
|
| 1096 |
|
|
TH2F*hexpCLs_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1097 |
|
|
plotTools->Area(hexpCLs_MHT_NLO_gq, MGluino, MSquarkL, NLOExpExclCL);
|
| 1098 |
|
|
//
|
| 1099 |
|
|
TH2F*hexp1p_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1100 |
|
|
plotToolsHT->Area(hexp1p_HT_NLO_gq, MGluino, MSquarkL, NLOExpExclCLp1sigma);
|
| 1101 |
|
|
TH2F*hexp1p_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1102 |
|
|
plotTools->Area(hexp1p_MHT_NLO_gq, MGluino, MSquarkL, NLOExpExclCLp1sigma);
|
| 1103 |
|
|
TH2F*hexp1m_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1104 |
|
|
plotToolsHT->Area(hexp1m_HT_NLO_gq, MGluino, MSquarkL, NLOExpExclCLm1sigma);
|
| 1105 |
|
|
TH2F*hexp1m_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1106 |
|
|
plotTools->Area(hexp1m_MHT_NLO_gq, MGluino, MSquarkL, NLOExpExclCLm1sigma);
|
| 1107 |
|
|
|
| 1108 |
|
|
TH2F*hobsCLs_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1109 |
|
|
plotToolsHT->Area(hobsCLs_HT_NLO_gq, MGluino, MSquarkL, NLOObsExclCL);
|
| 1110 |
|
|
TH2F*hobsCLs_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1111 |
|
|
plotTools->Area(hobsCLs_MHT_NLO_gq, MGluino, MSquarkL, NLOObsExclCL);
|
| 1112 |
|
|
TH2F*hobsCLs_HT_LO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1113 |
|
|
plotToolsHT->Area(hobsCLs_HT_LO_gq, MGluino, MSquarkL, ObsExclCL);
|
| 1114 |
|
|
TH2F*hobsCLs_MHT_LO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1115 |
|
|
plotTools->Area(hobsCLs_MHT_LO_gq, MGluino, MSquarkL, ObsExclCL);
|
| 1116 |
|
|
//
|
| 1117 |
|
|
//TH2 * hexpFC_LO_gq = plotToolsHT->BinWiseOr(hexpFC_HT_LO_gq, hexpFC_MHT_LO_gq);
|
| 1118 |
|
|
TH2 * hexpCLs_LO_gq = plotToolsHT->BinWiseOr(hexpCLs_HT_LO_gq, hexpCLs_MHT_LO_gq);
|
| 1119 |
|
|
TH2 * hexpCLs_NLO_gq = plotToolsHT->BinWiseOr(hexpCLs_HT_NLO_gq, hexpCLs_MHT_NLO_gq);
|
| 1120 |
|
|
TH2 * hobsCLs_NLO_gq = plotToolsHT->BinWiseOr(hobsCLs_HT_NLO_gq, hobsCLs_MHT_NLO_gq);
|
| 1121 |
|
|
TH2 * hobsCLs_LO_gq = plotToolsHT->BinWiseOr(hobsCLs_HT_LO_gq, hobsCLs_MHT_LO_gq);
|
| 1122 |
|
|
TH2 * hexp1p_NLO_gq = plotToolsHT->BinWiseOr(hexp1p_HT_NLO_gq, hexp1p_MHT_NLO_gq);
|
| 1123 |
|
|
TH2 * hexp1m_NLO_gq = plotToolsHT->BinWiseOr(hexp1m_HT_NLO_gq, hexp1m_MHT_NLO_gq);
|
| 1124 |
|
|
//
|
| 1125 |
|
|
//TGraph * gexpFC_LO_gq = plotToolsHT->GetContour(hexpFC_LO_gq,3,0);
|
| 1126 |
|
|
//gexpFC_LO_gq->SetLineStyle(4); gexpFC_LO_gq->SetLineColor(3);gexpFC_LO_gq->SetLineWidth(2);
|
| 1127 |
|
|
TGraph * gexpCLs_LO_gq = plotToolsHT->GetContour(hexpCLs_LO_gq,3,0);
|
| 1128 |
|
|
gexpCLs_LO_gq->SetLineStyle(4); gexpCLs_LO_gq->SetLineColor(1);gexpCLs_LO_gq->SetLineWidth(2);
|
| 1129 |
|
|
TGraph * gexpCLs_NLO_gq = plotToolsHT->GetContour(hexpCLs_NLO_gq,3,0);
|
| 1130 |
|
|
gexpCLs_NLO_gq->SetLineStyle(2); gexpCLs_NLO_gq->SetLineColor(1);gexpCLs_NLO_gq->SetLineWidth(2);
|
| 1131 |
|
|
TGraph * gobsCLs_NLO_gq = plotToolsHT->GetContour(hobsCLs_NLO_gq,3,0);
|
| 1132 |
|
|
gobsCLs_NLO_gq->SetLineStyle(1); gobsCLs_NLO_gq->SetLineColor(1);gobsCLs_NLO_gq->SetLineWidth(2);
|
| 1133 |
|
|
TGraph * gobsCLs_LO_gq = plotToolsHT->GetContour(hobsCLs_LO_gq,3,0);
|
| 1134 |
|
|
gobsCLs_LO_gq->SetLineStyle(1); gobsCLs_LO_gq->SetLineColor(2);gobsCLs_LO_gq->SetLineWidth(2);
|
| 1135 |
|
|
TGraph * gexp1p_NLO_gq = plotToolsHT->GetContour(hexp1p_NLO_gq,3,0);
|
| 1136 |
|
|
TGraph * gexp1m_NLO_gq = plotToolsHT->GetContour(hexp1m_NLO_gq,3,0);
|
| 1137 |
|
|
TGraph * ch1 = (TGraph*)gobsCLs_NLO_gq->Clone();
|
| 1138 |
|
|
//Smooth2D(gexpFC_LO_gq, 11);
|
| 1139 |
|
|
//Smooth2D(gexpCLs_LO_gq, 11);
|
| 1140 |
|
|
Smooth2D(gexpCLs_NLO_gq, 11);
|
| 1141 |
|
|
Smooth2D(gobsCLs_NLO_gq, 11);
|
| 1142 |
|
|
Smooth2D(gobsCLs_LO_gq, 11);
|
| 1143 |
|
|
Smooth2D(gexp1p_NLO_gq, 11);
|
| 1144 |
|
|
Smooth2D(gexp1m_NLO_gq, 11);
|
| 1145 |
|
|
//Close2D(gobsCLs_NLO_gq);
|
| 1146 |
|
|
TGraph * g1Sigma_gq = MakeBand(gexp1p_NLO_gq, gexp1m_NLO_gq);
|
| 1147 |
|
|
g1Sigma_gq->SetFillColor(5);
|
| 1148 |
|
|
h_gq->GetYaxis()->SetTitleOffset(1.5);
|
| 1149 |
|
|
h_gq->GetXaxis()->SetTitleOffset(0.98);
|
| 1150 |
|
|
h_gq->Draw("");
|
| 1151 |
|
|
TGraph * sqDEZ = sq_DEZ(); sqDEZ->Draw("f");
|
| 1152 |
|
|
TGraph * sqCDF = sq_CDF(); sqCDF->Draw("f");
|
| 1153 |
|
|
TGraph * sqTEV = sq_TEV(); sqTEV->Draw("f");
|
| 1154 |
|
|
TGraph * glLEP = gl_LEP(); glLEP->Draw("f");
|
| 1155 |
|
|
TGraph * noSol = glsq_NoSol(); noSol->Draw("f"); noSol->Draw("l");
|
| 1156 |
|
|
|
| 1157 |
|
|
ms.DrawLatex(10,1013,"L_{int} = 36 pb^{-1}, #sqrt{s} = 7 TeV");
|
| 1158 |
|
|
ms.DrawLatex(660,1013,"tan #beta=10, #mu>0, A_{0}=0");
|
| 1159 |
|
|
g1Sigma_gq->Draw("f");
|
| 1160 |
|
|
//ch1->Draw("l");
|
| 1161 |
|
|
// gexpFC_LO_gq->Draw("l");
|
| 1162 |
|
|
//gexpCLs_LO_gq->Draw("l");
|
| 1163 |
|
|
gexpCLs_NLO_gq->Draw("l");
|
| 1164 |
|
|
gobsCLs_NLO_gq->Draw("l");
|
| 1165 |
|
|
gobsCLs_LO_gq->Draw("l");
|
| 1166 |
|
|
TGraph * expL_gq = (TGraph*)gexpCLs_NLO_gq->Clone(); expL_gq->SetFillStyle(g1Sigma_gq->GetFillStyle());expL_gq->SetFillColor(g1Sigma_gq->GetFillColor());
|
| 1167 |
|
|
TLegend* legexp_gq = new TLegend(0.2,0.62,0.52,0.93,NULL,"brNDC");
|
| 1168 |
|
|
legexp_gq->SetFillColor(0);legexp_gq->SetShadowColor(0);//legexp_gq->SetFillStyle(4000);
|
| 1169 |
|
|
legexp_gq->SetTextFont(42);legexp_gq->SetTextSize(0.025);legexp_gq->SetBorderSize(1);
|
| 1170 |
|
|
//TEV_sg_cdf.SetLineColor(1);
|
| 1171 |
|
|
legexp_gq->SetHeader("CMS L_{int} = 36 pb^{-1}, #sqrt{s} = 7 TeV");
|
| 1172 |
|
|
if (gobsCLs_NLO_gq) legexp_gq->AddEntry(gobsCLs_NLO_gq,"Observed, NLO","l");
|
| 1173 |
|
|
if (gobsCLs_LO_gq) legexp_gq->AddEntry(gobsCLs_LO_gq,"Observed, LO","l");
|
| 1174 |
|
|
if (expL_gq) legexp_gq->AddEntry(expL_gq, "Expected #pm 1#sigma, NLO","lf");
|
| 1175 |
|
|
//if (gexpCLs_LO_gq) legexp_gq->AddEntry(gexpCLs_LO_gq, "Expected, LO","l");
|
| 1176 |
|
|
legexp_gq->AddEntry(sqTEV,"Tevatron RunI","f");
|
| 1177 |
|
|
legexp_gq->AddEntry(sqCDF,"CDF RunII","f");
|
| 1178 |
|
|
legexp_gq->AddEntry(sqDEZ,"D0 RunII","f");
|
| 1179 |
|
|
legexp_gq->AddEntry(glLEP,"LEP2","f"); //NOT FOR tb=50!
|
| 1180 |
|
|
legexp_gq->Draw();
|
| 1181 |
|
|
//TGraph * glWHT = gl_WHT(); glWHT->Draw("f");
|
| 1182 |
|
|
//b.DrawLatex( 800, 50,"LEP #tilde{q}#tilde{q}");
|
| 1183 |
|
|
//b.DrawLatex( 50, 1000 ,"Tevatron RunI");
|
| 1184 |
|
|
//b.DrawLatex( 200, 1100,"CDF RunII");
|
| 1185 |
|
|
//b.DrawLatex( 240, 1200,"Dzero RunII");
|
| 1186 |
|
|
|
| 1187 |
|
|
TPaveText tpave(0.6,0.33,0.75,0.4,"NDC");
|
| 1188 |
|
|
tpave.AddText("no CMSSM");tpave.AddText("solution");
|
| 1189 |
|
|
tpave.SetBorderSize(0);tpave.SetFillColor(0);tpave.SetTextFont(42);
|
| 1190 |
|
|
tpave.Draw();
|
| 1191 |
|
|
//
|
| 1192 |
|
|
gPad->RedrawAxis();
|
| 1193 |
|
|
c1->SaveAs("results/Exclusion_mGl_mSql_tb10.pdf");
|
| 1194 |
|
|
//====================================================================================================================
|
| 1195 |
|
|
|
| 1196 |
|
|
// TestContours in gluino - squark mass
|
| 1197 |
|
|
TH2F*texcl_gq = new TH2F("texcl_gq",";m_{#tilde{g}} (GeV); m_{#tilde{q}} (GeV); 95% CL Observed Exclusion",
|
| 1198 |
|
|
50,200,1200,60,200,1400);
|
| 1199 |
|
|
TH2F*area_gq=(TH2F*)texcl_gq->Clone();
|
| 1200 |
|
|
plotTools->Area(area_gq, MGluino, MSquarkL, NLOObsExclCL);
|
| 1201 |
|
|
std::vector<TGraph*> contours_qg = plotTools->GetContours(hobsCLs_NLO_gq,3);
|
| 1202 |
|
|
area_gq->Draw("colz");
|
| 1203 |
|
|
//hobsCLs_MHT_NLO_gq->Draw("colz,same");
|
| 1204 |
|
|
int col_gl=kBlue-10;
|
| 1205 |
|
|
for (std::vector<TGraph*>::iterator cont=contours_qg.begin(); cont!=contours_qg.end(); ++cont){
|
| 1206 |
|
|
if (! *cont) continue;
|
| 1207 |
|
|
double x, y;
|
| 1208 |
|
|
(*cont)->GetPoint(0, x, y);
|
| 1209 |
|
|
(*cont)->SetLineColor(col_gl);(*cont)->SetLineWidth(3);
|
| 1210 |
|
|
(*cont)->Draw("l");
|
| 1211 |
|
|
TLatex l; l.SetTextSize(0.04); l.SetTextColor(col_gl++);
|
| 1212 |
|
|
char val[20];
|
| 1213 |
|
|
sprintf(val,"%d",(int)(cont-contours_qg.begin()));
|
| 1214 |
|
|
l.DrawLatex(x,y,val);
|
| 1215 |
|
|
if (cont-contours_qg.begin()>3) break;
|
| 1216 |
|
|
}
|
| 1217 |
|
|
c1->SaveAs("results/ExclusionTestContours_mGl_mSq_tb10.pdf");
|
| 1218 |
|
|
|
| 1219 |
|
|
//====================================================================================================================
|
| 1220 |
|
|
// Exclusion in gluino - squark mass ----------------------------------------
|
| 1221 |
|
|
//
|
| 1222 |
|
|
TH2F*hexpNoSCLs_HT_NLO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1223 |
|
|
plotToolsHT->Area(hexpNoSCLs_HT_NLO_gq, MGluino, MSquarkL, NLOExpNoSigExclCL);
|
| 1224 |
|
|
TH2F*hexpNoSCLs_MHT_NLO_gq=(TH2F*)hExcl_gq->Clone();
|
| 1225 |
|
|
plotTools->Area(hexpNoSCLs_MHT_NLO_gq, MGluino, MSquarkL, NLOExpNoSigExclCL);
|
| 1226 |
|
|
//
|
| 1227 |
|
|
TH2 * hexpNoSCLs_NLO_gq = plotToolsHT->BinWiseOr(hexpNoSCLs_HT_NLO_gq, hexpNoSCLs_MHT_NLO_gq);
|
| 1228 |
|
|
TGraph * gexpNoSCLs_NLO_gq = plotToolsHT->GetContour(hexpNoSCLs_NLO_gq,3,0);
|
| 1229 |
|
|
gexpNoSCLs_NLO_gq->SetLineStyle(2); gexpNoSCLs_NLO_gq->SetLineColor(1);gexpNoSCLs_NLO_gq->SetLineWidth(2);
|
| 1230 |
|
|
TGraph * gexpNoS1p_NLO_gq = plotToolsHT->ModifyExpSigmaY(gexp1p_NLO_gq, gexpCLs_NLO_gq, gexpNoSCLs_NLO_gq);
|
| 1231 |
|
|
TGraph * gexpNoS1m_NLO_gq = plotToolsHT->ModifyExpSigmaY(gexp1m_NLO_gq, gexpCLs_NLO_gq, gexpNoSCLs_NLO_gq);
|
| 1232 |
|
|
TGraph * ch2 = (TGraph*)gobsCLs_NLO_gq->Clone();
|
| 1233 |
|
|
h_gq->Draw("");
|
| 1234 |
|
|
Smooth(gexpNoS1p_NLO_gq, 11);
|
| 1235 |
|
|
Smooth(gexpNoS1m_NLO_gq, 11);
|
| 1236 |
|
|
Smooth2D(gexpNoSCLs_NLO_gq, 11);
|
| 1237 |
|
|
TGraph * g1SigmaNoS_gq = MakeBand(gexpNoS1p_NLO_gq, gexpNoS1m_NLO_gq);
|
| 1238 |
|
|
g1SigmaNoS_gq->SetFillColor(5);
|
| 1239 |
|
|
sqDEZ->Draw("f");
|
| 1240 |
|
|
sqCDF->Draw("f");
|
| 1241 |
|
|
sqTEV->Draw("f");
|
| 1242 |
|
|
glLEP->Draw("f");
|
| 1243 |
|
|
noSol->Draw("f"); noSol->Draw("l");
|
| 1244 |
|
|
//ms.DrawLatex(10,1013,"L_{int} = 36/pb, #sqrt{s} = 7 TeV");
|
| 1245 |
|
|
ms.DrawLatex(650,940,"tan #beta=10, #mu>0, A_{0}=0");
|
| 1246 |
|
|
g1SigmaNoS_gq->Draw("f");
|
| 1247 |
|
|
//gexpNoS1p_NLO_gq->Draw("l");
|
| 1248 |
|
|
//gexpNoS1m_NLO_gq->Draw("l");
|
| 1249 |
|
|
gexpNoSCLs_NLO_gq->Draw("l");
|
| 1250 |
|
|
//gexpCLs_NLO_gq->Draw("l");
|
| 1251 |
|
|
gobsCLs_NLO_gq->Draw("l");
|
| 1252 |
|
|
gobsCLs_LO_gq->Draw("l");
|
| 1253 |
|
|
legexp_gq->Draw();
|
| 1254 |
|
|
tpave.Draw();
|
| 1255 |
|
|
gPad->RedrawAxis();
|
| 1256 |
|
|
c1->SaveAs("results/Exclusion_mGl_mSql_tb10_NoSigHypPseudoData.pdf");
|
| 1257 |
|
|
//====================================================================================================================
|
| 1258 |
|
|
|
| 1259 |
|
|
|
| 1260 |
|
|
|
| 1261 |
|
|
|
| 1262 |
|
|
|
| 1263 |
|
|
|
| 1264 |
|
|
TH2F*hexcl_qg = new TH2F("hexcl_qg",";m_{#tilde{q}} [GeV]; m_{#tilde{g}} [GeV]; 95% CL Observed Exclusion",
|
| 1265 |
|
|
60,200,1400,50,200,1000);
|
| 1266 |
|
|
TGraph * gexpexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL, 3,0, 2,2);
|
| 1267 |
|
|
TGraph * gexpexcl_qg_LO = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL, 3,0, 2,4);
|
| 1268 |
|
|
TGraph * gobsexcl_qg = plotTools->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL, 3,0, 2,1);
|
| 1269 |
|
|
TGraph * gMCMCexpexcl_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCExpExclusion,3,0, 3,2);
|
| 1270 |
|
|
TGraph * gMCMCobsexcl_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,MCMCObsExclusion,3,0, 3,1);
|
| 1271 |
|
|
TGraph * gCLsExpExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCL, 3,0, 1,2);
|
| 1272 |
|
|
TGraph * gCLsExpExclHT_qg_LO = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,ExpExclCL, 3,0, 1,4);
|
| 1273 |
|
|
TGraph * gCLsObsExclHT_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOObsExclCL, 3,0, 1,1);
|
| 1274 |
|
|
TGraph * gCLsExpExclHTm1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLm1sigma,3,0, 5,2);
|
| 1275 |
|
|
TGraph * gCLsExpExclHTp1_qg = plotToolsHT->GetContour(hexcl_qg,MSquarkL,MGluino,NLOExpExclCLp1sigma,3,0, 5,2);
|
| 1276 |
|
|
Smooth(gexpexcl_qg, 11);
|
| 1277 |
|
|
Smooth(gexpexcl_qg_LO, 11);
|
| 1278 |
|
|
Smooth(gobsexcl_qg, 11);
|
| 1279 |
|
|
Smooth(gCLsExpExclHT_qg, 11);
|
| 1280 |
|
|
Smooth(gCLsExpExclHT_qg_LO, 11);
|
| 1281 |
|
|
Smooth(gCLsObsExclHT_qg, 11);
|
| 1282 |
|
|
Smooth(gCLsExpExclHTm1_qg, 11);
|
| 1283 |
|
|
Smooth(gCLsExpExclHTp1_qg, 11);
|
| 1284 |
|
|
TGraph * gCLsExp1Sigma_qg = MakeBand(gCLsExpExclHTp1_qg, gCLsExpExclHTm1_qg);
|
| 1285 |
|
|
hexcl_qg->Draw("");
|
| 1286 |
|
|
//
|
| 1287 |
|
|
if (gCLsExp1Sigma_qg) gCLsExp1Sigma_qg->Draw("lf");
|
| 1288 |
|
|
if (gobsexcl_qg) gobsexcl_qg->Draw("l");
|
| 1289 |
|
|
if (gCLsObsExclHT_qg) gCLsObsExclHT_qg->Draw("l");
|
| 1290 |
|
|
if (gexpexcl_qg) gexpexcl_qg->Draw("l");
|
| 1291 |
|
|
if (gCLsExpExclHT_qg) gCLsExpExclHT_qg->Draw("l");
|
| 1292 |
|
|
if (gexpexcl_qg_LO) gexpexcl_qg_LO->Draw("l");
|
| 1293 |
|
|
if (gCLsExpExclHT_qg_LO) gCLsExpExclHT_qg_LO->Draw("l");
|
| 1294 |
|
|
|
| 1295 |
|
|
TLegend * leg_qg = new TLegend(0.45,0.78,0.85,0.89);
|
| 1296 |
|
|
leg_qg->SetBorderSize(0);leg_qg->SetFillColor(0);
|
| 1297 |
|
|
TGraph * expLeg_qg = (TGraph*)gCLsExpExclHT_qg->Clone();expLeg_qg->SetFillStyle(gCLsExp1Sigma_qg->GetFillStyle());expLeg_qg->SetFillColor(gCLsExp1Sigma_qg->GetFillColor());
|
| 1298 |
|
|
if (gobsexcl_qg) leg_qg->AddEntry(gobsexcl_qg,"NLO Observed (MHT, CLs)","l");
|
| 1299 |
|
|
if (gexpexcl_qg) leg_qg->AddEntry(gexpexcl_qg,"NLO Expected (MHT, CLs)","l");
|
| 1300 |
|
|
if (gexpexcl_qg_LO) leg_qg->AddEntry(gexpexcl_qg_LO,"LO Expected (MHT, CLs)","l");
|
| 1301 |
|
|
if (gCLsObsExclHT_qg) leg_qg->AddEntry(gCLsObsExclHT_qg,"NLO Observed (HT, CLs)","l");
|
| 1302 |
|
|
if (gCLsExpExclHT_qg) leg_qg->AddEntry(expLeg,"NLO Expected (HT, CLs)","lf");
|
| 1303 |
|
|
if (gCLsExpExclHT_qg_LO) leg_qg->AddEntry(gCLsExpExclHT_qg_LO,"LO Expected (HT, CLs)","lf");
|
| 1304 |
|
|
//if (gPLobsexcl_qg) leg_qg->AddEntry(gPLobsexcl_qg,"Observed (PL, RooStat)","l");
|
| 1305 |
|
|
//if (gPLexpexcl_qg) leg_qg->AddEntry(gPLexpexcl_qg,"Expected (PL, RooStat)","l");
|
| 1306 |
|
|
leg_qg->Draw();
|
| 1307 |
|
|
c1->SaveAs("results/Exclusion_mSql_mGl_tb10_details.pdf");
|
| 1308 |
|
|
//c1->SaveAs("plot_tb10.pdf");
|
| 1309 |
|
|
|
| 1310 |
|
|
c1->SetLogy(1);
|
| 1311 |
|
|
c1->SetLogx(1);
|
| 1312 |
|
|
TH2F*hUncertScan = new TH2F("hUncertScan",";relative signal uncertainty; number of signal events",
|
| 1313 |
|
|
10,0.1,0.6,10,0.9,200);
|
| 1314 |
|
|
//hUncertScan->GetXaxis()->SetNdivisions(505, false);
|
| 1315 |
|
|
//hUncertScan->SetMinimum(1);
|
| 1316 |
|
|
TGraph * gCLsExpUncert = new TGraph(0);
|
| 1317 |
|
|
TGraph * gCLsObsUncert = new TGraph(0);
|
| 1318 |
|
|
TGraph * gPLExpUncert = new TGraph(0);
|
| 1319 |
|
|
TGraph * gPLObsUncert = new TGraph(0);
|
| 1320 |
|
|
TGraph * gFCExpUncert = new TGraph(0);
|
| 1321 |
|
|
TGraph * gFCObsUncert = new TGraph(0);
|
| 1322 |
|
|
TGraph * gMCMCExpUncert= new TGraph(0);
|
| 1323 |
|
|
TGraph * gMCMCObsUncert= new TGraph(0);
|
| 1324 |
|
|
plotTools->Graph(gCLsExpUncert, SignalRelUncertainty, ExpNSignLimit,1);
|
| 1325 |
|
|
plotTools->Graph(gCLsObsUncert, SignalRelUncertainty, ObsNSignLimit,1);
|
| 1326 |
|
|
plotTools->Graph(gPLExpUncert, SignalRelUncertainty, PLExpNSignLimit,1);
|
| 1327 |
|
|
plotTools->Graph(gPLObsUncert, SignalRelUncertainty, PLObsNSignLimit, 1);
|
| 1328 |
|
|
plotTools->Graph(gFCExpUncert, SignalRelUncertainty, FCExpNSignLimit,1);
|
| 1329 |
|
|
plotTools->Graph(gFCObsUncert, SignalRelUncertainty, FCObsNSignLimit, 1);
|
| 1330 |
|
|
plotTools->Graph(gMCMCExpUncert,SignalRelUncertainty, MCMCExpNSignLimit,1);
|
| 1331 |
|
|
plotTools->Graph(gMCMCObsUncert,SignalRelUncertainty, MCMCObsNSignLimit, 1);
|
| 1332 |
|
|
gCLsExpUncert->SetLineStyle(8);
|
| 1333 |
|
|
gPLObsUncert->SetLineColor(2);
|
| 1334 |
|
|
gPLExpUncert->SetLineColor(2); gPLExpUncert->SetLineStyle(8);
|
| 1335 |
|
|
gFCObsUncert->SetLineColor(3);
|
| 1336 |
|
|
gFCExpUncert->SetLineColor(3); gFCExpUncert->SetLineStyle(8);
|
| 1337 |
|
|
gMCMCObsUncert->SetLineColor(4);
|
| 1338 |
|
|
gMCMCExpUncert->SetLineColor(4); gMCMCExpUncert->SetLineStyle(8);
|
| 1339 |
|
|
hUncertScan->Draw("h");
|
| 1340 |
|
|
gCLsExpUncert->Draw("l");
|
| 1341 |
|
|
gCLsObsUncert->Draw("l");
|
| 1342 |
|
|
//gPLObsUncert->Draw("l");
|
| 1343 |
|
|
gPLExpUncert->Draw("l");
|
| 1344 |
|
|
//gFCObsUncert->Draw("l");
|
| 1345 |
|
|
//gFCExpUncert->Draw("l");
|
| 1346 |
|
|
//gMCMCObsUncert->Draw("l");
|
| 1347 |
|
|
//gMCMCExpUncert->Draw("l");
|
| 1348 |
|
|
TLatex l; l.SetTextSize(0.03); l.SetTextColor(1);
|
| 1349 |
|
|
l.DrawLatex(0.025,25,"excluded");
|
| 1350 |
|
|
TMarker m; m.SetMarkerStyle(29); m.SetMarkerColor(2);
|
| 1351 |
|
|
m.DrawMarker(0.128, 28.5);l.DrawLatex(0.14,28.5,"LM1");
|
| 1352 |
|
|
m.DrawMarker(0.147, 77.4);l.DrawLatex(0.16,77.4,"LM0");
|
| 1353 |
|
|
|
| 1354 |
|
|
TLegend * leg_S = new TLegend(0.5,0.14,0.85,0.38);
|
| 1355 |
|
|
leg_S->SetBorderSize(0);leg_S->SetFillColor(0);
|
| 1356 |
|
|
leg_S->SetHeader("95% CL upper limits");
|
| 1357 |
|
|
leg_S->AddEntry(gCLsExpUncert,"Expected (CLs, TLimit)","l");
|
| 1358 |
|
|
leg_S->AddEntry(gCLsObsUncert,"Observed (CLs, TLimit)","l");
|
| 1359 |
|
|
//leg_S->AddEntry(gPLExpUncert,"Expected (PL, RooStat)","l");
|
| 1360 |
|
|
//leg_S->AddEntry(gPLObsUncert,"Observed (PL, RooStat)","l");
|
| 1361 |
|
|
//leg_S->AddEntry(gFCObsUncert,"Observed (FC, RooStat)","l");
|
| 1362 |
|
|
//leg_S->AddEntry(gFCExpUncert,"Expected (FC, RooStat)","l");
|
| 1363 |
|
|
//leg_S->AddEntry(gMCMCObsUncert,"Observed (MCMC, RooStat)","l");
|
| 1364 |
|
|
//leg_S->AddEntry(gMCMCExpUncert,"Expected (MCMC, RooStat)","l");
|
| 1365 |
|
|
leg_S->Draw();
|
| 1366 |
|
|
c1->SaveAs("results/UncertaintyScan.pdf");
|
| 1367 |
|
|
|
| 1368 |
|
|
|
| 1369 |
|
|
|
| 1370 |
|
|
PrintTools<SusyScan> * printTools = new PrintTools<SusyScan>(genpoints->GetScan());
|
| 1371 |
|
|
//PlotTools<SusyScan> * plotToolsHT = new PlotTools<SusyScan>(genpointsHT->GetScan());
|
| 1372 |
|
|
|
| 1373 |
|
|
//printTools->AddVar( Mzero );
|
| 1374 |
|
|
//printTools->AddVar( Mhalf );
|
| 1375 |
|
|
//printTools->AddVar( MChi1 );
|
| 1376 |
|
|
//printTools->AddMatch( new PrintTools<SusyScan>::equal(Mzero, 250, 5) );
|
| 1377 |
|
|
//printTools->AddMatch( new PrintTools<SusyScan>::equal(Mhalf, 250, 50) );
|
| 1378 |
|
|
//printTools->AddMatch( new PrintTools<SusyScan>::deltaR(Mhalf, 250, 50) );
|
| 1379 |
|
|
//printTools->Print( Mzero );
|
| 1380 |
|
|
|
| 1381 |
|
|
printTools->AddMatch( new PrintTools<SusyScan>::Xbinning(10, 10) );
|
| 1382 |
|
|
//printTools->Print( CLsObsNLO, "Mzero = ", "; Mhalf = " );
|
| 1383 |
|
|
|
| 1384 |
|
|
|
| 1385 |
|
|
|
| 1386 |
|
|
|
| 1387 |
|
|
//theApp.Run();
|
| 1388 |
|
|
}
|
| 1389 |
|
|
|
| 1390 |
|
|
|
| 1391 |
|
|
int main(int argc, char** argv)
|
| 1392 |
|
|
{
|
| 1393 |
|
|
return plot(argc, argv);
|
| 1394 |
|
|
}
|
