|
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
|
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/ParameterSet/interface/FileInPath.h"
#include "DQM/TrackingMonitor/interface/TrackFoldedOccupancyClient.h"
//-----------------------------------------------------------------------------------
TrackFoldedOccupancyClient::TrackFoldedOccupancyClient(edm::ParameterSet const& iConfig)
//-----------------------------------------------------------------------------------
{
edm::LogInfo("TrackFoldedOccupancyClient") << "TrackFoldedOccupancyClient::Deleting TrackFoldedOccupancyClient ";
TopFolder_ = iConfig.getParameter<std::string>("FolderName");
quality_ = iConfig.getParameter<std::string>("TrackQuality");
algoName_ = iConfig.getParameter<std::string>("AlgoName");
state_ = iConfig.getParameter<std::string>("MeasurementState");
histTag_ = (state_ == "default") ? algoName_ : state_ + "_" + algoName_;
conf_ = iConfig;
}
//-----------------------------------------------------------------------------------
TrackFoldedOccupancyClient::~TrackFoldedOccupancyClient()
//-----------------------------------------------------------------------------------
{
edm::LogInfo("TrackFoldedOccupancyClient") << "TrackFoldedOccupancyClient::Deleting TrackFoldedOccupancyClient ";
}
//-----------------------------------------------------------------------------------
void TrackFoldedOccupancyClient::beginJob(void)
//-----------------------------------------------------------------------------------
{
edm::LogInfo("TrackFoldedOccupancyClient") << "TrackFoldedOccupancyClient::beginJob done";
}
//-----------------------------------------------------------------------------------
void TrackFoldedOccupancyClient::beginRun(edm::Run const& run, edm::EventSetup const& eSetup)
//-----------------------------------------------------------------------------------
{
edm::LogInfo("TrackFoldedOccupancyClient") << "TrackFoldedOccupancyClient:: Begining of Run";
}
//-----------------------------------------------------------------------------------
void TrackFoldedOccupancyClient::bookMEs(DQMStore::IBooker& ibooker)
//-----------------------------------------------------------------------------------
{
ibooker.setCurrentFolder(TopFolder_ + "/" + quality_ + "/GeneralProperties/");
int Phi2DBin = conf_.getParameter<int>("Phi2DBin");
int Eta2DBin = conf_.getParameter<int>("Eta2DBin");
double EtaMin = conf_.getParameter<double>("EtaMin");
double EtaMax = conf_.getParameter<double>("EtaMax");
double PhiMin = conf_.getParameter<double>("PhiMin");
double PhiMax = conf_.getParameter<double>("PhiMax");
// use the AlgoName and Quality Name
std::string histname = "TkEtaPhi_RelativeDifference_byFoldingmap_" + histTag_;
TkEtaPhi_RelativeDifference_byFoldingmap =
ibooker.book2D(histname, histname, Eta2DBin, EtaMin, EtaMax, Phi2DBin, PhiMin, PhiMax);
TkEtaPhi_RelativeDifference_byFoldingmap->setAxisTitle("Track #eta", 1);
TkEtaPhi_RelativeDifference_byFoldingmap->setAxisTitle("Track #phi", 2);
histname = "TkEtaPhi_RelativeDifference_byFoldingmap_op_" + histTag_;
TkEtaPhi_RelativeDifference_byFoldingmap_op =
ibooker.book2D(histname, histname, Eta2DBin, EtaMin, EtaMax, Phi2DBin, PhiMin, PhiMax);
TkEtaPhi_RelativeDifference_byFoldingmap_op->setAxisTitle("Track #eta", 1);
TkEtaPhi_RelativeDifference_byFoldingmap_op->setAxisTitle("Track #phi", 2);
histname = "TkEtaPhi_Ratio_byFoldingmap_" + histTag_;
TkEtaPhi_Ratio_byFoldingmap = ibooker.book2D(histname, histname, Eta2DBin, EtaMin, EtaMax, Phi2DBin, PhiMin, PhiMax);
TkEtaPhi_Ratio_byFoldingmap->setAxisTitle("Track #eta", 1);
TkEtaPhi_Ratio_byFoldingmap->setAxisTitle("Track #phi", 2);
histname = "TkEtaPhi_Ratio_byFoldingmap_op_" + histTag_;
TkEtaPhi_Ratio_byFoldingmap_op =
ibooker.book2D(histname, histname, Eta2DBin, EtaMin, EtaMax, Phi2DBin, PhiMin, PhiMax);
TkEtaPhi_Ratio_byFoldingmap_op->setAxisTitle("Track #eta", 1);
TkEtaPhi_Ratio_byFoldingmap_op->setAxisTitle("Track #phi", 2);
}
//-----------------------------------------------------------------------------------
void TrackFoldedOccupancyClient::dqmEndJob(DQMStore::IBooker& ibooker, DQMStore::IGetter& igetter)
//-----------------------------------------------------------------------------------
{
edm::LogInfo("TrackFoldedOccupancyClient") << "TrackFoldedOccupancyClient::endLuminosityBlock";
bookMEs(ibooker);
std::string inFolder = TopFolder_ + "/" + quality_ + "/GeneralProperties/";
std::string hname;
hname = "TrackEtaPhi_";
MonitorElement* TrackEtaPhi = igetter.get(inFolder + hname + histTag_);
hname = "TrackEtaPhiInverted_";
MonitorElement* TrackEtaPhiInverted = igetter.get(inFolder + hname + histTag_);
hname = "TrackEtaPhiInvertedoutofphase_";
MonitorElement* TrackEtaPhiInvertedoutofphase = igetter.get(inFolder + hname + histTag_);
if (TrackEtaPhi == nullptr || TrackEtaPhiInverted == nullptr || TrackEtaPhiInvertedoutofphase == nullptr) {
edm::LogWarning("TrackFoldedOccupancyClient") << "MEs needed for this module not found. Skipping.";
return;
}
TkEtaPhi_Ratio_byFoldingmap->divide(TrackEtaPhi, TrackEtaPhiInverted, 1., 1., "");
TkEtaPhi_Ratio_byFoldingmap_op->divide(TrackEtaPhi, TrackEtaPhiInvertedoutofphase, 1., 1., "");
int nx = TrackEtaPhi->getNbinsX();
int ny = TrackEtaPhi->getNbinsY();
for (int ii = 1; ii <= nx; ii++) {
for (int jj = 1; jj <= ny; jj++) {
double Sum1 = TrackEtaPhi->getBinContent(ii, jj) + TrackEtaPhiInverted->getBinContent(ii, jj);
double Sum2 = TrackEtaPhi->getBinContent(ii, jj) + TrackEtaPhiInvertedoutofphase->getBinContent(ii, jj);
double Sub1 = TrackEtaPhi->getBinContent(ii, jj) - TrackEtaPhiInverted->getBinContent(ii, jj);
double Sub2 = TrackEtaPhi->getBinContent(ii, jj) - TrackEtaPhiInvertedoutofphase->getBinContent(ii, jj);
if (Sum1 == 0 || Sum2 == 0) {
TkEtaPhi_RelativeDifference_byFoldingmap->setBinContent(ii, jj, 1);
TkEtaPhi_RelativeDifference_byFoldingmap_op->setBinContent(ii, jj, 1);
} else {
double ratio1 = Sub1 / Sum1;
double ratio2 = Sub2 / Sum2;
TkEtaPhi_RelativeDifference_byFoldingmap->setBinContent(ii, jj, ratio1);
TkEtaPhi_RelativeDifference_byFoldingmap_op->setBinContent(ii, jj, ratio2);
}
}
}
}
DEFINE_FWK_MODULE(TrackFoldedOccupancyClient);
|