|
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
|
#include "DataFormats/FTLRecHit/interface/FTLRecHit.h"
#include <cassert>
#include <cmath>
#include <limits>
namespace {
constexpr float timereso_max = 10000;
}
FTLRecHit::FTLRecHit()
: FTLRecHit(DetId(), 0, 0, -1.f, -1.f, -1.f, -1.f, -1.f, std::numeric_limits<unsigned char>::max()) {}
FTLRecHit::FTLRecHit(const DetId& id,
uint8_t row,
uint8_t column,
float energy,
float time,
float timeError,
float position,
float positionError,
uint32_t flagBits)
: id_(id),
energy_(energy),
time_(time),
timeError_(timeError),
position_(position),
positionError_(positionError),
row_(row),
column_(column),
flagBits_(flagBits) {}
FTLRecHit::FTLRecHit(
const DetId& id, float energy, float time, float timeError, float position, float positionError, uint32_t flagBits)
: FTLRecHit(id, 0, 0, energy, time, timeError, position, positionError, flagBits) {}
bool FTLRecHit::isTimeValid() const {
if (timeError() < 0)
return false;
else
return true;
}
bool FTLRecHit::isTimeErrorValid() const {
if (!isTimeValid())
return false;
if (timeError() >= timereso_max)
return false;
return true;
}
/// check if one of the flags in a set is true
bool FTLRecHit::checkFlags(const std::vector<int>& flagsvec) const {
for (std::vector<int>::const_iterator flagPtr = flagsvec.begin(); flagPtr != flagsvec.end();
++flagPtr) { // check if one of the flags is up
if (checkFlag(*flagPtr))
return true;
}
return false;
}
std::ostream& operator<<(std::ostream& s, const FTLRecHit& hit) {
if (hit.detid().det() == DetId::Forward && hit.detid().subdetId() == FastTime)
return s << MTDDetId(hit.detid()) << ": " << hit.energy() << " GeV, " << hit.time() << " ns";
else
return s << "FTLRecHit undefined subdetector";
}
|