|
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
|
#ifndef TrackingRecHit_DTRecSegment2D_h
#define TrackingRecHit_DTRecSegment2D_h
/** \class DTRecSegment2D
*
* Base class for 2-parameters segments measuring position and direction in X
* projection.
*
* Implements the AbstractDetMeasurement part of the interface
* for 2D RecHits in terms of localPosition() and localPositionError() and
* Direction. This segment is measuring the position and the direction in just
* one projection, the "X". Typical use case is a segment reconstructed only in
* X projection.
* To be used as base class for all 2D positional-directional segments.
* The coordinate measured is assumend to be the local "x" and "dx/dz"
*
* 2D Segments for the muon barrel system.
* 2D means that this segment has information about position and direction in
* one projection (r-phi or r-theta/zeta).
*
* \author Stefano Lacaprara - INFN Legnaro <stefano.lacaprara@pd.infn.it>
* \author Riccardo Bellan - INFN TO <riccardo.bellan@cern.ch>
*
*/
/* Base Class Headers */
#include "DataFormats/TrackingRecHit/interface/RecSegment.h"
/* Collaborating Class Declarations */
#include "DataFormats/GeometrySurface/interface/LocalError.h"
#include "DataFormats/GeometryVector/interface/LocalVector.h"
#include "DataFormats/GeometryVector/interface/LocalPoint.h"
#include "DataFormats/DTRecHit/interface/DTRecHit1D.h"
/* C++ Headers */
#include <iosfwd>
/* Fwd declaration */
class DTSegmentUpdator;
/* ====================================================================== */
/* Class DTRecSegment2D Interface */
class DTRecSegment2D : public RecSegment {
public:
/// Constructor
/// empty c'tor needed by POOL (I guess)
DTRecSegment2D() : theChi2(0.0), theT0(0.), theVdrift(0.) {}
/// c'tor from hits
DTRecSegment2D(DetId id, const std::vector<DTRecHit1D>& hits);
/// complete constructor
DTRecSegment2D(DetId id,
LocalPoint& position,
LocalVector& direction,
AlgebraicSymMatrix& covMatrix,
double chi2,
std::vector<DTRecHit1D>& hits1D);
/// Destructor
~DTRecSegment2D() override;
/* Operations */
DTRecSegment2D* clone() const override { return new DTRecSegment2D(*this); }
/// the vector of parameters (dx/dz,x)
AlgebraicVector parameters() const override { return param(localPosition(), localDirection()); }
// The parameter error matrix
AlgebraicSymMatrix parametersError() const override;
/** return the projection matrix, which must project a parameter vector,
* whose components are (q/p, dx/dz, dy/dz, x, y), into the vector returned
* by parameters() */
AlgebraicMatrix projectionMatrix() const override { return theProjectionMatrix; }
/// return 2. The dimension of the matrix
int dimension() const override { return 2; }
/// local position in SL frame
LocalPoint localPosition() const override { return thePosition; }
/// local position error in SL frame
LocalError localPositionError() const override;
/// the local direction in SL frame
LocalVector localDirection() const override { return theDirection; }
/// the local direction error (xx,xy,yy) in SL frame: only xx is not 0.
LocalError localDirectionError() const override;
/// the chi2 of the fit
double chi2() const override { return theChi2; }
/// return the DOF of the segment
int degreesOfFreedom() const override;
// Access to component RecHits (if any)
std::vector<const TrackingRecHit*> recHits() const override;
// Non-const access to component RecHits (if any)
std::vector<TrackingRecHit*> recHits() override;
/// Access to specific components
std::vector<DTRecHit1D> specificRecHits() const;
/// the Covariance Matrix
AlgebraicSymMatrix covMatrix() const { return theCovMatrix; }
/// Get the segment t0 (if recomputed, 0 is returned otherwise)
double t0() const { return theT0; }
bool ist0Valid() const { return (theT0 > -998.) ? true : false; }
/// Get the vDirft as computed by the algo for the computation of the segment t0
/// (if recomputed, 0 is returned otherwise)
double vDrift() const { return theVdrift; }
protected:
friend class DTSegmentUpdator;
void setPosition(const LocalPoint& pos);
void setDirection(const LocalVector& dir);
void setCovMatrix(const AlgebraicSymMatrix& cov);
void setChi2(const double& chi2);
void update(std::vector<DTRecHit1D>& updatedRecHits);
void setT0(const double& t0);
void setVdrift(const double& vdrift);
LocalPoint thePosition; // in SL frame
LocalVector theDirection; // in SL frame
/// mat[0][0]=sigma (dx/dz)
/// mat[1][1]=sigma (x)
/// mat[0][1]=cov(dx/dz,x)
AlgebraicSymMatrix theCovMatrix; // the covariance matrix
double theChi2; // chi2 of the fit
double theT0; // T0 as coming from the fit
double theVdrift; // vDrift as coming from the fit
std::vector<DTRecHit1D> theHits; // the hits with defined R/L
private:
static const AlgebraicMatrix theProjectionMatrix;
AlgebraicVector param(const LocalPoint& lp, const LocalVector& lv) const {
AlgebraicVector result(2);
result[1] = lp.x();
result[0] = lv.x() / lv.z();
return result;
}
};
std::ostream& operator<<(std::ostream& os, const DTRecSegment2D& seg);
#endif // TrackingRecHit_DTRecSegment2D_h
|