#include "TFile.h"
#include "TCanvas.h"
#include "TH1.h"
#include "TMath.h"


void patTau_test()
{
//--- open ROOT file
  TFile inputFile("patTau_Histograms_100GeV_scaled_gen.root");

//--- descend into directory containing histograms 
//    written by TFileService
  TDirectory* inputDir = (TDirectory*)inputFile.Get("analyzePatTau");

  // the following is valid using PF2PAT, were the TauJets are already "cleaned" and good
  TH1* hGenTauJetPt = (TH1*)inputDir->Get("GenTauPt");
  TH1* hGenTauJetEta = (TH1*)inputDir->Get("GenTauEta");
  TH1* hGenTauJetPhi = (TH1*)inputDir->Get("GenTauPhi");

  TH1* hTauJetPt = (TH1*)inputDir->Get("TauJetPt");
  TH1* hTauJetEta = (TH1*)inputDir->Get("TauJetEta");
  TH1* hTauJetPhi = (TH1*)inputDir->Get("TauJetPhi");

  hTauJetPt->Sumw2();
  hTauJetPt->Divide(hTauJetPt,hGenTauJetPt, 1., 1., "B");

  hTauJetEta->Sumw2();
  hTauJetEta->Divide(hTauJetEta,hGenTauJetEta, 1., 1., "B");

  hTauJetPhi->Sumw2();
  hTauJetPhi->Divide(hTauJetPhi,hGenTauJetPhi, 1., 1., "B");

  TCanvas * canvas = new TCanvas("plotPatTauEff", "plotPatTauEff", 800, 600);
  canvas->SetFillColor(10);
  canvas->cd(1);
  hTauJetPt->Draw("p");

//  TString outputFileName = TString("tauIdEff_100GeV_pt").Append(".gif");
//  canvas->Print(outputFileName.Data());

  inputFile.Close();
}