FastJetMedianBkg/src_JetMedianTree/PrintTowers.h

0001 #ifndef PRINTTOWERS___H
0002 #define PRINTTOWERS___H
0003
0004 #include <fun4all/SubsysReco.h>
0005 #include <jetbase/Jet.h>
0006
0007 #include <memory>
0008 #include <string>
0009 #include <utility>  // std::pair, std::make_pair
0010
0011 #include <array>
0012 #include <vector>
0013
0014 class Jet;
0015 class PHCompositeNode;
0016 class JetEvalStack;
0017 class TTree;
0018 class TH2D;
0019 class TH1;
0020 class JetInput;
0021
0022 class PrintTowers : public SubsysReco
0023 {
0024  public:
0025   PrintTowers
0026     (
0027         std::vector<std::string> jet_names
0028       , std::vector<std::pair<std::string,std::string>> matchlist = {}
0029       , std::vector<std::pair<Jet::SRC,std::string>> tower_sources = {}
0030       , const float min_jet_pt = 5.
0031       , const float m_tow_tow_pt = 0.01
0032     /* , const float        _jet_R          = 0.4 */
0033     /* , const std::string& truthjetname   = "AntiKt_Truth_r04" */
0034     /* , const std::string& sub1jetname   = "AntiKt_Truth_r04" */
0035     /* , const float min_lead_truth_pt = 10. */
0036     );
0037
0038   virtual ~PrintTowers();
0039   float m_max_cent = 100.;
0040   int nmax_jetprint = 10.;
0041   float m_pt_min_match = 30.;
0042
0043   /* void use_initial_vertex(const bool b = true) {initial_vertex = b;} */
0044   int Init          (PHCompositeNode *topNode);
0045   int InitRun       (PHCompositeNode *topNode);
0046   int process_event (PHCompositeNode *topNode);
0047   int End           (PHCompositeNode *topNode);
0048   void clear_vectors();
0049   /* void add_input(JetInput *input) { m_inputs.push_back(input); } */
0050
0051  private:
0052
0053   const std::vector<std::string> m_jet_names;
0054   const std::vector<std::pair<std::string,std::string>> m_matchlist;
0055   const std::vector<std::pair<Jet::SRC,std::string>> m_SRC;
0056   const float m_min_jet_pt;
0057   const float m_min_tow_pt;
0058   int nevnt {0};
0059
0060 };
0061
0062 #endif  // JETRHOMEDIAN__H