|
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
38
39
40
41
42
43
44
45
46
47
48
|
#ifndef MatchJet_H
#define MatchJet_H
#include <vector>
#include <string>
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/Common/interface/RefToBaseVector.h"
#include "DataFormats/JetReco/interface/Jet.h"
#include "DQMOffline/RecoB/interface/CorrectJet.h"
/** \class MatchJet
*
* Match jets
*
*/
class MatchJet {
public:
MatchJet() {}
MatchJet(const edm::ParameterSet& pSet);
void setThreshold(const double& energy) { threshold = energy; }
/// match the collections
void matchCollections(const edm::RefToBaseVector<reco::Jet>& refJets,
const edm::RefToBaseVector<reco::Jet>& recJets,
const reco::JetCorrector* corrector);
/// Returns the matched "reference" jet
edm::RefToBase<reco::Jet> operator()(const edm::RefToBase<reco::Jet>& recJet) const;
private:
std::vector<int> refToRec, recToRef;
edm::RefToBaseVector<reco::Jet> refJets;
edm::RefToBaseVector<reco::Jet> recJets;
CorrectJet refJetCorrector;
CorrectJet recJetCorrector;
double maxChi2;
double sigmaDeltaR2;
double sigmaDeltaE2;
double threshold;
};
#endif
|