|
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
|
/*
* adapting to CTPPS pixel detector March 2017 - F.Ferro
*/
#include "EventFilter/CTPPSRawToDigi/interface/CTPPSPixelRawToDigi.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DataFormats/CTPPSDigi/interface/CTPPSPixelDigi.h"
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
#include "DataFormats/FEDRawData/interface/FEDRawData.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"
#include "EventFilter/CTPPSRawToDigi/interface/CTPPSPixelDataFormatter.h"
#include "FWCore/Framework/interface/ConsumesCollector.h"
using namespace std;
CTPPSPixelRawToDigi::CTPPSPixelRawToDigi(const edm::ParameterSet& conf)
: config_(conf),
eSummary_("CTPPSPixelDataFormatter", "[ctppsPixelRawToDigi]", edm::isDebugEnabled())
{
FEDRawDataCollection_ = consumes<FEDRawDataCollection>(config_.getParameter<edm::InputTag>("inputLabel"));
CTPPSPixelDAQMapping_ = esConsumes<CTPPSPixelDAQMapping, CTPPSPixelDAQMappingRcd>();
produces<edm::DetSetVector<CTPPSPixelDigi>>();
isRun3_ = config_.getParameter<bool>("isRun3");
includeErrors_ = config_.getParameter<bool>("includeErrors");
mappingLabel_ = config_.getParameter<std::string>("mappingLabel");
if (includeErrors_) {
produces<edm::DetSetVector<CTPPSPixelDataError>>();
}
}
CTPPSPixelRawToDigi::~CTPPSPixelRawToDigi() {
edm::LogInfo("CTPPSPixelRawToDigi") << " CTPPSPixelRawToDigi destructor!";
}
void CTPPSPixelRawToDigi::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
edm::ParameterSetDescription desc;
desc.add<bool>("isRun3", true);
desc.add<bool>("includeErrors", true);
desc.add<edm::InputTag>("inputLabel", edm::InputTag("rawDataCollector"));
desc.add<std::string>("mappingLabel", "RPix");
descriptions.add("ctppsPixelDigis", desc);
}
void CTPPSPixelRawToDigi::produce(edm::Event& ev, const edm::EventSetup& es) {
edm::Handle<FEDRawDataCollection> buffers;
ev.getByToken(FEDRawDataCollection_, buffers);
edm::ESHandle<CTPPSPixelDAQMapping> mapping;
bool data_exist = false;
for (int fed = FEDNumbering::MINCTPPSPixelsFEDID; fed <= FEDNumbering::MAXCTPPSPixelsFEDID; fed++) {
const FEDRawData& tempRawData = buffers->FEDData(fed);
if (tempRawData.size() != 0) {
data_exist = true;
break;
}
}
/// create product (digis & errors)
auto collection = std::make_unique<edm::DetSetVector<CTPPSPixelDigi>>();
auto errorcollection = std::make_unique<edm::DetSetVector<CTPPSPixelDataError>>();
if (data_exist) {
mapping = es.getHandle(CTPPSPixelDAQMapping_);
fedIds_ = mapping->fedIds();
CTPPSPixelDataFormatter formatter(mapping->ROCMapping, eSummary_);
formatter.setErrorStatus(includeErrors_);
bool errorsInEvent = false;
CTPPSPixelDataFormatter::DetErrors nodeterrors;
for (auto aFed = fedIds_.begin(); aFed != fedIds_.end(); ++aFed) {
int fedId = *aFed;
edm::LogInfo("CTPPSPixelRawToDigi") << " PRODUCE DIGI FOR FED: " << dec << fedId << endl;
CTPPSPixelDataFormatter::Errors errors;
/// get event data for this fed
const FEDRawData& fedRawData = buffers->FEDData(fedId);
formatter.interpretRawData(isRun3_, errorsInEvent, fedId, fedRawData, *collection, errors);
if (includeErrors_) {
for (auto const& is : errors) {
uint32_t errordetid = is.first;
/// errors given dummy detId must be sorted by Fed
if (errordetid == RPixErrorChecker::dummyDetId) {
nodeterrors.insert(nodeterrors.end(), errors[errordetid].begin(), errors[errordetid].end());
} else {
edm::DetSet<CTPPSPixelDataError>& errorDetSet = errorcollection->find_or_insert(errordetid);
errorDetSet.data.insert(errorDetSet.data.end(), is.second.begin(), is.second.end());
}
}
}
}
if (includeErrors_) {
errorcollection->find_or_insert(RPixErrorChecker::dummyDetId).data = nodeterrors;
}
if (errorsInEvent)
LogDebug("CTPPSPixelRawToDigi") << "Error words were stored in this event";
}
///send digis and errors back to framework
ev.put(std::move(collection));
if (includeErrors_) {
ev.put(std::move(errorcollection));
}
}
void CTPPSPixelRawToDigi::endStream() { eSummary_.printSummary(); }
DEFINE_FWK_MODULE(CTPPSPixelRawToDigi);
|