30 if (m_jetMinDR < 0.1) m_jetMinDR = 0.1;
40 desc.
add<
double>(
"JetMinPt",25.0);
41 desc.
add<
double>(
"JetMaxEta",2.4);
42 desc.
add<
double>(
"JetMinDR",0.2);
43 desc.
add<
double>(
"JetMaxDZ",0.2);
44 desc.
add<
int>(
"TriggerType",84);
58 typedef vector<T> TCollection;
72 const TCollection &
jets = *jetsHandle;
73 const size_t n_jets = jets.size();
77 size_t npairs = 0, nfail_dz = 0;
78 if (n_jets > 1)
for(
size_t j1 = 0; j1 < n_jets; ++j1)
83 for(
size_t j2 = j1+1; j2 < n_jets; ++j2)
87 float deta = jets[j1].eta() - jets[j2].eta();
89 float dr2 = dphi*dphi + deta*deta;
90 float dz = jets[j1].vz() - jets[j2].vz();
106 ref = TRef(jetsHandle, j1);
109 ref = TRef(jetsHandle, j2);
constexpr double deltaPhi(double phi1, double phi2)
T getParameter(std::string const &) const
bool getByToken(EDGetToken token, Handle< PROD > &result) const
std::string defaultModuleLabel()
void addObject(int id, const reco::RecoEcalCandidateRef &ref)
setters for L3 collections: (id=physics type, and Ref<C>)
Abs< T >::type abs(const T &t)
bool hltFilter(edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override
ParameterDescriptionBase * add(U const &iLabel, T const &value)
static void makeHLTFilterDescription(edm::ParameterSetDescription &desc)
void addCollectionTag(const edm::InputTag &collectionTag)
collectionTags
void add(std::string const &label, ParameterSetDescription const &psetDescription)
HLTJetPairDzMatchFilter(const edm::ParameterSet &)
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
edm::EDGetTokenT< std::vector< T > > m_jetToken
~HLTJetPairDzMatchFilter() override