SBSSimEvent.h

#ifndef __SBSSimEvent_h
#define __SBSSimEvent_h

//#include "gmn_dig_tree.h"
//#include "g4sbs_tree.h"

#include "gmn_tree_digitized.h"
#include "genrp_tree_digitized.h"

enum Exp_t    { kGEp, kGEnRP, kGMN, kSIDIS};

class TTree;

class SBSSimEvent {
 public:
  //SBSSimEvent(){};                 // Default constructor, for ROOT I/O
  //Now we want to initialize the ROOT tree the same way 
  //SBSSimEvent(TTree* tree, TString experiment="gmn");//, std::vector<TString> det_list);
  SBSSimEvent(TTree* tree, Exp_t experiment=kGMN);//, std::vector<TString> det_list);
  
  virtual ~SBSSimEvent(){};
  
  virtual Int_t    GetEntry(Long64_t entry);
  virtual void Clear( const Option_t* opt="" );
  virtual void Print( const Option_t* opt="" ) const;

  ULong64_t RunID;
  ULong64_t EvtID;

  //AJRP: redesign for simplicity/portability/maintainability:

  //We probably don't actually need these getter and setter functions, but it doesn't hurt to have them.
  //void SetExperiment( TString expname ){ fExperiment = expname; }
  //TString GetExperiment() const { return fExperiment; }
  
  //TString fExperiment; //string that tells us which simulated experiment we are doing
  
  void SetExperiment( Exp_t exp ){ fExperiment = exp; }
  Exp_t GetExperiment() const { return fExperiment; }
  
  Exp_t fExperiment;
  
  //Auto-generated ROOT Tree classes for each experiment ROOT tree; generated using TTree::MakeClass()
  
  //Later on, any time we want to analyze a g4sbs root file whose format has changed, we can just run TTree::MakeClass on that root file with the appropriate
  //class name, and copy the source and header files into SBS-offline,
  //recompile, and voila: compatibility guaranteed:
  gmn_tree_digitized *Tgmn;
  //gep_tree_digitized *Tgep;
  genrp_tree_digitized *Tgenrp;//EPAF: for now, genrp tree is thought as a complement of the GMN tree. 
  // we might keep it this way unless it induces crashes or significant slowdown!
  //sidis_tree_digitized *Tsidis;
  //  gen_tree_digitized *Tgen; //This actually seems like it wouldn't require anything different from gmn.

  
  
  ClassDef(SBSSimEvent, 1) // Simulated data for one event
};

#endif