|
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 DQM_BeamMonitor_AlcaBeamMonitor_h
#define DQM_BeamMonitor_AlcaBeamMonitor_h
/** \class AlcaBeamMonitor
* *
* \author Lorenzo Uplegger/FNAL
* modified by Simone Gennai INFN/Bicocca
*/
// C++
#include <map>
#include <array>
#include <vector>
#include <string>
#include <utility>
// CMS
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/ESGetToken.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/DQMOneEDAnalyzer.h"
#include "DataFormats/Provenance/interface/LuminosityBlockID.h"
#include "DataFormats/BeamSpot/interface/BeamSpot.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "CondFormats/DataRecord/interface/BeamSpotObjectsRcd.h"
#include "CondFormats/BeamSpotObjects/interface/BeamSpotObjects.h"
class BeamFitter;
class PVFitter;
namespace alcabeammonitor {
struct pvPosAndErr {
// Array of pairs: (value, error) for x, y, z
std::array<std::pair<double, double>, 3> data;
// Constructor initializes the array with values and errors from a reco::Vertex
pvPosAndErr(const reco::Vertex& vertex)
: data{{{vertex.x(), vertex.xError()}, {vertex.y(), vertex.yError()}, {vertex.z(), vertex.zError()}}} {}
// Accessor functions that return pairs (value, error) directly
std::pair<double, double> xWithError() const { return data[0]; }
std::pair<double, double> yWithError() const { return data[1]; }
std::pair<double, double> zWithError() const { return data[2]; }
};
struct BeamSpotInfo {
std::vector<std::vector<pvPosAndErr>> vertices_;
typedef std::map<std::string, reco::BeamSpot> BeamSpotContainer;
BeamSpotContainer beamSpotMap_;
};
} // namespace alcabeammonitor
class AlcaBeamMonitor : public DQMOneEDAnalyzer<edm::LuminosityBlockCache<alcabeammonitor::BeamSpotInfo>> {
public:
AlcaBeamMonitor(const edm::ParameterSet&);
static void fillDescriptions(edm::ConfigurationDescriptions&);
protected:
void bookHistograms(DQMStore::IBooker&, edm::Run const&, edm::EventSetup const&) override;
void analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) override;
std::shared_ptr<alcabeammonitor::BeamSpotInfo> globalBeginLuminosityBlock(
const edm::LuminosityBlock& iLumi, const edm::EventSetup& iSetup) const override;
void globalEndLuminosityBlock(const edm::LuminosityBlock& iLumi, const edm::EventSetup& iSetup) override;
void dqmEndRun(edm::Run const&, edm::EventSetup const&) override;
private:
//Typedefs
// BF,BS...
typedef std::map<std::string, reco::BeamSpot> BeamSpotContainer;
// x,y,z,sigmax(y,z)... [run,lumi] Histo name
typedef std::map<std::string, std::map<std::string, std::map<std::string, MonitorElement*>>> HistosContainer;
// x,y,z,sigmax(y,z)... [run,lumi] Histo name
typedef std::map<std::string, std::map<std::string, std::map<std::string, int>>> PositionContainer;
//Parameters
std::string monitorName_;
const edm::EDGetTokenT<reco::VertexCollection> primaryVertexLabel_;
const edm::EDGetTokenT<reco::TrackCollection> trackLabel_;
const edm::EDGetTokenT<reco::BeamSpot> scalerLabel_;
const edm::ESGetToken<BeamSpotObjects, BeamSpotObjectsRcd> beamSpotToken_;
bool perLSsaving_; //to avoid nanoDQMIO crashing, driven by DQMServices/Core/python/DQMStore_cfi.py
//Service variables
int numberOfValuesToSave_;
std::unique_ptr<BeamFitter> theBeamFitter_;
std::unique_ptr<PVFitter> thePVFitter_;
std::vector<int> processedLumis_;
// MonitorElements:
MonitorElement* hD0Phi0_;
MonitorElement* hDxyBS_;
//mutable MonitorElement* theValuesContainer_;
//Containers
HistosContainer histosMap_;
PositionContainer positionsMap_;
std::vector<std::string> varNamesV_; //x,y,z,sigmax(y,z)
std::multimap<std::string, std::string> histoByCategoryNames_; //run, lumi
};
#endif
|