isoCluster.h

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file isoCluster.h

// Tell emacs that this is a C++ source
//  -*- C++ -*-.
#ifndef ISOCLUSTER_H
#define ISOCLUSTER_H

#include <fun4all/SubsysReco.h>

#include <string>

#include <vector>

class PHCompositeNode;
class TFile;
class TTree;

class isoCluster : public SubsysReco
{
 public:

  isoCluster(const std::string &name = "isoCluster");

  ~isoCluster() override;

  int Init(PHCompositeNode *topNode) override;

  int InitRun(PHCompositeNode *topNode) override;

  int process_event(PHCompositeNode *topNode) override;

  int ResetEvent(PHCompositeNode *topNode) override;

  int EndRun(const int runnumber) override;

  int End(PHCompositeNode *topNode) override;

  int Reset(PHCompositeNode * /*topNode*/) override;

  void Print(const std::string &what = "ALL") const override;

  void setGenEvent(int eventGet)     {getEvent = eventGet;}

 private:

  TTree *T;
  //cluster information
  std::vector<float> m_clusterPhi;
  std::vector<float> m_clusterEta;
  std::vector<float> m_clusterE;
  std::vector<float> m_clusterPt;
  std::vector<float> m_clusterEtIso;
  std::vector<float> m_clusterChisq;
  
  //truth particle information
  std::vector<float> m_photonPhi;
  std::vector<float> m_photonEta;
  std::vector<float> m_photonE;
  std::vector<float> m_photonPt;

  TFile *fout;
  std::string outname;
  int getEvent;
};

#endif // ISOCLUSTER_H