103 cout <<
" InitStatistics() " << endl;
136 cout <<
"InitStatistics(): vectors resized " << endl;
139 main_titles[0] =
"Total Tracker Momentum, #Delta R = 0.3";
140 main_titles[1] =
"Total EM Cal Energy, #Delta R = 0.3";
141 main_titles[2] =
"Total Had Cal Energy, #Delta R = 0.3";
142 main_titles[3] =
"Total HO Cal Energy, #Delta R = 0.3";
143 main_titles[4] =
"Number of Tracker Tracks, #Delta R = 0.3";
144 main_titles[5] =
"Number of Jets around Muon, #Delta R = 0.3";
145 main_titles[6] =
"Tracker p_{T} within veto cone, #Delta R = 0.3";
146 main_titles[7] =
"EM E_{T} within veto cone, #Delta R = 0.3";
147 main_titles[8] =
"Had E_{T} within veto cone, #Delta R = 0.3";
148 main_titles[9] =
"HO E_{T} within veto cone, #Delta R = 0.3";
149 main_titles[10] =
"Average Momentum per Track, #Delta R = 0.3";
150 main_titles[11] =
"Weighted Energy, #Delta R = 0.3";
152 main_titles[12] =
"Total Tracker Momentum, #Delta R = 0.5";
153 main_titles[13] =
"Total EM Cal Energy, #Delta R = 0.5";
154 main_titles[14] =
"Total Had Cal Energy, #Delta R = 0.5";
155 main_titles[15] =
"Total HO Cal Energy, #Delta R = 0.5";
156 main_titles[16] =
"Number of Tracker Tracks, #Delta R = 0.5";
157 main_titles[17] =
"Number of Jets around Muon, #Delta R = 0.5";
158 main_titles[18] =
"Tracker p_{T} within veto cone, #Delta R = 0.5";
159 main_titles[19] =
"EM E_{T} within veto cone, #Delta R = 0.5";
160 main_titles[20] =
"Had E_{T} within veto cone, #Delta R = 0.5";
161 main_titles[21] =
"HO E_{T} within veto cone, #Delta R = 0.5";
162 main_titles[22] =
"Average Momentum per Track, #Delta R = 0.5";
163 main_titles[23] =
"Weighted Energy, #Delta R = 0.5";
165 main_titles[24] =
"Relative Detector-Based Isolation, #Delta R = 0.3";
166 main_titles[25] =
"Relative Detector-Based Isolation, #Delta R = 0.5";
169 main_titles[26] =
"Sum PF Charged Hadron Pt, #Delta R = 0.3";
170 main_titles[27] =
"Sum PF Neutral Hadron Pt, #Delta R = 0.3";
171 main_titles[28] =
"Sum PF Photon Et, #Delta R = 0.3";
172 main_titles[29] =
"Sum PF Neutral Hadron Pt (Higher Pt threshold), #Delta R = 0.3";
173 main_titles[30] =
"Sum PF Photon Et (Higher Pt threshold), #Delta R = 0.3";
174 main_titles[31] =
"Sum PF Charged Particles Pt not from PV (for Pu corrections), #Delta R = 0.3";
177 main_titles[32] =
"Sum PF Charged Hadron Pt, #Delta R = 0.4";
178 main_titles[33] =
"Sum PF Neutral Hadron Pt, #Delta R = 0.4";
179 main_titles[34] =
"Sum PF Photon Et, #Delta R = 0.4";
180 main_titles[35] =
"Sum PF Neutral Hadron Pt (Higher Pt threshold), #Delta R = 0.4";
181 main_titles[36] =
"Sum PF Photon Et (Higher Pt threshold), #Delta R = 0.4";
182 main_titles[37] =
"Sum PF Charged Particles Pt not from PV (for Pu corrections), #Delta R = 0.4";
184 main_titles[38] =
"Relative PF Isolation, #Delta R = 0.3";
185 main_titles[39] =
"Relative PF Isolation, #Delta R = 0.4";
187 main_titles[40] =
"Relative PF Isolation (Higher Pt threshold), #Delta R = 0.3";
188 main_titles[41] =
"Relative PF Isolation (Higher Pt threshold), #Delta R = 0.4";
190 main_titles[42] =
"Sum DR Isolation Profile for Charged Hadron, #Delta R = 0.4";
192 main_titles[43] =
"Sum DR Isolation Profile for Neutral Hadron, #Delta R = 0.4";
194 main_titles[44] =
"Sum DR Isolation Profile for Photon, #Delta R = 0.4";
196 main_titles[45] =
"Mean DR Isolation Profile for Charged Hadron, #Delta R = 0.4";
198 main_titles[46] =
"Mean DR Isolation Profile for Neutral Hadron, #Delta R = 0.4";
200 main_titles[47] =
"Mean DR Isolation Profile for Photon, #Delta R = 0.4";
203 cout <<
"InitStatistics(): main titles 1D DONE " << endl;
205 titles_2D[0] =
"Total Tracker Momentum, #Delta R = 0.3";
206 titles_2D[1] =
"Total EM Cal Energy, #Delta R = 0.3";
207 titles_2D[2] =
"Total Had Cal Energy, #Delta R = 0.3";
208 titles_2D[3] =
"Total HO Cal Energy, #Delta R = 0.3";
209 titles_2D[4] =
"Sum PF Charged Hadron Pt, #Delta R = 0.4";
210 titles_2D[5] =
"Sum PF Neutral Hadron Pt, #Delta R = 0.4";
211 titles_2D[6] =
"Sum PF Photon Et, #Delta R = 0.4";
212 titles_2D[7] =
"Sum PF Charged Pt Not from PV, #Delta R = 0.4";
213 titles_2D[8] =
"Relative Detector-Based Isolation, #Delta R = 0.4";
214 titles_2D[9] =
"Relative PF Isolation, #Delta R = 0.4";
216 main_titles_NVtxs[0] =
"Sum PF Neutral Hadron Pt, #DeltaR = 0.4 ( 20 < N_{Vtx} < 50)";
217 main_titles_NVtxs[1] =
"Sum PF Neutral Hadron Pt, #DeltaR = 0.4 (50 < N_{Vtx} < 80)";
224 cout <<
"InitStatistics(): main titles 2D DONE " << endl;
238 axis_titles[10] =
"#Sigma p_{T} / N_{Tracks} (GeV)";
239 axis_titles[11] =
"(1.5) X #Sigma E_{T}^{EM} + #Sigma E_{T}^{Had}";
251 axis_titles[22] =
"#Sigma p_{T} / N_{Tracks} (GeV)";
252 axis_titles[23] =
"(1.5) X #Sigma E_{T}^{EM} + #Sigma E_{T}^{Had}";
254 axis_titles[24] =
"(#Sigma Tk p_{T} + #Sigma ECAL p_{T} + #Sigma HCAL p_{T})/ Mu p_{T} (GeV)";
255 axis_titles[25] =
"(#Sigma Tk p_{T} + #Sigma ECAL p_{T} + #Sigma HCAL p_{T})/ Mu p_{T} (GeV)";
271 axis_titles[38] =
"(#Sigma PFCharged p_{T} + #Sigma PFNeutral p_{T} + #Sigma PFPhoton p_{T}) Mu p_{T} (GeV)";
272 axis_titles[39] =
"(#Sigma PFCharged p_{T} + #Sigma PFNeutral p_{T} + #Sigma PFPhoton p_{T}) Mu p_{T} (GeV)";
273 axis_titles[40] =
"(#Sigma PFCharged p_{T} + #Sigma PFNeutral p_{T} + #Sigma PFPhoton p_{T}) Mu p_{T} (GeV)";
274 axis_titles[41] =
"(#Sigma PFCharged p_{T} + #Sigma PFNeutral p_{T} + #Sigma PFPhoton p_{T}) Mu p_{T} (GeV)";
292 cout <<
"InitStatistics(): main titles 1D DONE " << endl;
296 names[0] =
"sumPt_R03";
297 names[1] =
"emEt_R03";
298 names[2] =
"hadEt_R03";
299 names[3] =
"hoEt_R03";
300 names[4] =
"nTracks_R03";
301 names[5] =
"nJets_R03";
302 names[6] =
"trackerVetoPt_R03";
303 names[7] =
"emVetoEt_R03";
304 names[8] =
"hadVetoEt_R03";
305 names[9] =
"hoVetoEt_R03";
306 names[10] =
"avgPt_R03";
307 names[11] =
"weightedEt_R03";
309 names[12] =
"sumPt_R05";
310 names[13] =
"emEt_R05";
311 names[14] =
"hadEt_R05";
312 names[15] =
"hoEt_R05";
313 names[16] =
"nTracks_R05";
314 names[17] =
"nJets_R05";
315 names[18] =
"trackerVetoPt_R05";
316 names[19] =
"emVetoEt_R05";
317 names[20] =
"hadVetoEt_R05";
318 names[21] =
"hoVetoEt_R05";
319 names[22] =
"avgPt_R05";
320 names[23] =
"weightedEt_R05";
322 names[24] =
"relDetIso_R03";
323 names[25] =
"relDetIso_R05";
325 names[26] =
"pfChargedPt_R03";
326 names[27] =
"pfNeutralPt_R03";
327 names[28] =
"pfPhotonPt_R03";
328 names[29] =
"pfNeutralPt_HT_R03";
329 names[30] =
"pfPhotonPt_HT_R03";
330 names[31] =
"pfChargedPt_PU_R03";
332 names[32] =
"pfChargedPt_R04";
333 names[33] =
"pfNeutralPt_R04";
334 names[34] =
"pfPhotonPt_R04";
335 names[35] =
"pfNeutralPt_HT_R04";
336 names[36] =
"pfPhotonPt_HT_R04";
337 names[37] =
"pfChargedPt_PU_R04";
339 names[38] =
"relPFIso_R03";
340 names[39] =
"relPFIso_R04";
342 names[40] =
"relPFIso_HT_R03";
343 names[41] =
"relPFIso_HT_R04";
345 names[42] =
"SumDR_PFCharged_R04";
346 names[43] =
"SumDR_PFNeutral_R04";
347 names[44] =
"SumDR_PFPhoton_R04";
349 names[45] =
"MeanDR_PFCharged_R04";
350 names[46] =
"MeanDR_PFNeutral_R04";
351 names[47] =
"MeanDR_PFPhoton_R04";
354 cout <<
"InitStatistics(): names 1D DONE " << endl;
369 cout <<
"InitStatistics(): names 2D DONE " << endl;
390 param[1][0] = (int)(20.0 / S_BIN_WIDTH);
393 param[2][0] = (int)(20.0 / S_BIN_WIDTH);
405 param[6][0] = (int)(40.0 / S_BIN_WIDTH);
411 param[8][0] = (int)(20.0 / S_BIN_WIDTH);
417 param[10][0] = (int)(15.0 / S_BIN_WIDTH);
420 param[11][0] = (int)(20.0 / S_BIN_WIDTH);
424 param[12][0] = (int)(20.0 / S_BIN_WIDTH);
427 param[13][0] = (int)(20.0 / S_BIN_WIDTH);
430 param[14][0] = (int)(20.0 / S_BIN_WIDTH);
442 param[18][0] = (int)(40.0 / S_BIN_WIDTH);
448 param[20][0] = (int)(20.0 / S_BIN_WIDTH);
454 param[22][0] = (int)(15.0 / S_BIN_WIDTH);
457 param[23][0] = (int)(20.0 / S_BIN_WIDTH);
468 param[26][0] = (int)(20.0 / S_BIN_WIDTH);
471 param[27][0] = (int)(20.0 / S_BIN_WIDTH);
474 param[28][0] = (int)(20.0 / S_BIN_WIDTH);
477 param[29][0] = (int)(20.0 / S_BIN_WIDTH);
480 param[30][0] = (int)(20.0 / S_BIN_WIDTH);
483 param[31][0] = (int)(20.0 / S_BIN_WIDTH);
487 param[32][0] = (int)(20.0 / S_BIN_WIDTH);
490 param[33][0] = (int)(20.0 / S_BIN_WIDTH);
493 param[34][0] = (int)(20.0 / S_BIN_WIDTH);
496 param[35][0] = (int)(20.0 / S_BIN_WIDTH);
499 param[36][0] = (int)(20.0 / S_BIN_WIDTH);
502 param[37][0] = (int)(20.0 / S_BIN_WIDTH);
594 cout <<
"InitStatistics(): DONE " << endl;
static const int NUM_VARS_NVTX
std::vector< std::string > main_titles
std::vector< std::string > names
std::vector< std::string > names_2D
static const int NUM_VARS_2D
std::vector< std::string > names_NVtxs
std::vector< int > isContinuous
std::vector< std::string > titles_2D
std::vector< std::string > axis_titles
std::vector< std::string > axis_titles_NVtxs
std::vector< std::vector< double > > param
static const int NUM_VARS
std::vector< std::string > main_titles_NVtxs