#include <TROOT.h>
#include <TH1F.h>
#include <TFile.h>
#include <TList.h>
#include <TTree.h>

//using namespace std;
void extract(){
TFile *f = TFile::Open("flav_Akt4EMPf.root");
f->GetName();
TTree *t1 = (TTree*)f->Get("bTag_AntiKt4EMPFlowJets");
//TFile* fout = new TFile("analysis.root","RECREATE");

vector<float> *jet_mv2c10,*mymc_decayVtx_x,*mymc_decayVtx_y;
float_t truth_PVx,truth_PVy;
   t1->SetBranchAddress("jet_mv2c10",&jet_mv2c10);
   t1->SetBranchAddress("mymc_decayVtx_x",&mymc_decayVtx_x);
   t1->SetBranchAddress("mymc_decayVtx_y",&mymc_decayVtx_y);
   t1->SetBranchAddress("truth_PVx",&truth_PVx);
   t1->SetBranchAddress("truth_PVy",&truth_PVy);

//TList *l   = new TList();
TH1F *c1   = new TH1F("jet_mv2c10","jet_mv2c10",100,-1,1);
TH1F *c2   = new TH1F("mymc_decayVtx_x","jet_mv2c10",100,-1,1);
TH1F *c3   = new TH1F("mymc_decayVtx_y","jet_mv2c10",100,-1,1);
TH1F *c4   = new TH1F("truth_PVx","jet_mv2c10",100,-1,1);
TH1F *c5   = new TH1F("truth_PVy","jet_mv2c10",100,-1,1);

 Long64_t nentries = t1->GetEntries();
   for (Long64_t i=0;i<nentries;i++) {
     t1->GetEntry(i);
     c1->Fill((*jet_mv2c10)[i]);
     c2->Fill((*mymc_decayVtx_x)[i]);
     c3->Fill((*mymc_decayVtx_y)[i]);
     c4->Fill(truth_PVx);
     c5->Fill(truth_PVy);
     
  }
//l->Add(c1);
//l->Add(c2);
//l->Add(c3);
//l->Add(c4);
//l->Add(c5);
//l->Write("analysis", TObject::kSingleKey);
TFile* fout = new TFile("analysis.root","RECREATE");
c1->Write();
c2->Write();
c3->Write();
c4->Write();
c5->Write();
fout->Close();
}