Belle II Software light-2406-ragdoll
VertexFitKFit Class Reference

VertexFitKFit is a derived class from KFitBase to perform vertex-constraint kinematical fit. More...

#include <VertexFitKFit.h>

Inheritance diagram for VertexFitKFit:
Collaboration diagram for VertexFitKFit:

Public Member Functions

 VertexFitKFit (void)
 Construct an object with no argument.
 
 ~VertexFitKFit (void)
 Destruct the object.
 
enum KFitError::ECode setInitialVertex (const HepPoint3D &v)
 Set an initial vertex point for the vertex-vertex constraint fit.
 
enum KFitError::ECode setInitialVertex (const B2Vector3D &v)
 Set an initial vertex point for the mass-vertex constraint fit.
 
enum KFitError::ECode setIpProfile (const HepPoint3D &ip, const CLHEP::HepSymMatrix &ipe)
 Set an IP-ellipsoid shape for the vertex constraint fit.
 
enum KFitError::ECode setIpTubeProfile (const CLHEP::HepLorentzVector &p, const HepPoint3D &x, const CLHEP::HepSymMatrix &e, const double q)
 Set a virtual IP-tube track for the vertex constraint fit.
 
enum KFitError::ECode setKnownVertex (const bool flag=true)
 Tell the object to perform a fit with vertex position fixed.
 
enum KFitError::ECode setCorrelationMode (const bool m)
 Tell the object to perform a fit with track correlations.
 
const HepPoint3D getVertex (const int flag=KFitConst::kAfterFit) const
 Get a vertex position.
 
const CLHEP::HepSymMatrix getVertexError (void) const
 Get a fitted vertex error matrix.
 
double getCHIsq (void) const override
 Get a chi-square of the fit.
 
double getCHIsqVertex (void) const
 Get a chi-square of the fit excluding IP-constraint part.
 
const CLHEP::HepMatrix getTrackVertexError (const int id) const
 Get a vertex error matrix of the track.
 
double getTrackCHIsq (const int id) const override
 Get a chi-square of the track.
 
double getTrackPartCHIsq (void) const
 Get a sum of the chi-square associated to the input tracks.
 
int getTrackPartNDF (void) const
 Get an NDF relevant to the getTrackPartCHIsq().
 
enum KFitError::ECode doFit (void)
 Perform a vertex-constraint fit.
 
enum KFitError::ECode updateMother (Particle *mother)
 Update mother particle.
 
enum KFitError::ECode addTrack (const KFitTrack &kp)
 Add a track to the fitter object.
 
enum KFitError::ECode addTrack (const CLHEP::HepLorentzVector &p, const HepPoint3D &x, const CLHEP::HepSymMatrix &e, const double q)
 Add a track to the fitter object with specifying its momentum, position, error matrix, and charge.
 
enum KFitError::ECode addParticle (const Particle *particle)
 Add a particle to the fitter.
 
virtual enum KFitError::ECode setCorrelation (const CLHEP::HepMatrix &c)
 Set a correlation matrix.
 
virtual enum KFitError::ECode setZeroCorrelation (void)
 Indicate no correlation between tracks.
 
enum KFitError::ECode setMagneticField (const double mf)
 Change a magnetic field from the default value KFitConst::kDefaultMagneticField.
 
enum KFitError::ECode getErrorCode (void) const
 Get a code of the last error.
 
int getTrackCount (void) const
 Get the number of added tracks.
 
virtual int getNDF (void) const
 Get an NDF of the fit.
 
double getMagneticField (void) const
 Get a magnetic field.
 
const CLHEP::HepLorentzVector getTrackMomentum (const int id) const
 Get a Lorentz vector of the track.
 
const HepPoint3D getTrackPosition (const int id) const
 Get a position of the track.
 
const CLHEP::HepSymMatrix getTrackError (const int id) const
 Get an error matrix of the track.
 
const KFitTrack getTrack (const int id) const
 Get a specified track object.
 
virtual const CLHEP::HepMatrix getCorrelation (const int id1, const int id2, const int flag=KFitConst::kAfterFit) const
 Get a correlation matrix between two tracks.
 
bool isFitted (void) const
 Return false if fit is not performed yet or performed fit is failed; otherwise true.
 

Protected Member Functions

const CLHEP::HepSymMatrix makeError1 (const CLHEP::HepLorentzVector &p, const CLHEP::HepMatrix &e) const
 Rebuild an error matrix from a Lorentz vector and an error matrix.
 
const CLHEP::HepMatrix makeError1 (const CLHEP::HepLorentzVector &p1, const CLHEP::HepLorentzVector &p2, const CLHEP::HepMatrix &e) const
 Rebuild an error matrix from a pair of Lorentz vectors and an error matrix.
 
const CLHEP::HepMatrix makeError2 (const CLHEP::HepLorentzVector &p, const CLHEP::HepMatrix &e) const
 Rebuild an error matrix from a Lorentz vector and an error matrix.
 
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.
 
const CLHEP::HepMatrix makeError3 (const CLHEP::HepLorentzVector &p1, const CLHEP::HepLorentzVector &p2, const CLHEP::HepMatrix &e, const bool is_fix_mass1, const bool is_fix_mass2) const
 Rebuild an error matrix from a pair of Lorentz vectors and an error matrix.
 
const CLHEP::HepMatrix makeError4 (const CLHEP::HepLorentzVector &p, const CLHEP::HepMatrix &e) const
 Rebuild an error matrix from a Lorentz vector and an error matrix.
 
virtual enum KFitError::ECode prepareCorrelation (void)
 Build a grand correlation matrix from input-track properties.
 
enum KFitError::ECode doFit1 (void)
 Perform a fit (used in MassFitKFit::doFit()).
 
enum KFitError::ECode doFit2 (void)
 Perform a fit (used in VertexFitKFit::doFit() and MassVertexFitKFit::doFit()).
 
bool isTrackIDInRange (const int id) const
 Check if the id is in the range.
 
bool isNonZeroEnergy (const CLHEP::HepLorentzVector &p) const
 Check if the energy is non-zero.
 

Protected Attributes

enum KFitError::ECode m_ErrorCode
 Error code.
 
bool m_FlagFitted
 Flag to indicate if the fit is performed and succeeded.
 
std::vector< KFitTrackm_Tracks
 Container of input tracks.
 
std::vector< CLHEP::HepMatrix > m_BeforeCorrelation
 Container of input correlation matrices.
 
CLHEP::HepSymMatrix m_V_al_0
 See J.Tanaka Ph.D (2001) p137 for definition.
 
CLHEP::HepMatrix m_al_0
 See J.Tanaka Ph.D (2001) p136 for definition.
 
CLHEP::HepMatrix m_al_1
 See J.Tanaka Ph.D (2001) p136 for definition.
 
CLHEP::HepMatrix m_al_a
 See J.Tanaka Ph.D (2001) p137 for definition.
 
CLHEP::HepMatrix m_property
 Container of charges and masses.
 
CLHEP::HepMatrix m_D
 See J.Tanaka Ph.D (2001) p137 for definition.
 
CLHEP::HepMatrix m_d
 See J.Tanaka Ph.D (2001) p137 for definition.
 
CLHEP::HepMatrix m_V_D
 See J.Tanaka Ph.D (2001) p138 for definition.
 
CLHEP::HepMatrix m_V_al_1
 See J.Tanaka Ph.D (2001) p138 for definition.
 
CLHEP::HepMatrix m_lam
 See J.Tanaka Ph.D (2001) p137 for definition.
 
CLHEP::HepMatrix m_E
 See J.Tanaka Ph.D (2001) p137 for definition.
 
CLHEP::HepMatrix m_V_E
 See J.Tanaka Ph.D (2001) p138 for definition.
 
CLHEP::HepMatrix m_lam0
 See J.Tanaka Ph.D (2001) p138 for definition.
 
CLHEP::HepMatrix m_v
 See J.Tanaka Ph.D (2001) p137 for definition.
 
CLHEP::HepMatrix m_v_a
 See J.Tanaka Ph.D (2001) p137 for definition.
 
CLHEP::HepMatrix m_V_Dt
 See J.Tanaka Ph.D (2001) p138 for definition.
 
CLHEP::HepMatrix m_Cov_v_al_1
 See J.Tanaka Ph.D (2001) p137 for definition.
 
int m_NDF
 NDF of the fit.
 
double m_CHIsq
 chi-square of the fit.
 
int m_TrackCount
 Number of tracks.
 
int m_NecessaryTrackCount
 Number needed tracks to perform fit.
 
bool m_FlagCorrelation
 Flag whether a correlation among tracks exists.
 
bool m_FlagOverIteration
 Flag whether the iteration count exceeds the limit.
 
double m_MagneticField
 Magnetic field.
 

Private Member Functions

enum KFitError::ECode doFit3 (void)
 Perform a standard vertex-constraint fit including IP-tube constraint.
 
enum KFitError::ECode doFit4 (void)
 Perform a IP-ellipsoid and vertex-constraint fit.
 
enum KFitError::ECode doFit5 (void)
 Perform a fixed-vertex-position fit mainly for slow pion.
 
enum KFitError::ECode prepareInputMatrix (void) override
 Build grand matrices for minimum search from input-track properties.
 
enum KFitError::ECode prepareInputSubMatrix (void) override
 Build sub-matrices for minimum search from input-track properties.
 
enum KFitError::ECode prepareOutputMatrix (void) override
 Build an output error matrix.
 
enum KFitError::ECode makeCoreMatrix (void) override
 Build matrices using the kinematical constraint.
 
enum KFitError::ECode calculateNDF (void) override
 Calculate an NDF of the fit.
 
enum KFitError::ECode appendTube (void)
 Add the virtual tube track to m_Tracks just before the internal minimization call.
 
enum KFitError::ECode deleteTube (void)
 Delete the virtual tube track to m_Tracks just after the internal minimization call.
 

Private Attributes

bool m_CorrelationMode
 Flag controlled by setCorrelationMode().
 
double m_EachCHIsq [KFitConst::kMaxTrackCount2]
 Container of chi-square's of the input tracks.
 
double m_CHIsqVertex
 chi-square of the fit excluding IP-constraint part.
 
HepPoint3D m_BeforeVertex
 Vertex position before the fit.
 
HepPoint3D m_AfterVertex
 Vertex position after the fit.
 
CLHEP::HepSymMatrix m_AfterVertexError
 Vertex error matrix after the fit.
 
std::vector< CLHEP::HepMatrix > m_AfterTrackVertexError
 Array of vertex error matrices after the fit.
 
bool m_FlagBeam
 Flag if to perform IP-ellipsoid constraint fit.
 
CLHEP::HepSymMatrix m_BeamError
 Error matrix modeling the IP ellipsoid.
 
bool m_FlagKnownVertex
 Flag controlled by setKnownVertex().
 
bool m_FlagTube
 Flag if to perform IP-tube constraint fit.
 
int m_iTrackTube
 ID of the virtual tube track in the m_Tracks.
 
KFitTrack m_TubeTrack
 Entity of the virtual IP-tube track.
 

Detailed Description

VertexFitKFit is a derived class from KFitBase to perform vertex-constraint kinematical fit.

Definition at line 34 of file VertexFitKFit.h.

Constructor & Destructor Documentation

◆ VertexFitKFit()

VertexFitKFit ( void  )

Construct an object with no argument.

Definition at line 26 of file VertexFitKFit.cc.

26 :
27 m_BeforeVertex(HepPoint3D(0, 0, 0)),
28 m_AfterVertexError(HepSymMatrix(3, 0)),
29 m_BeamError(HepSymMatrix(3, 0))
30{
31 m_CorrelationMode = false;
32 m_FlagFitted = false;
33 m_FlagKnownVertex = false;
34 m_FlagBeam = false;
35 m_FlagTube = false;
36 m_iTrackTube = -1;
37 m_CHIsqVertex = 0;
39 m_V_E = HepMatrix(3, 3, 0);
40 m_v = HepMatrix(3, 1, 0);
41 m_v_a = HepMatrix(3, 1, 0);
42
44}
int m_NecessaryTrackCount
Number needed tracks to perform fit.
Definition: KFitBase.h:303
CLHEP::HepMatrix m_v_a
See J.Tanaka Ph.D (2001) p137 for definition.
Definition: KFitBase.h:287
CLHEP::HepMatrix m_V_E
See J.Tanaka Ph.D (2001) p138 for definition.
Definition: KFitBase.h:281
bool m_FlagFitted
Flag to indicate if the fit is performed and succeeded.
Definition: KFitBase.h:245
CLHEP::HepMatrix m_v
See J.Tanaka Ph.D (2001) p137 for definition.
Definition: KFitBase.h:285
bool m_FlagTube
Flag if to perform IP-tube constraint fit.
bool m_FlagKnownVertex
Flag controlled by setKnownVertex().
CLHEP::HepSymMatrix m_BeamError
Error matrix modeling the IP ellipsoid.
double m_EachCHIsq[KFitConst::kMaxTrackCount2]
Container of chi-square's of the input tracks.
bool m_FlagBeam
Flag if to perform IP-ellipsoid constraint fit.
int m_iTrackTube
ID of the virtual tube track in the m_Tracks.
double m_CHIsqVertex
chi-square of the fit excluding IP-constraint part.
CLHEP::HepSymMatrix m_AfterVertexError
Vertex error matrix after the fit.
bool m_CorrelationMode
Flag controlled by setCorrelationMode().
HepPoint3D m_BeforeVertex
Vertex position before the fit.
static const int kMaxTrackCount2
Maximum track size (internal use)
Definition: KFitConst.h:42

Member Function Documentation

◆ addParticle()

enum KFitError::ECode addParticle ( const Particle particle)
inherited

Add a particle to the fitter.

The function gets track parameters from the Particle dataobject and calls addTrack().

Parameters
[in]particleParticle.
Returns
error code (zero if success)

Definition at line 59 of file KFitBase.cc.

60{
61 return addTrack(
62 ROOTToCLHEP::getHepLorentzVector(particle->get4Vector()),
63 ROOTToCLHEP::getPoint3D(particle->getVertex()),
64 ROOTToCLHEP::getHepSymMatrix(particle->getMomentumVertexErrorMatrix()),
65 particle->getCharge());
66}
ROOT::Math::XYZVector getVertex() const
Returns vertex position (POCA for charged, IP for neutral FS particles)
Definition: Particle.h:631
double getCharge(void) const
Returns particle charge.
Definition: Particle.cc:622
ROOT::Math::PxPyPzEVector get4Vector() const
Returns Lorentz vector.
Definition: Particle.h:547
TMatrixFSym getMomentumVertexErrorMatrix() const
Returns 7x7 error matrix.
Definition: Particle.cc:420
enum KFitError::ECode addTrack(const KFitTrack &kp)
Add a track to the fitter object.
Definition: KFitBase.cc:38

◆ addTrack() [1/2]

enum KFitError::ECode addTrack ( const CLHEP::HepLorentzVector &  p,
const HepPoint3D &  x,
const CLHEP::HepSymMatrix &  e,
const double  q 
)
inherited

Add a track to the fitter object with specifying its momentum, position, error matrix, and charge.

This function internally calls addTrack(const KFitTrack &kp).

Parameters
pLorentz vector of the track
xposition of the track
e(7x7) error matrix of the track
qcharge of the track
Returns
error code (zero if success)

Definition at line 47 of file KFitBase.cc.

47 {
48 if (e.num_row() != KFitConst::kNumber7)
49 {
51 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
52 return m_ErrorCode;
53 }
54
55 return this->addTrack(KFitTrack(p, x, e, q));
56}
enum KFitError::ECode m_ErrorCode
Error code.
Definition: KFitBase.h:243
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:72
@ kBadMatrixSize
Wrong correlation matrix size.
Definition: KFitError.h:49
KFitTrack is a container of the track information (Lorentz vector, position, and error matrix),...
Definition: KFitTrack.h:38
static const int kNumber7
Constant 7 to check matrix size (internal use)
Definition: KFitConst.h:32

◆ addTrack() [2/2]

enum KFitError::ECode addTrack ( const KFitTrack kp)
inherited

Add a track to the fitter object.

Parameters
kpobject of the track
Returns
error code (zero if success)

Definition at line 38 of file KFitBase.cc.

38 {
39 m_Tracks.push_back(p);
40 m_TrackCount = m_Tracks.size();
41
43}
std::vector< KFitTrack > m_Tracks
Container of input tracks.
Definition: KFitBase.h:249
int m_TrackCount
Number of tracks.
Definition: KFitBase.h:301

◆ appendTube()

enum KFitError::ECode appendTube ( void  )
private

Add the virtual tube track to m_Tracks just before the internal minimization call.

Returns
error code (zero if success)

Definition at line 867 of file VertexFitKFit.cc.

867 {
869
870 if (m_iTrackTube != -1)
871 {
872 char buf[1024];
873 sprintf(buf, "%s:%s(): internal error; duplicated appendTube() call?", __FILE__, __func__);
874 B2FATAL(buf);
875 }
876
877 m_Tracks.push_back(m_TubeTrack);
878 m_TrackCount = m_Tracks.size();
880
882}
KFitTrack m_TubeTrack
Entity of the virtual IP-tube track.

◆ calculateNDF()

enum KFitError::ECode calculateNDF ( void  )
overrideprivatevirtual

Calculate an NDF of the fit.

Returns
error code (zero if success)

Implements KFitBase.

Definition at line 856 of file VertexFitKFit.cc.

856 {
857 if (m_FlagBeam) m_NDF = 2 * m_TrackCount;
858 else if (m_FlagTube) m_NDF = 2 * (m_TrackCount - 1) - 1;
859 else if (m_FlagKnownVertex) m_NDF = 2 * m_TrackCount;
860 else m_NDF = 2 * m_TrackCount - 3;
861
863}
int m_NDF
NDF of the fit.
Definition: KFitBase.h:295

◆ deleteTube()

enum KFitError::ECode deleteTube ( void  )
private

Delete the virtual tube track to m_Tracks just after the internal minimization call.

Returns
error code (zero if success)

Definition at line 886 of file VertexFitKFit.cc.

886 {
888
889 if (m_iTrackTube == -1)
890 {
891 char buf[1024];
892 sprintf(buf, "%s:%s(): internal error; duplicated deleteTube() call?", __FILE__, __func__);
893 B2FATAL(buf);
894 }
895
896 m_Tracks.pop_back();
897 m_TrackCount = m_Tracks.size();
898 m_iTrackTube = -1;
899
901}

◆ doFit()

enum KFitError::ECode doFit ( void  )

Perform a vertex-constraint fit.

Returns
error code (zero if success)

Definition at line 215 of file VertexFitKFit.cc.

215 {
216 if (m_FlagTube) this->appendTube();
217
218 if (m_FlagBeam) m_ErrorCode = this->doFit4();
219 else if (m_FlagKnownVertex) m_ErrorCode = this->doFit5();
220 else if (m_CorrelationMode) m_ErrorCode = this->doFit2();
221 else
222 m_ErrorCode = this->doFit3();
223
224 const enum KFitError::ECode tmp_ErrorCode = m_ErrorCode;
225
226 if (m_FlagTube) this->deleteTube();
227
228 if (tmp_ErrorCode == KFitError::kNoError) m_FlagFitted = true;
229
230 return m_ErrorCode = tmp_ErrorCode;
231}
enum KFitError::ECode doFit2(void)
Perform a fit (used in VertexFitKFit::doFit() and MassVertexFitKFit::doFit()).
Definition: KFitBase.cc:578
ECode
ECode is a error code enumerate.
Definition: KFitError.h:34
enum KFitError::ECode doFit5(void)
Perform a fixed-vertex-position fit mainly for slow pion.
enum KFitError::ECode doFit4(void)
Perform a IP-ellipsoid and vertex-constraint fit.
enum KFitError::ECode deleteTube(void)
Delete the virtual tube track to m_Tracks just after the internal minimization call.
enum KFitError::ECode doFit3(void)
Perform a standard vertex-constraint fit including IP-tube constraint.
enum KFitError::ECode appendTube(void)
Add the virtual tube track to m_Tracks just before the internal minimization call.

◆ doFit1()

enum KFitError::ECode doFit1 ( void  )
protectedinherited

Perform a fit (used in MassFitKFit::doFit()).

Returns
error code (zero if success)

Definition at line 502 of file KFitBase.cc.

502 {
504
506 {
508 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
509 return m_ErrorCode;
510 }
511
514
515
516 double chisq = 0;
517 double tmp_chisq = KFitConst::kInitialCHIsq;
518 int err_inverse = 0;
519
520 HepMatrix tmp_al_1(m_al_1);
521 HepMatrix tmp_V_al_1(m_V_al_1);
522
523 m_al_a = m_al_0;
524 HepMatrix tmp_al_a(m_al_a);
525
526
527 for (int i = 0; i < KFitConst::kMaxIterationCount; i++)
528 {
530
531 m_V_D = (m_V_al_0.similarity(m_D)).inverse(err_inverse);
532 if (err_inverse != 0) {
534 return m_ErrorCode;
535 }
536
537 m_lam = m_V_D * (m_D * (m_al_0 - m_al_1) + m_d);
538 chisq = ((m_lam.T()) * (m_D * (m_al_0 - m_al_1) + m_d))(1, 1);
539 m_al_1 = m_al_0 - m_V_al_0 * (m_D.T()) * m_lam;
540 m_V_al_1 = m_V_al_0 - m_V_al_0 * (m_D.T()) * m_V_D * m_D * m_V_al_0;
541
542 if (tmp_chisq <= chisq) {
543 if (i == 0) {
545 return m_ErrorCode;
546 } else {
547 chisq = tmp_chisq;
548 m_al_1 = tmp_al_1;
549 m_al_a = tmp_al_a;
550 m_V_al_1 = tmp_V_al_1;
551 break;
552 }
553 } else {
554 tmp_chisq = chisq;
555 tmp_al_a = tmp_al_1;
556 tmp_al_1 = m_al_1;
557 tmp_V_al_1 = m_V_al_1;
558 if (i == KFitConst::kMaxIterationCount - 1) {
559 m_al_a = tmp_al_1;
560 m_FlagOverIteration = true;
561 }
562 }
563 }
564
566
568
569 m_CHIsq = chisq;
570
571 m_FlagFitted = true;
572
574}
virtual enum KFitError::ECode prepareInputMatrix(void)=0
Build grand matrices for minimum search from input-track properties.
virtual enum KFitError::ECode prepareOutputMatrix(void)=0
Build an output error matrix.
CLHEP::HepMatrix m_al_1
See J.Tanaka Ph.D (2001) p136 for definition.
Definition: KFitBase.h:259
CLHEP::HepMatrix m_lam
See J.Tanaka Ph.D (2001) p137 for definition.
Definition: KFitBase.h:276
bool m_FlagOverIteration
Flag whether the iteration count exceeds the limit.
Definition: KFitBase.h:308
CLHEP::HepMatrix m_V_al_1
See J.Tanaka Ph.D (2001) p138 for definition.
Definition: KFitBase.h:274
CLHEP::HepMatrix m_d
See J.Tanaka Ph.D (2001) p137 for definition.
Definition: KFitBase.h:268
CLHEP::HepMatrix m_al_a
See J.Tanaka Ph.D (2001) p137 for definition.
Definition: KFitBase.h:261
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
CLHEP::HepSymMatrix m_V_al_0
See J.Tanaka Ph.D (2001) p137 for definition.
Definition: KFitBase.h:255
virtual enum KFitError::ECode makeCoreMatrix(void)=0
Build matrices using the kinematical constraint.
double m_CHIsq
chi-square of the fit.
Definition: KFitBase.h:297
virtual enum KFitError::ECode calculateNDF(void)=0
Calculate an NDF of the fit.
CLHEP::HepMatrix m_al_0
See J.Tanaka Ph.D (2001) p136 for definition.
Definition: KFitBase.h:257
@ kCannotGetMatrixInverse
Cannot calculate matrix inverse (bad track property or internal error)
Definition: KFitError.h:58
@ kBadInitialCHIsq
Bad initial chi-square (internal error)
Definition: KFitError.h:53
@ kBadTrackSize
Track count too small to perform fit.
Definition: KFitError.h:47
static constexpr double kInitialCHIsq
Initial chi-square value (internal use)
Definition: KFitConst.h:48
static const int kMaxIterationCount
Maximum iteration step (internal use)
Definition: KFitConst.h:45

◆ doFit2()

enum KFitError::ECode doFit2 ( void  )
protectedinherited

Perform a fit (used in VertexFitKFit::doFit() and MassVertexFitKFit::doFit()).

Returns
error code (zero if success)

Definition at line 578 of file KFitBase.cc.

578 {
580
582 {
584 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
585 return m_ErrorCode;
586 }
587
590
591
592 double chisq = 0;
593 double tmp2_chisq = KFitConst::kInitialCHIsq;
594 int err_inverse = 0;
595
596 m_al_a = m_al_0;
597 HepMatrix tmp_al_a(m_al_a);
598
599 HepMatrix tmp_D(m_D), tmp_E(m_E);
600 HepMatrix tmp_V_D(m_V_D), tmp_V_E(m_V_E);
601 HepMatrix tmp_lam0(m_lam0), tmp_v_a(m_v_a);
602
603 HepMatrix tmp2_D(m_D), tmp2_E(m_E);
604 HepMatrix tmp2_V_D(m_V_D), tmp2_V_E(m_V_E);
605 HepMatrix tmp2_lam0(m_lam0), tmp2_v_a(m_v_a), tmp2_v(m_v_a);
606
607
608 for (int j = 0; j < KFitConst::kMaxIterationCount; j++) // j'th loop start
609 {
610
611 double tmp_chisq = KFitConst::kInitialCHIsq;
612
613 for (int i = 0; i < KFitConst::kMaxIterationCount; i++) { // i'th loop start
614
617
618 m_V_D = (m_V_al_0.similarity(m_D)).inverse(err_inverse);
619 if (err_inverse) {
621 return m_ErrorCode;
622 }
623
624 m_V_E = ((m_E.T()) * m_V_D * m_E).inverse(err_inverse);
625 if (err_inverse) {
627 return m_ErrorCode;
628 }
629 m_lam0 = m_V_D * (m_D * (m_al_0 - m_al_1) + m_d);
630 chisq = ((m_lam0.T()) * (m_D * (m_al_0 - m_al_1) + m_E * (m_v - m_v_a) + m_d))(1, 1);
631 m_v_a = m_v_a - m_V_E * (m_E.T()) * m_lam0;
632
633 if (tmp_chisq <= chisq) {
634 if (i == 0) {
636 return m_ErrorCode;
637 } else {
638 chisq = tmp_chisq;
639 m_v_a = tmp_v_a;
640 m_V_E = tmp_V_E;
641 m_V_D = tmp_V_D;
642 m_lam0 = tmp_lam0;
643 m_E = tmp_E;
644 m_D = tmp_D;
645 break;
646 }
647 } else {
648 tmp_chisq = chisq;
649 tmp_v_a = m_v_a;
650 tmp_V_E = m_V_E;
651 tmp_V_D = m_V_D;
652 tmp_lam0 = m_lam0;
653 tmp_E = m_E;
654 tmp_D = m_D;
655 if (i == KFitConst::kMaxIterationCount - 1) {
656 m_FlagOverIteration = true;
657 }
658 }
659 } // i'th loop over
660
661
662 m_al_a = m_al_1;
663 m_lam = m_lam0 - m_V_D * m_E * m_V_E * (m_E.T()) * m_lam0;
664 m_al_1 = m_al_0 - m_V_al_0 * (m_D.T()) * m_lam;
665
666 if (j == 0) {
667
668 tmp2_chisq = chisq;
669 tmp2_v_a = m_v_a;
670 tmp2_v = m_v;
671 tmp2_V_E = m_V_E;
672 tmp2_V_D = m_V_D;
673 tmp2_lam0 = m_lam0;
674 tmp2_E = m_E;
675 tmp2_D = m_D;
676 tmp_al_a = m_al_a;
677
678 } else {
679
680 if (tmp2_chisq <= chisq) {
681 chisq = tmp2_chisq;
682 m_v_a = tmp2_v_a;
683 m_v = tmp2_v;
684 m_V_E = tmp2_V_E;
685 m_V_D = tmp2_V_D;
686 m_lam0 = tmp2_lam0;
687 m_E = tmp2_E;
688 m_D = tmp2_D;
689 m_al_a = tmp_al_a;
690 break;
691 } else {
692 tmp2_chisq = chisq;
693 tmp2_v_a = m_v_a;
694 tmp2_v = m_v;
695 tmp2_V_E = m_V_E;
696 tmp2_V_D = m_V_D;
697 tmp2_lam0 = m_lam0;
698 tmp2_E = m_E;
699 tmp2_D = m_D;
700 tmp_al_a = m_al_a;
701 if (j == KFitConst::kMaxIterationCount - 1) {
702 m_FlagOverIteration = true;
703 }
704 }
705 }
706 } // j'th loop over
707
708
710
711 m_lam = m_lam0 - m_V_D * m_E * m_V_E * (m_E.T()) * m_lam0;
712 m_al_1 = m_al_0 - m_V_al_0 * (m_D.T()) * m_lam;
713 m_V_Dt = m_V_D - m_V_D * m_E * m_V_E * (m_E.T()) * m_V_D;
714 m_V_al_1 = m_V_al_0 - m_V_al_0 * (m_D.T()) * m_V_Dt * m_D * m_V_al_0;
715 m_Cov_v_al_1 = -m_V_E * (m_E.T()) * m_V_D * m_D * m_V_al_0;
716
718
719 m_CHIsq = chisq;
720
721 m_FlagFitted = true;
722
724}
CLHEP::HepMatrix m_V_Dt
See J.Tanaka Ph.D (2001) p138 for definition.
Definition: KFitBase.h:289
CLHEP::HepMatrix m_E
See J.Tanaka Ph.D (2001) p137 for definition.
Definition: KFitBase.h:279
virtual enum KFitError::ECode prepareInputSubMatrix(void)=0
Build sub-matrices for minimum search from input-track properties.
CLHEP::HepMatrix m_lam0
See J.Tanaka Ph.D (2001) p138 for definition.
Definition: KFitBase.h:283
CLHEP::HepMatrix m_Cov_v_al_1
See J.Tanaka Ph.D (2001) p137 for definition.
Definition: KFitBase.h:291

◆ doFit3()

enum KFitError::ECode doFit3 ( void  )
private

Perform a standard vertex-constraint fit including IP-tube constraint.

Returns
error code (zero if success)

Definition at line 235 of file VertexFitKFit.cc.

235 {
236 // use small Matrix --> No Correlation
238
240 {
242 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
243 return m_ErrorCode;
244 }
245
248
249
250 double chisq = 0;
251 double tmp2_chisq = KFitConst::kInitialCHIsq;
252 int err_inverse = 0;
253
254 m_al_a = m_al_0;
255 HepMatrix tmp_al_a(m_al_a);
256
257 HepMatrix tmp_D(m_D), tmp_E(m_E);
258 HepMatrix tmp_V_D(m_V_D), tmp_V_E(m_V_E);
259 HepMatrix tmp_lam0(m_lam0), tmp_v_a(m_v_a);
260
261 HepMatrix tmp2_D(m_D), tmp2_E(m_E);
262 HepMatrix tmp2_V_D(m_V_D), tmp2_V_E(m_V_E);
263 HepMatrix tmp2_lam0(m_lam0), tmp2_v_a(m_v_a);
264
265 std::vector<double> tmp_each_chisq(m_TrackCount);
266 std::vector<double> tmp2_each_chisq(m_TrackCount);
267
268 for (int j = 0; j < KFitConst::kMaxIterationCount; j++) // j'th loop start
269 {
270
271 double tmp_chisq = KFitConst::kInitialCHIsq;
272
273 for (int i = 0; i < KFitConst::kMaxIterationCount; i++) { // i'th loop start
274
277
278 HepMatrix tV_Ein(3, 3, 0);
279 chisq = 0;
280
281 for (int k = 0; k < m_TrackCount; k++) { // k'th loop start
282
283 HepMatrix tD = m_D.sub(2 * k + 1, 2 * (k + 1), KFitConst::kNumber6 * k + 1, KFitConst::kNumber6 * (k + 1)); // 2x6
284 HepMatrix tV_D = ((m_V_al_0.sub(KFitConst::kNumber6 * k + 1,
285 (int)(KFitConst::kNumber6 * (k + 1)))).similarity(tD)).inverse(err_inverse); // 2x2
286 if (err_inverse) {
288 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
289 return m_ErrorCode;
290 }
291
292 m_V_D.sub(2 * k + 1, 2 * k + 1, tV_D);
293 HepMatrix tE = m_E.sub(2 * k + 1, 2 * (k + 1), 1, 3); // 2x3
294 tV_Ein += (tE.T()) * tV_D * tE; // 3x3
295 HepMatrix tDeltaAl = (m_al_0 - m_al_1).sub(KFitConst::kNumber6 * k + 1, KFitConst::kNumber6 * (k + 1), 1, 1); // 6x1
296 HepMatrix td = m_d.sub(2 * k + 1, 2 * (k + 1), 1, 1); // 2x1
297 HepMatrix tlam0 = tV_D * (tD * tDeltaAl + td); // 2x2x(2x6x6x1+2x1) = 2x1
298 m_lam0.sub(2 * k + 1, 1, tlam0);
299 m_EachCHIsq[k] = ((tlam0.T()) * (tD * tDeltaAl + tE * (m_v - m_v_a) + td))(1, 1); // 1x2x(2x6x6x1+2x3x3x1+2x1)
300 chisq += m_EachCHIsq[k];
301 } // k'th loop over
302
303 m_V_E = tV_Ein.inverse(err_inverse);
304 if (err_inverse) {
306 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
307 return m_ErrorCode;
308 }
309
310 m_v_a = m_v_a - m_V_E * (m_E.T()) * m_lam0;
311
312 if (tmp_chisq <= chisq) {
313 if (i == 0) {
315 } else {
316 for (int k = 0; k < m_TrackCount; k++) m_EachCHIsq[k] = tmp_each_chisq[k];
317 chisq = tmp_chisq;
318 m_v_a = tmp_v_a;
319 m_V_E = tmp_V_E;
320 m_V_D = tmp_V_D;
321 m_lam0 = tmp_lam0;
322 m_E = tmp_E;
323 m_D = tmp_D;
324 }
325 break;
326 } else {
327 for (int k = 0; k < m_TrackCount; k++) tmp_each_chisq[k] = m_EachCHIsq[k];
328 tmp_chisq = chisq;
329 tmp_v_a = m_v_a;
330 tmp_V_E = m_V_E;
331 tmp_V_D = m_V_D;
332 tmp_lam0 = m_lam0;
333 tmp_E = m_E;
334 tmp_D = m_D;
335 if (i == KFitConst::kMaxIterationCount - 1) {
336 m_FlagOverIteration = true;
337 }
338 }
339 } // i'th loop over
340
341 m_al_a = m_al_1;
342 m_lam = m_lam0 - m_V_D * m_E * m_V_E * (m_E.T()) * m_lam0;
343 m_al_1 = m_al_0 - m_V_al_0 * (m_D.T()) * m_lam;
344
345 if (j == 0) {
346
347 for (int k = 0; k < m_TrackCount; k++) tmp2_each_chisq[k] = m_EachCHIsq[k];
348 tmp2_chisq = chisq;
349 tmp2_v_a = m_v_a;
350 tmp2_V_E = m_V_E;
351 tmp2_V_D = m_V_D;
352 tmp2_lam0 = m_lam0;
353 tmp2_E = m_E;
354 tmp2_D = m_D;
355 tmp_al_a = m_al_a;
356
357 } else {
358
359 if (tmp2_chisq <= chisq) {
360 for (int k = 0; k < m_TrackCount; k++) m_EachCHIsq[k] = tmp2_each_chisq[k];
361 chisq = tmp2_chisq;
362 m_v_a = tmp2_v_a;
363 m_V_E = tmp2_V_E;
364 m_V_D = tmp2_V_D;
365 m_lam0 = tmp2_lam0;
366 m_E = tmp2_E;
367 m_D = tmp2_D;
368 m_al_a = tmp_al_a;
369 break;
370 } else {
371 for (int k = 0; k < m_TrackCount; k++) tmp2_each_chisq[k] = m_EachCHIsq[k];
372 tmp2_chisq = chisq;
373 tmp2_v_a = m_v_a;
374 tmp2_V_E = m_V_E;
375 tmp2_V_D = m_V_D;
376 tmp2_lam0 = m_lam0;
377 tmp2_E = m_E;
378 tmp2_D = m_D;
379 tmp_al_a = m_al_a;
380 }
381 }
382 } // j'th loop over
383
384
386
387 m_lam = m_lam0 - m_V_D * m_E * m_V_E * (m_E.T()) * m_lam0;
388 m_al_1 = m_al_0 - m_V_al_0 * (m_D.T()) * m_lam;
389 m_V_Dt = m_V_D - m_V_D * m_E * m_V_E * (m_E.T()) * m_V_D;
390 m_V_al_1 = m_V_al_0 - m_V_al_0 * (m_D.T()) * m_V_Dt * m_D * m_V_al_0;
391 m_Cov_v_al_1 = -m_V_E * (m_E.T()) * m_V_D * m_D * m_V_al_0;
392
394
395 m_CHIsq = chisq;
396
398}
enum KFitError::ECode prepareInputMatrix(void) override
Build grand matrices for minimum search from input-track properties.
enum KFitError::ECode calculateNDF(void) override
Calculate an NDF of the fit.
enum KFitError::ECode prepareInputSubMatrix(void) override
Build sub-matrices for minimum search from input-track properties.
enum KFitError::ECode makeCoreMatrix(void) override
Build matrices using the kinematical constraint.
enum KFitError::ECode prepareOutputMatrix(void) override
Build an output error matrix.
static const int kNumber6
Constant 6 to check matrix size (internal use)
Definition: KFitConst.h:30

◆ doFit4()

enum KFitError::ECode doFit4 ( void  )
private

Perform a IP-ellipsoid and vertex-constraint fit.

Returns
error code (zero if success)

Definition at line 402 of file VertexFitKFit.cc.

402 {
403 // included beam position constraint (only no correlation)
405
407 {
409 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
410 return m_ErrorCode;
411 }
412
415
416
417 double chisq = 0;
418 double tmp_chisq = KFitConst::kInitialCHIsq;
419 int err_inverse = 0;
420
421 m_al_a = m_al_0;
422 HepMatrix tmp_al_a(m_al_a);
423
424 HepMatrix tmp_D(m_D), tmp_E(m_E);
425 HepMatrix tmp_lam(m_lam);
426
427 // vertex
428 m_v[0][0] = m_BeforeVertex.x();
429 m_v[1][0] = m_BeforeVertex.y();
430 m_v[2][0] = m_BeforeVertex.z();
431
432 std::vector<double> tmp_each_chisq(m_TrackCount);
433 double tmp_vertex_chisq = 1.e+30; // An init-value is not needed but the C++ compiler requires the init-value.
434
435 // to avoid overestimation of vertex-z error.
436 bool it_flag = false;
437
438 for (int i = 0; i < KFitConst::kMaxIterationCount ; i++)
439 {
440
442
443 chisq = 0;
444
445 HepMatrix tV_Dtin = m_V_al_0.similarity(m_D) + m_BeamError.similarity(m_E);
446 HepMatrix tV_Dt = tV_Dtin.inverse(err_inverse);
447 if (err_inverse) {
449 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
450 return m_ErrorCode;
451 }
452 m_lam = tV_Dt * (m_D * (m_al_0 - m_al_1) + m_E * (m_v - m_v_a) + m_d); // (2*nTrk)x1
453 for (int k = 0; k < m_TrackCount; k++) {
454 HepMatrix tD = m_D.sub(2 * k + 1, 2 * (k + 1), KFitConst::kNumber6 * k + 1, KFitConst::kNumber6 * (k + 1)); // 2x6
455 HepMatrix tDeltaAl = (m_al_0 - m_al_1).sub(KFitConst::kNumber6 * k + 1, KFitConst::kNumber6 * (k + 1), 1, 1); // 6x1
456 HepMatrix td = m_d.sub(2 * k + 1, 2 * (k + 1), 1, 1); // 2x1
457 HepMatrix tE = m_E.sub(2 * k + 1, 2 * (k + 1), 1, 3); // 2x3
458 chisq += ((m_lam.sub(2 * k + 1, 2 * (k + 1), 1, 1).T()) * (tD * tDeltaAl + tE * (m_v - m_v_a) + td))(1,
459 1); // 1x2x(2x6x6x1+2x3x3x1+2x1)
460 m_EachCHIsq[k] = (m_lam.sub(2 * k + 1, 2 * (k + 1), 1, 1).T() * tD * m_V_al_0.sub(KFitConst::kNumber6 * k + 1,
461 KFitConst::kNumber6 * (k + 1)) * (tD.T()) * m_lam.sub(2 * k + 1, 2 * (k + 1), 1, 1))(1, 1);
462 }
463
464 m_CHIsqVertex = (m_lam.T() * m_E * m_BeamError * (m_E.T()) * m_lam)(1, 1);
465 m_al_a = m_al_1;
466 m_v_a = m_v - m_BeamError * (m_E.T()) * m_lam;
467 m_al_1 = m_al_0 - m_V_al_0 * (m_D.T()) * m_lam;
468
469 if (tmp_chisq <= chisq && it_flag) {
470 if (i == 0) {
472 } else {
473 for (int k = 0; k < m_TrackCount; k++) m_EachCHIsq[k] = tmp_each_chisq[k];
474 m_CHIsqVertex = tmp_vertex_chisq;
475 chisq = tmp_chisq;
476 m_lam = tmp_lam;
477 m_E = tmp_E;
478 m_D = tmp_D;
479 m_al_a = tmp_al_a;
480 }
481 break;
482 } else {
483 if (tmp_chisq <= chisq) it_flag = true;
484 for (int k = 0; k < m_TrackCount; k++) tmp_each_chisq[k] = m_EachCHIsq[k];
485 tmp_vertex_chisq = m_CHIsqVertex;
486 tmp_chisq = chisq;
487 tmp_lam = m_lam;
488 tmp_E = m_E;
489 tmp_D = m_D;
490 tmp_al_a = m_al_a;
491 if (i == KFitConst::kMaxIterationCount - 1) {
492 m_FlagOverIteration = true;
493 }
494 }
495 }
496
497
499
500 m_al_1 = m_al_0 - m_V_al_0 * (m_D.T()) * m_lam;
501 m_v_a = m_v - m_BeamError * (m_E.T()) * m_lam;
502 HepMatrix tV_Dtin = m_V_al_0.similarity(m_D) + m_BeamError.similarity(m_E);
503 m_V_Dt = tV_Dtin.inverse(err_inverse);
504 if (err_inverse)
505 {
507 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
508 return m_ErrorCode;
509 }
510
511 m_V_al_1 = m_V_al_0 - m_V_al_0 * (m_D.T()) * m_V_Dt * m_D * m_V_al_0;
513 // m_V_v is m_V_E
514 // --> need to replace m_V_E for my implementation.
516
518
519 m_CHIsq = chisq;
520
522}

◆ doFit5()

enum KFitError::ECode doFit5 ( void  )
private

Perform a fixed-vertex-position fit mainly for slow pion.

Returns
error code (zero if success)

Definition at line 526 of file VertexFitKFit.cc.

526 {
527 // known vertex --> do not find vertex. (only no correlation)
529
531 {
533 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
534 return m_ErrorCode;
535 }
536
539
540
541 double chisq = 0;
542 double tmp_chisq = KFitConst::kInitialCHIsq;
543 int err_inverse = 0;
544
545 m_al_a = m_al_0;
546 HepMatrix tmp_al_a(m_al_a);
547
548 HepMatrix tmp_al_0(m_al_1);
549 HepMatrix tmp_V_al_0(m_V_al_1);
550
551 std::vector<double> tmp_each_chisq(m_TrackCount);
552
553 for (int i = 0; i < KFitConst::kMaxIterationCount; i++)
554 {
555
557
558 chisq = 0;
559
560 for (int k = 0; k < m_TrackCount; k++) {
561 HepMatrix tD = m_D.sub(2 * k + 1, 2 * (k + 1), KFitConst::kNumber6 * k + 1, KFitConst::kNumber6 * (k + 1)); // 2x6
562 HepMatrix tV_D = ((m_V_al_0.sub(KFitConst::kNumber6 * k + 1,
563 (int)(KFitConst::kNumber6 * (k + 1)))).similarity(tD)).inverse(err_inverse); // 2x2
564 if (err_inverse) {
566 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
567 return m_ErrorCode;
568 }
569 m_V_D.sub(2 * k + 1, 2 * k + 1, tV_D);
570
571 HepMatrix tDeltaAl = (m_al_0 - m_al_1).sub(KFitConst::kNumber6 * k + 1, KFitConst::kNumber6 * (k + 1), 1, 1); // 6x1
572 HepMatrix td = m_d.sub(2 * k + 1, 2 * (k + 1), 1, 1); // 2x1
573 HepMatrix tlam = tV_D * (tD * tDeltaAl + td); // 2x2x(2x6x6x1+2x1) = 2x1
574 m_lam.sub(2 * k + 1, 1, tlam);
575 m_EachCHIsq[k] = ((tlam.T()) * (tD * tDeltaAl + td))(1, 1); // 1x2x(2x6x6x1+2x1)
576 chisq += m_EachCHIsq[k];
577 }
578
579 m_al_a = m_al_1;
580 m_al_1 = m_al_0 - m_V_al_0 * (m_D.T()) * m_lam;
581 m_V_al_1 = m_V_al_0 - m_V_al_0 * (m_D.T()) * m_V_D * m_D * m_V_al_0;
582
583 if (tmp_chisq <= chisq) {
584 if (i == 0) {
586 } else {
587 for (int k = 0; k < m_TrackCount; k++) m_EachCHIsq[k] = tmp_each_chisq[k];
588 chisq = tmp_chisq;
589 m_al_1 = tmp_al_0;
590 m_V_al_1 = tmp_V_al_0;
591 m_al_a = tmp_al_a;
592 }
593 break;
594 } else {
595 for (int k = 0; k < m_TrackCount; k++) tmp_each_chisq[k] = m_EachCHIsq[k];
596 tmp_chisq = chisq;
597 tmp_al_0 = m_al_1;
598 tmp_V_al_0 = m_V_al_1;
599 tmp_al_a = m_al_a;
600 if (i == KFitConst::kMaxIterationCount - 1) {
601 m_FlagOverIteration = true;
602 }
603 }
604 }
605
606
608
610
611 m_CHIsq = chisq;
612
614}

◆ getCHIsq()

double getCHIsq ( void  ) const
overridevirtual

Get a chi-square of the fit.

Returns
chi-square of the fit

Reimplemented from KFitBase.

Definition at line 139 of file VertexFitKFit.cc.

140{
141 return KFitBase::getCHIsq();
142}
virtual double getCHIsq(void) const
Get a chi-square of the fit.
Definition: KFitBase.cc:121

◆ getCHIsqVertex()

double getCHIsqVertex ( void  ) const

Get a chi-square of the fit excluding IP-constraint part.

Returns
chi-square of the fit excluding IP-constraint part.

Definition at line 146 of file VertexFitKFit.cc.

147{
148 // only m_FlagBeam = 1
149 return m_CHIsqVertex;
150}

◆ getCorrelation()

const HepMatrix getCorrelation ( const int  id1,
const int  id2,
const int  flag = KFitConst::kAfterFit 
) const
virtualinherited

Get a correlation matrix between two tracks.

Parameters
id1first track id
id2second track id
flagKFitConst::kBeforeFit or KFitConst::kAfterFit
Returns
(7x7) correlation matrix

Reimplemented in FourCFitKFit, MassFitKFit, MassFourCFitKFit, MassPointingVertexFitKFit, MassVertexFitKFit, and RecoilMassKFit.

Definition at line 183 of file KFitBase.cc.

184{
185 if (flag == KFitConst::kAfterFit && !isFitted()) return HepMatrix(KFitConst::kNumber7, KFitConst::kNumber7, 0);
186 if (!isTrackIDInRange(id1)) return HepMatrix(KFitConst::kNumber7, KFitConst::kNumber7, 0);
187 if (!isTrackIDInRange(id2)) return HepMatrix(KFitConst::kNumber7, KFitConst::kNumber7, 0);
188
189 switch (flag) {
191 return makeError1(
192 getTrackMomentum(id1),
193 getTrackMomentum(id2),
194 m_V_al_1.sub(KFitConst::kNumber6 * id1 + 1, KFitConst::kNumber6 * (id1 + 1), KFitConst::kNumber6 * id2 + 1,
195 KFitConst::kNumber6 * (id2 + 1))
196 );
197
198 default:
199 if (id1 == id2) {
200
201 return static_cast<HepMatrix>(m_Tracks[id1].getError(KFitConst::kBeforeFit));
202
203 } else {
204 const int idx1 = id1 < id2 ? id1 : id2, idx2 = id1 < id2 ? id2 : id1;
205
206 int index = 0;
207
208 for (int i = 0; i < idx1; i++) index += m_TrackCount - 1 - i;
209 index -= idx1 + 1;
210 index += idx2;
211 if (id1 == idx1)
212 return m_BeforeCorrelation[index + idx2];
213 else
214 return m_BeforeCorrelation[index + idx2].T();
215 }
216 }
217}
const CLHEP::HepSymMatrix makeError1(const CLHEP::HepLorentzVector &p, const CLHEP::HepMatrix &e) const
Rebuild an error matrix from a Lorentz vector and an error matrix.
Definition: KFitBase.cc:221
const CLHEP::HepLorentzVector getTrackMomentum(const int id) const
Get a Lorentz vector of the track.
Definition: KFitBase.cc:154
bool isFitted(void) const
Return false if fit is not performed yet or performed fit is failed; otherwise true.
Definition: KFitBase.cc:728
bool isTrackIDInRange(const int id) const
Check if the id is in the range.
Definition: KFitBase.cc:739
std::vector< CLHEP::HepMatrix > m_BeforeCorrelation
Container of input correlation matrices.
Definition: KFitBase.h:251
static const int kAfterFit
Input parameter to specify after-fit when setting/getting a track attribute.
Definition: KFitConst.h:37
static const int kBeforeFit
Input parameter to specify before-fit when setting/getting a track attribute.
Definition: KFitConst.h:35

◆ getErrorCode()

enum KFitError::ECode getErrorCode ( void  ) const
inherited

Get a code of the last error.

Returns
the last error code

Definition at line 101 of file KFitBase.cc.

101 {
102 return m_ErrorCode;
103}

◆ getMagneticField()

double getMagneticField ( void  ) const
inherited

Get a magnetic field.

Returns
magnetic field

Definition at line 128 of file KFitBase.cc.

129{
130 return m_MagneticField;
131}
double m_MagneticField
Magnetic field.
Definition: KFitBase.h:311

◆ getNDF()

int getNDF ( void  ) const
virtualinherited

Get an NDF of the fit.

Returns
NDF of the fit

Definition at line 114 of file KFitBase.cc.

115{
116 return m_NDF;
117}

◆ getTrack()

const KFitTrack getTrack ( const int  id) const
inherited

Get a specified track object.

Parameters
idtrack id
Returns
specified track object

Definition at line 175 of file KFitBase.cc.

176{
177 if (!isTrackIDInRange(id)) return KFitTrack();
178 return m_Tracks[id];
179}

◆ getTrackCHIsq()

double getTrackCHIsq ( const int  id) const
overridevirtual

Get a chi-square of the track.

Parameters
idtrack id
Returns
chi-square of the track

Reimplemented from KFitBase.

Definition at line 162 of file VertexFitKFit.cc.

163{
164 if (!isTrackIDInRange(id)) return -1;
165
167 return KFitBase::getTrackCHIsq(id);
168 }
169
170 return m_EachCHIsq[id];
171}
virtual double getTrackCHIsq(const int id) const
Get a chi-square of the track.
Definition: KFitBase.cc:135

◆ getTrackCount()

int getTrackCount ( void  ) const
inherited

Get the number of added tracks.

Returns
the number of added tracks

Definition at line 107 of file KFitBase.cc.

108{
109 return m_TrackCount;
110}

◆ getTrackError()

const HepSymMatrix getTrackError ( const int  id) const
inherited

Get an error matrix of the track.

Parameters
idtrack id
Returns
error matrix of the track

Definition at line 168 of file KFitBase.cc.

169{
170 if (!isTrackIDInRange(id)) return HepSymMatrix(KFitConst::kNumber7, 0);
171 return m_Tracks[id].getError();
172}

◆ getTrackMomentum()

const HepLorentzVector getTrackMomentum ( const int  id) const
inherited

Get a Lorentz vector of the track.

Parameters
idtrack id
Returns
Lorentz vector of the track

Definition at line 154 of file KFitBase.cc.

155{
156 if (!isTrackIDInRange(id)) return HepLorentzVector();
157 return m_Tracks[id].getMomentum();
158}

◆ getTrackPartCHIsq()

double getTrackPartCHIsq ( void  ) const

Get a sum of the chi-square associated to the input tracks.

The return value should be the same as the one from getCHIsqVertex().

Returns
sum of the chi-square associated to the input tracks

Definition at line 175 of file VertexFitKFit.cc.

176{
178 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
179 return -1;
180 }
181
182 if (m_TrackCount == 0) {
183 KFitError::displayError(__FILE__, __LINE__, __func__, KFitError::kBadTrackSize);
184 return -1;
185 }
186
187 double chisq = 0.0;
188 for (int i = 0; i < m_TrackCount; i++) {
189 const double i_chisq = this->getTrackCHIsq(i);
190 chisq += i_chisq;
191 }
192
193 return chisq;
194}
double getTrackCHIsq(const int id) const override
Get a chi-square of the track.

◆ getTrackPartNDF()

int getTrackPartNDF ( void  ) const

Get an NDF relevant to the getTrackPartCHIsq().

Returns
NDF relevant to the getTrackPartCHIsq()

Definition at line 198 of file VertexFitKFit.cc.

199{
201 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
202 return 0;
203 }
204
205 if (m_TrackCount == 0) {
206 KFitError::displayError(__FILE__, __LINE__, __func__, KFitError::kBadTrackSize);
207 return 0;
208 }
209
210 return m_TrackCount * 2 - 2;
211}

◆ getTrackPosition()

const HepPoint3D getTrackPosition ( const int  id) const
inherited

Get a position of the track.

Parameters
idtrack id
Returns
position of the track

Definition at line 161 of file KFitBase.cc.

162{
163 if (!isTrackIDInRange(id)) return HepPoint3D();
164 return m_Tracks[id].getPosition();
165}

◆ getTrackVertexError()

const HepMatrix getTrackVertexError ( const int  id) const

Get a vertex error matrix of the track.

Parameters
idtrack id
Returns
vertex error matrix

Definition at line 154 of file VertexFitKFit.cc.

155{
156 if (!isTrackIDInRange(id)) return HepMatrix(3, KFitConst::kNumber7, 0);
157 return m_AfterTrackVertexError[id];
158}
std::vector< CLHEP::HepMatrix > m_AfterTrackVertexError
Array of vertex error matrices after the fit.

◆ getVertex()

const HepPoint3D getVertex ( const int  flag = KFitConst::kAfterFit) const

Get a vertex position.

Parameters
flagKFitConst::kBeforeFit or KFitConst::kAfterFit
Returns
vertex position

Definition at line 113 of file VertexFitKFit.cc.

114{
115 if (flag == KFitConst::kAfterFit && !isFitted()) return HepPoint3D();
116
117 switch (flag) {
119 return m_BeforeVertex;
120
122 return m_AfterVertex;
123
124 default:
125 KFitError::displayError(__FILE__, __LINE__, __func__, KFitError::kOutOfRange);
126 return HepPoint3D();
127 }
128}
@ kOutOfRange
Specified track-id out of range.
Definition: KFitError.h:42
HepPoint3D m_AfterVertex
Vertex position after the fit.

◆ getVertexError()

const HepSymMatrix getVertexError ( void  ) const

Get a fitted vertex error matrix.

Returns
vertex error matrix

Definition at line 132 of file VertexFitKFit.cc.

133{
134 return m_AfterVertexError;
135}

◆ isFitted()

bool isFitted ( void  ) const
inherited

Return false if fit is not performed yet or performed fit is failed; otherwise true.

Returns
see description

Definition at line 728 of file KFitBase.cc.

729{
730 if (m_FlagFitted) return true;
731
732 KFitError::displayError(__FILE__, __LINE__, __func__, KFitError::kNotFittedYet);
733
734 return false;
735}
@ kNotFittedYet
Not fitted yet.
Definition: KFitError.h:39

◆ isNonZeroEnergy()

bool isNonZeroEnergy ( const CLHEP::HepLorentzVector &  p) const
protectedinherited

Check if the energy is non-zero.

Parameters
pLorentz vector
Returns
true for non-zero energy, false for otherwise

Definition at line 750 of file KFitBase.cc.

751{
752 if (p.t() != 0) return true;
753
754 KFitError::displayError(__FILE__, __LINE__, __func__, KFitError::kDivisionByZero);
755
756 return false;
757}
@ kDivisionByZero
Division by zero (bad track property or internal error)
Definition: KFitError.h:56

◆ isTrackIDInRange()

bool isTrackIDInRange ( const int  id) const
protectedinherited

Check if the id is in the range.

Parameters
idtrack id
Returns
true if the id is in the range, false otherwise

Definition at line 739 of file KFitBase.cc.

740{
741 if (0 <= id && id < m_TrackCount) return true;
742
743 KFitError::displayError(__FILE__, __LINE__, __func__, KFitError::kOutOfRange);
744
745 return false;
746}

◆ makeCoreMatrix()

enum KFitError::ECode makeCoreMatrix ( void  )
overrideprivatevirtual

Build matrices using the kinematical constraint.

Returns
error code (zero if success)

Implements KFitBase.

Definition at line 762 of file VertexFitKFit.cc.

762 {
763 // vertex fit
764 for (int i = 0; i < m_TrackCount; i++)
765 {
766 double S, U;
767 double sininv;
768
769 double px = m_al_1[i * KFitConst::kNumber6 + 0][0];
770 double py = m_al_1[i * KFitConst::kNumber6 + 1][0];
771 double pz = m_al_1[i * KFitConst::kNumber6 + 2][0];
772 double x = m_al_1[i * KFitConst::kNumber6 + 3][0];
773 double y = m_al_1[i * KFitConst::kNumber6 + 4][0];
774 double z = m_al_1[i * KFitConst::kNumber6 + 5][0];
775 double a = m_property[i][2];
776
777 double pt = sqrt(px * px + py * py);
778 if (pt == 0) {
780 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
781 return m_ErrorCode;
782 }
783
784 double invPt = 1. / pt;
785 double invPt2 = invPt * invPt;
786 double dlx = m_v_a[0][0] - x;
787 double dly = m_v_a[1][0] - y;
788 double dlz = m_v_a[2][0] - z;
789 double a1 = -dlx * py + dly * px;
790 double a2 = dlx * px + dly * py;
791 double r2d2 = dlx * dlx + dly * dly;
792 double Rx = dlx - 2.*px * a2 * invPt2;
793 double Ry = dly - 2.*py * a2 * invPt2;
794
795 if (a != 0) { // charged
796
797 double B = a * a2 * invPt2;
798 if (fabs(B) > 1) {
800 B2DEBUG(10, "KFitError: Cannot calculate arcsin");
801 //KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
802 return m_ErrorCode;
803 }
804 // sin^(-1)(B)
805 sininv = asin(B);
806 double tmp0 = 1.0 - B * B;
807 if (tmp0 == 0) {
809 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
810 return m_ErrorCode;
811 }
812 // 1/sqrt(1-B^2)
813 double sqrtag = 1.0 / sqrt(tmp0);
814 S = sqrtag * invPt2;
815 U = dlz - pz * sininv / a;
816
817 } else { // neutral
818
819 sininv = 0.0;
820 S = invPt2;
821 U = dlz - pz * a2 * invPt2;
822 }
823
824 // d
825 m_d[i * 2 + 0][0] = a1 - 0.5 * a * r2d2;
826 m_d[i * 2 + 1][0] = U * pt;
827
828 // D
829 m_D[i * 2 + 0][i * KFitConst::kNumber6 + 0] = dly;
830 m_D[i * 2 + 0][i * KFitConst::kNumber6 + 1] = -dlx;
831 m_D[i * 2 + 0][i * KFitConst::kNumber6 + 2] = 0.0;
832 m_D[i * 2 + 0][i * KFitConst::kNumber6 + 3] = py + a * dlx;
833 m_D[i * 2 + 0][i * KFitConst::kNumber6 + 4] = -px + a * dly;
834 m_D[i * 2 + 0][i * KFitConst::kNumber6 + 5] = 0.0;
835 m_D[i * 2 + 1][i * KFitConst::kNumber6 + 0] = -pz * pt * S * Rx + U * px * invPt;
836 m_D[i * 2 + 1][i * KFitConst::kNumber6 + 1] = -pz * pt * S * Ry + U * py * invPt;
837 m_D[i * 2 + 1][i * KFitConst::kNumber6 + 2] = a != 0 ? -sininv * pt / a : -a2 * invPt;
838 m_D[i * 2 + 1][i * KFitConst::kNumber6 + 3] = px * pz * pt * S;
839 m_D[i * 2 + 1][i * KFitConst::kNumber6 + 4] = py * pz * pt * S;
840 m_D[i * 2 + 1][i * KFitConst::kNumber6 + 5] = -pt;
841
842 // E
843 m_E[i * 2 + 0][0] = -py - a * dlx;
844 m_E[i * 2 + 0][1] = px - a * dly;
845 m_E[i * 2 + 0][2] = 0.0;
846 m_E[i * 2 + 1][0] = -px * pz * pt * S;
847 m_E[i * 2 + 1][1] = -py * pz * pt * S;
848 m_E[i * 2 + 1][2] = pt;
849 }
850
852}
CLHEP::HepMatrix m_property
Container of charges and masses.
Definition: KFitBase.h:263
@ kCannotGetARCSIN
Cannot get arcsin (bad track property or internal error)
Definition: KFitError.h:60

◆ makeError1() [1/2]

const HepSymMatrix makeError1 ( const CLHEP::HepLorentzVector &  p,
const CLHEP::HepMatrix &  e 
) const
protectedinherited

Rebuild an error matrix from a Lorentz vector and an error matrix.

Parameters
pLorentz vector
e(6x6) error matrix
Returns
(7x7) error matrix

Definition at line 221 of file KFitBase.cc.

222{
223 // self track
224 // Error(6x6,e) ==> Error(7x7,output(hsm)) using Momentum(p).
225
226 if (!isNonZeroEnergy(p)) return HepSymMatrix(KFitConst::kNumber7, 0);
227
228 HepSymMatrix hsm(KFitConst::kNumber7, 0);
229
230 for (int i = 0; i < 3; i++) for (int j = i; j < 3; j++) {
231 hsm[i][j] = e[i][j];
232 hsm[4 + i][4 + j] = e[3 + i][3 + j];
233 }
234 for (int i = 0; i < 3; i++) for (int j = 0; j < 3; j++) {
235 hsm[i][4 + j] = e[i][3 + j];
236 }
237
238 const double invE = 1 / p.t();
239 hsm[0][3] = (p.x() * hsm[0][0] + p.y() * hsm[0][1] + p.z() * hsm[0][2]) * invE;
240 hsm[1][3] = (p.x() * hsm[0][1] + p.y() * hsm[1][1] + p.z() * hsm[1][2]) * invE;
241 hsm[2][3] = (p.x() * hsm[0][2] + p.y() * hsm[1][2] + p.z() * hsm[2][2]) * invE;
242 hsm[3][3] = (p.x() * p.x() * hsm[0][0] + p.y() * p.y() * hsm[1][1] + p.z() * p.z() * hsm[2][2]
243 + 2.0 * p.x() * p.y() * hsm[0][1]
244 + 2.0 * p.x() * p.z() * hsm[0][2]
245 + 2.0 * p.y() * p.z() * hsm[1][2]) * invE * invE;
246 hsm[3][4] = (p.x() * hsm[0][4] + p.y() * hsm[1][4] + p.z() * hsm[2][4]) * invE;
247 hsm[3][5] = (p.x() * hsm[0][5] + p.y() * hsm[1][5] + p.z() * hsm[2][5]) * invE;
248 hsm[3][6] = (p.x() * hsm[0][6] + p.y() * hsm[1][6] + p.z() * hsm[2][6]) * invE;
249
250 return hsm;
251}
bool isNonZeroEnergy(const CLHEP::HepLorentzVector &p) const
Check if the energy is non-zero.
Definition: KFitBase.cc:750

◆ makeError1() [2/2]

const HepMatrix makeError1 ( const CLHEP::HepLorentzVector &  p1,
const CLHEP::HepLorentzVector &  p2,
const CLHEP::HepMatrix &  e 
) const
protectedinherited

Rebuild an error matrix from a pair of Lorentz vectors and an error matrix.

Parameters
p1first Lorentz vector
p2second Lorentz vector
e(6x6) error matrix
Returns
(7x7) error matrix

Definition at line 255 of file KFitBase.cc.

256{
257 // track and track
258 // Error(6x6,e) ==> Error(7x7,output(hm)) using Momentum(p1&p2).
259
260 if (!isNonZeroEnergy(p1)) return HepSymMatrix(KFitConst::kNumber7, 0);
261 if (!isNonZeroEnergy(p2)) return HepSymMatrix(KFitConst::kNumber7, 0);
262
264
265 for (int i = 0; i < 3; i++) for (int j = 0; j < 3; j++) {
266 hm[i][j] = e[i][j];
267 hm[4 + i][4 + j] = e[3 + i][3 + j];
268 hm[4 + i][j] = e[3 + i][j];
269 hm[i][4 + j] = e[i][3 + j];
270 }
271
272 const double invE1 = 1 / p1.t();
273 const double invE2 = 1 / p2.t();
274 hm[0][3] = (p2.x() * hm[0][0] + p2.y() * hm[0][1] + p2.z() * hm[0][2]) * invE2;
275 hm[1][3] = (p2.x() * hm[1][0] + p2.y() * hm[1][1] + p2.z() * hm[1][2]) * invE2;
276 hm[2][3] = (p2.x() * hm[2][0] + p2.y() * hm[2][1] + p2.z() * hm[2][2]) * invE2;
277 hm[4][3] = (p2.x() * hm[4][0] + p2.y() * hm[4][1] + p2.z() * hm[4][2]) * invE2;
278 hm[5][3] = (p2.x() * hm[5][0] + p2.y() * hm[5][1] + p2.z() * hm[5][2]) * invE2;
279 hm[6][3] = (p2.x() * hm[6][0] + p2.y() * hm[6][1] + p2.z() * hm[6][2]) * invE2;
280 hm[3][3] = (p1.x() * p2.x() * hm[0][0] + p1.y() * p2.y() * hm[1][1] + p1.z() * p2.z() * hm[2][2] +
281 p1.x() * p2.y() * hm[0][1] + p2.x() * p1.y() * hm[1][0] +
282 p1.x() * p2.z() * hm[0][2] + p2.x() * p1.z() * hm[2][0] +
283 p1.y() * p2.z() * hm[1][2] + p2.y() * p1.z() * hm[2][1]) * invE1 * invE2;
284 hm[3][0] = (p1.x() * hm[0][0] + p1.y() * hm[1][0] + p1.z() * hm[2][0]) * invE1;
285 hm[3][1] = (p1.x() * hm[0][1] + p1.y() * hm[1][1] + p1.z() * hm[2][1]) * invE1;
286 hm[3][2] = (p1.x() * hm[0][2] + p1.y() * hm[1][2] + p1.z() * hm[2][2]) * invE1;
287 hm[3][4] = (p1.x() * hm[0][4] + p1.y() * hm[1][4] + p1.z() * hm[2][4]) * invE1;
288 hm[3][5] = (p1.x() * hm[0][5] + p1.y() * hm[1][5] + p1.z() * hm[2][5]) * invE1;
289 hm[3][6] = (p1.x() * hm[0][6] + p1.y() * hm[1][6] + p1.z() * hm[2][6]) * invE1;
290
291 return hm;
292}

◆ makeError2()

const HepMatrix makeError2 ( const CLHEP::HepLorentzVector &  p,
const CLHEP::HepMatrix &  e 
) const
protectedinherited

Rebuild an error matrix from a Lorentz vector and an error matrix.

Parameters
pLorentz vector
e(3x6) error matrix
Returns
(3x7) error matrix

Definition at line 296 of file KFitBase.cc.

297{
298 // vertex and track
299 // Error(3x6,e) ==> Error(3x7,output(hm)) using Momentum(p).
300
301 if (!isNonZeroEnergy(p)) return HepSymMatrix(KFitConst::kNumber7, 0);
302
303 HepMatrix hm(3, KFitConst::kNumber7, 0);
304
305 for (int i = 0; i < 3; i++) for (int j = 0; j < 3; j++) {
306 hm[i][j] = e[i][j];
307 hm[i][4 + j] = e[i][3 + j];
308 }
309
310 const double invE = 1 / p.t();
311 hm[0][3] = (p.x() * hm[0][0] + p.y() * hm[0][1] + p.z() * hm[0][2]) * invE;
312 hm[1][3] = (p.x() * hm[1][0] + p.y() * hm[1][1] + p.z() * hm[1][2]) * invE;
313 hm[2][3] = (p.x() * hm[2][0] + p.y() * hm[2][1] + p.z() * hm[2][2]) * invE;
314
315 return hm;
316}

◆ makeError3() [1/2]

const HepSymMatrix makeError3 ( const CLHEP::HepLorentzVector &  p,
const CLHEP::HepMatrix &  e,
const bool  is_fix_mass 
) const
protectedinherited

Rebuild an error matrix from a Lorentz vector and an error matrix.

Parameters
pLorentz vector
e(7x7) error matrix
is_fix_masstrue to recalculate energy term from other parameters, false to do nothing
Returns
(7x7) error matrix

Definition at line 320 of file KFitBase.cc.

321{
322 // self track
323 // Error(7x7,e) ==> Error(7x7,output(hsm)) using Momentum(p).
324 // is_fix_mass = 1 : Energy term is recalculated from the other parameters.
325 // is_fix_mass = 0 : hsm = e.
326
327 if (!isNonZeroEnergy(p)) return HepSymMatrix(KFitConst::kNumber7, 0);
328
329 if (!is_fix_mass) {
330 HepSymMatrix hsm(KFitConst::kNumber7, 0);
331 for (int i = 0; i < 7; i++) for (int j = i; j < 7; j++) {
332 hsm[i][j] = e[i][j];
333 }
334 return hsm;
335 }
336
337 HepSymMatrix hsm(KFitConst::kNumber7, 0);
338
339 for (int i = 0; i < 7; i++) {
340 if (i != 3)
341 for (int j = i; j < 7; j++) hsm[i][j] = e[i][j];
342 }
343
344 double invE = 1 / p.t();
345 hsm[0][3] = (p.x() * hsm[0][0] + p.y() * hsm[0][1] + p.z() * hsm[0][2]) * invE;
346 hsm[1][3] = (p.x() * hsm[0][1] + p.y() * hsm[1][1] + p.z() * hsm[1][2]) * invE;
347 hsm[2][3] = (p.x() * hsm[0][2] + p.y() * hsm[1][2] + p.z() * hsm[2][2]) * invE;
348 hsm[3][3] = (p.x() * p.x() * hsm[0][0] + p.y() * p.y() * hsm[1][1] + p.z() * p.z() * hsm[2][2]
349 + 2.0 * p.x() * p.y() * hsm[0][1]
350 + 2.0 * p.x() * p.z() * hsm[0][2]
351 + 2.0 * p.y() * p.z() * hsm[1][2]) * invE * invE;
352 hsm[3][4] = (p.x() * hsm[0][4] + p.y() * hsm[1][4] + p.z() * hsm[2][4]) * invE;
353 hsm[3][5] = (p.x() * hsm[0][5] + p.y() * hsm[1][5] + p.z() * hsm[2][5]) * invE;
354 hsm[3][6] = (p.x() * hsm[0][6] + p.y() * hsm[1][6] + p.z() * hsm[2][6]) * invE;
355
356 return hsm;
357}

◆ makeError3() [2/2]

const HepMatrix makeError3 ( const CLHEP::HepLorentzVector &  p1,
const CLHEP::HepLorentzVector &  p2,
const CLHEP::HepMatrix &  e,
const bool  is_fix_mass1,
const bool  is_fix_mass2 
) const
protectedinherited

Rebuild an error matrix from a pair of Lorentz vectors and an error matrix.

Parameters
p1first Lorentz vector
p2second Lorentz vector
e(7x7) error matrix
is_fix_mass1true to recalculate energy term from other parameters, false to do nothing
is_fix_mass2true to recalculate energy term from other parameters, false to do nothing
Returns
(7x7) error matrix

Definition at line 361 of file KFitBase.cc.

364{
365 // track and track
366 // Error(7x7,e) ==> Error(7x7,output(hm)) using Momentum(p1&p2).
367 // is_fix_mass = 1 : Energy term is recalculated from the other parameters.
368 // is_fix_mass = 0 : not.
369
370 if (is_fix_mass1 && is_fix_mass2) {
371 if (!isNonZeroEnergy(p1)) return HepSymMatrix(KFitConst::kNumber7, 0);
372 if (!isNonZeroEnergy(p2)) return HepSymMatrix(KFitConst::kNumber7, 0);
373
374 HepMatrix hm(e);
375
376 const double invE1 = 1 / p1.t();
377 const double invE2 = 1 / p2.t();
378 hm[0][3] = (p2.x() * hm[0][0] + p2.y() * hm[0][1] + p2.z() * hm[0][2]) * invE2;
379 hm[1][3] = (p2.x() * hm[1][0] + p2.y() * hm[1][1] + p2.z() * hm[1][2]) * invE2;
380 hm[2][3] = (p2.x() * hm[2][0] + p2.y() * hm[2][1] + p2.z() * hm[2][2]) * invE2;
381 hm[4][3] = (p2.x() * hm[4][0] + p2.y() * hm[4][1] + p2.z() * hm[4][2]) * invE2;
382 hm[5][3] = (p2.x() * hm[5][0] + p2.y() * hm[5][1] + p2.z() * hm[5][2]) * invE2;
383 hm[6][3] = (p2.x() * hm[6][0] + p2.y() * hm[6][1] + p2.z() * hm[6][2]) * invE2;
384 hm[3][0] = (p1.x() * hm[0][0] + p1.y() * hm[1][0] + p1.z() * hm[2][0]) * invE1;
385 hm[3][1] = (p1.x() * hm[0][1] + p1.y() * hm[1][1] + p1.z() * hm[2][1]) * invE1;
386 hm[3][2] = (p1.x() * hm[0][2] + p1.y() * hm[1][2] + p1.z() * hm[2][2]) * invE1;
387 hm[3][3] = (p1.x() * p2.x() * hm[0][0] + p1.y() * p2.y() * hm[1][1] + p1.z() * p2.z() * hm[2][2] +
388 p1.x() * p2.y() * hm[0][1] + p2.x() * p1.y() * hm[1][0] +
389 p1.x() * p2.z() * hm[0][2] + p2.x() * p1.z() * hm[2][0] +
390 p1.y() * p2.z() * hm[1][2] + p2.y() * p1.z() * hm[2][1]) * invE1 * invE2;
391 hm[3][4] = (p1.x() * hm[0][4] + p1.y() * hm[1][4] + p1.z() * hm[2][4]) * invE1;
392 hm[3][5] = (p1.x() * hm[0][5] + p1.y() * hm[1][5] + p1.z() * hm[2][5]) * invE1;
393 hm[3][6] = (p1.x() * hm[0][6] + p1.y() * hm[1][6] + p1.z() * hm[2][6]) * invE1;
394
395 return hm;
396 }
397
398
399 if (is_fix_mass1 && !is_fix_mass2) {
400 if (!isNonZeroEnergy(p1)) return HepSymMatrix(KFitConst::kNumber7, 0);
401
402 HepMatrix hm(e);
403
404 const double invE1 = 1 / p1.t();
405 hm[3][0] = (p1.x() * hm[0][0] + p1.y() * hm[1][0] + p1.z() * hm[2][0]) * invE1;
406 hm[3][1] = (p1.x() * hm[0][1] + p1.y() * hm[1][1] + p1.z() * hm[2][1]) * invE1;
407 hm[3][2] = (p1.x() * hm[0][2] + p1.y() * hm[1][2] + p1.z() * hm[2][2]) * invE1;
408 hm[3][3] = (p1.x() * hm[0][3] + p1.y() * hm[1][3] + p1.z() * hm[2][3]) * invE1;
409 hm[3][4] = (p1.x() * hm[0][4] + p1.y() * hm[1][4] + p1.z() * hm[2][4]) * invE1;
410 hm[3][5] = (p1.x() * hm[0][5] + p1.y() * hm[1][5] + p1.z() * hm[2][5]) * invE1;
411 hm[3][6] = (p1.x() * hm[0][6] + p1.y() * hm[1][6] + p1.z() * hm[2][6]) * invE1;
412
413 return hm;
414 }
415
416
417 if (!is_fix_mass1 && is_fix_mass2) {
418 if (!isNonZeroEnergy(p2)) return HepSymMatrix(KFitConst::kNumber7, 0);
419
420 HepMatrix hm(e);
421
422 const double invE2 = 1 / p2.t();
423 hm[0][3] = (p2.x() * hm[0][0] + p2.y() * hm[0][1] + p2.z() * hm[0][2]) * invE2;
424 hm[1][3] = (p2.x() * hm[1][0] + p2.y() * hm[1][1] + p2.z() * hm[1][2]) * invE2;
425 hm[2][3] = (p2.x() * hm[2][0] + p2.y() * hm[2][1] + p2.z() * hm[2][2]) * invE2;
426 hm[3][3] = (p2.x() * hm[3][0] + p2.y() * hm[3][1] + p2.z() * hm[3][2]) * invE2;
427 hm[4][3] = (p2.x() * hm[4][0] + p2.y() * hm[4][1] + p2.z() * hm[4][2]) * invE2;
428 hm[5][3] = (p2.x() * hm[5][0] + p2.y() * hm[5][1] + p2.z() * hm[5][2]) * invE2;
429 hm[6][3] = (p2.x() * hm[6][0] + p2.y() * hm[6][1] + p2.z() * hm[6][2]) * invE2;
430
431 return hm;
432 }
433
434 return e;
435}

◆ makeError4()

const HepMatrix makeError4 ( const CLHEP::HepLorentzVector &  p,
const CLHEP::HepMatrix &  e 
) const
protectedinherited

Rebuild an error matrix from a Lorentz vector and an error matrix.

Parameters
pLorentz vector
e(3x7) error matrix
Returns
(7x7) error matrix

Definition at line 439 of file KFitBase.cc.

440{
441 // vertex and track
442 // Error(3x7,e) ==> Error(3x7,output(hm)) using Momentum(p).
443 // Energy term is recalculated from the other parameters.
444
445 if (!isNonZeroEnergy(p)) return HepSymMatrix(KFitConst::kNumber7, 0);
446
447 HepMatrix hm(e);
448
449 const double invE = 1 / p.t();
450 hm[0][3] = (p.x() * hm[0][0] + p.y() * hm[0][1] + p.z() * hm[0][2]) * invE;
451 hm[1][3] = (p.x() * hm[1][0] + p.y() * hm[1][1] + p.z() * hm[1][2]) * invE;
452 hm[2][3] = (p.x() * hm[2][0] + p.y() * hm[2][1] + p.z() * hm[2][2]) * invE;
453
454 return hm;
455}

◆ prepareCorrelation()

enum KFitError::ECode prepareCorrelation ( void  )
protectedvirtualinherited

Build a grand correlation matrix from input-track properties.

Returns
error code (zero if success)

Reimplemented in FourCFitKFit, MassFitKFit, MassFourCFitKFit, MassPointingVertexFitKFit, MassVertexFitKFit, and RecoilMassKFit.

Definition at line 459 of file KFitBase.cc.

459 {
460 if (m_BeforeCorrelation.size() != (double)m_TrackCount * ((double)m_TrackCount - 1)*.5)
461 {
463 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
464 return m_ErrorCode;
465 }
466
467 HepMatrix tmp_hm(KFitConst::kNumber6, KFitConst::kNumber6, 0);
468 int row = 0, col = 0;
469
470 for (auto& hm : m_BeforeCorrelation)
471 {
472 row++;
473 if (row == m_TrackCount) {
474 col++;
475 row = col + 1;
476 }
477
478 // 7x7 --> 6x6
479 for (int i = 0; i < 3; i++) for (int j = 0; j < 3; j++) {
480 tmp_hm[i][j] = hm[i][j];
481 tmp_hm[3 + i][3 + j] = hm[4 + i][4 + j];
482 tmp_hm[3 + i][j] = hm[4 + i][j];
483 tmp_hm[i][3 + j] = hm[i][4 + j];
484 }
485
486 int ii = 0, jj = 0;
487 for (int i = KFitConst::kNumber6 * row; i < KFitConst::kNumber6 * (row + 1); i++) {
488 for (int j = KFitConst::kNumber6 * col; j < KFitConst::kNumber6 * (col + 1); j++) {
489 m_V_al_0[i][j] = tmp_hm[ii][jj];
490 jj++;
491 }
492 jj = 0;
493 ii++;
494 }
495 }
496
498}
@ kBadCorrelationSize
Wrong correlation matrix size (internal error)
Definition: KFitError.h:51

◆ prepareInputMatrix()

enum KFitError::ECode prepareInputMatrix ( void  )
overrideprivatevirtual

Build grand matrices for minimum search from input-track properties.

Returns
error code (zero if success)

Implements KFitBase.

Definition at line 618 of file VertexFitKFit.cc.

618 {
620 {
623 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
624 return m_ErrorCode;
625 }
626 } else
627 {
630 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
631 return m_ErrorCode;
632 }
633 }
634
635
636 int index = 0;
637 HepMatrix tmp_al_0(KFitConst::kNumber6 * m_TrackCount, 1, 0);
638 HepSymMatrix tmp_V_al_0(KFitConst::kNumber6 * m_TrackCount, 0);
639 HepMatrix tmp_property(m_TrackCount, 3, 0);
640
641
642 for (auto& track : m_Tracks)
643 {
644 // momentum x,y,z and position x,y,z
645 for (int j = 0; j < KFitConst::kNumber6; j++)
646 tmp_al_0[index * KFitConst::kNumber6 + j][0] = track.getFitParameter(j, KFitConst::kBeforeFit);
647 // these error
648 tmp_V_al_0.sub(index * KFitConst::kNumber6 + 1, track.getFitError(KFitConst::kBeforeFit));
649 // charge , mass , a
650 tmp_property[index][0] = track.getCharge();
651 tmp_property[index][1] = track.getMass();
652 const double c = Belle2::Const::speedOfLight * 1e-4;
653 tmp_property[index][2] = -c * m_MagneticField * track.getCharge();
654 index++;
655 }
656
657 // error between tarck and track
658 m_V_al_0 = tmp_V_al_0;
660 {
661 if (m_FlagCorrelation) {
662 this->prepareCorrelation();
664 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
665 return m_ErrorCode;
666 }
667 }
668 }
669
670 // vertex
671 m_v_a[0][0] = m_BeforeVertex.x();
672 m_v_a[1][0] = m_BeforeVertex.y();
673 m_v_a[2][0] = m_BeforeVertex.z();
674
675 // set member matrix
676 m_al_0 = tmp_al_0;
677 m_al_1 = m_al_0;
678 m_property = tmp_property;
679
680 // define size of matrix
683 m_E = m_V_al_1.sub(1, m_TrackCount * 2, 1, 3);
684 m_d = m_V_al_1.sub(1, m_TrackCount * 2, 1, 1);
685 m_V_D = m_V_al_1.sub(1, m_TrackCount * 2, 1, m_TrackCount * 2);
686 m_lam = m_V_al_1.sub(1, m_TrackCount * 2, 1, 1);
687 m_lam0 = m_V_al_1.sub(1, m_TrackCount * 2, 1, 1);
688 m_V_Dt = m_V_al_1.sub(1, m_TrackCount * 2, 1, m_TrackCount * 2);
690
692}
static const double speedOfLight
[cm/ns]
Definition: Const.h:695
bool m_FlagCorrelation
Flag whether a correlation among tracks exists.
Definition: KFitBase.h:306
virtual enum KFitError::ECode prepareCorrelation(void)
Build a grand correlation matrix from input-track properties.
Definition: KFitBase.cc:459
static const int kMaxTrackCount
Maximum track size.
Definition: KFitConst.h:40

◆ prepareInputSubMatrix()

enum KFitError::ECode prepareInputSubMatrix ( void  )
overrideprivatevirtual

Build sub-matrices for minimum search from input-track properties.

Returns
error code (zero if success)

Implements KFitBase.

Definition at line 696 of file VertexFitKFit.cc.

696 {
697 // vertex
698 for (int i = 0; i < 3; i++)
699 {
700 m_v[i][0] = m_v_a[i][0];
701 }
703}

◆ prepareOutputMatrix()

enum KFitError::ECode prepareOutputMatrix ( void  )
overrideprivatevirtual

Build an output error matrix.

Returns
error code (zero if success)

Implements KFitBase.

Definition at line 707 of file VertexFitKFit.cc.

707 {
708 Hep3Vector h3v;
709 unsigned index = 0;
710
711 for (auto& pdata : m_Tracks)
712 {
713 // tracks
714 // momentum
715 h3v.setX(m_al_1[index * KFitConst::kNumber6 + 0][0]);
716 h3v.setY(m_al_1[index * KFitConst::kNumber6 + 1][0]);
717 h3v.setZ(m_al_1[index * KFitConst::kNumber6 + 2][0]);
718 pdata.setMomentum(HepLorentzVector(h3v, sqrt(h3v.mag2() + pdata.getMass()*pdata.getMass())), KFitConst::kAfterFit);
719 // position
720 pdata.setPosition(HepPoint3D(
721 m_al_1[index * KFitConst::kNumber6 + 3][0],
722 m_al_1[index * KFitConst::kNumber6 + 4][0],
724 // error of the tracks
725 pdata.setError(makeError1(pdata.getMomentum(),
726 m_V_al_1.sub(
727 index * KFitConst::kNumber6 + 1,
728 (index + 1)*KFitConst::kNumber6,
729 index * KFitConst::kNumber6 + 1,
730 (index + 1)*KFitConst::kNumber6)),
732 if (m_ErrorCode != KFitError::kNoError) break;
733 index++;
734 }
735
736 // vertex
737 m_AfterVertex.setX(m_v_a[0][0]);
738 m_AfterVertex.setY(m_v_a[1][0]);
739 m_AfterVertex.setZ(m_v_a[2][0]);
740
741 // error of the vertex
742 for (int i = 0; i < 3; i++) for (int j = i; j < 3; j++)
743 {
744 m_AfterVertexError[i][j] = m_V_E[i][j];
745 }
746
747 // error between vertex and tracks
748 for (int i = 0; i < m_TrackCount; i++)
749 {
750 HepMatrix hm(3, KFitConst::kNumber6, 0);
751 for (int j = 0; j < 3; j++) for (int k = 0; k < KFitConst::kNumber6; k++) {
752 hm[j][k] = m_Cov_v_al_1[j][KFitConst::kNumber6 * i + k];
753 }
754 m_AfterTrackVertexError.push_back(makeError2(m_Tracks[i].getMomentum(), hm));
755 }
756
758}
const CLHEP::HepMatrix makeError2(const CLHEP::HepLorentzVector &p, const CLHEP::HepMatrix &e) const
Rebuild an error matrix from a Lorentz vector and an error matrix.
Definition: KFitBase.cc:296

◆ setCorrelation()

enum KFitError::ECode setCorrelation ( const CLHEP::HepMatrix &  c)
virtualinherited

Set a correlation matrix.

Not intended for end user's use.

Parameters
c(7x7) correlation matrix
Returns
error code (zero if success)

Reimplemented in FourCFitKFit, MassFitKFit, MassFourCFitKFit, MassPointingVertexFitKFit, MassVertexFitKFit, and RecoilMassKFit.

Definition at line 70 of file KFitBase.cc.

70 {
71 if (e.num_row() != KFitConst::kNumber7)
72 {
74 KFitError::displayError(__FILE__, __LINE__, __func__, m_ErrorCode);
75 return m_ErrorCode;
76 }
77 m_BeforeCorrelation.push_back(e);
78 m_FlagCorrelation = true;
79
81}

◆ setCorrelationMode()

enum KFitError::ECode setCorrelationMode ( const bool  m)

Tell the object to perform a fit with track correlations.

Parameters
mtrue for with correlation, false for otherwise
Returns
error code (zero if success)

Definition at line 105 of file VertexFitKFit.cc.

105 {
107
109}

◆ setInitialVertex() [1/2]

enum KFitError::ECode setInitialVertex ( const B2Vector3D v)

Set an initial vertex point for the mass-vertex constraint fit.

Parameters
vinitial vertex point
Returns
error code (zero if success)

Definition at line 56 of file VertexFitKFit.cc.

57{
58 m_BeforeVertex = HepPoint3D(v.X(), v.Y(), v.Z());
60 return m_ErrorCode;
61}
DataType Z() const
access variable Z (= .at(2) without boundary check)
Definition: B2Vector3.h:435
DataType X() const
access variable X (= .at(0) without boundary check)
Definition: B2Vector3.h:431
DataType Y() const
access variable Y (= .at(1) without boundary check)
Definition: B2Vector3.h:433

◆ setInitialVertex() [2/2]

enum KFitError::ECode setInitialVertex ( const HepPoint3D &  v)

Set an initial vertex point for the vertex-vertex constraint fit.

Parameters
vinitial vertex point
Returns
error code (zero if success)

Definition at line 50 of file VertexFitKFit.cc.

50 {
52
54}

◆ setIpProfile()

enum KFitError::ECode setIpProfile ( const HepPoint3D &  ip,
const CLHEP::HepSymMatrix &  ipe 
)

Set an IP-ellipsoid shape for the vertex constraint fit.

Parameters
ipIP position
ipeerror matrix of the IP
Returns
error code (zero if success)

Definition at line 64 of file VertexFitKFit.cc.

64 {
65 if (m_FlagTube)
66 {
67 char buf[1024];
68 sprintf(buf, "%s:%s(): already constrained to IPtube", __FILE__, __func__);
69 B2FATAL(buf);
70 }
71
72 m_FlagBeam = true;
73 m_BeforeVertex = ip;
74 m_BeamError = ipe;
75
77}

◆ setIpTubeProfile()

enum KFitError::ECode setIpTubeProfile ( const CLHEP::HepLorentzVector &  p,
const HepPoint3D &  x,
const CLHEP::HepSymMatrix &  e,
const double  q 
)

Set a virtual IP-tube track for the vertex constraint fit.

Parameters
pLorentz vector of the virtual IP-tube track
xIP position
eerror matrix of IP-tube track and IP position
qcharge of the virtual IP-tube track
Returns
error code (zero if success)

Definition at line 81 of file VertexFitKFit.cc.

81 {
82 if (m_FlagBeam)
83 {
84 char buf[1024];
85 sprintf(buf, "%s:%s(): already constrained to IP", __FILE__, __func__);
86 B2FATAL(buf);
87 }
88
89 m_FlagTube = true;
90 m_TubeTrack = KFitTrack(p, x, e, q);
91
93}

◆ setKnownVertex()

enum KFitError::ECode setKnownVertex ( const bool  flag = true)

Tell the object to perform a fit with vertex position fixed.

Parameters
flagtrue for fixed vertex, false for otherwise
Returns
error code (zero if success)

Definition at line 97 of file VertexFitKFit.cc.

97 {
98 m_FlagKnownVertex = flag;
99
101}

◆ setMagneticField()

enum KFitError::ECode setMagneticField ( const double  mf)
inherited

Change a magnetic field from the default value KFitConst::kDefaultMagneticField.

Parameters
mfmagnetic field to set
Returns
error code (zero if success)

Definition at line 93 of file KFitBase.cc.

93 {
94 m_MagneticField = mf;
95
97}

◆ setZeroCorrelation()

enum KFitError::ECode setZeroCorrelation ( void  )
virtualinherited

Indicate no correlation between tracks.

Not intended for end user's use.

Returns
error code (zero if success)

Reimplemented in FourCFitKFit, MassFitKFit, MassFourCFitKFit, MassPointingVertexFitKFit, MassVertexFitKFit, and RecoilMassKFit.

Definition at line 85 of file KFitBase.cc.

85 {
86 HepMatrix zero(KFitConst::kNumber7, KFitConst::kNumber7, 0);
87
88 return this->setCorrelation(zero);
89}
virtual enum KFitError::ECode setCorrelation(const CLHEP::HepMatrix &c)
Set a correlation matrix.
Definition: KFitBase.cc:70

◆ updateMother()

enum KFitError::ECode updateMother ( Particle mother)

Update mother particle.

Parameters
[in]motherMother particle.

Definition at line 903 of file VertexFitKFit.cc.

904{
905 MakeMotherKFit kmm;
907 unsigned n = getTrackCount();
908 for (unsigned i = 0; i < n; ++i) {
910 getTrack(i).getCharge());
912 for (unsigned j = i + 1; j < n; ++j) {
914 }
915 }
916 kmm.setVertex(getVertex());
918 m_ErrorCode = kmm.doMake();
920 return m_ErrorCode;
921 double chi2 = getCHIsq();
922 int ndf = getNDF();
923 double prob = TMath::Prob(chi2, ndf);
924 //
925 bool haschi2 = mother->hasExtraInfo("chiSquared");
926 if (haschi2) {
927 mother->setExtraInfo("chiSquared", chi2);
928 mother->setExtraInfo("ndf", ndf);
929 } else {
930 mother->addExtraInfo("chiSquared", chi2);
931 mother->addExtraInfo("ndf", ndf);
932 }
933
934 mother->updateMomentum(
935 CLHEPToROOT::getLorentzVector(kmm.getMotherMomentum()),
936 CLHEPToROOT::getXYZVector(kmm.getMotherPosition()),
937 CLHEPToROOT::getTMatrixFSym(kmm.getMotherError()),
938 prob);
940 return m_ErrorCode;
941}
void setExtraInfo(const std::string &name, double value)
Sets the user-defined data of given name to the given value.
Definition: Particle.cc:1317
bool hasExtraInfo(const std::string &name) const
Return whether the extra info with the given name is set.
Definition: Particle.cc:1266
void addExtraInfo(const std::string &name, double value)
Sets the user-defined data of given name to the given value.
Definition: Particle.cc:1336
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:386
const CLHEP::HepSymMatrix getTrackError(const int id) const
Get an error matrix of the track.
Definition: KFitBase.cc:168
const HepPoint3D getTrackPosition(const int id) const
Get a position of the track.
Definition: KFitBase.cc:161
virtual int getNDF(void) const
Get an NDF of the fit.
Definition: KFitBase.cc:114
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
const KFitTrack getTrack(const int id) const
Get a specified track object.
Definition: KFitBase.cc:175
int getTrackCount(void) const
Get the number of added tracks.
Definition: KFitBase.cc:107
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.
double getCHIsq(void) const override
Get a chi-square of the fit.
const CLHEP::HepSymMatrix getVertexError(void) const
Get a fitted vertex error matrix.
const HepPoint3D getVertex(const int flag=KFitConst::kAfterFit) const
Get a vertex position.
const CLHEP::HepMatrix getTrackVertexError(const int id) const
Get a vertex error matrix of the track.

Member Data Documentation

◆ m_AfterTrackVertexError

std::vector<CLHEP::HepMatrix> m_AfterTrackVertexError
private

Array of vertex error matrices after the fit.

Definition at line 174 of file VertexFitKFit.h.

◆ m_AfterVertex

HepPoint3D m_AfterVertex
private

Vertex position after the fit.

Definition at line 170 of file VertexFitKFit.h.

◆ m_AfterVertexError

CLHEP::HepSymMatrix m_AfterVertexError
private

Vertex error matrix after the fit.

Definition at line 172 of file VertexFitKFit.h.

◆ m_al_0

CLHEP::HepMatrix m_al_0
protectedinherited

See J.Tanaka Ph.D (2001) p136 for definition.

Definition at line 257 of file KFitBase.h.

◆ m_al_1

CLHEP::HepMatrix m_al_1
protectedinherited

See J.Tanaka Ph.D (2001) p136 for definition.

Definition at line 259 of file KFitBase.h.

◆ m_al_a

CLHEP::HepMatrix m_al_a
protectedinherited

See J.Tanaka Ph.D (2001) p137 for definition.

Definition at line 261 of file KFitBase.h.

◆ m_BeamError

CLHEP::HepSymMatrix m_BeamError
private

Error matrix modeling the IP ellipsoid.

Definition at line 179 of file VertexFitKFit.h.

◆ m_BeforeCorrelation

std::vector<CLHEP::HepMatrix> m_BeforeCorrelation
protectedinherited

Container of input correlation matrices.

Definition at line 251 of file KFitBase.h.

◆ m_BeforeVertex

HepPoint3D m_BeforeVertex
private

Vertex position before the fit.

Definition at line 167 of file VertexFitKFit.h.

◆ m_CHIsq

double m_CHIsq
protectedinherited

chi-square of the fit.

Definition at line 297 of file KFitBase.h.

◆ m_CHIsqVertex

double m_CHIsqVertex
private

chi-square of the fit excluding IP-constraint part.

Definition at line 164 of file VertexFitKFit.h.

◆ m_CorrelationMode

bool m_CorrelationMode
private

Flag controlled by setCorrelationMode().

Definition at line 160 of file VertexFitKFit.h.

◆ m_Cov_v_al_1

CLHEP::HepMatrix m_Cov_v_al_1
protectedinherited

See J.Tanaka Ph.D (2001) p137 for definition.

Definition at line 291 of file KFitBase.h.

◆ m_D

CLHEP::HepMatrix m_D
protectedinherited

See J.Tanaka Ph.D (2001) p137 for definition.

Definition at line 266 of file KFitBase.h.

◆ m_d

CLHEP::HepMatrix m_d
protectedinherited

See J.Tanaka Ph.D (2001) p137 for definition.

Definition at line 268 of file KFitBase.h.

◆ m_E

CLHEP::HepMatrix m_E
protectedinherited

See J.Tanaka Ph.D (2001) p137 for definition.

Definition at line 279 of file KFitBase.h.

◆ m_EachCHIsq

double m_EachCHIsq[KFitConst::kMaxTrackCount2]
private

Container of chi-square's of the input tracks.

Definition at line 162 of file VertexFitKFit.h.

◆ m_ErrorCode

enum KFitError::ECode m_ErrorCode
protectedinherited

Error code.

Definition at line 243 of file KFitBase.h.

◆ m_FlagBeam

bool m_FlagBeam
private

Flag if to perform IP-ellipsoid constraint fit.

Definition at line 177 of file VertexFitKFit.h.

◆ m_FlagCorrelation

bool m_FlagCorrelation
protectedinherited

Flag whether a correlation among tracks exists.

Definition at line 306 of file KFitBase.h.

◆ m_FlagFitted

bool m_FlagFitted
protectedinherited

Flag to indicate if the fit is performed and succeeded.

Definition at line 245 of file KFitBase.h.

◆ m_FlagKnownVertex

bool m_FlagKnownVertex
private

Flag controlled by setKnownVertex().

Definition at line 182 of file VertexFitKFit.h.

◆ m_FlagOverIteration

bool m_FlagOverIteration
protectedinherited

Flag whether the iteration count exceeds the limit.

Definition at line 308 of file KFitBase.h.

◆ m_FlagTube

bool m_FlagTube
private

Flag if to perform IP-tube constraint fit.

Definition at line 185 of file VertexFitKFit.h.

◆ m_iTrackTube

int m_iTrackTube
private

ID of the virtual tube track in the m_Tracks.

Definition at line 187 of file VertexFitKFit.h.

◆ m_lam

CLHEP::HepMatrix m_lam
protectedinherited

See J.Tanaka Ph.D (2001) p137 for definition.

Definition at line 276 of file KFitBase.h.

◆ m_lam0

CLHEP::HepMatrix m_lam0
protectedinherited

See J.Tanaka Ph.D (2001) p138 for definition.

Definition at line 283 of file KFitBase.h.

◆ m_MagneticField

double m_MagneticField
protectedinherited

Magnetic field.

Definition at line 311 of file KFitBase.h.

◆ m_NDF

int m_NDF
protectedinherited

NDF of the fit.

Definition at line 295 of file KFitBase.h.

◆ m_NecessaryTrackCount

int m_NecessaryTrackCount
protectedinherited

Number needed tracks to perform fit.

Definition at line 303 of file KFitBase.h.

◆ m_property

CLHEP::HepMatrix m_property
protectedinherited

Container of charges and masses.

Definition at line 263 of file KFitBase.h.

◆ m_TrackCount

int m_TrackCount
protectedinherited

Number of tracks.

Definition at line 301 of file KFitBase.h.

◆ m_Tracks

std::vector<KFitTrack> m_Tracks
protectedinherited

Container of input tracks.

Definition at line 249 of file KFitBase.h.

◆ m_TubeTrack

KFitTrack m_TubeTrack
private

Entity of the virtual IP-tube track.

Definition at line 189 of file VertexFitKFit.h.

◆ m_v

CLHEP::HepMatrix m_v
protectedinherited

See J.Tanaka Ph.D (2001) p137 for definition.

Definition at line 285 of file KFitBase.h.

◆ m_v_a

CLHEP::HepMatrix m_v_a
protectedinherited

See J.Tanaka Ph.D (2001) p137 for definition.

Definition at line 287 of file KFitBase.h.

◆ m_V_al_0

CLHEP::HepSymMatrix m_V_al_0
protectedinherited

See J.Tanaka Ph.D (2001) p137 for definition.

Definition at line 255 of file KFitBase.h.

◆ m_V_al_1

CLHEP::HepMatrix m_V_al_1
protectedinherited

See J.Tanaka Ph.D (2001) p138 for definition.

Definition at line 274 of file KFitBase.h.

◆ m_V_D

CLHEP::HepMatrix m_V_D
protectedinherited

See J.Tanaka Ph.D (2001) p138 for definition.

Definition at line 271 of file KFitBase.h.

◆ m_V_Dt

CLHEP::HepMatrix m_V_Dt
protectedinherited

See J.Tanaka Ph.D (2001) p138 for definition.

Definition at line 289 of file KFitBase.h.

◆ m_V_E

CLHEP::HepMatrix m_V_E
protectedinherited

See J.Tanaka Ph.D (2001) p138 for definition.

Definition at line 281 of file KFitBase.h.


The documentation for this class was generated from the following files: