|
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
|
#ifndef EventFilter_RPCRawToDigi_debugDigisPrintout_h
#define EventFilter_RPCRawToDigi_debugDigisPrintout_h
#include <algorithm>
#include <string>
#include <sstream>
#include <vector>
#include "DataFormats/RPCDigi/interface/RPCDigiCollection.h"
namespace rpcrawtodigi {
class DebugDigisPrintout {
struct MyDigi {
uint32_t det;
int strip;
int bx;
bool operator==(const MyDigi& o) const { return (det == o.det && strip == o.strip && bx == o.bx); }
bool operator<(const MyDigi& o) const {
if (this->det < o.det)
return true;
if (this->det > o.det)
return false;
if (this->strip < o.strip)
return true;
return false;
}
};
public:
std::string operator()(const RPCDigiCollection* digis) {
std::ostringstream str;
str << "DebugDigisPrintout:";
if (!digis)
return str.str();
typedef DigiContainerIterator<RPCDetId, RPCDigi> DigiRangeIterator;
std::vector<MyDigi> myDigis;
int nDet = 0;
int nDigisAll = 0;
for (DigiRangeIterator it = digis->begin(); it != digis->end(); it++) {
nDet++;
RPCDetId rpcDetId = (*it).first;
uint32_t rawDetId = rpcDetId.rawId();
RPCDigiCollection::Range range = digis->get(rpcDetId);
for (std::vector<RPCDigi>::const_iterator id = range.first; id != range.second; id++) {
nDigisAll++;
const RPCDigi& digi = (*id);
MyDigi myDigi = {rawDetId, digi.strip(), digi.bx()};
if (myDigis.end() == std::find(myDigis.begin(), myDigis.end(), myDigi))
myDigis.push_back(myDigi);
}
}
std::sort(myDigis.begin(), myDigis.end());
str << " dets: " << nDet << " allDigis: " << nDigisAll << " unigueDigis: " << myDigis.size() << std::endl;
for (std::vector<MyDigi>::const_iterator it = myDigis.begin(); it != myDigis.end(); ++it)
str << "debugDIGI: " << it->det << ", " << it->strip << ", " << it->bx << std::endl;
return str.str();
}
};
} // namespace rpcrawtodigi
#endif
|