|
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
|
#include "CondCore/Utilities/interface/PayloadInspectorModule.h"
#include "CondCore/Utilities/interface/PayloadInspector.h"
#include "CondCore/CondDB/interface/Time.h"
//#include "DataFormats/EcalDetId/interface/EBDetId.h"
//#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "DataFormats/HcalDetId/interface/HcalDetId.h"
//#include "Geometry/HcalCommonData/interface/HcalTopologyMode.h"
#include "Geometry/CaloTopology/interface/HcalTopology.h"
#include "CondCore/HcalPlugins/interface/HcalObjRepresent.h"
// the data format of the condition to be inspected
#include "CondFormats/HcalObjects/interface/HcalChannelQuality.h" //or ChannelStatus.h???
#include "TH2F.h"
#include "TCanvas.h"
#include "TLine.h"
#include "TStyle.h"
#include "TLatex.h"
#include "TPave.h"
#include "TPaveStats.h"
#include <string>
#include <fstream>
namespace {
class HcalChannelStatusContainer : public HcalObjRepresent::HcalDataContainer<HcalChannelQuality, HcalChannelStatus> {
public:
HcalChannelStatusContainer(std::shared_ptr<HcalChannelQuality> payload, unsigned int run)
: HcalObjRepresent::HcalDataContainer<HcalChannelQuality, HcalChannelStatus>(payload, run) {}
float getValue(const HcalChannelStatus* chan) override { return chan->getValue() / 32770; }
};
/******************************************
2d plot of ECAL ChannelStatusRatios of 1 IOV
******************************************/
class HcalChannelQualityPlot : public cond::payloadInspector::PlotImage<HcalChannelQuality> {
public:
HcalChannelQualityPlot()
: cond::payloadInspector::PlotImage<HcalChannelQuality>("HCAL ChannelStatus Ratios - map ") {
setSingleIov(true);
}
bool fill(const std::vector<std::tuple<cond::Time_t, cond::Hash> >& iovs) override {
auto iov = iovs.front();
std::shared_ptr<HcalChannelQuality> payload = fetchPayload(std::get<1>(iov));
if (payload.get()) {
HcalChannelStatusContainer* objContainer = new HcalChannelStatusContainer(payload, std::get<0>(iov));
std::string ImageName(m_imageFileName);
objContainer->getCanvasAll()->SaveAs(ImageName.c_str());
return true;
} else
return false;
} // fill method
};
class HcalChannelQualityChange : public cond::payloadInspector::PlotImage<HcalChannelQuality> {
public:
HcalChannelQualityChange()
: cond::payloadInspector::PlotImage<HcalChannelQuality>("HCAL ChannelStatus Ratios - map ") {
setSingleIov(false);
}
bool fill(const std::vector<std::tuple<cond::Time_t, cond::Hash> >& iovs) override {
auto iov1 = iovs.front();
auto iov2 = iovs.back();
std::shared_ptr<HcalChannelQuality> payload1 = fetchPayload(std::get<1>(iov1));
std::shared_ptr<HcalChannelQuality> payload2 = fetchPayload(std::get<1>(iov2));
if (payload1.get() && payload2.get()) {
HcalChannelStatusContainer* objContainer1 = new HcalChannelStatusContainer(payload1, std::get<0>(iov1));
HcalChannelStatusContainer* objContainer2 = new HcalChannelStatusContainer(payload2, std::get<0>(iov2));
objContainer2->Subtract(objContainer1);
//
// std::map< std::pair< std::string, int >, TH2F* > depths = objContainer1->GetDepths();
//
//
// TODO: How do I display this?
//
//
std::string ImageName(m_imageFileName);
objContainer2->getCanvasAll()->SaveAs(ImageName.c_str());
return true;
} else
return false;
} // fill method
};
} // namespace
// Register the classes as boost python plugin
PAYLOAD_INSPECTOR_MODULE(HcalChannelQuality) {
PAYLOAD_INSPECTOR_CLASS(HcalChannelQualityPlot);
PAYLOAD_INSPECTOR_CLASS(HcalChannelQualityChange);
}
|