#define XSelector_cxx
// The class definition in XSelector.h has been generated automatically
// by the ROOT utility TTree::MakeSelector(). This class is derived
// from the ROOT class TSelector. For more information on the TSelector
// framework see $ROOTSYS/README/README.SELECTOR or the ROOT User Manual.


// The following methods are defined in this file:
//    Begin():        called every time a loop on the tree starts,
//                    a convenient place to create your histograms.
//    SlaveBegin():   called after Begin(), when on PROOF called only on the
//                    slave servers.
//    Process():      called for each event, in this function you decide what
//                    to read and fill your histograms.
//    SlaveTerminate: called at the end of the loop on the tree, when on PROOF
//                    called only on the slave servers.
//    Terminate():    called at the end of the loop on the tree,
//                    a convenient place to draw/fit your histograms.
//
// To use this file, try the following session on your Tree T:
//
// root> T->Process("XSelector.C")
// root> T->Process("XSelector.C","some options")
// root> T->Process("XSelector.C+")
//


#include "XSelector.h"
#include <TH2.h>
#include <TStyle.h>

void XSelector::Begin(TTree *fChain)
{
	// The Begin() function is called at the start of the query.
	// When running with PROOF Begin() is only called on the client.
	// The tree argument is deprecated (on PROOF 0 is passed).

	TString option = GetOption();
}

void XSelector::SlaveBegin(TTree * /*tree*/)
{
	// The SlaveBegin() function is called after the Begin() function.
	// When running with PROOF SlaveBegin() is called on each slave server.
	// The tree argument is deprecated (on PROOF 0 is passed).

	hPID = new TH2D("hPID", "PID", 1000, 0, 30000, 1000, 0, 30000); 
	GetOutputList()->Add(hPID); 
	hSega5 = new TH1D("hSega5", "Sega 5", 1500, 0, 3000); 
	GetOutputList()->Add(hSega5); 

	TString option = GetOption();
}

Bool_t XSelector::Process(Long64_t entry)
{
	// The Process() function is called for each entry in the tree (or possibly
	// keyed object in the case of PROOF) to be processed. The entry argument
	// specifies which entry in the currently loaded tree is to be processed.
	// When processing keyed objects with PROOF, the object is already loaded
	// and is available via the fObject pointer.
	//
	// This function should contain the \"body\" of the analysis. It can contain
	// simple or elaborate selection criteria, run algorithms on the data
	// of the event and typically fill histograms.
	//
	// The processing can be stopped by calling Abort().
	//
	// Use fStatus to set the return value of TTree::Process().
	//
	// The return value is currently not used.

	fReader.SetEntry(entry);

	//	if(blob53->IsInside(*tac_i2scorr, *pin01_ecal))
	//	if((bool)!blob53)
	//	if(!canplot)
		hPID->Fill(*tac_i2scorr, *pin01_ecal); 

	if(sega_ecal[5] > 0) {
		hSega5->Fill(sega_ecal[5]); 
	}
	return kTRUE;
}

void XSelector::SlaveTerminate()
{
	// The SlaveTerminate() function is called after all entries or objects
	// have been processed. When running with PROOF SlaveTerminate() is called
	// on each slave server.

}

void XSelector::Terminate()
{
	// The Terminate() function is the last function to be called during
	// a query. It always runs on the client, it can be used to present
	// the results graphically or save the results to file.

	//	hPID->Draw("colz");
	TFile *f = new TFile("temp.root", "RECREATE");
	fOutput->FindObject("hPID")->Write(); 	
	fOutput->FindObject("hSega5")->Write(); 

	f->Close(); 


}