|
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
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
|
#include "DQM/L1TMonitor/interface/L1TCaloLayer1Summary.h"
L1TCaloLayer1Summary::L1TCaloLayer1Summary(const edm::ParameterSet& iConfig)
: caloLayer1CICADAScoreToken_(
consumes<l1t::CICADABxCollection>(iConfig.getParameter<edm::InputTag>("caloLayer1CICADAScore"))),
gtCICADAScoreToken_(consumes<l1t::CICADABxCollection>(iConfig.getParameter<edm::InputTag>("gtCICADAScore"))),
simCICADAScoreToken_(consumes<l1t::CICADABxCollection>(iConfig.getParameter<edm::InputTag>("simCICADAScore"))),
caloLayer1RegionsToken_(
consumes<L1CaloRegionCollection>(iConfig.getParameter<edm::InputTag>("caloLayer1Regions"))),
simRegionsToken_(consumes<L1CaloRegionCollection>(iConfig.getParameter<edm::InputTag>("simRegions"))),
fedRawData_(consumes<FEDRawDataCollection>(iConfig.getParameter<edm::InputTag>("fedRawDataLabel"))),
histFolder_(iConfig.getParameter<std::string>("histFolder")) {}
// ------------ method called for each event ------------
void L1TCaloLayer1Summary::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
edm::Handle<FEDRawDataCollection> fedRawDataCollection;
iEvent.getByToken(fedRawData_, fedRawDataCollection);
if (fedRawDataCollection.isValid()) {
for (int iFed = FEDNumbering::MINRCTFEDID + 4; iFed < FEDNumbering::MAXRCTFEDID; iFed += 2) {
const FEDRawData& fedRawData = fedRawDataCollection->FEDData(iFed);
if (fedRawData.size() == 0) {
continue;
}
const uint64_t* fedRawDataArray = (const uint64_t*)fedRawData.data();
UCTDAQRawData daqData(fedRawDataArray);
if (daqData.nAMCs() == 7) {
UCTAMCRawData amcSlot7(daqData.amcPayload(3));
if (amcSlot7.amcNo() == 7) {
histoSlot7MinusDaqBxid->Fill(amcSlot7.BXID() - daqData.BXID());
}
}
}
}
L1CaloRegionCollection caloLayer1Regions = iEvent.get(caloLayer1RegionsToken_);
L1CaloRegionCollection simRegions = iEvent.get(simRegionsToken_);
int nRegions = caloLayer1Regions.size();
unsigned int maxEtaIdx = 0;
for (int iRegion = 0; iRegion < nRegions; iRegion++) {
if (maxEtaIdx < caloLayer1Regions[iRegion].gctEta()) {
maxEtaIdx = caloLayer1Regions[iRegion].gctEta();
}
}
int matrixSize = maxEtaIdx + 1;
bool foundMatrix[2][matrixSize][matrixSize];
int etMatrix[2][matrixSize][matrixSize];
for (int i = 0; i < 2; i++) {
for (int j = 0; j < matrixSize; j++) {
for (int k = 0; k < matrixSize; k++) {
foundMatrix[i][j][k] = false;
etMatrix[i][j][k] = 0;
}
}
}
for (int iRegion = 0; iRegion < nRegions; iRegion++) {
L1CaloRegion cRegion = caloLayer1Regions[iRegion];
L1CaloRegion sRegion = simRegions[iRegion];
foundMatrix[0][cRegion.gctEta()][cRegion.gctPhi()] = true;
etMatrix[0][cRegion.gctEta()][cRegion.gctPhi()] = cRegion.et();
foundMatrix[1][sRegion.gctEta()][sRegion.gctPhi()] = true;
etMatrix[1][sRegion.gctEta()][sRegion.gctPhi()] = sRegion.et();
}
int iRegion = 0;
for (int iEta = 0; iEta < matrixSize; iEta++) {
for (int iPhi = 0; iPhi < matrixSize; iPhi++) {
if (foundMatrix[0][iEta][iPhi] && foundMatrix[1][iEta][iPhi]) {
histoCaloRegions->Fill(iRegion, etMatrix[0][iEta][iPhi]);
histoSimRegions->Fill(iRegion, etMatrix[1][iEta][iPhi]);
histoCaloMinusSimRegions->Fill(iRegion, etMatrix[0][iEta][iPhi] - etMatrix[1][iEta][iPhi]);
iRegion++;
}
}
}
auto caloCICADAScores = iEvent.get(caloLayer1CICADAScoreToken_);
const auto& gtCICADAScores = iEvent.get(gtCICADAScoreToken_);
auto simCICADAScores = iEvent.get(simCICADAScoreToken_);
if (caloCICADAScores.size() > 0) {
histoCaloLayer1CICADAScore->Fill(caloCICADAScores[0]);
if (gtCICADAScores.size() > 0) {
histoGtCICADAScore->Fill(gtCICADAScores.at(0, 0));
histoCaloMinusGt->Fill(caloCICADAScores[0] - gtCICADAScores.at(0, 0));
}
if (simCICADAScores.size() > 0) {
histoSimCICADAScore->Fill(simCICADAScores[0]);
histoCaloMinusSim->Fill(caloCICADAScores[0] - simCICADAScores[0]);
}
}
}
void L1TCaloLayer1Summary::bookHistograms(DQMStore::IBooker& ibooker, edm::Run const&, edm::EventSetup const&) {
ibooker.setCurrentFolder(histFolder_);
histoSlot7MinusDaqBxid = ibooker.book1D("slot7BXID", "Slot 7- DAQ BXID", 50, -20, 20);
ibooker.setCurrentFolder(histFolder_ + "/CICADAScore");
histoCaloLayer1CICADAScore = ibooker.book1D("caloLayer1CICADAScore", "CaloLayer1 CICADAScore", 50, 0, 200);
histoGtCICADAScore = ibooker.book1D("gtCICADAScore", "GT CICADAScore at BX0", 50, 0, 200);
histoCaloMinusGt = ibooker.book1D("caloMinusGtCICADAScore", "CaloLayer1 - GT CICADAScore at BX0", 50, -50, 50);
histoSimCICADAScore =
ibooker.book1D("simCaloLayer1CICADAScore", "simCaloLayer1 CICADAScore (input: DAQ regions)", 50, 0, 200);
histoCaloMinusSim = ibooker.book1D(
"caloMinusSimCICADAScore", "CaloLayer1 - simCaloLayer1 (input: DAQ regions) CICADAScore", 50, -50, 50);
ibooker.setCurrentFolder(histFolder_ + "/Regions");
histoCaloMinusSimRegions =
ibooker.book2D("caloMinusSumRegions",
"CaloLayer1 - simCaloLayer1 (input: DAQ trigger primatives) Regions;Region;ET Difference",
252,
-0.5,
252.5,
100,
-400,
400);
histoCaloRegions = ibooker.book2D("caloLayer1Regions", "CaloLayer1 Regions;Region;ET", 252, -0.5, 252.5, 100, 0, 800);
histoSimRegions = ibooker.book2D("simCaloLayer1Regions",
"simCaloLayer1 Regions (input: DAQ trigger primatives);Region;ET",
252,
-0.5,
252.5,
100,
0,
800);
}
// ------------ method fills 'descriptions' with the allowed parameters for the module ------------
void L1TCaloLayer1Summary::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
// l1tCaloLayer1Summary
edm::ParameterSetDescription desc;
desc.add<edm::InputTag>("caloLayer1CICADAScore", edm::InputTag("caloLayer1Digis", "CICADAScore"));
desc.add<edm::InputTag>("gtCICADAScore", edm::InputTag("gtTestcrateStage2Digis", "CICADAScore"));
desc.add<edm::InputTag>("simCICADAScore", edm::InputTag("simCaloStage2Layer1Summary", "CICADAScore"));
desc.add<edm::InputTag>("caloLayer1Regions", edm::InputTag("caloLayer1Digis"));
desc.add<edm::InputTag>("simRegions", edm::InputTag("simCaloStage2Layer1Digis"));
desc.add<edm::InputTag>("fedRawDataLabel", edm::InputTag("rawDataCollector"));
desc.add<std::string>("histFolder", "L1T/L1TCaloLayer1Summary");
descriptions.add("l1tCaloLayer1Summary", desc);
}
|