//______________________________________________________________________________
//
// PRD Plot Macro
//                                         
//
// Petra Haefner <haefner@fnal.gov>
// Generated 2006-03-06 
// Modified  2006-03-06
//______________________________________________________________________________


int plots_prd()
{ 

   gStyle->SetOptTitle(0);
   gStyle->SetOptStat(0); 
   gStyle->SetOptFit(0);
   gStyle->SetLineWidth(0.8);
   gStyle->SetFuncWidth(1.0);
   gStyle->SetFrameLineWidth(1.0);
   gStyle->SetHistLineWidth(1.0);
   gStyle->SetTitleOffset(1.5,"Y");       // axis description
   gStyle->SetEndErrorSize(1);
     
   // save histos in this file:
   TFile* fplots = new TFile("plots_prd.root", 
                                      "RECREATE", 
                                      "PRD plots");
    
//______________________________________________________________________________
//
//  psgn minus pbkg plots  
//______________________________________________________________________________
    
   TFile* f=new TFile("psgnminuspbkg.root","READ");
  
   TCanvas* PsgnMPbkgC = new TCanvas("PsgnMPbkgC","psgn minus pbkg",400,600);
   PsgnMPbkgC->Range(0,0,100,100);
   PsgnMPbkgC->SetFillColor(10);
   PsgnMPbkgC->SetBorderSize(0);
// PsgnMPbkgC->SetRightMargin(0.1);
   PsgnMPbkgC->Divide(1,2,0.1,0.);
// PsgnMPbkgC->Divide(1,2,0.,0);
             
   TH1F* hist_em       = (TH1F*)gROOT->FindObject("hpsgnminuspbkgem");
   TH1F* hist_em_njpm  = (TH1F*)gROOT->FindObject("hpsgnminuspbkgemnotjpm");
   TH1F* hist_em_wjjjj = (TH1F*)gROOT->FindObject("hpsgnminuspbkgem_wjjjj");
   TH1F* hist_mu       = (TH1F*)gROOT->FindObject("hpsgnminuspbkgmu");
   TH1F* hist_mu_njpm  = (TH1F*)gROOT->FindObject("hpsgnminuspbkgmunotjpm");
   TH1F* hist_mu_wjjjj = (TH1F*)gROOT->FindObject("hpsgnminuspbkgmu_wjjjj");
    
        
//  em plot 
//------------------------------------------------------------------------------
  
   PsgnMPbkgC->cd(1);
   gPad->SetTicks(1,1);
   gPad->SetLeftMargin(0.20);
   gPad->SetBottomMargin(0.020); // To see 0 on the Y axis of the top plot
  
   hist_em       -> GetYaxis()->SetRangeUser(0.0,0.15); 
   hist_em       -> GetYaxis()->SetNdivisions(508);
   hist_em       -> GetXaxis()->SetLabelOffset(99);
   hist_em       -> Draw();
   hist_em_njpm  -> SetFillStyle(3004);
   hist_em_njpm  -> SetFillColor(1);
   hist_em_njpm  -> Draw("same");
   hist_em_wjjjj -> SetLineStyle(2);
   hist_em_wjjjj -> Draw("same");
   gPad->RedrawAxis();
  
   hist_em       -> SetXTitle("log_{10}(P_{sgn} / P_{bkg})");
   hist_em       -> SetYTitle("dN/d(log_{10}(P_{sgn} / P_{bkg}))");

//  mu plot 
//------------------------------------------------------------------------------
  
   PsgnMPbkgC->cd(2);
   gPad->SetTicks(1,1);
   gPad->SetBottomMargin(0.15);
   gPad->SetLeftMargin(0.20);
  
  
   hist_mu       -> GetYaxis()->SetRangeUser(0.0,0.15); 
   hist_mu       -> GetYaxis()->SetTitleOffset(99.);
   hist_mu       -> GetYaxis()->SetNdivisions(508);
   hist_mu       -> Draw();
   hist_mu_njpm  -> SetFillStyle(3004);
   hist_mu_njpm  -> SetFillColor(1);
   hist_mu_njpm  -> Draw("same");
   hist_mu_wjjjj -> SetLineStyle(2);
   hist_mu_wjjjj -> Draw("same");
  
   gPad->RedrawAxis();
  
   hist_mu       -> SetXTitle("log_{10}(P_{sgn} / P_{bkg})");
   hist_mu       -> SetYTitle("dN/d(log_{10}(P_{sgn} / P_{bkg}))");

// save histograms:
   fplots->Write();
   return 0;
}