|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
|
#ifndef MatchLessByDPt_h_
#define MatchLessByDPt_h_
/** Provides the less operator for two pairs of matched objects based
* on deltaPt. */
#include "FWCore/ParameterSet/interface/ParameterSet.h"
// #include "CommonTools/UtilAlgos/interface/DeltaR.h"
namespace reco {
template <typename C1, typename C2>
struct MatchLessByDPt {
public:
MatchLessByDPt(const edm::ParameterSet& cfg, const C1& c1, const C2& c2)
: // deltaR_(cfg),
c1_(c1),
c2_(c2) {}
bool operator()(const std::pair<size_t, size_t>& p1, const std::pair<size_t, size_t>& p2) const {
typedef typename C1::value_type T1;
typedef typename C2::value_type T2;
const T1& p1_1 = c1_[p1.first];
const T2& p1_2 = c2_[p1.second];
const T1& p2_1 = c1_[p2.first];
const T2& p2_2 = c2_[p2.second];
if (fabs(p1_1.pt() - p1_2.pt()) / p1_2.pt() < fabs(p2_1.pt() - p2_2.pt()) / p2_2.pt())
return true;
return false;
}
private:
// DeltaR deltaR_;
const C1& c1_;
const C2& c2_;
};
} // namespace reco
#endif
|