|
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
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
|
#ifndef BaseTrackerRecHit_H
#define BaseTrackerRecHit_H
#include "DataFormats/TrackerRecHit2D/interface/trackerHitRTTI.h"
#include "DataFormats/TrackingRecHit/interface/TrackingRecHitGlobalState.h"
#include "DataFormats/GeometryCommonDetAlgo/interface/ErrorFrameTransformer.h"
#include "Geometry/CommonDetUnit/interface/TrackerGeomDet.h"
#include "DataFormats/GeometrySurface/interface/Surface.h"
//#define DO_INTERNAL_CHECKS_BTR
//#define VI_DEBUG
class OmniClusterRef;
class BaseTrackerRecHit : public TrackingRecHit {
public:
BaseTrackerRecHit() : qualWord_(0) {}
// fake TTRH interface
BaseTrackerRecHit const* hit() const final { return this; }
~BaseTrackerRecHit() override {}
// no position (as in persistent)
BaseTrackerRecHit(DetId id, trackerHitRTTI::RTTI rt) : TrackingRecHit(id, (unsigned int)(rt)), qualWord_(0) {}
BaseTrackerRecHit(const GeomDet& idet, trackerHitRTTI::RTTI rt)
: TrackingRecHit(idet, (unsigned int)(rt)), qualWord_(0) {}
BaseTrackerRecHit(const LocalPoint& p, const LocalError& e, GeomDet const& idet, trackerHitRTTI::RTTI rt)
: TrackingRecHit(idet, (unsigned int)(rt)), pos_(p), err_(e), qualWord_(0) {
auto trackerDet = static_cast<TrackerGeomDet const*>(det());
LocalError lape = trackerDet->localAlignmentError();
if (lape.valid())
err_ = LocalError(err_.xx() + lape.xx(), err_.xy() + lape.xy(), err_.yy() + lape.yy());
}
trackerHitRTTI::RTTI rtti() const { return trackerHitRTTI::rtti(*this); }
bool isSingle() const { return trackerHitRTTI::isSingle(*this); }
bool isMatched() const { return trackerHitRTTI::isMatched(*this); }
bool isProjected() const { return trackerHitRTTI::isProjected(*this); }
bool isProjMono() const { return trackerHitRTTI::isProjMono(*this); }
bool isProjStereo() const { return trackerHitRTTI::isProjStereo(*this); }
bool isMulti() const { return trackerHitRTTI::isMulti(*this); }
virtual bool isPixel() const { return false; }
virtual bool isPhase2() const { return false; }
// used by trackMerger (to be improved)
virtual OmniClusterRef const& firstClusterRef() const = 0;
// verify that hits can share clusters...
inline bool sameDetModule(TrackingRecHit const& hit) const;
bool hasPositionAndError() const override;
LocalPoint localPosition() const override {
check();
return pos_;
}
LocalError localPositionError() const override {
check();
return err_;
}
const LocalPoint& localPositionFast() const {
check();
return pos_;
}
const LocalError& localPositionErrorFast() const {
check();
return err_;
}
// to be specialized for 1D and 2D
void getKfComponents(KfComponentsHolder& holder) const override = 0;
int dimension() const override = 0;
void getKfComponents1D(KfComponentsHolder& holder) const;
void getKfComponents2D(KfComponentsHolder& holder) const;
// global coordinates
// Extension of the TrackingRecHit interface
const Surface* surface() const final { return &(det()->surface()); }
GlobalPoint globalPosition() const final { return surface()->toGlobal(localPosition()); }
GlobalError globalPositionError() const final {
return ErrorFrameTransformer().transform(localPositionError(), *surface());
}
float errorGlobalR() const final { return std::sqrt(globalPositionError().rerr(globalPosition())); }
float errorGlobalZ() const final { return std::sqrt(globalPositionError().czz()); }
float errorGlobalRPhi() const final {
return globalPosition().perp() * sqrt(globalPositionError().phierr(globalPosition()));
}
// once cache removed will obsolete the above
TrackingRecHitGlobalState globalState() const {
GlobalError globalError = ErrorFrameTransformer::transform(localPositionError(), *surface());
GlobalPoint gp = globalPosition();
float r = gp.perp();
float errorRPhi = r * std::sqrt(float(globalError.phierr(gp)));
float errorR = std::sqrt(float(globalError.rerr(gp)));
float errorZ = std::sqrt(float(globalError.czz()));
return (TrackingRecHitGlobalState){gp.basicVector(), r, gp.barePhi(), errorR, errorZ, errorRPhi};
}
/// cluster probability, overloaded by pixel rechits.
virtual float clusterProbability() const { return 1.f; }
public:
// obsolete (for what tracker is concerned...) interface
AlgebraicVector parameters() const override;
AlgebraicSymMatrix parametersError() const override;
AlgebraicMatrix projectionMatrix() const override;
private:
#ifdef VI_DEBUG
void check() const { assert(det()); }
#elif defined(DO_INTERNAL_CHECKS_BTR)
void check() const;
#else
static void check() {}
#endif
protected:
LocalPoint pos_;
LocalError err_;
protected:
//this comes for free (padding)
unsigned int qualWord_;
};
bool BaseTrackerRecHit::sameDetModule(TrackingRecHit const& hit) const {
unsigned int myid = geographicalId().rawId();
unsigned int mysubd = myid >> (DetId::kSubdetOffset);
unsigned int id = hit.geographicalId().rawId();
unsigned int subd = id >> (DetId::kSubdetOffset);
if (mysubd != subd)
return false;
//Protection against invalid hits
if (!hit.isValid())
return false;
const unsigned int limdet = 10; // TIB=11
if (mysubd > limdet) { // strip
// mask glue and stereo
myid |= 3;
id |= 3;
}
return id == myid;
}
// Comparison operators
inline bool operator<(const BaseTrackerRecHit& one, const BaseTrackerRecHit& other) {
return (one.geographicalId() < other.geographicalId());
}
#endif // BaseTrackerRecHit_H
|