|
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
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
|
#ifndef FakeBeamMonitor_H
#define FakeBeamMonitor_H
/** \class FakeBeamMonitor
* *
* \author Francesco Brivio - Milano-Bicocca
* \decription: Same behaviour of BeamMonitor except the
* beamspot is randomly generated and not fitted from tracks/PVs
*
*/
// C++
#include <string>
#include <memory>
// CMS
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/Common/interface/TriggerResults.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/DQMOneEDAnalyzer.h"
#include "CondCore/DBOutputService/interface/OnlineDBOutputService.h"
#include <fstream>
#include "TRandom3.h"
//
// class declaration
//
class FakeBeamMonitor : public DQMOneEDAnalyzer<edm::one::WatchLuminosityBlocks> {
public:
FakeBeamMonitor(const edm::ParameterSet&);
~FakeBeamMonitor() override;
protected:
//BeginRun
void dqmBeginRun(edm::Run const&, edm::EventSetup const&) override;
void bookHistograms(DQMStore::IBooker& i, const edm::Run& r, const edm::EventSetup& c) override;
void analyze(const edm::Event& e, const edm::EventSetup& c) override;
void beginLuminosityBlock(const edm::LuminosityBlock& lumiSeg, const edm::EventSetup& context) override;
void endLuminosityBlock(const edm::LuminosityBlock& lumiSeg, const edm::EventSetup& c) override;
// EndRun
void dqmEndRun(const edm::Run& r, const edm::EventSetup& c) override;
private:
void FitAndFill(const edm::LuminosityBlock& lumiSeg, int&, int&, int&);
void RestartFitting();
void scrollTH1(TH1*, std::time_t);
bool testScroll(std::time_t&, std::time_t&);
void formatFitTime(char*, const std::time_t&);
std::string getGMTstring(const std::time_t&);
const int dxBin_;
const double dxMin_;
const double dxMax_;
const int vxBin_;
const double vxMin_;
const double vxMax_;
const int phiBin_;
const double phiMin_;
const double phiMax_;
const int dzBin_;
const double dzMin_;
const double dzMax_;
std::string monitorName_;
std::string recordName_; // output BeamSpotOnline Record name
edm::Service<cond::service::OnlineDBOutputService> onlineDbService_; // DB service
int DBloggerReturn_; // DB logger return value
int fitNLumi_;
int fitPVNLumi_;
int resetFitNLumi_;
int resetPVNLumi_;
int intervalInSec_;
bool debug_;
bool onlineMode_;
int countEvt_; //counter
int countLumi_; //counter
int beginLumiOfBSFit_;
int endLumiOfBSFit_;
int beginLumiOfPVFit_;
int endLumiOfPVFit_;
int lastlumi_; // previous LS processed
int nextlumi_; // next LS of Fit
std::time_t refBStime[2];
std::time_t refPVtime[2];
unsigned int nthBSTrk_;
int nFitElements_;
int nFits_;
double deltaSigCut_;
unsigned int min_Ntrks_;
double maxZ_;
unsigned int minNrVertices_;
double minVtxNdf_;
double minVtxWgt_;
bool resetHistos_;
bool StartAverage_;
int firstAverageFit_;
int countGapLumi_;
bool processed_;
bool useLockRecords_;
// ----------member data ---------------------------
// std::vector<BSTrkParameters> fBSvector;
reco::BeamSpot refBS;
reco::BeamSpot preBS;
// MonitorElements:
MonitorElement* h_nTrk_lumi;
MonitorElement* h_nVtx_lumi;
MonitorElement* h_nVtx_lumi_all;
MonitorElement* h_d0_phi0;
MonitorElement* h_trk_z0;
MonitorElement* h_vx_vy;
MonitorElement* h_vx_dz;
MonitorElement* h_vy_dz;
MonitorElement* h_trkPt;
MonitorElement* h_trkVz;
MonitorElement* fitResults;
MonitorElement* h_x0;
MonitorElement* h_y0;
MonitorElement* h_z0;
MonitorElement* h_sigmaX0;
MonitorElement* h_sigmaY0;
MonitorElement* h_sigmaZ0;
MonitorElement* h_nVtx;
MonitorElement* h_nVtx_st;
MonitorElement* h_PVx[2];
MonitorElement* h_PVy[2];
MonitorElement* h_PVz[2];
MonitorElement* h_PVxz;
MonitorElement* h_PVyz;
MonitorElement* pvResults;
std::vector<MonitorElement*> hs;
// The histo of the primary vertex for DQM gui
std::map<int, std::vector<float> > mapPVx, mapPVy, mapPVz;
std::map<int, std::vector<int> > mapNPV;
//keep track of beginLuminosity block and time
std::map<int, int> mapBeginBSLS, mapBeginPVLS;
std::map<int, std::time_t> mapBeginBSTime, mapBeginPVTime;
//these maps are needed to keep track of size of theBSvector and pvStore
std::map<int, std::size_t> mapLSBSTrkSize;
std::map<int, size_t> mapLSPVStoreSize;
//to correct the cutFlot Table
std::map<int, TH1F> mapLSCF;
// Summary:
Float_t reportSummary_;
Float_t summarySum_;
Float_t summaryContent_[3];
MonitorElement* reportSummary;
MonitorElement* reportSummaryContents[3];
MonitorElement* reportSummaryMap;
MonitorElement* cutFlowTable;
// variables for beam fit
//
std::time_t tmpTime;
std::time_t startTime;
std::time_t refTime;
edm::TimeValue_t ftimestamp;
// Current run
int frun;
int lastNZbin; // last non zero bin of time histos
TRandom3* rndm_;
};
#endif
// Local Variables:
// show-trailing-whitespace: t
// truncate-lines: t
// End:
|