#include <TROOT.h>
#include <TChain.h>
#include <TFile.h>
#include <TSelector.h>
#include <TTreeReader.h>
#include <TTreeReaderValue.h>
#include <TTreeReaderArray.h>
#include <TClass.h>
#include <TH2.h>

void nhat6() {
	TH1F* h_pxpy = new TH1F("h_pxpy", "PPAC", 50, 0, 100000);
	TH1F* h_pt = new TH1F("h_pt", "tPPAC", 1, 0.5, 2);

	// Open file root
	//TFile* outF = new TFile("PPAC1", "RECREATE");
	//TTree* outTree = new TTree("Tree1", "ppac");
	TFile* inF = new TFile("run01.root");
	TTree* inTree = (TTree*)inF->Get("AnalysisxTree");
	//inTree->SetBranchAddress(&AnalysisxTree)
	Double_t NeEvent_PPAC[16], NeEvent_tPPAC[16];
	
	
	inTree->SetBranchAddress("NeEvent_PPAC[16]", NeEvent_PPAC);
	inTree->SetBranchAddress("NeEvent_tPPAC[16]",NeEvent_tPPAC);
	Long64_t nentries = inTree->GetEntries();
	for (Long64_t  entry = 0; entry<inTree->GetEntries(); entry++) {
		inTree->GetEntry(entry);
			h_pxpy->Fill(NeEvent_PPAC[entry]);
			h_pt->Fill(NeEvent_tPPAC[entry]);
	}
	
	TCanvas* c1 = new TCanvas();
	
	h_pxpy ->Draw();
	//c1 ->SaveAs("nhat6_pxpy.png");

	//c1->Update();
	/*h_pt->Draw();
	//c1->SaveAs("nhat6_pt.png");
	outF->Write();
	outF->Close();
	inF->Close();*/
}