|
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
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
|
#include "DQMOffline/PFTau/plugins/PFJetDQMAnalyzer.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "DataFormats/JetReco/interface/JetCollection.h"
#include "DataFormats/JetReco/interface/PFJet.h"
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidateFwd.h"
#include "DQMServices/Core/interface/DQMStore.h"
//
// -- Constructor
//
PFJetDQMAnalyzer::PFJetDQMAnalyzer(const edm::ParameterSet ¶meterSet)
{
pSet_ = parameterSet;
inputLabel_ = pSet_.getParameter<edm::InputTag>("InputCollection");
matchLabel_ = pSet_.getParameter<edm::InputTag>("MatchCollection");
benchmarkLabel_ = pSet_.getParameter<std::string>("BenchmarkLabel");
pfJetMonitor_.setParameters(parameterSet); // set parameters for booking histograms and validating jet
myJet_ = consumes<edm::View<reco::Jet>>(inputLabel_);
myMatchedJet_ = consumes<edm::View<reco::Jet>>(matchLabel_);
std::string folder = benchmarkLabel_;
subsystemname_ = "ParticleFlow";
eventInfoFolder_ = subsystemname_ + "/" + folder;
nBadEvents_ = 0;
}
//
// -- BookHistograms
//
void PFJetDQMAnalyzer::bookHistograms(DQMStore::IBooker &ibooker,
edm::Run const & /* iRun */,
edm::EventSetup const & /* iSetup */) {
ibooker.setCurrentFolder(eventInfoFolder_);
edm::LogInfo("PFJetDQMAnalyzer") << " PFJetDQMAnalyzer::bookHistograms "
<< "Histogram Folder path set to " << eventInfoFolder_;
pfJetMonitor_.setup(ibooker, pSet_);
}
//
// -- Analyze
//
void PFJetDQMAnalyzer::analyze(edm::Event const &iEvent, edm::EventSetup const &iSetup) {
edm::Handle<edm::View<reco::Jet>> jetCollection;
iEvent.getByToken(myJet_, jetCollection);
edm::Handle<edm::View<reco::Jet>> matchedJetCollection;
iEvent.getByToken(myMatchedJet_, matchedJetCollection);
float maxRes = 0.0;
float minRes = 99.99;
float jetpT = 0.0;
if (jetCollection.isValid() && matchedJetCollection.isValid()) {
pfJetMonitor_.fill(*jetCollection,
*matchedJetCollection,
minRes,
maxRes,
jetpT,
pSet_); // match collections and fill pt eta phi and charge histos for
// candidate jet, fill delta_x_VS_y histos for matched couples,
// book and fill delta_frac_VS_frac histos for matched couples
}
}
#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(PFJetDQMAnalyzer);
|