packages/MvtxPrototype1/MvtxQAHisto.h

0001 #ifndef __MvtxQAHisto_H__
0002 #define __MvtxQAHisto_H__
0003
0004 #include <fun4all/SubsysReco.h>
0005 #include <string>
0006
0007 //Forward declerations
0008 class PHCompositeNode;
0009 class TrkrClusterContainer;
0010 class TrkrHitSetContainer;
0011 class TFile;
0012 class TH1F;
0013 class TH2F;
0014
0015 //Brief: basic ntuple and histogram creation for sim evaluation
0016 class MvtxQAHisto: public SubsysReco
0017 {
0018  public:
0019   //Default constructor
0020   MvtxQAHisto(const std::string &name="MvtxQAHisto");
0021
0022   //Initialization, called for initialization
0023   int Init(PHCompositeNode *);
0024
0025   //Process Event, called for each event
0026   int process_event(PHCompositeNode *);
0027
0028   //End, write and close files
0029   int End(PHCompositeNode *);
0030
0031   //Change output filename
0032   void set_filename(const char* file)
0033   { if(file) _outfile = file; }
0034
0035  private:
0036   //output filename
0037   std::string _outfile;
0038
0039   //Event counter
0040   int _event;
0041
0042   //Get all the nodes
0043   void GetNodes(PHCompositeNode *);
0044
0045   //flags
0046   unsigned int _flags;
0047
0048     float beam_x[4];
0049     float beam_y[4];
0050
0051     TH2F *h2d_hit[4];
0052     TH2F *h2d_hit_beam[4];
0053     TH2F *h2d_hit_trk[4];
0054     TH1F *h1d_hit_per_evt[4];
0055
0056     TH2F *h2d_clus[4];
0057     TH2F *h2d_clus_beam[4];
0058     TH1F *h1d_clus_per_evt[4]; //cluster in a chip per event
0059     TH1F *h1d_clus_size_x[4];
0060     TH1F *h1d_clus_size_z[4];
0061
0062     TrkrClusterContainer *cluscon;
0063     TrkrHitSetContainer *hitsetcon;
0064
0065     TH1F *h1d_trk_finder_x;
0066     TH1F *h1d_trk_finder_z;
0067     TH2F *h2d_trk_finder;
0068
0069     TH1F *h1d_clus_associated;
0070     TH1F *h1d_clus_eff;
0071
0072
0073 };
0074
0075 #endif //* __MvtxQAHisto_H__ *//