45 TMatrixT<Double_t> returnMat(1, 1);
52 returnMat[0][0] = rC[6][0];
60 TMatrixT<Double_t> returnCov(1, 1);
63 returnCov[0][0] = rawCov[6][6];
73 assert(rC.GetNrows() == 7);
75 TVector3 wire1(rC[0][0], rC[1][0], rC[2][0]);
76 TVector3 wire2(rC[3][0], rC[4][0], rC[5][0]);
77 TVector3 wiredirection = wire1 - wire2;
79 TVector3 vaxis = plane.
getV();
80 wiredirection.SetMag(1.);
83 if (fabs(TMath::Abs(wiredirection.Dot(vaxis)) - 1) > 1
e-3) {
85 std::cout <<
"GFWireHitPolicy: plane not valid!!" << std::endl;
93 assert(x.GetNrows() == 7);
94 TVector3 wire1(x[0][0], x[1][0], x[2][0]);
95 TVector3 wire2(x[3][0], x[4][0], x[5][0]);
102 TVector3 poca, poca_onwire, dirInPoca;
107 distance = TMath::Sqrt(fabs(
108 ((wire1 - poca).Mag2() * (wire2 - wire1).Mag2() - pow((wire1 - poca).Dot(wire2 - wire1), 2)) /
109 (wire2 - wire1).Mag2()));
113 GFException exc(
"distance poca-wire > maxdistance", __LINE__, __FILE__);
120 TVector3 fromwiretoextr = poca - poca_onwire;
121 fromwiretoextr.SetMag(1.);
123 TVector3 wiredirection = wire2 - wire1;
124 wiredirection.SetMag(1.);
127 if (fabs(fromwiretoextr * wiredirection) > 1
e-3) {
128 GFException exc(
"fromwiretoextr*wiredirection > 1e-3", __LINE__, __FILE__);
139 TVector3 O = poca_onwire;
virtual void extrapolateToLine(const TVector3 &point1, const TVector3 &point2, TVector3 &poca, TVector3 &normVec, TVector3 &poca_onwire)
This method extrapolates to the point of closest approach to a line.
TVector3 dist(const TVector3 &point) const
const GFDetPlane & detPlane(GFAbsRecoHit *, GFAbsTrackRep *)
Get detector plane Calls GFAbsTrackRep::extrapolateToLine for POCA.
TMatrixT< Double_t > getRawHitCoord() const
Get raw hit coordinates.
Base Class for genfit track representations. Defines interface for track parameterizations.
void checkPlane(GFAbsRecoHit *, const GFDetPlane &)
Check if the detector plane is valid.
TMatrixT< Double_t > hitCoord(GFAbsRecoHit *, const GFDetPlane &)
Hit coordinates in detector plane.
Detector simulation of raw signals on wires.
Exception class for error handling in GENFIT (provides storage for diagnostic information) ...
TMatrixT< Double_t > hitCov(GFAbsRecoHit *, const GFDetPlane &)
Hit covariances in detector plane.
TMatrixT< Double_t > getRawHitCov() const
Get raw hit covariances.
static const std::string fPolicyName