//////////////////////////////////////////////////////////
// This class has been automatically generated on
// Fri May 29 10:04:35 2020 by ROOT version 6.20/00
// from TTree TreeMaster/TreeMaster
// found on file: Track_AGATA_Only_1MEvent_0.5MeVThread00.root
//////////////////////////////////////////////////////////

#ifndef MySelector_h
#define MySelector_h

#include <TROOT.h>
#include <TChain.h>
#include <TFile.h>
#include <TSelector.h>
#include <TTreeReader.h>
#include <TTreeReaderValue.h>
#include <TTreeReaderArray.h>

#include <TH2D.h>
#include <TH3D.h>
#include <TF1.h>
#include <TF2.h>
#include <TF3.h>

#include <TVector3.h>
#include <TFile.h>

// Headers needed by this particular selector


class MySelector : public TSelector {
public :
   TTreeReader     fReader;  //!the tree reader
   TTree          *fChain = 0;   //!pointer to the analyzed TTree or TChain

   TFile *fOutputFile = nullptr;

   TH1D *hCosThetaTheo = nullptr;
   TH1F *h_trackE = nullptr;

   // Readers to access the data (delete the ones you do not need).
   TTreeReaderValue<Int_t> nbHits = {fReader, "nbHits"};
   TTreeReaderArray<Float_t> hitE = {fReader, "hitE"};
   TTreeReaderArray<Float_t> hitX = {fReader, "hitX"};
   TTreeReaderArray<Float_t> hitY = {fReader, "hitY"};
   TTreeReaderArray<Float_t> hitZ = {fReader, "hitZ"};
   TTreeReaderArray<Float_t> hitGX = {fReader, "hitGX"};
   TTreeReaderArray<Float_t> hitGY = {fReader, "hitGY"};
   TTreeReaderArray<Float_t> hitGZ = {fReader, "hitGZ"};
   TTreeReaderArray<Int_t> hitId = {fReader, "hitId"};
   TTreeReaderArray<Int_t> hitSg = {fReader, "hitSg"};
   TTreeReaderValue<Int_t> nbCores = {fReader, "nbCores"};
   TTreeReaderArray<Int_t> nbHitsperCry = {fReader, "nbHitsperCry"};
   TTreeReaderArray<Int_t> coreId = {fReader, "coreId"};
   TTreeReaderArray<Float_t> coreE0 = {fReader, "coreE0"};
   TTreeReaderArray<Float_t> coreE1 = {fReader, "coreE1"};
   TTreeReaderArray<Float_t> coreT0 = {fReader, "coreT0"};
   TTreeReaderArray<ULong64_t> coreTS = {fReader, "coreTS"};
   TTreeReaderValue<Int_t> nbAdd = {fReader, "nbAdd"};
   TTreeReaderArray<Int_t> AddId = {fReader, "AddId"};
   TTreeReaderArray<Float_t> AddE = {fReader, "AddE"};
   TTreeReaderArray<Float_t> AddX = {fReader, "AddX"};
   TTreeReaderArray<Float_t> AddY = {fReader, "AddY"};
   TTreeReaderArray<Float_t> AddZ = {fReader, "AddZ"};
   TTreeReaderArray<ULong64_t> AddTS = {fReader, "AddTS"};
   TTreeReaderValue<ULong64_t> TSHit = {fReader, "TSHit"};
   TTreeReaderValue<Int_t> nbTrack = {fReader, "nbTrack"};
   TTreeReaderArray<Float_t> trackE = {fReader, "trackE"};
   TTreeReaderArray<Float_t> trackX1 = {fReader, "trackX1"};
   TTreeReaderArray<Float_t> trackY1 = {fReader, "trackY1"};
   TTreeReaderArray<Float_t> trackZ1 = {fReader, "trackZ1"};
   TTreeReaderArray<Float_t> trackX2 = {fReader, "trackX2"};
   TTreeReaderArray<Float_t> trackY2 = {fReader, "trackY2"};
   TTreeReaderArray<Float_t> trackZ2 = {fReader, "trackZ2"};
   TTreeReaderArray<Int_t> trackType = {fReader, "trackType"};
   TTreeReaderArray<Float_t> trackFOM = {fReader, "trackFOM"};
   TTreeReaderArray<Float_t> trackT = {fReader, "trackT"};
   TTreeReaderArray<Int_t> trackCrystalID = {fReader, "trackCrystalID"};
   TTreeReaderValue<ULong64_t> TStrack = {fReader, "TStrack"};

   static const int MaxSize = 50;
   std::array <Float_t, MaxSize> trackE1;
   std::array <Float_t, MaxSize> trackE2;

   MySelector(TTree * /*tree*/ =0) { }
   virtual ~MySelector() { }
   virtual Int_t   Version() const { return 2; }
   virtual void    Begin(TTree *tree);
   virtual void    SlaveBegin(TTree *tree);
   virtual void    Init(TTree *tree);
   virtual Bool_t  Notify();
   virtual Bool_t  Process(Long64_t entry);
   virtual Int_t   GetEntry(Long64_t entry, Int_t getall = 0) { return fChain ? fChain->GetTree()->GetEntry(entry, getall) : 0; }
   virtual void    SetOption(const char *option) { fOption = option; }
   virtual void    SetObject(TObject *obj) { fObject = obj; }
   virtual void    SetInputList(TList *input) { fInput = input; }
   virtual TList  *GetOutputList() const { return fOutput; }
   virtual void    SlaveTerminate();
   virtual void    Terminate();

   Double_t CalcCosTheta(Double_t X, Double_t Y, Double_t Z, TVector3 &T1, TVector3 &Int2);

   ClassDef(MySelector,0);

};

#endif

#ifdef MySelector_cxx
void MySelector::Init(TTree *tree)
{
   // The Init() function is called when the selector needs to initialize
   // a new tree or chain. Typically here the reader is initialized.
   // It is normally not necessary to make changes to the generated
   // code, but the routine can be extended by the user if needed.
   // Init() will be called many times when running on PROOF
   // (once per file to be processed).

   fReader.SetTree(tree);
}

Bool_t MySelector::Notify()
{
   // The Notify() function is called when a new file is opened. This
   // can be either for a new TTree in a TChain or when when a new TTree
   // is started when using PROOF. It is normally not necessary to make changes
   // to the generated code, but the routine can be extended by the
   // user if needed. The return value is currently not used.

   return kTRUE;
}


#endif // #ifdef MySelector_cxx