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
|
/**
*
* \author G. Franzoni
*
*/
#include "CaloOnlineTools/EcalTools/plugins/EcalHexDisplay.h"
EcalHexDisplay::EcalHexDisplay(const edm::ParameterSet& ps)
: verbosity_(ps.getUntrackedParameter<int>("verbosity", 1)),
beg_fed_id_(ps.getUntrackedParameter<int>("beg_fed_id", 0)),
end_fed_id_(ps.getUntrackedParameter<int>("end_fed_id", 654)),
first_event_(ps.getUntrackedParameter<int>("first_event", 1)),
last_event_(ps.getUntrackedParameter<int>("last_event", 9999999)),
event_(0),
writeDcc_(ps.getUntrackedParameter<bool>("writeDCC", false)),
filename_(ps.getUntrackedParameter<std::string>("filename", "dump.bin")),
fedRawDataCollectionToken_(
consumes<FEDRawDataCollection>(ps.getParameter<edm::InputTag>("fedRawDataCollectionTag"))) {}
void EcalHexDisplay::analyze(const edm::Event& e, const edm::EventSetup& c) {
event_++;
if (event_ < first_event_ || last_event_ < event_)
return;
const edm::Handle<FEDRawDataCollection>& rawdata = e.getHandle(fedRawDataCollectionToken_);
std::ofstream dumpFile(filename_.c_str(), std::ios::app);
for (int id = 0; id <= FEDNumbering::MAXFEDID; ++id) {
if (id < beg_fed_id_ || end_fed_id_ < id)
continue;
const FEDRawData& data = rawdata->FEDData(id);
if (data.size() > 4) {
edm::LogVerbatim("EcalTools") << "\n\n\n[EcalHexDumperModule] Event: " << std::dec << event_ << " fed_id: " << id
<< " size_fed: " << data.size() << "\n";
if ((data.size() % 16) != 0) {
edm::LogVerbatim("EcalTools") << "***********************************************";
edm::LogVerbatim("EcalTools") << "Fed size in bits not multiple of 64, strange.";
edm::LogVerbatim("EcalTools") << "***********************************************";
}
int length = data.size();
const unsigned long* pData = (reinterpret_cast<unsigned long*>(const_cast<unsigned char*>(data.data())));
edm::LogVerbatim("EcalTools") << std::setfill('0');
for (int words = 0; words < length / 4; (words += 2)) {
edm::LogVerbatim("EcalTools") << std::setw(8) << std::hex << pData[words + 1] << " " << std::setw(8) << std::hex
<< pData[words];
}
edm::LogVerbatim("EcalTools") << "\n";
if (beg_fed_id_ <= id && id <= end_fed_id_ && writeDcc_) {
dumpFile.write(reinterpret_cast<const char*>(pData), length);
}
}
}
dumpFile.close();
if (!writeDcc_)
remove(filename_.c_str());
}
|