Generate and assign the contained variables.
29{
30 if (not ptrFacetRelation) return false;
31
32 const CDCFacet* fromFacet = ptrFacetRelation->first;
33 const CDCFacet* toFacet = ptrFacetRelation->second;
34
35 const UncertainParameterLine2D& fromFitLine = fromFacet->
getFitLine();
36 const UncertainParameterLine2D& toFitLine = toFacet->
getFitLine();
37
40
43
46
48
51
53
54 var<
named(
"from_middle_cos_delta")>() = fromMiddleCos;
55 var<
named(
"to_middle_cos_delta")>() = toMiddleCos;
56
59 {
60 int nSteps = 0;
61 UncertainParameterLine2D fitLine =
FacetFitter::fit(*fromFacet, *toFacet, nSteps);
62 double s = fitLine->
lengthOnCurve(frontWirePos2D, backWirePos2D);
70 }
71
72 {
73 int nSteps = 1;
74 UncertainParameterLine2D fitLine =
FacetFitter::fit(*fromFacet, *toFacet, nSteps);
75 double s = fitLine->
lengthOnCurve(frontWirePos2D, backWirePos2D);
80 }
81
82 {
84 double s = fitLine->
lengthOnCurve(frontWirePos2D, backWirePos2D);
89 }
90
91
92 {
93 using namespace NLineParameterIndices;
94 {
96 if (not std::isfinite(phi0_var)) {
97 B2INFO("from addr " << *fromFacet);
100 }
101 }
102 {
104 if (not std::isfinite(phi0_var)) {
105 B2INFO("to addr " << *toFacet);
108 }
109 }
110
111 var<
named(
"phi0_from_sigma")>() = std::sqrt(fromCov(c_Phi0, c_Phi0));
112 var<
named(
"phi0_to_sigma")>() = std::sqrt(toCov(c_Phi0, c_Phi0));
113 var<
named(
"phi0_ref_sigma")>() = std::sqrt(fromCov(c_Phi0, c_Phi0) + toCov(c_Phi0, c_Phi0));
117 std::sqrt((toCov(c_Phi0, c_Phi0) + fromCov(c_Phi0, c_Phi0))));
118
119 LineParameters avgPar;
120 LineCovariance avgCov;
122
123 LineParameters meanPar = (fromPar + toPar) / 2.0;
125
126 LineParameters relAvgPar = avgPar - meanPar;
128
131 std::fabs(relAvgPar(c_Phi0) / std::sqrt(avgCov(c_Phi0, c_Phi0)));
132 var<
named(
"phi0_comb_diff")>() = relAvgPar(c_Phi0);
133 var<
named(
"phi0_comb_sigma")>() = std::sqrt(avgCov(c_Phi0, c_Phi0));
134 }
135
136
137 {
138 CDCObservations2D observations2D(EFitPos::c_RecoPos, EFitVariance::c_Unit);
139 observations2D.append(*fromFacet);
140 observations2D.append(*toFacet);
141
145 }
146
147 {
148 CDCObservations2D observations2D(EFitPos::c_RecoPos, EFitVariance::c_DriftLength);
149 observations2D.append(*fromFacet);
150 observations2D.append(*toFacet);
151
155 }
156
157 {
158 CDCObservations2D observations2D(EFitPos::c_RecoPos, EFitVariance::c_Pseudo);
159 observations2D.append(*fromFacet);
160 observations2D.append(*toFacet);
161
165 }
166
167 {
168 CDCObservations2D observations2D(EFitPos::c_RecoPos, EFitVariance::c_Proper);
169 observations2D.append(*fromFacet);
170 observations2D.append(*toFacet);
171
175 }
176
177 return true;
178}
TrackingUtilities::CDCTrajectory2D fit(const CDCObservations2D &observations2D) const
Fits a collection of observation drift circles.
static const CDCKarimakiFitter & getFitter()
Static getter for a general fitter instance with Karimakis method.
static double fit(const TrackingUtilities::CDCFacet &facet, int nSteps=100)
Fits a proper line to facet and returns the chi2.
ParameterLine2D getStartToEndLine() const
Getter for the tangential line from the first to the third hit.
const UncertainParameterLine2D & getFitLine() const
Getter for the contained line fit information.
ParameterLine2D getStartToMiddleLine() const
Getter for the tangential line from the first to the second hit.
ParameterLine2D getMiddleToEndLine() const
Getter for the tangential line from the second to the third hit.
const CDCWireHit & getEndWireHit() const
Getter for the hit wire of the third oriented wire hit.
const CDCWireHit & getStartWireHit() const
Getter for the hit wire of the first oriented wire hit.
double getChi2() const
Getter for the chi2 value of the circle fit.
double getCurvature() const
Getter for the curvature as seen from the xy projection.
const ROOT::Math::XYVector & getRefPos2D() const
The two dimensional reference position (z=0) of the underlying wire.
const Vector2D & support() const
Gives the support vector of the line.
const Vector2D & tangential() const
Gives the tangential vector of the line.
double lengthOnCurve(const Vector2D &from, const Vector2D &to) const
Denotes the length on the line between the two points.
const LineCovariance & lineCovariance() const
Getter for the whole covariance matrix of the line parameters.
double chi2() const
Getter for the chi square value of the line fit.
LineParameters lineParameters() const
Getter for the three perigee parameters in the order defined by EPerigeeParameter....
static constexpr int named(const char *name)
static Vector2D average(const Vector2D &one, const Vector2D &two)
Constructs the average of two vectors.
double cosWith(const Vector2D &rhs) const
double phi() const
Gives the azimuth angle being the angle to the x axes ( range -M_PI to M_PI )
double angleWith(const Vector2D &rhs) const
The angle between this and rhs.
@ c_Phi0
Constant to address the azimuth angle of the direction of flight.
static double average(const double angle1, double angle2)
Combines two angular values to the one that lies half way between them on the short arc.
static double normalised(const double angle)
Normalise an angle to lie in the range from [-pi, pi].
static double average(const LineUtil::ParameterVector &fromPar, const LineUtil::CovarianceMatrix &fromCov, const LineUtil::ParameterVector &toPar, const LineUtil::CovarianceMatrix &toCov, LineUtil::ParameterVector &avgPar, LineUtil::CovarianceMatrix &avgCov)
Calculates the weighted average between two line parameter sets with their respective covariance matr...
1.13.2