#include <TH2.h>
#include <TStyle.h>
#include <TCanvas.h>
#include <TF1.h>
#include <TLine.h>
#include <TMath.h>
#include <TFile.h>

#include <iostream>
  
void plot_profile_likelihood()
{
  Int_t colors[2]={4, 2};
  gStyle->SetPalette(2, colors);
  
  const double math_pi = TMath::Pi();

  TFile* f = new TFile("./profile_likelihood.root");
 
  TH2F* h_bp_dg = (TH2F*)gDirectory->Get("h_bp_dg3");


  TMarker* marker_sm = new TMarker(0.02, 0.096, 8);
  marker_sm->SetMarkerColor(kBlack);
  marker_sm->SetMarkerSize(1.0);

  TLine* line_sm1 = new TLine( 0.02, 0.096-0.039, 0.02, 0.096+0.039 );
  line_sm1->SetLineWidth(2);
  line_sm1->SetLineColor(kBlack);
  line_sm1->SetLineStyle(1);

  TLine* line_sm2 = new TLine( 0.02-0.005, 0.096, 0.02+0.005, 0.096 );
  line_sm2->SetLineWidth(2);
  line_sm2->SetLineColor(kBlack);
  line_sm2->SetLineStyle(1);

  TLine* line = new TLine( math_pi/4.0, -0.7, math_pi/4.0, 0.7 );
  line->SetLineWidth(1);
  line->SetLineColor(kBlack);
  line->SetLineStyle(2);
 
  TLine* line2 = new TLine( -math_pi/2.0, 0.0, math_pi/2.0, 0.0 );
  line2->SetLineWidth(1);
  line2->SetLineColor(kBlack);
  line2->SetLineStyle(2);
 

  TCanvas* c0 = new TCanvas( "c0", "c0", 0, 0, 600, 600 );
  c0->Divide(1,1);
  c0->cd(1);
  h_bp_dg->SetXTitle("#beta_{s} (rad)     ");
  h_bp_dg->SetYTitle("#Delta#Gamma (ps^{-1})     ");
  h_bp_dg->SetLineWidth(2);
  h_bp_dg->SetLineStyle(1);
  h_bp_dg->SetContourLevel(0, 2.30);
  h_bp_dg->SetContourLevel(1, 5.99);
  h_bp_dg->Draw("cont1");
  //h_bp_dg->Draw("lego2");

  marker_sm->Draw("same");
  line->Draw("same");
  line2->Draw("same");
  line_sm1->Draw("same");
  line_sm2->Draw("same");



  TLatex* latex = new TLatex(-math_pi/2.0+0.05, 0.72, "CDF Run II Preliminary        L = 1.35 fb^{-1}");
  latex->SetTextSize(0.045);
  latex->Draw("same");

  TLatex* latex2 = new TLatex(-1.0, 0.55, "2#Deltalog(L) = 2.30");
  latex2->Draw("same");
  TLatex* latex3 = new TLatex(-1.0, 0.45, "2#Deltalog(L) = 5.99");
  latex3->Draw("same");
  TLatex* latex4 = new TLatex(-1.0, 0.35, "SM prediction");
  latex4->Draw("same");

  TLine* line_leg1 = new TLine( -1.3, 0.58, -1.1, 0.58 );
  line_leg1->SetLineWidth(2);
  line_leg1->SetLineColor(kBlue);
  line_leg1->SetLineStyle(1);

  TLine* line_leg2 = new TLine( -1.3, 0.48, -1.1, 0.48 );
  line_leg2->SetLineWidth(2);
  line_leg2->SetLineColor(kRed);
  line_leg2->SetLineStyle(1);

  TLine* line_leg3 = new TLine( -1.3, 0.38, -1.1, 0.38 );
  line_leg3->SetLineWidth(2);
  line_leg3->SetLineColor(kBlack);
  line_leg3->SetLineStyle(1);
  
  TMarker* marker_leg3 = new TMarker(-1.2, 0.38, 8);
  marker_leg3->SetMarkerColor(kBlack);
  marker_leg3->SetMarkerSize(1.0);

  line_leg1->Draw("same");
  line_leg2->Draw("same");
  line_leg3->Draw("same");
  marker_leg3->Draw("same");

}