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


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

void BKGSelection::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();
}

void BKGSelection::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).
    std::cout<<"slave begin\n";

    h_dts          = new TH1D("dts", "#DeltaT BKG in Protodune.", 10000, 0, 50);
    GetOutputList()->Add(h_dts);
    h_nevt         = new TH1I("nevents", "NEvents in #DeltaT Background", 1, 0, 1);
    GetOutputList()->Add(h_nevt);
    
    std::cout<<"Hists made\n";

    events.clear();

    std::cout<<"events cleared\n";

    TString option = GetOption();

}

Bool_t BKGSelection::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.

    std::cout<<"Getting Entry from reader\n";
    fReader.SetLocalEntry(entry);
    std::cout<<"Got entryfrom reader\n";
        /* const char *current_file_name = ((TChain*)(BKGSelection::fChain))->GetFile()->GetName(); */
        /* std::cout<<"File Name for error: "<<current_file_name<<"\n"; */
    /* std::cout<<"entry: "<<entry<<"\n"; */

    // run, subrun, event
    std::cout<<"1\n";
    events.emplace(*run, *subRun, *event);
    std::cout<<"2\n";
    h_dts->Fill(*deltat);
    std::cout<<"3\n";
    h_nevt->Fill(1, events.size());
    std::cout<<"4\n";
    return kTRUE;
}

void BKGSelection::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 BKGSelection::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.
    std::cout<<"Have Reached Terminate\n";
    TFile* f = new TFile("./PDune_DT_Ana_BKG.root","RECREATE");
    TIter next(GetOutputList());
    while(TH1D* obj = (TH1D*)next())
    {
        obj->Write();
    }
    f->Write();


}