|
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
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
|
#ifndef DQMServices_Core_DQMOneEDAnalyzer_h
#define DQMServices_Core_DQMOneEDAnalyzer_h
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/Run.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Framework/interface/one/EDProducer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/Histograms/interface/DQMToken.h"
/**
* A "one" module base class that can only produce per-run histograms. This
* allows to easily wrap non-thread-safe code and is ok'ish performance-wise,
* since it only blocks concurrent runs, not concurrent lumis.
* It can be combined with edm::LuminosityBlockCache to watch per-lumi things,
* and fill per-run histograms with the results.
*/
template <typename... Args>
class DQMOneEDAnalyzer
: public edm::one::EDProducer<edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args...> {
public:
typedef dqm::reco::DQMStore DQMStore;
typedef dqm::reco::MonitorElement MonitorElement;
virtual bool getCanSaveByLumi() { return true; }
// framework calls in the order of invocation
DQMOneEDAnalyzer() {
// for whatever reason we need the explicit `template` keyword here.
runToken_ = this->template produces<DQMToken, edm::Transition::EndRun>("DQMGenerationRecoRun");
}
void beginRun(edm::Run const& run, edm::EventSetup const& setup) final {
// if we run booking multiple times because there are multiple runs in a
// job, this is needed to make sure all existing MEs are in a valid state
// before the booking code runs.
edm::Service<DQMStore>()->initLumi(run.run(), /* lumi */ 0, this->moduleDescription().id());
edm::Service<DQMStore>()->enterLumi(run.run(), /* lumi */ 0, this->moduleDescription().id());
dqmBeginRun(run, setup);
edm::Service<DQMStore>()->bookTransaction(
[this, &run, &setup](DQMStore::IBooker& booker) {
booker.cd();
this->bookHistograms(booker, run, setup);
},
this->moduleDescription().id(),
this->getCanSaveByLumi());
edm::Service<DQMStore>()->initLumi(run.run(), /* lumi */ 0, this->moduleDescription().id());
edm::Service<DQMStore>()->enterLumi(run.run(), /* lumi */ 0, this->moduleDescription().id());
}
void accumulate(edm::Event const& event, edm::EventSetup const& setup) override {
auto& lumi = event.getLuminosityBlock();
edm::Service<dqm::legacy::DQMStore>()->enterLumi(
lumi.run(), lumi.luminosityBlock(), this->moduleDescription().id());
analyze(event, setup);
edm::Service<dqm::legacy::DQMStore>()->leaveLumi(
lumi.run(), lumi.luminosityBlock(), this->moduleDescription().id());
}
void endRunProduce(edm::Run& run, edm::EventSetup const& setup) final {
dqmEndRun(run, setup);
edm::Service<DQMStore>()->leaveLumi(run.run(), /* lumi */ 0, this->moduleDescription().id());
run.emplace<DQMToken>(runToken_);
}
// Subsystems could safely override this, but any changes to MEs would not be
// noticeable since the product was made already.
void endRun(edm::Run const&, edm::EventSetup const&) final {}
protected:
// methods to be implemented by the user, in order of invocation
virtual void dqmBeginRun(edm::Run const&, edm::EventSetup const&) {}
virtual void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) = 0;
virtual void analyze(edm::Event const&, edm::EventSetup const&) {}
virtual void dqmEndRun(edm::Run const&, edm::EventSetup const&) {}
edm::EDPutTokenT<DQMToken> runToken_;
};
/**
* A "one" module base class that can also watch lumi transitions and produce
* per-lumi MEs. This should be used carefully, since it will block concurrent
* lumisections in the entire job!
* Combining with edm::LuminosityBlockCache is pointless and will not work
* properly, due to the ordering of global/produce transitions.
*/
template <typename... Args>
class DQMOneLumiEDAnalyzer
: public DQMOneEDAnalyzer<edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, Args...> {
public:
bool getCanSaveByLumi() override { return true; }
// framework calls in the order of invocation
DQMOneLumiEDAnalyzer() {
// for whatever reason we need the explicit `template` keyword here.
lumiToken_ = this->template produces<DQMToken, edm::Transition::EndLuminosityBlock>("DQMGenerationRecoLumi");
}
void beginLuminosityBlock(edm::LuminosityBlock const& lumi, edm::EventSetup const& setup) final {
edm::Service<dqm::legacy::DQMStore>()->enterLumi(
lumi.run(), lumi.luminosityBlock(), this->moduleDescription().id());
dqmBeginLuminosityBlock(lumi, setup);
}
void accumulate(edm::Event const& event, edm::EventSetup const& setup) override { this->analyze(event, setup); }
void endLuminosityBlockProduce(edm::LuminosityBlock& lumi, edm::EventSetup const& setup) final {
dqmEndLuminosityBlock(lumi, setup);
// fully qualified name required for... reasons.
edm::Service<dqm::legacy::DQMStore>()->leaveLumi(
lumi.run(), lumi.luminosityBlock(), this->moduleDescription().id());
lumi.emplace(lumiToken_);
}
// Subsystems could safely override this, but any changes to MEs would not be
// noticeable since the product was made already.
void endLuminosityBlock(edm::LuminosityBlock const&, edm::EventSetup const&) final {}
protected:
// methods to be implemented by the user, in order of invocation
virtual void dqmBeginLuminosityBlock(edm::LuminosityBlock const&, edm::EventSetup const&) {}
virtual void dqmEndLuminosityBlock(edm::LuminosityBlock const&, edm::EventSetup const&) {}
edm::EDPutTokenT<DQMToken> lumiToken_;
};
#endif // DQMServices_Core_DQMOneEDAnalyzer_h
|