#include <iostream>
#include <fstream>
#include <cstdio>
#include "TFile.h"
#include "TMath.h"
#include "TH1.h"
#include "TCanvas.h"
#include "TF1.h"
#include "TGraph.h"
#include "TMultiGraph.h"
#include "TAxis.h"
#include "TLegend.h"

void
plot_graph() {

     TCanvas *c1 = new TCanvas("ewb", "ewb",281,93,700,500);
     c1->SetLogy();
     c1->SetGrid();
     TFile *theFile = new TFile("p_value_cuts_1600_max_order.root");
 
     TMultiGraph *p_value_chi_squared = (TMultiGraph*)theFile->Get("lnlike_cuts"); 
     
     TH1 *h1 = p_value_chi_squared->GetHistogram();
     
     h1->Rebin();
     
     std::vector<std::string> label;
     
     label.push_back("6.3x10^{-5}");
     label.push_back("1.0x10^{-4}");
     label.push_back("4.7x10^{-4}");
     label.push_back("0.0010");
     label.push_back("0.0027");
     label.push_back("0.0100");
     label.push_back("0.0460");
     label.push_back("0.0500");
     label.push_back("0.1000");
     label.push_back("0.3170");
     label.push_back("0.5000");
     label.push_back("0.7000");    
     label.push_back("0.9000"); 
     
     //p_value_chi_squared_cut1->Print();
     if (p_value_chi_squared != NULL) // check we got something useful
       p_value_chi_squared->Draw("ALP");
       
     Int_t i;
     for (i=1;i<=13;i++) {
       p_value_chi_squared->GetXaxis()->LabelsOption("d");
       int bin = p_value_chi_squared->GetXaxis()->FindBin(i);
       p_value_chi_squared->GetXaxis()->SetBinLabel(bin,label[i-1].c_str());
       p_value_chi_squared->GetYaxis()->SetRangeUser(0.00001,1);
     }
         
     p_value_chi_squared->SetTitle("1-Dim");
         
     TAxis *ay = p_value_chi_squared->GetYaxis();
     ay->LabelsOption("d");               
     ay->SetTitle("fraction of survivors");
     
     TAxis *ax = p_value_chi_squared->GetXaxis();
     ax->LabelsOption("d");               
     ax->SetTitle("cut on p-value");
     
     // lets draw it
     
     if (p_value_chi_squared != NULL) // check we got something useful
       
       p_value_chi_squared->Draw("ALP"); 

}