dd/ddf/EICG4dRICHTree_8h_source.html

//-----------------------------------------------

// EICG4dRICH event tree module

//-----------------------------------------------


#ifndef DRICHTREE_H

#define DRICHTREE_H


#include <Rtypes.h>

#include <TString.h>

#include <fun4all/SubsysReco.h>

#include <Geant4/G4String.hh>

#include <Geant4/G4ThreeVector.hh>


// class Fun4AllHistoManager; //---

class PHCompositeNode;

class TFile;

class TTree;

// class TH1; //---

// class TH2; //---


class EICG4dRICHTree : public SubsysReco

{

 public:

  // constructor

  EICG4dRICHTree(const std::string &name = "EICG4dRICHTree",

                 const std::string &fname = "EICG4dRICHTree.root");


  // destructor

  virtual ~EICG4dRICHTree();


  // SubsysReco processing methods:

  int Init(PHCompositeNode *);

  int process_event(PHCompositeNode *);

  int End(PHCompositeNode *);


 private:

  std::string m_outfileN;  // output file name

  TFile *m_outfile;        // output file

  TTree *m_tree;           // the EICG4dRICH tree


  // TH1 *m_histo; //---

  // Fun4AllHistoManager *m_hm; // histogram manager //---


  // data accessors

  void getHits(PHCompositeNode *topNode);

  // void getHEPMCTruth(PHCompositeNode *topNode); //---


  void vectorToArray(G4ThreeVector vec, Double_t *arr);

  void initTrees();

  // void resetVars(); //---


  //------------------------------

  // tree variables

  //------------------------------


  Int_t evnum;

  Int_t trackID;

  char hitType[128];

  char hitSubtype[128];

  Int_t petal, psst, pdg;

  char particleName[128];

  char process[256];

  Int_t parentID;

  Double_t hitPos[3];  // [xyz]

  Double_t hitP[3];

  Double_t hitPdir[3];

  Double_t hitVtxPos[3];

  Double_t hitVtxPdir[3];

  Double_t deltaT;

  Double_t edep;

};


#endif  // DRICHTREE_H