QAG4Decayer.h

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

#include <fun4all/SubsysReco.h>
#include <g4main/PHG4Particle.h>
#include <g4main/PHG4TruthInfoContainer.h>
#include <g4main/PHG4VtxPoint.h>

#include <g4eval/SvtxEvalStack.h>

#include <globalvertex/SvtxVertex.h>
#include <globalvertex/SvtxVertexMap.h>

#include <trackbase_historic/ActsTransformations.h>
#include <trackbase_historic/SvtxTrack.h>
#include <trackbase_historic/SvtxTrackMap.h>
#include <trackbase_historic/TrackSeed.h>

#include <TBranch.h>
#include <TFile.h>
#include <TH1.h>
#include <TH2.h>
#include <TTree.h>

#include <string>

class PHCompositeNode;
class PHG4TruthInfoContainer;
class PHG4Particle;
class PHG4VtxPoint;

class QAG4Decayer : public SubsysReco
{
 public:
  //            QAG4Decayer();
  QAG4Decayer(const std::string &name = "QAG4Decayer");

  virtual ~QAG4Decayer();

  int Init(PHCompositeNode *topNode) override;
  //            int Init() override;
  int process_event(PHCompositeNode *topNode) override;
  int End(PHCompositeNode *topNode) override;
  std::vector<int> Channel(int pdgid, std::vector<int> Daughter);

 private:
  PHG4TruthInfoContainer *m_truth_info = nullptr;
  // SvtxEvalStack *_svtxevalstack = NULL;

  TFile *fout = nullptr;
  TTree *QATree = nullptr;
  int EvtID = -1;
  float LifeTime = -1;
  /*
    TH1D * MassHis;
    TH1D * NPartHis;
    TH1D * BR1DHis;
   */
  int NParticles = -1;
  std::vector<int> PVDaughtersPDGID;

  bool m_write_nTuple;
  //              bool m_write_QAHists;
  bool m_SaveFiles;
};