|
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
|
/**
*
* \author G. Franzoni
*
*/
#include <FWCore/Framework/interface/one/EDAnalyzer.h>
#include <DataFormats/Common/interface/Handle.h>
#include <FWCore/Framework/interface/Event.h>
#include <FWCore/Framework/interface/MakerMacros.h>
#include "FWCore/Utilities/interface/InputTag.h"
#include <DataFormats/FEDRawData/interface/FEDRawData.h>
#include <DataFormats/FEDRawData/interface/FEDNumbering.h>
#include <DataFormats/FEDRawData/interface/FEDRawDataCollection.h>
#include <iostream>
#include <vector>
#include <string>
#include <stdio.h>
#include <fstream>
#include <iomanip>
class EcalHexDumperModule : public edm::one::EDAnalyzer<> {
public:
EcalHexDumperModule(const edm::ParameterSet& ps)
: fedRawDataCollectionTag_(ps.getParameter<edm::InputTag>("fedRawDataCollectionTag")) {
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");
}
protected:
int verbosity_;
bool writeDcc_;
int beg_fed_id_;
int end_fed_id_;
int first_event_;
int last_event_;
std::string filename_;
int event_;
void analyze(const edm::Event& e, const edm::EventSetup& c);
private:
edm::InputTag fedRawDataCollectionTag_;
};
void EcalHexDumperModule::analyze(const edm::Event& e, const edm::EventSetup& c) {
event_++;
if (event_ < first_event_ || last_event_ < event_)
return;
edm::Handle<FEDRawDataCollection> rawdata;
e.getByLabel(fedRawDataCollectionTag_, rawdata);
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) {
std::cout << "\n\n\n[EcalHexDumperModule] Event: " << std::dec << event_ << " fed_id: " << id
<< " size_fed: " << data.size() << "\n"
<< std::endl;
if ((data.size() % 16) != 0) {
std::cout << "***********************************************" << std::endl;
std::cout << "Fed size in bits not multiple of 64, strange." << std::endl;
std::cout << "***********************************************" << std::endl;
}
int length = data.size();
const unsigned long* pData = (reinterpret_cast<unsigned long*>(const_cast<unsigned char*>(data.data())));
std::cout << std::setfill('0');
for (int words = 0; words < length / 4; (words += 2)) {
std::cout << std::setw(8) << std::hex << pData[words + 1] << " ";
std::cout << std::setw(8) << std::hex << pData[words] << std::endl;
}
std::cout << "\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());
}
DEFINE_FWK_MODULE(EcalHexDumperModule);
|