Belle II Software development
BKLMTrackFitter.cc
1/**************************************************************************
2 * basf2 (Belle II Analysis Software Framework) *
3 * Author: The Belle II Collaboration *
4 * *
5 * See git log for contributors and copyright holders. *
6 * This file is licensed under LGPL-3.0, see LICENSE.md. *
7 **************************************************************************/
8
9/* Own header. */
10#include <klm/bklm/modules/bklmTracking/BKLMTrackFitter.h>
11
12/* KLM headers. */
13#include <klm/bklm/geometry/GeometryPar.h>
14#include <klm/bklm/geometry/Module.h>
15
16/* Basf2 headers. */
17#include <framework/logging/Logger.h>
18
19/* CLHEP headers. */
20#include <CLHEP/Matrix/DiagMatrix.h>
21#include <CLHEP/Matrix/Matrix.h>
22
23/* C++ headers. */
24#include <cfloat>
25
26using namespace CLHEP;
27using namespace Belle2;
28using namespace Belle2::bklm;
29
31enum { VX = 0, VY = 1, VZ = 2 };
32
34enum { AY = 0, BY = 1, AZ = 2, BZ = 3 };
35
37enum { MY = 0, MZ = 1 };
38
40BKLMTrackFitter::BKLMTrackFitter():
41 m_Valid(false),
42 m_Good(false),
43 m_Chi2(0.0),
44 m_NumHit(0),
45 m_globalFit(false),
46 m_SectorPar(4, 0),
47 m_SectorErr(4, 0),
48 m_GlobalPar(4, 0),
49 m_GlobalErr(4, 0),
50 m_GeoPar(nullptr)
51{
52}
53
55BKLMTrackFitter::~BKLMTrackFitter()
56{
57}
58
60double BKLMTrackFitter::fit(std::list<KLMHit2d* >& listHitSector)
61{
62
63 // We can only do fit if there are at least two hits
64 if (listHitSector.size() < 2) {
65 m_Valid = false;
66 m_Good = false;
67 return (false);
68 }
69
70 HepVector eta(2, 0); // ( a, b ) in y = a + bx
71 HepSymMatrix error(2, 0);
72 HepVector gloEta(2, 0); // ( a, b ) in y = a + bx in global system
73 HepSymMatrix gloError(2, 0);
74 m_Chi2 = 0;
75
76 // Create temporary vector and matrix... so size can be set.
77 HepVector sectorPar(4, 0);
78 HepSymMatrix sectorErr(4, 0);
79 HepVector globalPar(4, 0);
80 HepSymMatrix globalErr(4, 0);
81
82 if (m_globalFit) {
83 m_Chi2 = fit1dTrack(listHitSector, eta, error, VY, VX);
84 globalPar.sub(1, eta);
85 globalErr.sub(1, error);
86 } else {
87 m_Chi2 = fit1dSectorTrack(listHitSector, eta, error, VY, VX);
88 sectorPar.sub(1, eta);
89 sectorErr.sub(1, error);
90 }
91
92 if (m_globalFit) {
93 m_Chi2 += fit1dTrack(listHitSector, eta, error, VY, VZ);
94 globalPar.sub(3, eta);
95 globalErr.sub(3, error);
96 } else {
97 m_Chi2 += fit1dSectorTrack(listHitSector, eta, error, VZ, VX);
98 sectorPar.sub(3, eta);
99 sectorErr.sub(3, error);
100 }
101
102 if (!m_globalFit) {
103 //transfer to the global system, choose two arbitrary points on track within in the sector jpionts on this track
104 const Belle2::bklm::Module* refMod = m_GeoPar->findModule((*listHitSector.begin())->getSection(),
105 (*listHitSector.begin())->getSector(), 1);
106
107 Hep3Vector p1(0, 0, 0); Hep3Vector p2(0, 0, 0);
108 double x1 = 5; // sector localX
109 double x2 = 10;
110 double y1 = sectorPar[0] + sectorPar[1] * x1;
111 double y2 = sectorPar[0] + sectorPar[1] * x2;
112 double z1 = sectorPar[2] + sectorPar[3] * x1;
113 double z2 = sectorPar[2] + sectorPar[3] * x2;
114 p1.setX(x1); p1.setY(y1); p1.setZ(z1);
115 p2.setX(x2); p2.setY(y2); p2.setZ(z2);
116 Hep3Vector gl1 = refMod->localToGlobal(p1);
117 Hep3Vector gl2 = refMod->localToGlobal(p2);
118
119 //transfer the trackParameters to global system y = p0 + p1 * x and z= p2 + p3*x
120 if (gl2[0] != gl1[0]) {
121 globalPar[1] = (gl2[1] - gl1[1]) / (gl2[0] - gl1[0]);
122 globalPar[0] = gl1[1] - globalPar[1] * gl1[0];
123 globalPar[3] = (gl2[2] - gl1[2]) / (gl2[0] - gl1[0]);
124 globalPar[2] = gl1[2] - globalPar[3] * gl1[0];
125 globalErr = sectorErr;
126 } else {
127 globalPar[1] = DBL_MAX;
128 globalPar[0] = DBL_MAX;
129 globalPar[3] = DBL_MAX;
130 globalPar[2] = DBL_MAX;
131 globalErr = sectorErr; //?
132 }
133 } else { //transfer to the local system, can not do. One global track might go through two different sectors.
134 }
135
136 m_Chi2 /= 2.0;
137
138 m_SectorPar = sectorPar;
139 m_SectorErr = sectorErr;
140 m_GlobalPar = globalPar;
141 m_GlobalErr = globalErr;
142
143 m_Valid = true;
144 m_Good = false;
145 m_NumHit = listHitSector.size();
146 if (m_Chi2 <= 20.0) {
147 m_Good = true;
148 }
149
150 return (m_Chi2);
151
152}
153
155double BKLMTrackFitter::distanceToHit(KLMHit2d* hit,
156 double& error,
157 double& sigma)
158{
159
160 double x, y, z, dy, dz;
161
162 if (!m_Valid) {
163 error = DBL_MAX;
164 sigma = DBL_MAX;
165 return DBL_MAX;
166 }
167
168 m_GeoPar = GeometryPar::instance();
169 const Belle2::bklm::Module* refMod = m_GeoPar->findModule(hit->getSection(), hit->getSector(), 1);
170 const Belle2::bklm::Module* corMod = m_GeoPar->findModule(hit->getSection(), hit->getSector(), hit->getLayer());
171
172 CLHEP::Hep3Vector globalPos(hit->getPositionX(), hit->getPositionY(), hit->getPositionZ());
173
174 //+++ local coordinates of the hit
175 CLHEP::Hep3Vector local = refMod->globalToLocal(globalPos);
176
177 x = local[0] ;
178
179 y = m_SectorPar[ AY ] + x * m_SectorPar[ BY ];
180 z = m_SectorPar[ AZ ] + x * m_SectorPar[ BZ ];
181
182 dy = y - local[1];
183 dz = z - local[2];
184
185 double distance = sqrt(dy * dy + dz * dz);
186
187 // Error is composed of four parts: error due to tracking, y and z;
188 // and error in hit, y and z. We know the latter two, got to find
189 // the first two. We could calculate this from simple equations or
190 // using matrices. I choose the later because it is extendable.
191 HepMatrix errors(2, 2, 0); // Matrix for errors
192 HepMatrix A(2, 4, 0); // Matrix for derivatives
193
194 // Derivatives of y (z) = a + bx with respect to a and b.
195 A[ MY ][ AY ] = 1.0;
196 A[ MY ][ BY ] = x;
197 A[ MZ ][ AZ ] = 1.0;
198 A[ MZ ][ BZ ] = x;
199
200 errors = A * m_SectorErr * A.T();
201
202 double hit_localPhiErr = corMod->getPhiStripWidth() / sqrt(12);
203 double hit_localZErr = corMod->getZStripWidth() / sqrt(12);
204
205 if (hit->inRPC()) {
206 //+++ scale localErr based on measured-in-Belle resolution
207 int nStrips = hit->getPhiStripMax() - hit->getPhiStripMin() + 1;
208 double dn = nStrips - 1.5;
209 double factor = std::pow((0.9 + 0.4 * dn * dn), 1.5) * 0.60;//measured-in-Belle resolution
210 hit_localPhiErr = hit_localPhiErr * sqrt(factor);
211
212 nStrips = hit->getZStripMax() - hit->getZStripMin() + 1;
213 dn = nStrips - 1.5;
214 factor = std::pow((0.9 + 0.4 * dn * dn), 1.5) * 0.55;//measured-in-Belle resolution
215 hit_localZErr = hit_localZErr * sqrt(factor);
216 }
217
218 error = sqrt(errors[ MY ][ MY ] +
219 errors[ MZ ][ MZ ] +
220 pow(hit_localPhiErr, 2) +
221 pow(hit_localZErr, 2));
222
223 if (error != 0.0) {
224 sigma = distance / error;
225 } else {
226 sigma = DBL_MAX;
227 }
228
229 return (distance);
230}
231
233double BKLMTrackFitter::globalDistanceToHit(KLMHit2d* hit,
234 double& error,
235 double& sigma)
236{
237
238 if (!m_Valid) {
239 error = DBL_MAX;
240 sigma = DBL_MAX;
241 return DBL_MAX;
242 }
243
244 //in global fit, we have two functions y = a + b*x and y = c + d*z
245 double x_mea = hit->getPositionX();
246 double y_mea = hit->getPositionY();
247 double z_mea = hit->getPositionZ();
248
249 double x_pre = (y_mea - m_GlobalPar[ AY ]) / m_GlobalPar[ BY ]; //y_mea has uncertainties actually
250 double z_pre = (y_mea - m_GlobalPar[ AZ ]) / m_GlobalPar[ BZ ];
251
252 double dx = x_pre - x_mea;
253 double dz = z_pre - z_mea;
254
255 double distance = sqrt(dx * dx + dz * dz);
256
257 // Error is composed of four parts: error due to tracking;
258 // and error in hit, y, x or y, z.
259 HepMatrix errors(2, 2, 0); // Matrix for errors
260 HepMatrix A(2, 4, 0); // Matrix for derivatives
261
262 // Derivatives of x (z) = y/b - a/b with respect to a and b.
263 A[ MY ][ AY ] = -1. / m_GlobalPar[BY];
264 A[ MY ][ BY ] = -1 * (y_mea - m_GlobalPar[ AY ]) / (m_GlobalPar[ BY ] * m_GlobalPar[ BY ]);
265 A[ MZ ][ AZ ] = -1. / m_GlobalPar[ BZ ];
266 A[ MZ ][ BZ ] = -1 * (y_mea - m_GlobalPar[ AZ ]) / (m_GlobalPar[ BZ ] * m_GlobalPar[ BZ ]);
267
268 //error from trackPar is inclueded, error from y_mea is not included
269 errors = A * m_GlobalErr * A.T();
270
271 //here get the resolustion of a hit, repeated several times, ugly. should we store this in KLMHit2d object ?
272 const Belle2::bklm::Module* corMod = m_GeoPar->findModule(hit->getSection(), hit->getSector(), hit->getLayer());
273 double hit_localPhiErr = corMod->getPhiStripWidth() / sqrt(12);
274 double hit_localZErr = corMod->getZStripWidth() / sqrt(12);
275
276 if (hit->inRPC()) {
277 //+++ scale localErr based on measured-in-Belle resolution
278 int nStrips = hit->getPhiStripMax() - hit->getPhiStripMin() + 1;
279 double dn = nStrips - 1.5;
280 double factor = std::pow((0.9 + 0.4 * dn * dn), 1.5) * 0.60;//measured-in-Belle resolution
281 hit_localPhiErr = hit_localPhiErr * sqrt(factor);
282
283 nStrips = hit->getZStripMax() - hit->getZStripMin() + 1;
284 dn = nStrips - 1.5;
285 factor = std::pow((0.9 + 0.4 * dn * dn), 1.5) * 0.55;//measured-in-Belle resolution
286 hit_localZErr = hit_localZErr * sqrt(factor);
287 }
288
289 Hep3Vector globalOrigin = corMod->getGlobalOrigin();
290 double sinphi = globalOrigin[1] / globalOrigin.mag();
291 double cosphi = globalOrigin[0] / globalOrigin.mag();
292
293 HepMatrix globalHitErr(3, 3, 0);
294 globalHitErr[0][0] = pow(hit_localPhiErr * sinphi, 2); //x
295 globalHitErr[0][1] = (hit_localPhiErr * sinphi) * (hit_localPhiErr * cosphi);
296 globalHitErr[0][2] = 0;
297 globalHitErr[1][1] = pow(hit_localPhiErr * cosphi, 2);;
298 globalHitErr[1][0] = (hit_localPhiErr * sinphi) * (hit_localPhiErr * cosphi);
299 globalHitErr[1][2] = 0;
300 globalHitErr[2][2] = pow(hit_localZErr, 2);
301 globalHitErr[2][0] = 0;
302 globalHitErr[2][1] = 0;
303
304 //HepMatrix B(2, 2, 0); // Matrix for derivatives
305
306 // Derivatives of x (z) = y/b - a/b with respect to y.
307 //B[ MY ][ AY ] = 1./m_GlobalPar[BY];
308 //B[ MZ ][ BY ] = 1./m_GlobalPar[BZ] ;
309 //double errors_b[0] = B[ MY ][ AY ]*globalHitErr[1][1];
310 //double errors_b[1] = B[ MZ ][ BY ]*globalHitErr[1][1];
311
312 HepMatrix error_mea(2, 2, 0);
313 error = sqrt(errors[ MY ][ MY ] +
314 errors[ MZ ][ MZ ] +
315 //errors_b[0] + errors_b[1] + //error of prediced point due to error of inPos (y_mea)
316 globalHitErr[0][0] +
317 globalHitErr[1][1] + //y_mea error is correlated to the error of predicted point, but we didn't consider that part in errors
318 globalHitErr[2][2]); //we ignore y_mea error?
319
320 if (error != 0.0) {
321 sigma = distance / error;
322 } else {
323 sigma = DBL_MAX;
324 }
325
326 return (distance);
327}
328
329
331double BKLMTrackFitter::fit1dSectorTrack(std::list< KLMHit2d* > hitList,
332 HepVector& eta,
333 HepSymMatrix& error,
334 int depDir, int indDir)
335{
336
337// Fit d = a + bi, where d is dependent direction and i is independent
338 Hep3Vector localHitPos;
339 HepMatrix localHitErr(3, 3, 0);
340
341 double indPos = 0;
342 double depPos = 0;
343
344 // Matrix based solution
345
346 int noPoints = hitList.size();
347
348 // Derivative of y = a + bx, with respect to a and b evaluated at x.
349 HepMatrix A(noPoints, 2, 0);
350
351 // Measured data points.
352 HepVector y(noPoints, 0);
353
354 // Inverse of covariance (error) matrix, also known as the weight matrix.
355 // In plain English: V_y_inverse_nn = 1 / (error of nth measurement)^2
356 HepDiagMatrix V_y_inverse(noPoints, 0);
357 //HepSymMatrix V_y_inverse(noPoints, 0);
358
359 // Error or correlation matrix for coefficients (2x2 matrices)
360 HepSymMatrix V_A, V_A_inverse;
361
362 int section = (*hitList.begin())->getSection();
363 int sector = (*hitList.begin())->getSector();
364
365 m_GeoPar = GeometryPar::instance();
366 const Belle2::bklm::Module* refMod = m_GeoPar->findModule(section, sector, 1);
367
368 int n = 0;
369 for (KLMHit2d* hit : hitList) {
370
371 if (hit->getSection() != section || hit->getSector() != sector) {
372 continue;
373 }
374
375 // m_GeoPar = GeometryPar::instance();
376 //const Belle2::bklm::Module* refMod = m_GeoPar->findModule(hit->getSection(), hit->getSector(), 1);
377 const Belle2::bklm::Module* corMod = m_GeoPar->findModule(hit->getSection(), hit->getSector(), hit->getLayer());
378
379 CLHEP::Hep3Vector globalPos;
380 globalPos[0] = hit->getPositionX();
381 globalPos[1] = hit->getPositionY();
382 globalPos[2] = hit->getPositionZ();
383
384 localHitPos = refMod->globalToLocal(globalPos);
385 double hit_localPhiErr = corMod->getPhiStripWidth() / sqrt(12);
386 double hit_localZErr = corMod->getZStripWidth() / sqrt(12);
387
388 if (hit->inRPC()) {
389 //+++ scale localErr based on measured-in-Belle resolution
390 int nStrips = hit->getPhiStripMax() - hit->getPhiStripMin() + 1;
391 double dn = nStrips - 1.5;
392 double factor = std::pow((0.9 + 0.4 * dn * dn), 1.5) * 0.60;//measured-in-Belle resolution
393 hit_localPhiErr = hit_localPhiErr * sqrt(factor);
394
395 nStrips = hit->getZStripMax() - hit->getZStripMin() + 1;
396 dn = nStrips - 1.5;
397 factor = std::pow((0.9 + 0.4 * dn * dn), 1.5) * 0.55;//measured-in-Belle resolution
398 hit_localZErr = hit_localZErr * sqrt(factor);
399 }
400
401 localHitErr[0][0] = 0.0; //x
402 localHitErr[0][1] = 0;
403 localHitErr[0][2] = 0;
404 localHitErr[1][1] = hit_localPhiErr;
405 localHitErr[1][0] = 0;
406 localHitErr[1][2] = 0;
407 localHitErr[2][2] = hit_localZErr;
408 localHitErr[2][0] = 0;
409 localHitErr[2][1] = 0;
410
411 switch (indDir) {
412
413 case VX:
414 indPos = localHitPos.x();
415 break;
416
417 case VY:
418 indPos = localHitPos.y();
419 break;
420
421 case VZ:
422 indPos = localHitPos.z();
423 break;
424
425 default:
426 B2DEBUG(20, "error in klm_trackSectorFit: illegal direction");
427
428 }
429
430 switch (depDir) {
431
432 case VX:
433 depPos = localHitPos.x();
434 break;
435
436 case VY:
437 depPos = localHitPos.y();
438 break;
439
440 case VZ:
441 depPos = localHitPos.z();
442 break;
443
444 default:
445 B2DEBUG(20, "error in klm_trackSectorFit: illegal direction");
446
447 }
448
449
450 A[ n ][ 0 ] = 1;
451 A[ n ][ 1 ] = indPos;
452
453 y[ n ] = depPos;
454
455 if (localHitErr[ depDir ][ depDir ] > 0.0) {
456 V_y_inverse[ n ][ n ] = 1.0 / localHitErr[ depDir ][ depDir ];
457 } else {
458 V_y_inverse[ n ][ n ] = DBL_MAX;
459 }
460 ++n;//n is the index of measured points/hits
461 }
462
463 V_A_inverse = V_y_inverse.similarityT(A);
464
465 int ierr = 0;
466 V_A = V_A_inverse.inverse(ierr);
467
468 eta = V_A * A.T() * V_y_inverse * y;
469 error = V_A;
470
471// Calculate residuals:
472 HepMatrix residual = y - A * eta;
473 HepMatrix chisqr = residual.T() * V_y_inverse * residual;
474
475 return (chisqr.trace());
476
477}
478
480double BKLMTrackFitter::fit1dTrack(std::list< KLMHit2d* > hitList,
481 HepVector& eta,
482 HepSymMatrix& error,
483 int depDir, int indDir)
484{
485// Fit d = a + bi, where d is dependent direction and i is independent
486// in global system we assume y = a + b*x and y = c + d*z different from local fit
487
488 HepMatrix globalHitErr(3, 3, 0);
489
490 double indPos = 0;
491 double depPos = 0;
492
493 // Matrix based solution
494 int noPoints = hitList.size();
495
496 // Derivative of y = a + bx, with respect to a and b evaluated at x.
497 HepMatrix A(noPoints, 2, 0);
498
499 // Measured data points.
500 HepVector y(noPoints, 0);
501
502 // Inverse of covariance (error) matrix, also known as the weight matrix.
503 // In plain English: V_y_inverse_nn = 1 / (error of nth measurement)^2
504 HepDiagMatrix V_y_inverse(noPoints, 0);
505
506 // Error or correlation matrix for coefficients (2x2 matrices)
507 HepSymMatrix V_A, V_A_inverse;
508
509 m_GeoPar = GeometryPar::instance();
510 const Belle2::bklm::Module* corMod;
511
512 int n = 0;
513 for (KLMHit2d* hit : hitList) {
514
515 // m_GeoPar = GeometryPar::instance();
516 //const Belle2::bklm::Module* refMod = m_GeoPar->findModule(hit->getSection(), hit->getSector(), 1);
517 corMod = m_GeoPar->findModule(hit->getSection(), hit->getSector(), hit->getLayer());
518
519 CLHEP::Hep3Vector globalPos;
520 globalPos[0] = hit->getPositionX();
521 globalPos[1] = hit->getPositionY();
522 globalPos[2] = hit->getPositionZ();
523
524 //localHitPos = refMod->globalToLocal(globalPos);
525 double hit_localPhiErr = corMod->getPhiStripWidth() / sqrt(12);
526 double hit_localZErr = corMod->getZStripWidth() / sqrt(12);
527
528 if (hit->inRPC()) {
529 //+++ scale localErr based on measured-in-Belle resolution
530 int nStrips = hit->getPhiStripMax() - hit->getPhiStripMin() + 1;
531 double dn = nStrips - 1.5;
532 double factor = std::pow((0.9 + 0.4 * dn * dn), 1.5) * 0.60;//measured-in-Belle resolution
533 hit_localPhiErr = hit_localPhiErr * sqrt(factor);
534
535 nStrips = hit->getZStripMax() - hit->getZStripMin() + 1;
536 dn = nStrips - 1.5;
537 factor = std::pow((0.9 + 0.4 * dn * dn), 1.5) * 0.55;//measured-in-Belle resolution
538 hit_localZErr = hit_localZErr * sqrt(factor);
539 }
540
541 Hep3Vector globalOrigin = corMod->getGlobalOrigin();
542 double sinphi = globalOrigin[1] / globalOrigin.mag();
543 double cosphi = globalOrigin[0] / globalOrigin.mag();
544
545 globalHitErr[0][0] = pow(hit_localPhiErr * sinphi, 2); //x
546 globalHitErr[0][1] = (hit_localPhiErr * sinphi) * (hit_localPhiErr * cosphi);
547 globalHitErr[0][2] = 0;
548 globalHitErr[1][1] = pow(hit_localPhiErr * cosphi, 2);;
549 globalHitErr[1][0] = (hit_localPhiErr * sinphi) * (hit_localPhiErr * cosphi);
550 globalHitErr[1][2] = 0;
551 globalHitErr[2][2] = pow(hit_localZErr, 2);;
552 globalHitErr[2][0] = 0;
553 globalHitErr[2][1] = 0;
554
555 switch (indDir) {
556
557 case VX:
558 indPos = globalPos.x();
559 break;
560
561 case VY:
562 indPos = globalPos.y();
563 break;
564
565 case VZ:
566 indPos = globalPos.z();
567 break;
568
569 default:
570 B2DEBUG(20, "error in bklm trackFit: illegal direction");
571
572 }
573
574 switch (depDir) {
575
576 case VX:
577 depPos = globalPos.x();
578 break;
579
580 case VY:
581 depPos = globalPos.y();
582 break;
583
584 case VZ:
585 depPos = globalPos.z();
586 break;
587
588 default:
589 B2DEBUG(20, "error in bklm trackFit: illegal direction");
590
591 }
592
593
594 A[ n ][ 0 ] = 1;
595 A[ n ][ 1 ] = indPos;
596
597 y[ n ] = depPos;
598
599 double error_raw = globalHitErr[indDir][indDir] + globalHitErr[depDir][depDir];
600 //double correlate_ = 2.0*globalHitErr[indDir][depDir]; //?
601 //double weight= error_raw - correlate_;
602 double weight = error_raw;
603 if (weight > 0) {
604 V_y_inverse[ n ][ n ] = 1.0 / weight;
605 } else {
606 V_y_inverse[ n ][ n ] = DBL_MAX;
607 }
608 ++n;//n is the index of measured points/hits
609 }
610
611 //HepMatrix AT = A.T();
612 //HepMatrix tmp = AT*V_y_inverse;
613 //V_A_inverse = tmp*A;
614 V_A_inverse = V_y_inverse.similarityT(A);
615
616 int ierr = 0;
617 V_A = V_A_inverse.inverse(ierr);
618
619 eta = V_A * A.T() * V_y_inverse * y;
620 error = V_A;
621
622// Calculate residuals:
623 HepMatrix residual = y - A * eta;
624 HepMatrix chisqr = residual.T() * V_y_inverse * residual;
625
626 return (chisqr.trace());
627
628}
Belle2::BKLMTrackFitter::fit
double fit(std::list< KLMHit2d * > &listTrackPoint)
do fit and returns chi square of the fit.
Definition: BKLMTrackFitter.cc:60
Belle2::BKLMTrackFitter::m_Chi2
float m_Chi2
Chi square of fit.
Definition: BKLMTrackFitter.h:133
Belle2::BKLMTrackFitter::fit1dTrack
double fit1dTrack(std::list< KLMHit2d * > hitList, CLHEP::HepVector &eta, CLHEP::HepSymMatrix &error, int depDir, int indDir)
do fit in the global system
Definition: BKLMTrackFitter.cc:480
Belle2::BKLMTrackFitter::m_NumHit
int m_NumHit
the number of hits on this track
Definition: BKLMTrackFitter.h:136
Belle2::BKLMTrackFitter::globalDistanceToHit
double globalDistanceToHit(KLMHit2d *hit, double &error, double &sigma)
Distance from track to a hit in the global system.
Definition: BKLMTrackFitter.cc:233
Belle2::BKLMTrackFitter::BKLMTrackFitter
BKLMTrackFitter()
Default constructor.
Definition: BKLMTrackFitter.cc:40
Belle2::BKLMTrackFitter::distanceToHit
double distanceToHit(KLMHit2d *hit, double &error, double &sigma)
Distance from track to a hit in the plane of the module.
Definition: BKLMTrackFitter.cc:155
Belle2::BKLMTrackFitter::m_GlobalErr
CLHEP::HepSymMatrix m_GlobalErr
track params errors in global system
Definition: BKLMTrackFitter.h:151
Belle2::BKLMTrackFitter::m_GlobalPar
CLHEP::HepVector m_GlobalPar
track params in global system
Definition: BKLMTrackFitter.h:148
Belle2::BKLMTrackFitter::m_GeoPar
bklm::GeometryPar * m_GeoPar
pointer to GeometryPar singleton
Definition: BKLMTrackFitter.h:154
Belle2::BKLMTrackFitter::fit1dSectorTrack
double fit1dSectorTrack(std::list< KLMHit2d * > hitList, CLHEP::HepVector &eta, CLHEP::HepSymMatrix &error, int depDir, int indDir)
do fit in the y-x plane or z-x plane
Definition: BKLMTrackFitter.cc:331
Belle2::BKLMTrackFitter::m_SectorPar
CLHEP::HepVector m_SectorPar
track params in the sector local system
Definition: BKLMTrackFitter.h:142
Belle2::BKLMTrackFitter::m_SectorErr
CLHEP::HepSymMatrix m_SectorErr
track params errors in the sector local system
Definition: BKLMTrackFitter.h:145
Belle2::BKLMTrackFitter::m_Valid
bool m_Valid
Is fit valid.
Definition: BKLMTrackFitter.h:127
Belle2::BKLMTrackFitter::m_globalFit
bool m_globalFit
do fit in the local system or global system false: local sys; true: global sys.
Definition: BKLMTrackFitter.h:139
Belle2::BKLMTrackFitter::m_Good
bool m_Good
Is fit good.
Definition: BKLMTrackFitter.h:130
Belle2::KLMHit2d
KLM 2d hit.
Definition: KLMHit2d.h:33
Belle2::bklm::GeometryPar::findModule
const Module * findModule(int section, int sector, int layer) const
Get the pointer to the definition of a module.
Definition: GeometryPar.cc:721
Belle2::bklm::GeometryPar::instance
static GeometryPar * instance(void)
Static method to get a reference to the singleton GeometryPar instance.
Definition: GeometryPar.cc:27
Belle2::bklm::Module
Define the geometry of a BKLM module Each sector [octant] contains Modules.
Definition: Module.h:76
Belle2::bklm::Module::globalToLocal
const CLHEP::Hep3Vector globalToLocal(const CLHEP::Hep3Vector &v, bool reco=false) const
Transform space-point within this module from global to local coordinates.
Definition: Module.cc:339
Belle2::bklm::Module::getPhiStripWidth
double getPhiStripWidth() const
Get phi-strip width.
Definition: Module.h:137
Belle2::bklm::Module::getGlobalOrigin
const CLHEP::Hep3Vector getGlobalOrigin() const
Return the position (in global coordinates) of this module's sensitive-volume origin.
Definition: Module.h:285
Belle2::bklm::Module::localToGlobal
const CLHEP::Hep3Vector localToGlobal(const CLHEP::Hep3Vector &v, bool reco=false) const
Transform space-point within this module from local to global coordinates.
Definition: Module.cc:326
Belle2::bklm::Module::getZStripWidth
double getZStripWidth() const
Get z-strip width.
Definition: Module.h:155
Belle2::sqrt
double sqrt(double a)
sqrt for double
Definition: beamHelpers.h:28
Belle2
Abstract base class for different kinds of events.
Definition: MillepedeAlgorithm.h:17