|
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
|
#ifndef DQM_TRACKERREMAPPER_PHASE1PIXELMAPS_H
#define DQM_TRACKERREMAPPER_PHASE1PIXELMAPS_H
#include "TH2Poly.h"
#include "TGraph.h"
#include "TH1.h"
#include "TH2.h"
#include "TStyle.h"
#include "TCanvas.h"
#include <fmt/printf.h>
#include <fstream>
#include <boost/tokenizer.hpp>
#include <boost/range/adaptor/indexed.hpp>
#include "FWCore/ParameterSet/interface/FileInPath.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "CalibTracker/StandaloneTrackerTopology/interface/StandaloneTrackerTopology.h"
#ifndef PHASE1PIXELMAP_STANDALONE
#define LOGDEBUG(x) LogDebug(x)
#else
#define LOGDEBUG(x) std::cout << x << ": "
#endif
using indexedCorners = std::map<unsigned int, std::pair<std::vector<float>, std::vector<float>>>;
/*--------------------------------------------------------------------
/ Ancillary class to build pixel phase-1 tracker maps
/--------------------------------------------------------------------*/
class Phase1PixelMaps {
public:
Phase1PixelMaps(const char* option)
: m_option{option},
m_isBooked{std::make_pair(false, false)},
m_knownNames{{}},
m_trackerTopo{StandaloneTrackerTopology::fromTrackerParametersXMLFile(
edm::FileInPath("Geometry/TrackerCommonData/data/PhaseI/trackerParameters.xml").fullPath())} {
// set the rescale to true by default
m_autorescale = true;
// store the file in path for the corners (BPIX)
for (unsigned int i = 1; i <= 4; i++) {
m_cornersBPIX.push_back(edm::FileInPath(Form("DQM/SiStripMonitorClient/data/Geometry/vertices_barrel_%i", i)));
}
// store the file in path for the corners (BPIX)
for (int j : {-3, -2, -1, 1, 2, 3}) {
m_cornersFPIX.push_back(edm::FileInPath(Form("DQM/SiStripMonitorClient/data/Geometry/vertices_forward_%i", j)));
}
}
~Phase1PixelMaps() = default;
// set of no rescale
inline void setNoRescale() { m_autorescale = false; }
// set option, but only if not already set
void resetOption(const char* option);
// book them all
void book(const std::string& currentHistoName, const char* what, const char* zaxis);
// booking methods
void bookBarrelHistograms(const std::string& currentHistoName, const char* what, const char* zaxis);
void bookForwardHistograms(const std::string& currentHistoName, const char* what, const char* zaxis);
// fill them all
void fill(const std::string& currentHistoName, unsigned int id, double value);
// filling methods
void fillBarrelBin(const std::string& currentHistoName, unsigned int id, double value);
void fillForwardBin(const std::string& currentHistoName, unsigned int id, double value);
// beautification methods
void beautifyAllHistograms();
void setBarrelScale(const std::string& currentHistoName, std::pair<float, float> extrema);
void setForwardScale(const std::string& currentHistoName, std::pair<float, float> extrema);
// drawing methos
void drawBarrelMaps(const std::string& currentHistoName, TCanvas& canvas, const char* drawOption = nullptr);
void drawForwardMaps(const std::string& currentHistoName, TCanvas& canvas, const char* drawOption = nullptr);
void drawSummaryMaps(const std::string& currentHistoName, TCanvas& canvas, const char* drawOption = nullptr);
private:
Option_t* m_option;
bool m_autorescale;
std::pair<bool, bool> m_isBooked;
std::vector<std::string> m_knownNames;
TrackerTopology m_trackerTopo;
std::map<uint32_t, std::shared_ptr<TGraph>> bins, binsSummary;
std::map<std::string, std::vector<std::shared_ptr<TH2Poly>>> pxbTh2PolyBarrel;
std::map<std::string, std::shared_ptr<TH2Poly>> pxbTh2PolyBarrelSummary;
std::map<std::string, std::vector<std::shared_ptr<TH2Poly>>> pxfTh2PolyForward;
std::map<std::string, std::shared_ptr<TH2Poly>> pxfTh2PolyForwardSummary;
std::vector<edm::FileInPath> m_cornersBPIX;
std::vector<edm::FileInPath> m_cornersFPIX;
const indexedCorners retrieveCorners(const std::vector<edm::FileInPath>& cornerFiles, const unsigned int reads);
// called by book histograms
void bookBarrelBins(const std::string& currentHistoName);
void bookForwardBins(const std::string& currentHistoName);
// graphics
void makeNicePlotStyle(TH1* hist);
void adjustCanvasMargins(TVirtualPad* pad, float top, float bottom, float left, float right);
void rescaleAllBarrel(const std::string& currentHistoName);
void rescaleAllForward(const std::string& currentHistoName);
};
#endif
|