analysis.h

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

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

#include <fun4all/SubsysReco.h>

#include <string>
#include <TF1.h>
#include<string.h>
#include<TProfile.h>
#include<TH1.h>
#include<TH2.h>
#include <TFile.h>
#include <TTree.h>


class Fun4AllHistoManager;

class PHCompositeNode;
class TF1;
class TH1D;
class TH2D;
class TTree;

class analysis : public SubsysReco
{
 public:

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

  virtual ~analysis();

  void SetOutputFileName(std::string outfilename)
  {
    m_outfilename = outfilename;
  }

  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:
  //---------------------------------------------
  //define the histogram manager
  //---------------------------------------------
  Fun4AllHistoManager *
  get_HistoManager();

  //-------------------------------------------------------------------------------
  //define some vectors in which we store jet information
  //-------------------------------------------------------------------------------

  std::vector<float> m_tjet_pt;
  std::vector<float> m_tjet_phi;
  std::vector<float> m_tjet_eta;
  std::vector<float> m_rjet_pt;
  std::vector<float> m_rjet_phi;
  std::vector<float> m_rjet_eta;
  std::vector<float> m_dr;
  int m_centrality; 
  int m_impactparam;
  std::string m_outfilename;

  float m_tjet_pt2;


  
  TFile* outfile;

};

#endif // ANALYSIS_H