Belle II Software development
RecoilMassKFit.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#include <analysis/VertexFitting/KFit/RecoilMassKFit.h>
10
11#include <analysis/dataobjects/Particle.h>
12#include <analysis/VertexFitting/KFit/MakeMotherKFit.h>
13#include <analysis/utility/CLHEPToROOT.h>
14#include <framework/gearbox/Const.h>
15
16#include <cstdio>
17
18#include <TMath.h>
19
20using namespace std;
21using namespace Belle2;
22using namespace Belle2::analysis;
23using namespace CLHEP;
24using namespace ROOT::Math;
25
27{
28 m_FlagFitted = false;
31 m_FlagAtDecayPoint = true;
33 m_d = HepMatrix(1, 1, 0);
34 m_V_D = HepMatrix(1, 1, 0);
35 m_lam = HepMatrix(1, 1, 0);
36 m_AfterVertexError = HepSymMatrix(3, 0);
37 m_recoilMass = -1.0;
38 m_FourMomentum = PxPyPzEVector();
39}
40
41
43
44
46RecoilMassKFit::setVertex(const HepPoint3D& v) {
48
50}
51
52
54RecoilMassKFit::setVertexError(const HepSymMatrix& e) {
55 if (e.num_row() != 3)
56 {
58 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
59 return m_ErrorCode;
60 }
61
64
66}
67
70 m_recoilMass = m;
71
73}
74
75
77RecoilMassKFit::setFourMomentum(const PxPyPzEVector& m) {
79
81}
82
83
90
91
94 m_IsFixMass.push_back(true);
95
97}
98
99
102 m_IsFixMass.push_back(false);
103
105}
106
107
110 if (e.num_row() != 3 || e.num_col() != KFitConst::kNumber7)
111 {
113 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
114 return m_ErrorCode;
115 }
116
117 m_BeforeTrackVertexError.push_back(e);
120
122}
123
124
127 HepMatrix zero(3, KFitConst::kNumber7, 0);
128
129 return this->setTrackVertexError(zero);
130}
131
132
135 return KFitBase::setCorrelation(m);
136}
137
138
143
144
145const HepPoint3D
146RecoilMassKFit::getVertex(const int flag) const
147{
148 if (flag == KFitConst::kAfterFit && !isFitted()) return HepPoint3D();
149
150 switch (flag) {
152 return m_BeforeVertex;
153
155 return m_AfterVertex;
156
157 default:
158 KFitError::displayError(__FILE__, __LINE__, __func__, KFitError::kOutOfRange);
159 return HepPoint3D();
160 }
161}
162
163
164const HepSymMatrix
166{
167 if (flag == KFitConst::kAfterFit && !isFitted()) return HepSymMatrix(3, 0);
168
169 if (flag == KFitConst::kBeforeFit)
170 return m_BeforeVertexError;
172 return m_AfterVertexError;
173 else {
174 KFitError::displayError(__FILE__, __LINE__, __func__, KFitError::kOutOfRange);
175 return HepSymMatrix(3, 0);
176 }
177}
178
179
180bool
185
186
187bool
192
193
194double
196{
197 return m_CHIsq;
198}
199
200
201const HepMatrix
202RecoilMassKFit::getTrackVertexError(const int id, const int flag) const
203{
204 if (flag == KFitConst::kAfterFit && !isFitted()) return HepMatrix(3, KFitConst::kNumber7, 0);
205 if (!isTrackIDInRange(id)) return HepMatrix(3, KFitConst::kNumber7, 0);
206
207 if (flag == KFitConst::kBeforeFit)
208 return m_BeforeTrackVertexError[id];
210 return m_AfterTrackVertexError[id];
211 else {
212 KFitError::displayError(__FILE__, __LINE__, __func__, KFitError::kOutOfRange);
213 return HepMatrix(3, KFitConst::kNumber7, 0);
214 }
215}
216
217
218double
220{
221 if (!isFitted()) return -1;
222 if (!isTrackIDInRange(id)) return -1;
223
224 if (m_IsFixMass[id]) {
225
226 HepMatrix da(m_Tracks[id].getFitParameter(KFitConst::kBeforeFit) - m_Tracks[id].getFitParameter(KFitConst::kAfterFit));
227 int err_inverse = 0;
228 const double chisq = (da.T() * (m_Tracks[id].getFitError(KFitConst::kBeforeFit).inverse(err_inverse)) * da)[0][0];
229
230 if (err_inverse) {
231 KFitError::displayError(__FILE__, __LINE__, __func__, KFitError::kCannotGetMatrixInverse);
232 return -1;
233 }
234
235 return chisq;
236
237 } else {
238
239 HepMatrix da(m_Tracks[id].getMomPos(KFitConst::kBeforeFit) - m_Tracks[id].getMomPos(KFitConst::kAfterFit));
240 int err_inverse = 0;
241 const double chisq = (da.T() * (m_Tracks[id].getError(KFitConst::kBeforeFit).inverse(err_inverse)) * da)[0][0];
242
243 if (err_inverse) {
244 KFitError::displayError(__FILE__, __LINE__, __func__, KFitError::kCannotGetMatrixInverse);
245 return -1;
246 }
247
248 return chisq;
249
250 }
251}
252
253
254const HepMatrix
255RecoilMassKFit::getCorrelation(const int id1, const int id2, const int flag) const
256{
257 if (flag == KFitConst::kAfterFit && !isFitted()) return HepMatrix(KFitConst::kNumber7, KFitConst::kNumber7, 0);
258 if (!isTrackIDInRange(id1)) return HepMatrix(KFitConst::kNumber7, KFitConst::kNumber7, 0);
259 if (!isTrackIDInRange(id2)) return HepMatrix(KFitConst::kNumber7, KFitConst::kNumber7, 0);
260
261 switch (flag) {
263 return KFitBase::getCorrelation(id1, id2, flag);
264
266 return makeError3(
267 this->getTrackMomentum(id1),
268 this->getTrackMomentum(id2),
269 m_V_al_1.sub(KFitConst::kNumber7 * id1 + 1, KFitConst::kNumber7 * (id1 + 1), KFitConst::kNumber7 * id2 + 1,
270 KFitConst::kNumber7 * (id2 + 1)),
271 m_IsFixMass[id1],
272 m_IsFixMass[id2]);
273
274 default:
275 KFitError::displayError(__FILE__, __LINE__, __func__, KFitError::kOutOfRange);
276 return HepMatrix(KFitConst::kNumber7, KFitConst::kNumber7, 0);
277 }
278}
279
280
285
286
290 {
292 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
293 return m_ErrorCode;
294 }
295
296
297 if (m_IsFixMass.size() == 0)
298 {
299 // If no fix_mass flag at all,
300 // all tracks are considered to be fixed at mass.
301 for (int i = 0; i < m_TrackCount; i++) this->fixMass();
302 } else if (m_IsFixMass.size() != (unsigned int)m_TrackCount)
303 {
305 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
306 return m_ErrorCode;
307 }
308
309
311 {
312 int index = 0;
313 m_al_0 = HepMatrix(KFitConst::kNumber7 * m_TrackCount, 1, 0);
314 m_property = HepMatrix(m_TrackCount, 3, 0);
315 m_V_al_0 = HepSymMatrix(KFitConst::kNumber7 * m_TrackCount, 0);
316
317 for (auto& track : m_Tracks) {
318 // momentum x,y,z and position x,y,z
319 m_al_0[index * KFitConst::kNumber7 + 0][0] = track.getMomentum(KFitConst::kBeforeFit).x();
320 m_al_0[index * KFitConst::kNumber7 + 1][0] = track.getMomentum(KFitConst::kBeforeFit).y();
321 m_al_0[index * KFitConst::kNumber7 + 2][0] = track.getMomentum(KFitConst::kBeforeFit).z();
322 m_al_0[index * KFitConst::kNumber7 + 3][0] = track.getMomentum(KFitConst::kBeforeFit).t();
323 m_al_0[index * KFitConst::kNumber7 + 4][0] = track.getPosition(KFitConst::kBeforeFit).x();
324 m_al_0[index * KFitConst::kNumber7 + 5][0] = track.getPosition(KFitConst::kBeforeFit).y();
325 m_al_0[index * KFitConst::kNumber7 + 6][0] = track.getPosition(KFitConst::kBeforeFit).z();
326 // these error
327 m_V_al_0.sub(index * KFitConst::kNumber7 + 1, track.getError(KFitConst::kBeforeFit));
328 // charge, mass, a
329 m_property[index][0] = track.getCharge();
330 m_property[index][1] = track.getMass();
331 const double c = Const::speedOfLight * 1e-4;
332 m_property[index][2] = -c * m_MagneticField * track.getCharge();
333 index++;
334 }
335
336 // error between track and track
337 if (m_FlagCorrelation) {
338 this->prepareCorrelation();
340 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
341 return m_ErrorCode;
342 }
343 }
344
345 // set member matrix
346 m_al_1 = m_al_0;
347
348 // define size of matrix
351
352 } else
353 {
354 //TODO: Not Implemented
356 }
357
359}
360
361
364 char buf[1024];
365 sprintf(buf, "%s:%s(): internal error; this function should never be called", __FILE__, __func__);
366 B2FATAL(buf);
367
368 /* NEVER REACHEd HERE */
370}
371
372
375 if (m_BeforeCorrelation.size() != static_cast<unsigned int>(m_TrackCount * (m_TrackCount - 1) / 2))
376 {
378 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
379 return m_ErrorCode;
380 }
381
382 int row = 0, col = 0;
383
384 for (auto& hm : m_BeforeCorrelation)
385 {
386 // counter
387 row++;
388 if (row == m_TrackCount) {
389 col++;
390 row = col + 1;
391 }
392
393 int ii = 0, jj = 0;
394 for (int i = KFitConst::kNumber7 * row; i < KFitConst::kNumber7 * (row + 1); i++) {
395 for (int j = KFitConst::kNumber7 * col; j < KFitConst::kNumber7 * (col + 1); j++) {
396 m_V_al_0[i][j] = hm[ii][jj];
397 jj++;
398 }
399 jj = 0;
400 ii++;
401 }
402 }
403
405 {
406 //TODO: Not Implemented
408
409 // ...error of vertex
411
412 // ...error matrix between vertex and tracks
414 if (m_BeforeTrackVertexError.size() != (unsigned int)m_TrackCount) {
416 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
417 return m_ErrorCode;
418 }
419
420 int i = 0;
421 for (auto& hm : m_BeforeTrackVertexError) {
422 for (int j = 0; j < 3; j++) for (int k = 0; k < KFitConst::kNumber7; k++) {
424 }
425 i++;
426 }
427 }
428 }
429
431}
432
433
436 Hep3Vector h3v;
437 int index = 0;
438 for (auto& pdata : m_Tracks)
439 {
440 // tracks
441 // momentum
442 h3v.setX(m_al_1[index * KFitConst::kNumber7 + 0][0]);
443 h3v.setY(m_al_1[index * KFitConst::kNumber7 + 1][0]);
444 h3v.setZ(m_al_1[index * KFitConst::kNumber7 + 2][0]);
445 if (m_IsFixMass[index])
446 pdata.setMomentum(HepLorentzVector(h3v, sqrt(h3v.mag2() + pdata.getMass()*pdata.getMass())), KFitConst::kAfterFit);
447 else
448 pdata.setMomentum(HepLorentzVector(h3v, m_al_1[index * KFitConst::kNumber7 + 3][0]), KFitConst::kAfterFit);
449 // position
450 pdata.setPosition(HepPoint3D(
451 m_al_1[index * KFitConst::kNumber7 + 4][0],
452 m_al_1[index * KFitConst::kNumber7 + 5][0],
454 // error of the tracks
455 pdata.setError(this->makeError3(pdata.getMomentum(),
456 m_V_al_1.sub(
457 index * KFitConst::kNumber7 + 1,
458 (index + 1)*KFitConst::kNumber7,
459 index * KFitConst::kNumber7 + 1,
460 (index + 1)*KFitConst::kNumber7), m_IsFixMass[index]),
462 if (m_ErrorCode != KFitError::kNoError) break;
463 index++;
464 }
465
467 {
468 //TODO: Not Implemented
470 } else
471 {
472 // not fit
474 }
475
477}
478
479
483 {
484
485 HepMatrix al_1_prime(m_al_1);
486 HepMatrix Sum_al_1(4, 1, 0);
487 std::vector<double> energy(m_TrackCount);
488 double a;
489
490 for (int i = 0; i < m_TrackCount; i++) {
491 a = m_property[i][2];
492 if (!m_FlagAtDecayPoint) a = 0.;
493 al_1_prime[i * KFitConst::kNumber7 + 0][0] -= a * (m_BeforeVertex.y() - al_1_prime[i * KFitConst::kNumber7 + 5][0]);
494 al_1_prime[i * KFitConst::kNumber7 + 1][0] += a * (m_BeforeVertex.x() - al_1_prime[i * KFitConst::kNumber7 + 4][0]);
495 energy[i] = sqrt(al_1_prime[i * KFitConst::kNumber7 + 0][0] * al_1_prime[i * KFitConst::kNumber7 + 0][0] +
496 al_1_prime[i * KFitConst::kNumber7 + 1][0] * al_1_prime[i * KFitConst::kNumber7 + 1][0] +
497 al_1_prime[i * KFitConst::kNumber7 + 2][0] * al_1_prime[i * KFitConst::kNumber7 + 2][0] +
498 m_property[i][1] * m_property[i][1]);
499 if (m_IsFixMass[i])
500 Sum_al_1[3][0] += energy[i];
501 else
502 Sum_al_1[3][0] += al_1_prime[i * KFitConst::kNumber7 + 3][0];
503 }
504
505 for (int i = 0; i < m_TrackCount; i++) {
506 for (int j = 0; j < 3; j++) Sum_al_1[j][0] += al_1_prime[i * KFitConst::kNumber7 + j][0];
507 }
508
509 Sum_al_1[0][0] -= m_FourMomentum.Px();
510 Sum_al_1[1][0] -= m_FourMomentum.Py();
511 Sum_al_1[2][0] -= m_FourMomentum.Pz();
512 Sum_al_1[3][0] -= m_FourMomentum.E();
513
514 m_d[0][0] =
515 + Sum_al_1[3][0] * Sum_al_1[3][0] - Sum_al_1[0][0] * Sum_al_1[0][0]
516 - Sum_al_1[1][0] * Sum_al_1[1][0] - Sum_al_1[2][0] * Sum_al_1[2][0]
518
519 for (int i = 0; i < m_TrackCount; i++) {
520 if (energy[i] == 0) {
522 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
523 break;
524 }
525
526 a = m_property[i][2];
527 if (!m_FlagAtDecayPoint) a = 0.;
528
529 if (m_IsFixMass[i]) {
530 double invE = 1. / energy[i];
531 m_D[0][i * KFitConst::kNumber7 + 0] = 2.*(Sum_al_1[3][0] * al_1_prime[i * KFitConst::kNumber7 + 0][0] * invE - Sum_al_1[0][0]);
532 m_D[0][i * KFitConst::kNumber7 + 1] = 2.*(Sum_al_1[3][0] * al_1_prime[i * KFitConst::kNumber7 + 1][0] * invE - Sum_al_1[1][0]);
533 m_D[0][i * KFitConst::kNumber7 + 2] = 2.*(Sum_al_1[3][0] * al_1_prime[i * KFitConst::kNumber7 + 2][0] * invE - Sum_al_1[2][0]);
534 m_D[0][i * KFitConst::kNumber7 + 3] = 0.;
535 m_D[0][i * KFitConst::kNumber7 + 4] = -2.*(Sum_al_1[3][0] * al_1_prime[i * KFitConst::kNumber7 + 1][0] * invE - Sum_al_1[1][0]) * a;
536 m_D[0][i * KFitConst::kNumber7 + 5] = 2.*(Sum_al_1[3][0] * al_1_prime[i * KFitConst::kNumber7 + 0][0] * invE - Sum_al_1[0][0]) * a;
537 m_D[0][i * KFitConst::kNumber7 + 6] = 0.;
538 } else {
539 m_D[0][i * KFitConst::kNumber7 + 0] = -2.*Sum_al_1[0][0];
540 m_D[0][i * KFitConst::kNumber7 + 1] = -2.*Sum_al_1[1][0];
541 m_D[0][i * KFitConst::kNumber7 + 2] = -2.*Sum_al_1[2][0];
542 m_D[0][i * KFitConst::kNumber7 + 3] = 2.*Sum_al_1[3][0];
543 m_D[0][i * KFitConst::kNumber7 + 4] = 2.*Sum_al_1[1][0] * a;
544 m_D[0][i * KFitConst::kNumber7 + 5] = -2.*Sum_al_1[0][0] * a;
545 m_D[0][i * KFitConst::kNumber7 + 6] = 0.;
546 }
547 }
548
549 } else
550 {
551 //TODO: Not Implemented
553 }
554
556}
557
558
565
567{
568 MakeMotherKFit kmm;
570 unsigned n = getTrackCount();
571 for (unsigned i = 0; i < n; ++i) {
573 getTrack(i).getCharge());
576 for (unsigned j = i + 1; j < n; ++j) {
578 }
579 }
580 kmm.setVertex(getVertex());
583 m_ErrorCode = kmm.doMake();
585 return m_ErrorCode;
586 double chi2 = getCHIsq();
587 int ndf = getNDF();
588 double prob = TMath::Prob(chi2, ndf);
589 //
590 mother->writeExtraInfo("chiSquared", chi2);
591 mother->writeExtraInfo("ndf", ndf);
592
593 mother->updateMomentum(
594 CLHEPToROOT::getLorentzVector(kmm.getMotherMomentum()),
595 CLHEPToROOT::getXYZVector(kmm.getMotherPosition()),
596 CLHEPToROOT::getTMatrixFSym(kmm.getMotherError()),
597 prob);
599 return m_ErrorCode;
600}
static const double speedOfLight
[cm/ns]
Definition Const.h:695
Class to store reconstructed particles.
Definition Particle.h:76
void writeExtraInfo(const std::string &name, const double value)
Sets the user defined extraInfo.
Definition Particle.cc:1393
void updateMomentum(const ROOT::Math::PxPyPzEVector &p4, const ROOT::Math::XYZVector &vertex, const TMatrixFSym &errMatrix, double pValue)
Sets Lorentz vector, position, 7x7 error matrix and p-value.
Definition Particle.h:397
int m_NecessaryTrackCount
Number needed tracks to perform fit.
Definition KFitBase.h:303
double m_MagneticField
Magnetic field.
Definition KFitBase.h:311
CLHEP::HepMatrix m_al_1
See J.Tanaka Ph.D (2001) p136 for definition.
Definition KFitBase.h:259
virtual enum KFitError::ECode setCorrelation(const CLHEP::HepMatrix &c)
Set a correlation matrix.
Definition KFitBase.cc:70
const CLHEP::HepSymMatrix makeError3(const CLHEP::HepLorentzVector &p, const CLHEP::HepMatrix &e, const bool is_fix_mass) const
Rebuild an error matrix from a Lorentz vector and an error matrix.
Definition KFitBase.cc:320
const CLHEP::HepSymMatrix getTrackError(const int id) const
Get an error matrix of the track.
Definition KFitBase.cc:168
const CLHEP::HepLorentzVector getTrackMomentum(const int id) const
Get a Lorentz vector of the track.
Definition KFitBase.cc:154
CLHEP::HepMatrix m_lam
See J.Tanaka Ph.D (2001) p137 for definition.
Definition KFitBase.h:276
const HepPoint3D getTrackPosition(const int id) const
Get a position of the track.
Definition KFitBase.cc:161
CLHEP::HepMatrix m_property
Container of charges and masses.
Definition KFitBase.h:263
enum KFitError::ECode m_ErrorCode
Error code.
Definition KFitBase.h:243
virtual enum KFitError::ECode setZeroCorrelation(void)
Indicate no correlation between tracks.
Definition KFitBase.cc:85
CLHEP::HepMatrix m_V_al_1
See J.Tanaka Ph.D (2001) p138 for definition.
Definition KFitBase.h:274
virtual int getNDF(void) const
Get an NDF of the fit.
Definition KFitBase.cc:114
CLHEP::HepMatrix m_d
See J.Tanaka Ph.D (2001) p137 for definition.
Definition KFitBase.h:268
bool isFitted(void) const
Return false if fit is not performed yet or performed fit is failed; otherwise true.
Definition KFitBase.cc:728
CLHEP::HepMatrix m_D
See J.Tanaka Ph.D (2001) p137 for definition.
Definition KFitBase.h:266
CLHEP::HepMatrix m_V_D
See J.Tanaka Ph.D (2001) p138 for definition.
Definition KFitBase.h:271
bool isTrackIDInRange(const int id) const
Check if the id is in the range.
Definition KFitBase.cc:739
virtual const CLHEP::HepMatrix getCorrelation(const int id1, const int id2, const int flag=KFitConst::kAfterFit) const
Get a correlation matrix between two tracks.
Definition KFitBase.cc:183
bool m_FlagCorrelation
Flag whether a correlation among tracks exists.
Definition KFitBase.h:306
CLHEP::HepSymMatrix m_V_al_0
See J.Tanaka Ph.D (2001) p137 for definition.
Definition KFitBase.h:255
const KFitTrack getTrack(const int id) const
Get a specified track object.
Definition KFitBase.cc:175
std::vector< CLHEP::HepMatrix > m_BeforeCorrelation
Container of input correlation matrices.
Definition KFitBase.h:251
bool m_FlagFitted
Flag to indicate if the fit is performed and succeeded.
Definition KFitBase.h:245
double m_CHIsq
chi-square of the fit.
Definition KFitBase.h:297
int getTrackCount(void) const
Get the number of added tracks.
Definition KFitBase.cc:107
int m_NDF
NDF of the fit.
Definition KFitBase.h:295
std::vector< KFitTrack > m_Tracks
Container of input tracks.
Definition KFitBase.h:249
int m_TrackCount
Number of tracks.
Definition KFitBase.h:301
CLHEP::HepMatrix m_al_0
See J.Tanaka Ph.D (2001) p136 for definition.
Definition KFitBase.h:257
enum KFitError::ECode doFit1(void)
Perform a fit (used in MassFitKFit::doFit()).
Definition KFitBase.cc:502
static void displayError(const char *file, const int line, const char *func, const enum ECode code)
Display a description of error and its location.
Definition KFitError.h:71
ECode
ECode is a error code enumerate.
Definition KFitError.h:33
@ kUnimplemented
Unprepared.
Definition KFitError.h:43
@ kCannotGetMatrixInverse
Cannot calculate matrix inverse (bad track property or internal error)
Definition KFitError.h:57
@ kOutOfRange
Specified track-id out of range.
Definition KFitError.h:41
@ kDivisionByZero
Division by zero (bad track property or internal error)
Definition KFitError.h:55
@ kBadTrackSize
Track count too small to perform fit.
Definition KFitError.h:46
@ kBadMatrixSize
Wrong correlation matrix size.
Definition KFitError.h:48
@ kBadCorrelationSize
Wrong correlation matrix size (internal error)
Definition KFitError.h:50
MakeMotherKFit is a class to build mother particle from kinematically fitted daughters.
enum KFitError::ECode setVertex(const HepPoint3D &v)
Set a vertex position of the mother particle.
enum KFitError::ECode addTrack(const KFitTrack &kp)
Add a track to the make-mother object.
enum KFitError::ECode doMake(void)
Perform a reconstruction of mother particle.
const CLHEP::HepSymMatrix getMotherError(void) const
Get an error matrix of the mother particle.
enum KFitError::ECode setCorrelation(const CLHEP::HepMatrix &e)
Set a correlation matrix.
const HepPoint3D getMotherPosition(void) const
Get a position of the mother particle.
enum KFitError::ECode setVertexError(const CLHEP::HepSymMatrix &e)
Set a vertex error matrix of the mother particle.
enum KFitError::ECode setTrackVertexError(const CLHEP::HepMatrix &e)
Set a vertex error matrix of the child particle in the addTrack'ed order.
const CLHEP::HepLorentzVector getMotherMomentum(void) const
Get a Lorentz vector of the mother particle.
enum KFitError::ECode setMagneticField(const double mf)
Change a magnetic field from the default value KFitConst::kDefaultMagneticField.
enum KFitError::ECode setZeroCorrelation(void) override
Indicate no correlation between tracks.
enum KFitError::ECode setVertex(const HepPoint3D &v)
Set an initial vertex position for the four momentum-constraint fit.
bool m_FlagTrackVertexError
Flag to indicate if the vertex error matrix of the child particle is preset.
enum KFitError::ECode setFourMomentum(const ROOT::Math::PxPyPzEVector &m)
Set a recoil mass .
bool m_FlagAtDecayPoint
Flag controlled by setFlagAtDecayPoint().
bool getFlagAtDecayPoint(void) const
Get a flag if to constraint at the decay point in the four momentum-constraint fit.
enum KFitError::ECode prepareInputMatrix(void) override
Build grand matrices for minimum search from input-track properties.
bool getFlagFitWithVertex(void) const
Get a flag if the fit is allowed with moving the vertex position.
enum KFitError::ECode calculateNDF(void) override
Calculate an NDF of the fit.
double getCHIsq(void) const override
Get a chi-square of the fit.
enum KFitError::ECode setFlagAtDecayPoint(const bool flag)
Set a flag if to constraint at the decay point in the four momentum-constraint fit.
enum KFitError::ECode setTrackZeroVertexError(void)
Indicate no vertex uncertainty in the child particle in the addTrack'ed order.
std::vector< int > m_IsFixMass
Array of flags whether the track property is fixed at the mass.
enum KFitError::ECode updateMother(Particle *mother)
Update mother particle.
enum KFitError::ECode unfixMass(void)
Tell the object to unfix the last added track property at the invariant mass.
enum KFitError::ECode setRecoilMass(const double m)
Set an invariant mass for the four momentum-constraint fit.
enum KFitError::ECode prepareInputSubMatrix(void) override
Build sub-matrices for minimum search from input-track properties.
bool m_FlagFitIncludingVertex
Flag to indicate if the fit is allowed with moving the vertex position.
enum KFitError::ECode doFit(void)
Perform a recoil-mass constraint fit.
std::vector< CLHEP::HepMatrix > m_BeforeTrackVertexError
array of vertex error matrices before the fit.
enum KFitError::ECode makeCoreMatrix(void) override
Build matrices using the kinematical constraint.
enum KFitError::ECode setVertexError(const CLHEP::HepSymMatrix &e)
Set an initial vertex error matrix for the four momentum-constraint fit.
enum KFitError::ECode setTrackVertexError(const CLHEP::HepMatrix &e)
Set a vertex error matrix of the child particle in the addTrack'ed order.
enum KFitError::ECode prepareCorrelation(void) override
Build a grand correlation matrix from input-track properties.
HepPoint3D m_AfterVertex
Vertex position after the fit.
enum KFitError::ECode setCorrelation(const CLHEP::HepMatrix &m) override
Set a correlation matrix.
std::vector< CLHEP::HepMatrix > m_AfterTrackVertexError
array of vertex error matrices after the fit.
const CLHEP::HepMatrix getCorrelation(const int id1, const int id2, const int flag=KFitConst::kAfterFit) const override
Get a correlation matrix between two tracks.
double getTrackCHIsq(const int id) const override
Get a chi-square of the track.
const CLHEP::HepSymMatrix getVertexError(const int flag=KFitConst::kAfterFit) const
Get a vertex error matrix.
const HepPoint3D getVertex(const int flag=KFitConst::kAfterFit) const
Get a vertex position.
const CLHEP::HepMatrix getTrackVertexError(const int id, const int flag=KFitConst::kAfterFit) const
Get a vertex error matrix of the track.
CLHEP::HepSymMatrix m_AfterVertexError
Vertex error matrix after the fit.
enum KFitError::ECode prepareOutputMatrix(void) override
Build an output error matrix.
~RecoilMassKFit(void)
Destruct the object.
CLHEP::HepSymMatrix m_BeforeVertexError
Vertex error matrix before the fit.
HepPoint3D m_BeforeVertex
Vertex position before the fit.
enum KFitError::ECode fixMass(void)
Tell the object to fix the last added track property at the invariant mass.
ROOT::Math::PxPyPzEVector m_FourMomentum
Four Momentum.
RecoilMassKFit(void)
Construct an object with no argument.
double sqrt(double a)
sqrt for double
Definition beamHelpers.h:28
Abstract base class for different kinds of events.
STL namespace.
static const int kMaxTrackCount
Maximum track size.
Definition KFitConst.h:38
static const int kAfterFit
Input parameter to specify after-fit when setting/getting a track attribute.
Definition KFitConst.h:35
static const int kBeforeFit
Input parameter to specify before-fit when setting/getting a track attribute.
Definition KFitConst.h:33
static const int kNumber7
Constant 7 to check matrix size (internal use)
Definition KFitConst.h:30