15 #include "fastjet/PseudoJet.hh" 16 #include "fastjet/JetDefinition.hh" 17 #include "fastjet/ClusterSequence.hh" 18 #include "fastjet/Selector.hh" 37 : tokenLHEEvent_(consumes<
LHEEventProduct>(params.getParameter<
edm::InputTag>(
"src"))),
38 jetPtMin_(params.getParameter<double>(
"jetPtMin")),
39 jetdef_(antikt_algorithm, params.getParameter<double>(
"jetR")) {}
49 vector<PseudoJet> jetconsts;
50 jetconsts.reserve(10);
52 for (
size_t p = 0;
p < hepeup.
IDUP.size(); ++
p) {
53 if (hepeup.
ISTUP[
p] == 1) {
55 jetconsts.emplace_back(mom[0], mom[1], mom[2], mom[3]);
const lhef::HEPEUP & hepeup() const
bool getByToken(EDGetToken token, Handle< PROD > &result) const
unique_ptr< ClusterSequence > cs
edm::EDGetTokenT< LHEEventProduct > tokenLHEEvent_
#define DEFINE_FWK_MODULE(type)
LHEJetFilter(const edm::ParameterSet &)
std::vector< FiveVector > PUP
bool filter(edm::StreamID strid, edm::Event &evt, const edm::EventSetup ¶ms) const override