Generate and assign the contained variables.
28{
29 if (not ptrFacetRelation) return false;
30
31 const CDCFacet* fromFacet = ptrFacetRelation->first;
32 const CDCFacet* toFacet = ptrFacetRelation->second;
33
34 const UncertainParameterLine2D& fromFitLine = fromFacet->
getFitLine();
35 const UncertainParameterLine2D& toFitLine = toFacet->
getFitLine();
36
39
42
45
47
50
52
53 var<
named(
"from_middle_cos_delta")>() = fromMiddleCos;
54 var<
named(
"to_middle_cos_delta")>() = toMiddleCos;
55
58 {
59 int nSteps = 0;
60 UncertainParameterLine2D fitLine =
FacetFitter::fit(*fromFacet, *toFacet, nSteps);
61 double s = fitLine->
lengthOnCurve(frontWirePos2D, backWirePos2D);
69 }
70
71 {
72 int nSteps = 1;
73 UncertainParameterLine2D fitLine =
FacetFitter::fit(*fromFacet, *toFacet, nSteps);
74 double s = fitLine->
lengthOnCurve(frontWirePos2D, backWirePos2D);
79 }
80
81 {
83 double s = fitLine->
lengthOnCurve(frontWirePos2D, backWirePos2D);
88 }
89
90
91 {
92 using namespace NLineParameterIndices;
93 {
95 if (not std::isfinite(phi0_var)) {
96 B2INFO("from addr " << *fromFacet);
99 }
100 }
101 {
103 if (not std::isfinite(phi0_var)) {
104 B2INFO("to addr " << *toFacet);
107 }
108 }
109
110 var<
named(
"phi0_from_sigma")>() = std::sqrt(fromCov(c_Phi0, c_Phi0));
111 var<
named(
"phi0_to_sigma")>() = std::sqrt(toCov(c_Phi0, c_Phi0));
112 var<
named(
"phi0_ref_sigma")>() = std::sqrt(fromCov(c_Phi0, c_Phi0) + toCov(c_Phi0, c_Phi0));
116 std::sqrt((toCov(c_Phi0, c_Phi0) + fromCov(c_Phi0, c_Phi0))));
117
118 LineParameters avgPar;
119 LineCovariance avgCov;
121
122 LineParameters meanPar = (fromPar + toPar) / 2.0;
124
125 LineParameters relAvgPar = avgPar - meanPar;
127
130 std::fabs(relAvgPar(c_Phi0) / std::sqrt(avgCov(c_Phi0, c_Phi0)));
131 var<
named(
"phi0_comb_diff")>() = relAvgPar(c_Phi0);
132 var<
named(
"phi0_comb_sigma")>() = std::sqrt(avgCov(c_Phi0, c_Phi0));
133 }
134
135
136 {
137 CDCObservations2D observations2D(EFitPos::c_RecoPos, EFitVariance::c_Unit);
138 observations2D.append(*fromFacet);
139 observations2D.append(*toFacet);
140
144 }
145
146 {
147 CDCObservations2D observations2D(EFitPos::c_RecoPos, EFitVariance::c_DriftLength);
148 observations2D.append(*fromFacet);
149 observations2D.append(*toFacet);
150
154 }
155
156 {
157 CDCObservations2D observations2D(EFitPos::c_RecoPos, EFitVariance::c_Pseudo);
158 observations2D.append(*fromFacet);
159 observations2D.append(*toFacet);
160
164 }
165
166 {
167 CDCObservations2D observations2D(EFitPos::c_RecoPos, EFitVariance::c_Proper);
168 observations2D.append(*fromFacet);
169 observations2D.append(*toFacet);
170
174 }
175
176 return true;
177}
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.
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.
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 Vector2D & getRefPos2D() const
The two dimensional reference position (z=0) of the underlying wire.
static double fit(const CDCFacet &facet, int nSteps=100)
Fits a proper line to facet and returns the chi2.
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...