|
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
|
#ifndef ALIGNMENT_OFFLINEVALIDATION_OUTPVTRENDS_H_
#define ALIGNMENT_OFFLINEVALIDATION_OUTPVTRENDS_H_
#include "TArrow.h"
#include "TAxis.h"
#include "TCanvas.h"
#include "TF1.h"
#include "TFile.h"
#include "TGaxis.h"
#include "TGraph.h"
#include "TGraphAsymmErrors.h"
#include "TGraphErrors.h"
#include "TH1F.h"
#include "TH2F.h"
#include "TLatex.h"
#include "TLegend.h"
#include "TList.h"
#include "TMath.h"
#include "TObjArray.h"
#include "TObjString.h"
#include "TPad.h"
#include "TPaveText.h"
#include "TProcPool.h"
#include "TProfile.h"
#include "TROOT.h"
#include "TStyle.h"
#include "TSystem.h"
#include "TSystemDirectory.h"
#include "TSystemFile.h"
#include <TStopwatch.h>
#include <algorithm>
#include <bitset>
#include <fstream>
#include <functional>
#include <iostream>
#include <iterator>
#include <map>
#include <sstream>
#include <vector>
/*!
* \def basically the y-values of a TGraph
*/
typedef std::map<TString, std::vector<double> > alignmentTrend;
// auxilliary struct to store
// histogram features
struct unrolledHisto {
double m_y_min;
double m_y_max;
unsigned int m_n_bins;
std::vector<double> m_bincontents;
unrolledHisto() {
m_y_min = 0.;
m_y_max = 0.;
m_n_bins = 0;
m_bincontents.clear();
} // first constructor empty
unrolledHisto(const double &y_min,
const double &y_max,
const unsigned int &n_bins,
const std::vector<double> &bincontents) {
m_y_min = y_min;
m_y_max = y_max;
m_n_bins = n_bins;
m_bincontents = bincontents;
} //look, a constructor
double get_y_min() { return m_y_min; }
double get_y_max() { return m_y_max; }
unsigned int get_n_bins() { return m_n_bins; }
std::vector<double> get_bin_contents() { return m_bincontents; }
double get_integral() {
double ret(0.);
for (const auto &binc : m_bincontents) {
ret += binc;
}
return ret;
}
};
struct outPVtrends {
/*! \struct outPVtrends
* \brief Structure outPVtrends
* Contains the ensemble of all the alignmentTrends built by the functor
*
* @param m_index int, to keep track of which chunk of data has been processed
* @param m_runs std::vector, list of the run processed in this section
* @param m_dxyPhiMeans alignmentTrend of the mean values of the profile dxy vs phi
* @param m_dxyPhiChi2 alignmentTrend of chi2 of the linear fit per profile dxy vs phi
* @param m_dxyPhiKS alignmentTrend of Kolmogorow-Smirnov score of comparison of dxy vs phi profile with flat line
* @param m_dxyPhiHi alignmentTrend of the highest value of the profile dxy vs phi
* @param m_dxyPhiLo alignmentTrend of the lowest value of the profile dxy vs phi
* @param m_dxyEtaMeans alignmentTrend of the mean values of the profile dxy vs eta
* @param m_dxyEtaChi2 alignmentTrend of chi2 of the linear fit per profile dxy vs eta
* @param m_dxyEtaKS alignmentTrend of Kolmogorow-Smirnov score of comparison of dxy vs eta profile with flat line
* @param m_dxyEtaHi alignmentTrend of the highest value of the profile dxy vs eta
* @param m_dxyEtaLo alignmentTrend of the lowest value of the profile dxy vs eta
* @param m_dzPhiMeans alignmentTrend of the mean values of the profile dz vs phi
* @param m_dzPhiChi2 alignmentTrend of chi2 of the linear fit per profile dz vs phi
* @param m_dzPhiKS alignmentTrend of Kolmogorow-Smirnov score of comparison of dz vs phi profile with flat line
* @param m_dzPhiHi alignmentTrend of the highest value of the profile dz vs phi
* @param m_dzPhiLo alignmentTrend of the lowest value of the profile dz vs phi
* @param m_dzEtaMeans alignmentTrend of the mean values of the profile dz vs eta
* @param m_dzEtaChi2 alignmentTrend of chi2 of the linear fit per profile dz vs eta
* @param m_dzEtaKS alignmentTrend of Kolmogorow-Smirnov score of comparison of dz vs eta profile with flat line
* @param m_dzEtaHi alignmentTrend of the highest value of the profile dz vs eta
* @param m_dzEtaLo alignmentTrend of the lowest value of the profile dz vs eta
* @param m_dxyVect map of the unrolled histograms for dxy residuals
* @param m_dzVect map of the unrolled histograms for dz residulas
*/
// empty constructor
outPVtrends() { init(); }
int m_index;
std::vector<double> m_runs;
alignmentTrend m_dxyPhiMeans;
alignmentTrend m_dxyPhiChi2;
alignmentTrend m_dxyPhiKS;
alignmentTrend m_dxyPhiHi;
alignmentTrend m_dxyPhiLo;
alignmentTrend m_dxyEtaMeans;
alignmentTrend m_dxyEtaChi2;
alignmentTrend m_dxyEtaKS;
alignmentTrend m_dxyEtaHi;
alignmentTrend m_dxyEtaLo;
alignmentTrend m_dzPhiMeans;
alignmentTrend m_dzPhiChi2;
alignmentTrend m_dzPhiKS;
alignmentTrend m_dzPhiHi;
alignmentTrend m_dzPhiLo;
alignmentTrend m_dzEtaMeans;
alignmentTrend m_dzEtaChi2;
alignmentTrend m_dzEtaKS;
alignmentTrend m_dzEtaHi;
alignmentTrend m_dzEtaLo;
std::map<TString, std::vector<unrolledHisto> > m_dxyVect;
std::map<TString, std::vector<unrolledHisto> > m_dzVect;
void init() {
m_index = -1;
m_runs.clear();
m_dxyPhiMeans.clear();
m_dxyPhiChi2.clear();
m_dxyPhiKS.clear();
m_dxyPhiHi.clear();
m_dxyPhiLo.clear();
m_dxyEtaMeans.clear();
m_dxyEtaChi2.clear();
m_dxyEtaKS.clear();
m_dxyEtaHi.clear();
m_dxyEtaLo.clear();
m_dzPhiMeans.clear();
m_dzPhiChi2.clear();
m_dzPhiKS.clear();
m_dzPhiHi.clear();
m_dzPhiLo.clear();
m_dzEtaMeans.clear();
m_dzEtaChi2.clear();
m_dzEtaKS.clear();
m_dzEtaHi.clear();
m_dzEtaLo.clear();
m_dxyVect.clear();
m_dzVect.clear();
}
};
#if defined(__ROOTCLING__)
#pragma link C++ class std::map < TString, std::vector < double>> + ;
#pragma link C++ class std::map < TString, std::vector < unrolledHisto>> + ;
#pragma link C++ class outPVtrends + ;
#pragma link C++ class unrolledHisto + ;
#endif
#endif // ALIGNMENT_OFFLINEVALIDATION_OUTPVTRENDS_H_
|