|
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
|
#ifndef BTagPerformanceAnalyzerOnData_H
#define BTagPerformanceAnalyzerOnData_H
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "DataFormats/JetReco/interface/Jet.h"
#include "DQMOffline/RecoB/interface/AcceptJet.h"
#include "DQMOffline/RecoB/interface/JetTagPlotter.h"
#include "DQMOffline/RecoB/interface/TagCorrelationPlotter.h"
#include "DQMOffline/RecoB/interface/BaseTagInfoPlotter.h"
#include "SimDataFormats/GeneratorProducts/interface/GenEventInfoProduct.h"
#include "JetMETCorrections/JetCorrector/interface/JetCorrector.h"
/** \class BTagPerformanceAnalyzerOnData
*
* Top level steering routine for b tag performance analysis.
*
*/
class BTagPerformanceAnalyzerOnData : public DQMEDAnalyzer {
public:
explicit BTagPerformanceAnalyzerOnData(const edm::ParameterSet& pSet);
~BTagPerformanceAnalyzerOnData() override;
void analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) override;
private:
inline bool jetRefCompare(const edm::RefToBase<reco::Jet>& j1, const edm::RefToBase<reco::Jet>& j2) {
return j1.id() < j2.id() || (j1.id() == j2.id() && j1.key() < j2.key());
}
// Get histogram plotting options from configuration.
void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
EtaPtBin getEtaPtBin(const int& iEta, const int& iPt);
std::vector<std::string> tiDataFormatType;
AcceptJet jetSelector; // Decides if jet and parton satisfy kinematic cuts.
std::vector<double> etaRanges, ptRanges;
bool doJEC;
edm::InputTag slInfoTag;
std::vector<std::vector<std::unique_ptr<JetTagPlotter>>> binJetTagPlotters;
std::vector<std::vector<std::unique_ptr<TagCorrelationPlotter>>> binTagCorrelationPlotters;
std::vector<std::vector<std::unique_ptr<BaseTagInfoPlotter>>> binTagInfoPlotters;
std::vector<edm::InputTag> jetTagInputTags;
std::vector<std::pair<edm::InputTag, edm::InputTag>> tagCorrelationInputTags;
std::vector<std::vector<edm::InputTag>> tagInfoInputTags;
// Contains plots for each bin of rapidity and pt.
std::vector<edm::ParameterSet> moduleConfig;
//add consumes
edm::EDGetTokenT<reco::JetCorrector> jecMCToken;
edm::EDGetTokenT<reco::JetCorrector> jecDataToken;
edm::EDGetTokenT<GenEventInfoProduct> genToken;
edm::EDGetTokenT<reco::SoftLeptonTagInfoCollection> slInfoToken;
std::vector<edm::EDGetTokenT<reco::JetTagCollection>> jetTagToken;
std::vector<std::pair<edm::EDGetTokenT<reco::JetTagCollection>, edm::EDGetTokenT<reco::JetTagCollection>>>
tagCorrelationToken;
std::vector<std::vector<edm::EDGetTokenT<edm::View<reco::BaseTagInfo>>>> tagInfoToken;
};
#endif
|