HCal-analysis/HCALAnalysis/HCALAnalysis.h

0001 #ifndef __HCALANALYSIS_H__
0002 #define __HCALANALYSIS_H__
0003
0004 #include <fun4all/SubsysReco.h>
0005 #include <string>
0006
0007 class PHCompositeNode;
0008 class PHG4HitContainer;
0009 class TFile;
0010 class TNtuple;
0011
0012 class HCALAnalysis: public SubsysReco
0013 {
0014  public:
0015   HCALAnalysis(const char *outfile = "ntuple.root");
0016   virtual ~HCALAnalysis() {}
0017
0018   //  For this analysis we only use Init, process_event;
0019   int Init         (PHCompositeNode *topNode);
0020   int process_event(PHCompositeNode *topNode);
0021   int End          (PHCompositeNode *topNode);
0022
0023  protected:
0024
0025   const char *OutFileName;
0026
0027   void GetNodes(PHCompositeNode *topNode);
0028
0029   PHCompositeNode* _topNode;
0030   PHG4HitContainer* _hcalout_hit_container;
0031   PHG4HitContainer* _hcalin_hit_container;
0032   PHG4HitContainer* _cemc_hit_container;
0033   PHG4HitContainer* _hcalout_abs_hit_container;
0034   PHG4HitContainer* _hcalin_abs_hit_container;
0035   PHG4HitContainer* _cemc_abs_hit_container;
0036   PHG4HitContainer* _magnet_hit_container;
0037   PHG4HitContainer* _bh_hit_container;
0038   PHG4HitContainer* _bh_plus_hit_container;
0039   PHG4HitContainer* _bh_minus_hit_container;
0040   PHG4HitContainer* _cemc_electronics_hit_container;
0041   PHG4HitContainer* _hcalin_spt_hit_container;
0042   PHG4HitContainer* _svtx_hit_container;
0043   PHG4HitContainer* _svtx_support_hit_container;
0044
0045   TFile *outputfile;
0046   TNtuple *calenergy;
0047
0048 };
0049
0050 #endif /* __HCALANALYSIS_H__ */