| 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 |
}
|
