#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_integrated_likelihood_ratio()
{
  TFile* f = new TFile("./integrated_likelihood_ratio.root");

  TH1F* h_lr_integral_all  = (TH1F*)gDirectory->Get("h_lr_toy_all_integral");
  TH1F* h_lr_integral_syst = (TH1F*)gDirectory->Get("h_lr_toy_95_integral_16_1");

  TF1 *func = new TF1("func", "exp(-x/2.0)", 0.0, 100.0);
  func->SetNpx(10000);
  func->SetLineColor(kRed);
  func->SetLineWidth(2);
  

  TCanvas* c0 = new TCanvas( "c0", "c0", 0, 0, 800, 800 );
  c0->Divide(1,1);
  c0->cd(1);
  h_lr_integral_all->SetLineWidth(3);
  h_lr_integral_all->SetXTitle("Likelihood Ratio     ");
  h_lr_integral_all->SetYTitle("Int(LR, inf) / Int(0, inf)      ");
  h_lr_integral_all->Draw("");
  h_lr_integral_syst->Draw("same");
  func->Draw("same");

}