#define tsel_cxx
// The class definition in tsel.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 everytime 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("tsel.C")
// Root > T->Process("tsel.C","some options")
// Root > T->Process("tsel.C+")
//

#include "tsel.h"
#include <TH2.h>
#include <TStyle.h>
#include <iostream>

Int_t nentries=0;
Float_t ratio=0;

const Int_t maxmul=8;
Int_t mult, entries;
Float_t Ge[maxmul];

TTree *ggMat = new TTree("ggMat","g's in coincidence");
TBranch *b = ggMat->Branch("mult",&mult,"mult/I"); TBranch *b1 = ggMat->Branch("Ge",Ge,"Ge[mult]/F"); TBranch *b2 = ggMat->Branch("entries",&entries,"entries/I");

void tsel::Begin(TTree *tree)
{
   // 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();
	
   	TChain *t = (TChain*)gROOT->FindObject("t");
	nentries = Int_t(t->GetEntries()); //I should now have global access to nentries!
	printf("total events to be processed: %d \n \n",nentries);

}

void tsel::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).

   Init(tree);

   TString option = GetOption();

}

Bool_t tsel::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.
   // It can be passed to either TTree::GetEntry() or TBranch::GetEntry()
   // to read either all or the required parts of the data. 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.

   // WARNING when a selector is used with a TChain, you must use
   //  the pointer to the current TTree to call GetEntry(entry).
   //  The entry is always the local entry number in the current tree.
   //  Assuming that fChain is the pointer to the TChain being processed,
   //  use fChain->GetTree()->GetEntry(entry).

	ratio = (Float_t)100*entry/nentries;
	if ((entry%(nentries/10)==0)&&(entry>0)) {
		printf("Processed Events: %d %% (%d events out of %d) \n",int(TMath::Ceil(ratio)), (int)entry, nentries);
	}

	fChain->GetTree()->GetEntry(entry);
	
	if (mul>1) {
		entries=entry;
		mult=mul;
		for (Int_t j=0; j<mult; j++) Ge[j]=Eg[j];
		ggMat->Fill();
		}

   return kTRUE;
}

void tsel::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 tsel::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.
   TFile *f = new TFile("subtreehouse.root","recreate");
   ggMat->Write();
   f->Close();

}