File indexing completed on 2023-03-17 11:18:21
0001 #ifndef RECOJETS_JETALGORITHMS_SUBJETFILTERALGORITHM_H
0002 #define RECOJETS_JETALGORITHMS_SUBJETFILTERALGORITHM_H 1
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015 #include <vector>
0016
0017 #include "RecoJets/JetAlgorithms/interface/CompoundPseudoJet.h"
0018 #include "FWCore/Framework/interface/Event.h"
0019
0020 #include <fastjet/JetDefinition.hh>
0021 #include <fastjet/AreaDefinition.hh>
0022 #include <fastjet/PseudoJet.hh>
0023
0024 class SubjetFilterAlgorithm {
0025
0026
0027
0028 public:
0029 SubjetFilterAlgorithm(const std::string& moduleLabel,
0030 const std::string& jetAlgorithm,
0031 unsigned nFatMax,
0032 double rParam,
0033 double rFilt,
0034 double jetPtMin,
0035 double massDropCut,
0036 double asymmCut,
0037 bool asymmCutLater,
0038 bool doAreaFastjet,
0039 double ghostEtaMax,
0040 int activeAreaRepeats,
0041 double ghostArea,
0042 bool verbose);
0043 virtual ~SubjetFilterAlgorithm();
0044
0045
0046
0047
0048 public:
0049 void run(const std::vector<fastjet::PseudoJet>& inputs,
0050 std::vector<CompoundPseudoJet>& fatJets,
0051 const edm::EventSetup& iSetup);
0052
0053 std::string summary() const;
0054
0055
0056
0057
0058 private:
0059 std::string moduleLabel_;
0060 std::string jetAlgorithm_;
0061 unsigned nFatMax_;
0062 double rParam_;
0063 double rFilt_;
0064 double jetPtMin_;
0065 double massDropCut_;
0066 double asymmCut2_;
0067 bool asymmCutLater_;
0068 bool doAreaFastjet_;
0069 double ghostEtaMax_;
0070 int activeAreaRepeats_;
0071 double ghostArea_;
0072 bool verbose_;
0073
0074 unsigned nevents_;
0075 unsigned ntotal_;
0076 unsigned nfound_;
0077
0078 fastjet::JetDefinition* fjJetDef_;
0079 fastjet::AreaDefinition* fjAreaDef_;
0080 };
0081
0082 #endif