|
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
|
/*
* See header file for a description of this class.
*
* \author G. Mila - INFN Torino
*/
#include <DQMOffline/Muon/interface/MuonTrackResidualsTest.h>
// Framework
#include <FWCore/Framework/interface/Event.h>
#include "DataFormats/Common/interface/Handle.h"
#include <FWCore/Framework/interface/ESHandle.h>
#include <FWCore/Framework/interface/MakerMacros.h>
#include <FWCore/Framework/interface/EventSetup.h>
#include <FWCore/ParameterSet/interface/ParameterSet.h>
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/Run.h"
#include <iostream>
#include <cstdio>
#include <string>
#include <cmath>
#include "TF1.h"
using namespace edm;
using namespace std;
MuonTrackResidualsTest::MuonTrackResidualsTest(const edm::ParameterSet& ps) {
parameters = ps;
prescaleFactor = parameters.getUntrackedParameter<int>("diagnosticPrescale", 1);
GaussianCriterionName =
parameters.getUntrackedParameter<string>("resDistributionTestName", "ResidualsDistributionGaussianTest");
SigmaCriterionName = parameters.getUntrackedParameter<string>("sigmaTestName", "ResidualsSigmaInRange");
MeanCriterionName = parameters.getUntrackedParameter<string>("meanTestName", "ResidualsMeanInRange");
}
void MuonTrackResidualsTest::dqmEndRun(DQMStore::IBooker& ibooker,
DQMStore::IGetter& igetter,
edm::Run const&,
edm::EventSetup const&) {
////////////////////////////////////////////////////
//// BOOK NEW HISTOGRAMS
////////////////////////////////////////////////////
ibooker.setCurrentFolder("Muons/Tests/trackResidualsTest");
string histName, MeanHistoName, SigmaHistoName, MeanHistoTitle, SigmaHistoTitle;
vector<string> type;
type.push_back("eta");
type.push_back("theta");
type.push_back("phi");
for (unsigned int c = 0; c < type.size(); c++) {
MeanHistoName = "MeanTest_" + type[c];
SigmaHistoName = "SigmaTest_" + type[c];
MeanHistoTitle = "Mean of the #" + type[c] + " residuals distribution";
SigmaHistoTitle = "Sigma of the #" + type[c] + " residuals distribution";
histName = "Res_GlbSta_" + type[c];
histoNames[type[c]].push_back(histName);
histName = "Res_TkGlb_" + type[c];
histoNames[type[c]].push_back(histName);
histName = "Res_TkSta_" + type[c];
histoNames[type[c]].push_back(histName);
MeanHistos[type[c]] = ibooker.book1D(MeanHistoName.c_str(), MeanHistoTitle.c_str(), 3, 0.5, 3.5);
(MeanHistos[type[c]])->setBinLabel(1, "Res_StaGlb", 1);
(MeanHistos[type[c]])->setBinLabel(2, "Res_TkGlb", 1);
(MeanHistos[type[c]])->setBinLabel(3, "Res_TkSta", 1);
SigmaHistos[type[c]] = ibooker.book1D(SigmaHistoName.c_str(), SigmaHistoTitle.c_str(), 3, 0.5, 3.5);
(SigmaHistos[type[c]])->setBinLabel(1, "Res_StaGlb", 1);
(SigmaHistos[type[c]])->setBinLabel(2, "Res_TkGlb", 1);
(SigmaHistos[type[c]])->setBinLabel(3, "Res_TkSta", 1);
}
////////////////////////////////////////////////////
//// OPERATIONS WITH OTHER HISTOGRAMS
////////////////////////////////////////////////////
for (map<string, vector<string> >::const_iterator histo = histoNames.begin(); histo != histoNames.end(); histo++) {
for (unsigned int type = 0; type < (*histo).second.size(); type++) {
string path = "Muons/MuonRecoAnalyzer/" + (*histo).second[type];
MonitorElement* res_histo = igetter.get(path);
if (res_histo) {
// gaussian test
// const QReport *GaussianReport = res_histo->getQReport(GaussianCriterionName);
int BinNumber = type + 1;
float mean = (*res_histo).getMean(1);
float sigma = (*res_histo).getRMS(1);
MeanHistos.find((*histo).first)->second->setBinContent(BinNumber, mean);
SigmaHistos.find((*histo).first)->second->setBinContent(BinNumber, sigma);
}
}
}
}
|