singlePhotonClusterAna.h

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

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

#include <fun4all/SubsysReco.h>

#include <string>
#include <vector>

class PHCompositeNode;
class TTree;
class PHG4Particle;
class CaloEvalStack;
class TFile;
class Fun4AllHistoManager;

class singlePhotonClusterAna : public SubsysReco
{
 public:

  singlePhotonClusterAna(const std::string &name, const std::string &outName);

  ~singlePhotonClusterAna() 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;

 private:

  float getEta(PHG4Particle *particle);
    
  TTree *clusters_Towers;
  TTree *truth_photon;
  TTree *conversion;
  
  //CaloEvalStack *caloevalstack;
  
  //stuff for towers and clusters
  std::vector<float> m_eta_center;
  std::vector<float> m_phi_center;
  std::vector<float> m_tower_energy;
  std::vector<float> m_cluster_eta;
  std::vector<float> m_cluster_phi;
  std::vector<float> m_cluster_e;
  std::vector<float> m_cluster_ecore;
  std::vector<float> m_cluster_chi2;
  std::vector<float> m_cluster_prob;
  std::vector<float> m_cluster_nTowers;
  
  //stuff for truth photons
  std::vector<float> m_photon_E;
  std::vector<float> m_photon_eta;
  std::vector<float> m_photon_phi;

  //stuff for conversions
  std::vector<float> m_electron_E;
  std::vector<float> m_electron_eta;
  std::vector<float> m_electron_phi;
  std::vector<float> m_positron_E;
  std::vector<float> m_positron_eta;
  std::vector<float> m_positron_phi;
  std::vector<float> m_vtx_x;
  std::vector<float> m_vtx_y;
  std::vector<float> m_vtx_z;
  std::vector<float> m_vtx_t;
  std::vector<float> m_vtx_r;
  std::vector<bool> m_isConversionEvent;
  
  TFile *out;
  Fun4AllHistoManager *hm = nullptr;
  std::string Outfile;

  
};

#endif // SINGLEPHOTONCLUSTERANA_H