|
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
|
/**\class CSCALCTDigi
*
* Digi for ALCT trigger primitives.
*
*
* \author N. Terentiev, CMU
*/
#include "DataFormats/CSCDigi/interface/CSCALCTDigi.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <iomanip>
#include <iostream>
using namespace std;
namespace {
enum Pattern_Info { NUM_LAYERS = 6, ALCT_PATTERN_WIDTH = 5 };
CSCALCTDigi::WireContainer makeEmptyContainer() {
CSCALCTDigi::WireContainer ret;
ret.resize(NUM_LAYERS);
for (auto& p : ret) {
p.resize(ALCT_PATTERN_WIDTH);
}
return ret;
}
} // namespace
CSCALCTDigi::WireContainer const& CSCALCTDigi::emptyContainer() {
static WireContainer const s_container = makeEmptyContainer();
return s_container;
}
/// Constructors
CSCALCTDigi::CSCALCTDigi(const uint16_t valid,
const uint16_t quality,
const uint16_t accel,
const uint16_t patternb,
const uint16_t keywire,
const uint16_t bx,
const uint16_t trknmb,
const uint16_t hmt,
const Version version)
: valid_(valid),
quality_(quality),
accel_(accel),
patternb_(patternb),
keywire_(keywire),
bx_(bx),
trknmb_(trknmb),
hmt_(hmt),
version_(version) {}
/// Default
CSCALCTDigi::CSCALCTDigi() {
clear(); // set contents to zero
}
/// Clears this ALCT.
void CSCALCTDigi::clear() {
valid_ = 0;
quality_ = 0;
accel_ = 0;
patternb_ = 0;
keywire_ = 0;
bx_ = 0;
trknmb_ = 0;
fullbx_ = 0;
hmt_ = 0;
hits_.clear();
version_ = Version::Legacy;
}
uint16_t CSCALCTDigi::getHMT() const { return (isRun3() ? hmt_ : std::numeric_limits<uint16_t>::max()); }
void CSCALCTDigi::setHMT(const uint16_t h) { hmt_ = isRun3() ? h : std::numeric_limits<uint16_t>::max(); }
void CSCALCTDigi::setRun3(const bool isRun3) { version_ = isRun3 ? Version::Run3 : Version::Legacy; }
bool CSCALCTDigi::operator>(const CSCALCTDigi& rhs) const {
bool returnValue = false;
// Early ALCTs are always preferred to the ones found at later bx's.
if (getBX() < rhs.getBX()) {
returnValue = true;
}
if (getBX() != rhs.getBX()) {
return returnValue;
}
// The > operator then checks the quality of ALCTs.
// If two qualities are equal, the ALCT furthest from the beam axis
// (lowest eta, highest wire group number) is selected.
uint16_t quality1 = getQuality();
uint16_t quality2 = rhs.getQuality();
if (quality1 > quality2) {
returnValue = true;
} else if (quality1 == quality2 && getKeyWG() > rhs.getKeyWG()) {
returnValue = true;
}
return returnValue;
}
bool CSCALCTDigi::operator==(const CSCALCTDigi& rhs) const {
// Exact equality.
bool returnValue = false;
if (isValid() == rhs.isValid() && getQuality() == rhs.getQuality() && getAccelerator() == rhs.getAccelerator() &&
getCollisionB() == rhs.getCollisionB() && getKeyWG() == rhs.getKeyWG() && getBX() == rhs.getBX()) {
returnValue = true;
}
return returnValue;
}
bool CSCALCTDigi::operator!=(const CSCALCTDigi& rhs) const {
// True if == is false.
bool returnValue = true;
if ((*this) == rhs)
returnValue = false;
return returnValue;
}
/// Debug
void CSCALCTDigi::print() const {
if (isValid()) {
edm::LogVerbatim("CSCDigi") << "CSC ALCT #" << setw(1) << getTrknmb() << ": Valid = " << setw(1) << isValid()
<< " Quality = " << setw(2) << getQuality() << " Accel. = " << setw(1)
<< getAccelerator() << " PatternB = " << setw(1) << getCollisionB()
<< " Key wire group = " << setw(3) << getKeyWG() << " BX = " << setw(2) << getBX()
<< " Full BX = " << std::setw(1) << getFullBX();
} else {
edm::LogVerbatim("CSCDigi") << "Not a valid Anode LCT.";
}
}
std::ostream& operator<<(std::ostream& o, const CSCALCTDigi& digi) {
return o << "CSC ALCT #" << digi.getTrknmb() << ": Valid = " << digi.isValid() << " Quality = " << digi.getQuality()
<< " Accel. = " << digi.getAccelerator() << " PatternB = " << digi.getCollisionB()
<< " Key wire group = " << digi.getKeyWG() << " BX = " << digi.getBX();
}
|