analysis/EICAnalysis/Quarkonia2LeptonsMC.h

0001 #ifndef __Quarkonia2LeptonsMC_H__
0002 #define __Quarkonia2LeptonsMC_H__
0003
0004 #include <fun4all/SubsysReco.h>
0005 //#include <vector>
0006 //#include <cmath>
0007 //#include <string>
0008 //#include <limits.h>
0009 //#include <TNtuple.h>
0010
0011 class TNtuple;
0012 class TFile;
0013 //class TH2D;
0014 class TH1D;
0015
0016 class PHCompositeNode;
0017 namespace HepMC
0018 {
0019   class GenEvent;
0020 }
0021
0022 class Quarkonia2LeptonsMC : public SubsysReco
0023 {
0024
0025 public:
0026
0027   Quarkonia2LeptonsMC(std::string filename);
0028
0029   int
0030   Init(PHCompositeNode*);
0031   int
0032   process_event(PHCompositeNode*);
0033   int
0034   End(PHCompositeNode*);
0035
0036
0037   //! The embedding ID for the HepMC subevent to be analyzed.
0038   //! embedding ID for the event
0039   //! positive ID is the embedded event of interest, e.g. jetty event from pythia
0040   //! negative IDs are backgrounds, .e.g out of time pile up collisions
0041   //! Usually, ID = 0 means the primary Au+Au collision background
0042   int get_embedding_id() const { return _embedding_id; }
0043   //
0044   //! The embedding ID for the HepMC subevent to be analyzed.
0045   //! embedding ID for the event
0046   //! positive ID is the embedded event of interest, e.g. jetty event from pythia
0047   //! negative IDs are backgrounds, .e.g out of time pile up collisions
0048   //! Usually, ID = 0 means the primary Au+Au collision background
0049   void set_embedding_id(int id) { _embedding_id = id; }
0050 private:
0051
0052   bool _verbose;
0053
0054   int _ievent;
0055   int _total_pass;
0056
0057   std::string _foutname;
0058   TFile *_fout_root;
0059
0060   /* output tree and variables */
0061   TNtuple* _tree_quarkonia;
0062
0063   //! The embedding ID for the HepMC subevent to be analyzed.
0064   //! positive ID is the embedded event of interest, e.g. jetty event from pythia
0065   //! negative IDs are backgrounds, .e.g out of time pile up collisions
0066   //! Usually, ID = 0 means the primary Au+Au collision background
0067   int _embedding_id;
0068 };
0069
0070 #endif // __Quarkonia2LeptonsMC_H__