SQGeomAcc.h

Go to the documentation of this file.

#ifndef _SQ_GEOM_ACC__H_
#define _SQ_GEOM_ACC__H_
#include <map>
#include <fun4all/SubsysReco.h>
class SQHitVector;
 
 
class SQGeomAcc: public SubsysReco {
 public:
  typedef enum {
    UNDEF_MUON, 
    SINGLE,    //< Require one muon at either top or bottom half
    SINGLE_T,  //< Require one muon at top half
    SINGLE_B,  //< Require one muon at bottom half
    PAIR,      //< Require two muons at any halves
    PAIR_TBBT, //< Require two muons at T+B or B+T halves
    PAIR_TTBB, //< Require two muons at T+T or B+B halves
  } MuonMode_t;
  typedef enum {
    UNDEF_PLANE, 
    HODO,      //< Use the hodoscope acceptance
    CHAM,      //< Use the chamber acceptance
    HODO_CHAM, //< Use the hodoscope+chamber acceptance
  } PlaneMode_t;
 
 private:
  MuonMode_t   m_mode_muon;
  PlaneMode_t  m_mode_plane;
  unsigned int m_n_ele_h1_ex; //< N of H1T/B elements excluded
  SQHitVector* m_vec_hit;
 
 public:
  SQGeomAcc(const std::string& name = "SQGeomAcc");
  virtual ~SQGeomAcc();
  int Init(PHCompositeNode *topNode);
  int InitRun(PHCompositeNode *topNode);
  int process_event(PHCompositeNode *topNode);
  int End(PHCompositeNode *topNode);
 
  void SetMuonMode (MuonMode_t  mode) { m_mode_muon  = mode; }
  void SetPlaneMode(PlaneMode_t mode) { m_mode_plane = mode; }
  void SetNumOfH1EdgeElementsExcluded(const unsigned int num) { m_n_ele_h1_ex = num; }
 
 private:
  int  GetDetId(const std::string& det_name);
  bool FindDetIdSet(const std::vector<int>& vec_det_id_all, const std::vector<int>& vec_det_id_want);
};
 
#endif /* _SQ_GEOM_ACC__H_ */