|
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
|
#include <Alignment/HIPAlignmentAlgorithm/interface/HIPMonitorConfig.h>
HIPMonitorConfig::HIPMonitorConfig(const edm::ParameterSet& cfg)
: cfgMonitor(cfg.getParameter<edm::ParameterSet>("monitorConfig")),
outfilecore(cfgMonitor.getParameter<std::string>("outfile")),
maxEventsPerJob(cfgMonitor.getParameter<int>("maxEventsPerJob")),
fillTrackMonitoring(cfgMonitor.getParameter<bool>("fillTrackMonitoring")),
maxTracks(cfgMonitor.getParameter<int>("maxTracks")),
trackmonitorvars(maxEventsPerJob, maxTracks),
fillTrackHitMonitoring(cfgMonitor.getParameter<bool>("fillTrackHitMonitoring")),
maxHits(cfgMonitor.getParameter<int>("maxHits")),
hitmonitorvars(maxHits),
eventCounter(0),
hitCounter(0) {
outfile = cfg.getParameter<std::string>("outpath") + outfilecore;
}
HIPMonitorConfig::HIPMonitorConfig(const HIPMonitorConfig& other)
: cfgMonitor(other.cfgMonitor),
outfilecore(other.outfilecore),
maxEventsPerJob(other.maxEventsPerJob),
fillTrackMonitoring(other.fillTrackMonitoring),
maxTracks(other.maxTracks),
fillTrackHitMonitoring(other.fillTrackHitMonitoring),
maxHits(other.maxHits),
outfile(other.outfile),
eventCounter(other.eventCounter),
hitCounter(other.hitCounter) {}
bool HIPMonitorConfig::checkNevents() {
bool res = (maxEventsPerJob < 0 || maxEventsPerJob > eventCounter);
eventCounter++;
return res;
}
bool HIPMonitorConfig::checkNhits() {
bool res = (maxHits < 0 || maxHits > hitCounter);
hitCounter++;
return res;
}
void HIPTrackMonitorVariables::bookBranches() {
if (tree != nullptr) {
tree->Branch("Ntracks", &m_Ntracks);
tree->Branch("Nhits", &m_Nhits);
//tree->Branch("DataType", &m_datatype); This is done in HIPAlignmentAlgorithm
tree->Branch("nhPXB", &m_nhPXB);
tree->Branch("nhPXF", &m_nhPXF);
tree->Branch("nhTIB", &m_nhTIB);
tree->Branch("nhTOB", &m_nhTOB);
tree->Branch("nhTID", &m_nhTID);
tree->Branch("nhTEC", &m_nhTEC);
tree->Branch("Pt", &m_Pt);
tree->Branch("P", &m_P);
tree->Branch("Eta", &m_Eta);
tree->Branch("Phi", &m_Phi);
tree->Branch("Chi2n", &m_Chi2n);
tree->Branch("d0", &m_d0);
tree->Branch("dz", &m_dz);
tree->Branch("wt", &m_wt);
}
}
void HIPHitMonitorVariables::bookBranches() {
if (tree != nullptr) {
tree->Branch("Id", &m_detId, "Id/i");
tree->Branch("sinTheta", &m_sinTheta);
tree->Branch("impactAngle", &m_angle);
tree->Branch("wt", &m_hitwt);
tree->Branch("probPresent", &m_hasHitProb);
tree->Branch("probXY", &m_probXY);
tree->Branch("probQ", &m_probQ);
tree->Branch("qualityWord", &m_rawQualityWord);
}
}
void HIPTrackMonitorVariables::fill() {
if (tree == nullptr)
return;
if (maxNEvents >= 0 && nEvents >= maxNEvents)
return;
bool doFill = false;
m_Ntracks = m_Pt.size();
if (m_Ntracks == 0)
return;
if (maxTracksRcd < 0)
doFill = true;
else if (nTracks < maxTracksRcd) {
if ((nTracks + m_Ntracks) < maxTracksRcd)
doFill = true;
else {
int NewSize = maxTracksRcd - nTracks;
if ((int)m_Pt.size() < NewSize)
NewSize = m_Pt.size();
// Do not touch m_Ntracks, just resize these vectors
resizeVectors(NewSize);
doFill = true;
}
} else {
// Record only m_Ntracks
int tmpNtracks = m_Ntracks;
resetPerEvent();
m_Ntracks = tmpNtracks;
doFill = true;
}
if (doFill) {
tree->Fill();
nEvents++;
nTracks += m_Ntracks;
}
resetPerEvent();
}
void HIPHitMonitorVariables::fill() {
if (tree == nullptr)
return;
bool doFill = (maxHitsRcd < 0 || nHits < maxHitsRcd);
if (doFill) {
tree->Fill();
nHits++;
}
resetPerHit();
}
|