#include "TH2.h"
#include "TF1.h"
#include "TStyle.h"
#include "TCanvas.h"
#include "TPaveStats.h"
#include "TLine.h"
#include "TMath.h"
#include "TTree.h"
#include "TCutG.h"


int main()
{
	TCanvas *c1 = new TCanvas("c1","canvas");
	TFile * file = new TFile("list_ascii_time_correctedNewTreeRoot.root");
	TTree * tr = (TTree*)file -> Get("tr");
	ULong64_t Br_templong;
	long double Br_t100msec;
	ULong64_t Br_intpart;
	double Br_Ql, Br_Qs, psd;


	double temp;

	tr -> SetBranchAddress("Br_Ql",&Br_Ql);
	tr -> SetBranchAddress("psd",&psd);



	TH2D * hsep2 = new TH2D("hsep2","",65536,0,65536,1000,0,1);
	tr -> Draw("psd:Br_Ql >> hsep2","","colz");
	TCutG * mycutg = new TCutG("mycutg",7);
	mycutg->SetName("mycutg");
	mycutg->SetVarX("psd");
	mycutg->SetVarY("Br_Ql");
	mycutg->SetPoint(0,600,0.422889); //500
	mycutg->SetPoint(1,5953,0.393);
	mycutg->SetPoint(2,16000,0.393);
	mycutg->SetPoint(3,16000,0.46);
	mycutg->SetPoint(4,2207,0.48);  //500
	mycutg->SetPoint(5,600,0.48);  //500
	mycutg->SetPoint(6,600,0.422889);  //500
	mycutg->Print();
	//cutg->Write();
	tr -> Draw("psd:Br_Ql >> hsep2","mycutg","colz");
	c1->cd();
	//hsep2->Draw();

	cout << "Number of entries: " << tr -> GetEntries() << endl;
	//file->Close();

}