QAG4SimulationMicromegas.h

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

#ifndef QA_QAG4SIMULATIONMICROMEGAS_H
#define QA_QAG4SIMULATIONMICROMEGAS_H

#include <trackbase/ClusterErrorPara.h>
#include <trackbase/TrkrDefs.h>

#include <fun4all/SubsysReco.h>

#include <set>
#include <string>

class ActsGeometry;
class PHCompositeNode;
class PHG4CylinderGeomContainer;
class PHG4Hit;
class PHG4Particle;
class PHG4HitContainer;
class TrkrClusterContainer;
class TrkrHitSetContainer;
class TrkrClusterHitAssoc;
class TrkrHitTruthAssoc;
class PHG4TruthInfoContainer;
class SvtxEvalStack;

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

  int InitRun(PHCompositeNode* topNode) override;
  int process_event(PHCompositeNode* topNode) override;
 
 private:
  std::string get_histo_prefix() const;
  std::unique_ptr<SvtxEvalStack> m_svtxEvalStack;
  int load_nodes(PHCompositeNode*);

  void evaluate_hits();

  void evaluate_clusters();

  // get geant hits associated to a cluster
  using G4HitSet = std::set<PHG4Hit*>;
  G4HitSet find_g4hits(TrkrDefs::cluskey) const;

  bool m_initialized = false;

  PHG4CylinderGeomContainer* m_micromegas_geonode = nullptr;

  ActsGeometry* m_tGeometry = nullptr;

  TrkrClusterContainer* m_cluster_map = nullptr;

  TrkrHitSetContainer* m_hitsets = nullptr;

  TrkrClusterHitAssoc* m_cluster_hit_map = nullptr;

  TrkrHitTruthAssoc* m_hit_truth_map = nullptr;

  PHG4HitContainer* m_g4hits_micromegas = nullptr;

  PHG4TruthInfoContainer* m_truthContainer;
  /* it is filled at Init stage. It should not change for the full run */
  std::set<int> m_layers;
  ClusterErrorPara _ClusErrPara;
};

#endif