#define StrSelector_cxx
// The class definition in StrSelector.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("StrSelector.C")
// root> T->Process("StrSelector.C","some options")
// root> T->Process("StrSelector.C+")
//


#include "StrSelector.h"
#include <TH2.h>
#include <TStyle.h>
#include <TH1.h>
#include <iostream>


void StrSelector::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).
   htedep = new TH1F("htedep", "total Edep per crystal", 8, 0, 8);

   TString option = GetOption();
}

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

   TString option = GetOption();

}

Bool_t StrSelector::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);
    for(int i = 0,  n = StrSimHits_m_Edep.GetSize(); i < n; ++i){
        cout << "Z position  " << StrSimHits_m_Position[i].z() << endl;
     if (StrSimHits_m_Position[i].Z() >=0){
         htedep->Fill(1, StrSimHits_m_Edep[i]);
         /*if (StrSimHits_m_Position[i].X() > 20 && StrSimHits_m_Position[i].Y()< 10){
             htedep->Fill(0.01, StrSimHits_m_Edep[i]);
         }
         if (StrSimHits_m_Position[i].Y() > 20 && StrSimHits_m_Position[i].X() < 10){
             htedep->Fill(1, StrSimHits_m_Edep[i]);
         }
         if (StrSimHits_m_Position[i].X () < -20 && StrSimHits_m_Position[i].Y()< 10) {
             htedep->Fill(2, StrSimHits_m_Edep[i]);
         }
         if (StrSimHits_m_Position[i].Y() < -20 && StrSimHits_m_Position[i].X() < 10) {
             htedep->Fill(3, StrSimHits_m_Edep[i]);
         }*/
         
     }
        if (StrSimHits_m_Position[i].Z() < 0) {
            htedep->Fill(3, StrSimHits_m_Edep[i]);
         /*if (StrSimHits_m_Position[i].X() > 20){
             htedep->Fill(4, StrSimHits_m_Edep[i]);
         }
         if (StrSimHits_m_Position[i].Y() > 20){
             htedep->Fill(5, StrSimHits_m_Edep[i]);
         }
         if (StrSimHits_m_Position[i].X() < -20) {
             htedep->Fill(6, StrSimHits_m_Edep[i]);
         }
         if (StrSimHits_m_Position[i].Y() < -20) {
             htedep->Fill(7, StrSimHits_m_Edep[i]);
         }*/
        }
    }
//   std::cout << "now at event " << entry << std::endl;
//   std::cout << "Edep is " << *StrSimHits_m_Edep << std::endl;
	return kTRUE;
}

void StrSelector::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 StrSelector::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 = TFile::Open("results.root", "RECREATE");
    f->Delete("Edep;*");
   htedep->Write("Edep");
   f->Close();
}