|
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
|
#include "DQM/DataScouting/plugins/AlphaTVarAnalyzer.h"
#include "DataFormats/JetReco/interface/CaloJet.h"
#include "DataFormats/JetReco/interface/CaloJetCollection.h"
#include "DataFormats/EgammaCandidates/interface/Electron.h"
#include "DataFormats/EgammaCandidates/interface/ElectronFwd.h"
#include "DataFormats/MuonReco/interface/Muon.h"
#include "DataFormats/RecoCandidate/interface/RecoChargedCandidate.h"
#include <cmath>
// A simple constructor which takes as input only the name of the PF jet
// collection
AlphaTVarAnalyzer::AlphaTVarAnalyzer(const edm::ParameterSet &conf)
: ScoutingAnalyzerBase(conf),
m_jetCollectionTag(
conf.getUntrackedParameter<edm::InputTag>("jetCollectionName", edm::InputTag("hltCaloJetIDPassed"))),
m_alphaTVarCollectionTag(conf.getUntrackedParameter<edm::InputTag>("alphaTVarCollectionName")) {
// set Token(-s)
m_alphaTVarCollectionTagToken_ =
consumes<std::vector<double>>(conf.getUntrackedParameter<edm::InputTag>("alphaTVarCollectionName"));
}
AlphaTVarAnalyzer::~AlphaTVarAnalyzer() {}
// Function to book the Monitoring Elements.
void AlphaTVarAnalyzer::bookHistograms(DQMStore::IBooker &iBooker, edm::Run const &, edm::EventSetup const &) {
ScoutingAnalyzerBase::prepareBooking(iBooker);
// the full inclusive histograms
m_HTAlphaT = bookH2withSumw2(iBooker,
"HTvsAlphaT",
"H_{T} vs #alpha_{T} (All Events)",
400,
0.,
4000.,
50,
0.,
1.,
"H_{T} [GeV]",
"#alpha_{T}");
m_HTAlphaTg0p55 =
bookH1withSumw2(iBooker, "HTvsAlphaTg0p55", "H_{T} (#alpha_{T} > 0.55)", 400, 0., 4000., "H_{T} [GeV]");
m_HTAlphaTg0p60 =
bookH1withSumw2(iBooker, "HTvsAlphaTg0p60", "H_{T} (#alpha_{T} > 0.60)", 400, 0., 4000., "H_{T} [GeV]");
}
// Usual analyze method
void AlphaTVarAnalyzer::analyze(const edm::Event &iEvent, const edm::EventSetup &c) {
edm::Handle<std::vector<double>> alphaTvar_handle;
iEvent.getByToken(m_alphaTVarCollectionTagToken_, alphaTvar_handle);
if (alphaTvar_handle->size() > 1) {
const double AlphaT = alphaTvar_handle->at(0);
const double HT = alphaTvar_handle->at(1);
m_HTAlphaT->Fill(HT, AlphaT);
if (AlphaT > 0.55)
m_HTAlphaTg0p55->Fill(HT);
if (AlphaT > 0.60)
m_HTAlphaTg0p60->Fill(HT);
}
}
|