|
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
|
#ifndef BeamMonitorBx_H
#define BeamMonitorBx_H
/** \class BeamMonitorBx
* *
* \author Geng-yuan Jeng/UC Riverside
* Francisco Yumiceva/FNAL
*
*/
// C++
#include <string>
// CMS
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "RecoVertex/BeamSpotProducer/interface/BeamFitter.h"
#include <fstream>
//
// class declaration
//
class BeamMonitorBx
: public edm::one::EDAnalyzer<edm::one::SharedResources, edm::one::WatchRuns, edm::one::WatchLuminosityBlocks> {
public:
typedef dqm::legacy::MonitorElement MonitorElement;
typedef dqm::legacy::DQMStore DQMStore;
BeamMonitorBx(const edm::ParameterSet&);
~BeamMonitorBx() override;
typedef int BxNum;
typedef std::map<BxNum, reco::BeamSpot> BeamSpotMapBx;
protected:
// BeginJob
void beginJob() override;
// BeginRun
void beginRun(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 endRun(const edm::Run& r, const edm::EventSetup& c) override;
private:
void FitAndFill(const edm::LuminosityBlock& lumiSeg, int&, int&, int&);
void BookTables(int, std::map<std::string, std::string>&, std::string);
void BookTrendHistos(bool, int, std::map<std::string, std::string>&, std::string, const TString&, const TString&);
void FillTables(int, int, std::map<std::string, std::string>&, reco::BeamSpot&, std::string);
void FillTrendHistos(int, int, std::map<std::string, std::string>&, reco::BeamSpot&, const TString&);
void weight(BeamSpotMapBx&, const BeamSpotMapBx&);
void weight(double& mean, double& meanError, const double& val, const double& valError);
void formatFitTime(char*, const std::time_t&);
edm::ParameterSet parameters_;
std::string monitorName_;
edm::InputTag bsSrc_; // beam spot
int fitNLumi_;
int resetFitNLumi_;
bool debug_;
DQMStore* dbe_;
BeamFitter* theBeamFitter;
unsigned int countBx_;
int countEvt_; //counter
int countLumi_; //counter
int beginLumiOfBSFit_;
int endLumiOfBSFit_;
int lastlumi_; // previous LS processed
int nextlumi_; // next LS of Fit
int firstlumi_; // first LS with good fit
int countGoodFit_;
std::time_t refBStime[2];
bool resetHistos_;
bool processed_;
// ----------member data ---------------------------
BeamSpotMapBx fbspotMap; //for weighted beam spots of each bunch
std::map<std::string, std::string> varMap;
std::map<std::string, std::string> varMap1;
// MonitorElements:
std::map<TString, MonitorElement*> hs; // Tables
std::map<TString, MonitorElement*> hst; // Trending Histos
//Test
// MonitorElement * h_x0;
//
std::time_t tmpTime;
std::time_t refTime;
std::time_t startTime;
edm::TimeValue_t ftimestamp;
};
#endif
|