|
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
|
#include "EventFilter/CSCRawToDigi/interface/CSCRPCData.h"
#include "DataFormats/CSCDigi/interface/CSCRPCDigi.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <string>
#include <cstdio>
#include <strings.h> // for bzero
#include <cstring>
/** data format is
RPC0 Tbin 0 Pads[7:0]
RPC0 BXN Pads[15:8]
RPC0 Tbin 1 Pads[7:0]
...
RPC0 Tbin 6 Pads[7:0]
RPC3 BXN Pads[15:8]
*/
#ifdef LOCAL_UNPACK
bool CSCRPCData::debug = false;
#else
std::atomic<bool> CSCRPCData::debug{false};
#endif
CSCRPCData::CSCRPCData(int ntbins) : ntbins_(ntbins), size_(0) {
theData[0] = 0x6b04;
for (int i = 1; i < 257; ++i) {
// data format is bits 12-14 are RPC number, 0 to 3
int rpc = (i - 1) / 14;
theData[i] = rpc << 12;
// bits 8-11 of the first word of the pair is time bin
int tbin = ((i - 1) % 14) / 2;
theData[i] |= tbin << 8;
}
theData[257] = 0x6e04;
}
CSCRPCData::CSCRPCData(const unsigned short* buf, int length) : size_(length) {
//size_ = ntbins_*2*4+2;
// header & trailer word, + 4 RPCs per time bin, 2 lines per RPC
ntbins_ = (size_ - 2) / 8;
memcpy(theData, buf, size_ * 2);
}
void CSCRPCData::Print() const {
LogTrace("CSCRPCData|CSCRawToDigi") << "CSCRPCData.Print";
for (int line = 0; line < ((size_)); ++line) {
LogTrace("CSCRPCData|CSCRawToDigi") << std::hex << theData[line];
}
for (int linePair = 0; linePair < ((size_ - 2) / 2); ++linePair) {
// skip header word
int pos = linePair * 2 + 1;
// make the two pad words into one and see if it's empty
//int pad = theData[pos] & 0xff + ((theData[pos+1] & 0x3f) << 8);
int bxnnew = ((theData[pos + 1] >> 8) & 0x7);
int rpc = (theData[pos] >> 12) & 0x7;
int tbin = (theData[pos] >> 8) & 0xf;
int bxn = bxnnew;
LogTrace("CSCRPCData|CSCRawToDigi") << " RPC=" << rpc << " Tbin=" << tbin << " BXN=" << bxn;
}
}
std::vector<int> CSCRPCData::BXN() const {
std::vector<int> result;
for (int linePair = 0; linePair < ((size_ - 2) / 2); ++linePair) {
// skip header word
int pos = linePair * 2 + 1;
/// make the two pad words into one and see if it's empty
//int pad = theData[pos] & 0xff + ((theData[pos+1] & 0x3f) << 8);
int bxnnew = ((theData[pos + 1] >> 8) & 0x7);
//int bxnnew = (((theData[pos+1] >> 8) & 0x3 )<<2) | ((theData[pos+1]>>6)&0x3) ;
int rpc = (theData[pos] >> 12) & 0x7;
//int tbin = (theData[pos] >> 8) & 0xf;
//int bxn = bxnnew;
result.push_back(bxnnew);
result.push_back(rpc);
}
return result;
}
std::vector<CSCRPCDigi> CSCRPCData::digis() const {
std::vector<CSCRPCDigi> result;
int bxnold = 0;
int bxnnew = 0;
//int bxnewGreg;
for (int linePair = 0; linePair < ((size_ - 2) / 2); ++linePair) {
// skip header word
int pos = linePair * 2 + 1;
// LogTrace("RPC") << "+++ CSCRPCData " << std::hex << theData[pos]
// << " " << theData[pos+1];
if (debug)
LogTrace("CSCRPCData|CSCRawToDigi") << "+++ CSCRPCData " << std::hex << theData[pos] << " " << theData[pos + 1];
// make the two pad words into one and see if it's empty
int pad = (theData[pos] & 0xff) + ((theData[pos + 1] & 0xff) << 8);
//bxnnew = (((theData[pos+1] >> 8) & 0x3 )<<2) | ((theData[pos+1]>>6)&0x3) ;
bxnnew = ((theData[pos + 1] >> 8) & 0x7);
//LogTrace("RPC") << " " << "bxnnew" << " " << bxnnew;
//LogTrace("RPC") << " " << "bxnnewGreg" << " " << bxnewGreg;
if (linePair == 0)
bxnold = bxnnew;
if (bxnnew - bxnold > 1)
LogTrace("CSCRPCData|CSCRawToDigi")
<< "+++ CSCRPCData warning: RPC BXN is incrementing by more than 1 clock cycle";
bxnold = bxnnew;
if (pad != 0) {
if (debug)
LogTrace("CSCRPCData|CSCRawToDigi")
<< "+++ CSCRPCData Found a PAD =" << std::hex << pad << " " << theData[pos] << " + " << theData[pos + 1];
int rpc = (theData[pos] >> 12) & 0x7;
int tbin = (theData[pos] >> 8) & 0xf;
int bxn = bxnnew;
//LogTrace("RPC") << " rpc: " << rpc << " bxn: " << bxn << " tbin: " << tbin;
for (int i = 0; i < 16; ++i) {
// if the bit is set, make a digi
if ((pad >> i) & 1) {
result.push_back(CSCRPCDigi(rpc, i, bxn, tbin));
//LogTrace("RPC") << "digi-->" << " rpc: " << rpc << " i: " << i << " bxn: " << bxn << " tbin: " << tbin;
}
}
}
}
return result;
}
|