BOSS 7.1.1
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | Static Public Member Functions | Friends | List of all members
VertexFit Class Reference

#include <VertexFit.h>

+ Inheritance diagram for VertexFit:

Public Member Functions

 ~VertexFit ()
 
void init ()
 
void AddVertex (int number, VertexParameter vpar, std::vector< int > lis)
 
void AddVertex (int number, VertexParameter vpar, int n1, int n2)
 
void AddVertex (int number, VertexParameter vpar, int n1, int n2, int n3)
 
void AddVertex (int number, VertexParameter vpar, int n1, int n2, int n3, int n4)
 
void AddVertex (int number, VertexParameter vpar, int n1, int n2, int n3, int n4, int n5)
 
void AddVertex (int number, VertexParameter vpar, int n1, int n2, int n3, int n4, int n5, int n6)
 
void AddVertex (int number, VertexParameter vpar, int n1, int n2, int n3, int n4, int n5, int n6, int n7)
 
void AddVertex (int number, VertexParameter vpar, int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8)
 
void AddVertex (int number, VertexParameter vpar, int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8, int n9)
 
void AddVertex (int number, VertexParameter vpar, int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8, int n9, int n10)
 
void AddVertex (int number, VertexParameter vpar, int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8, int n9, int n10, int n11)
 
void AddVertex (int number, VertexParameter vpar, int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8, int n9, int n10, int n11, int n12)
 
void AddBeamFit (int number, VertexParameter vpar, int n)
 
void BuildVirtualParticle (int number)
 
void setIterNumber (const int niter=10)
 
void setChisqCut (const double chicut=1000, const double chiter=1.0e-3)
 
void setMagCorrFactor (const double factor=1.000)
 
double factor () const
 
bool Fit ()
 
bool Fit (int n)
 
bool BeamFit (int n)
 
void Swim (int n)
 
double chisq () const
 
double chisq (int n) const
 
bool pull (int n, int itk, HepVector &p)
 
HepLorentzVector pfit (int n) const
 
HepPoint3D xfit (int n) const
 
HepVector w (int n) const
 
HepSymMatrix Ew (int n) const
 
WTrackParameter wtrk (int n) const
 
HepVector cpu () const
 
HepPoint3D vx (int n) const
 
HepVector Vx (int n) const
 
HepSymMatrix Evx (int n) const
 
double errorVx (int n, int i) const
 
VertexParameter vpar (int n) const
 
WTrackParameter wVirtualTrack (int n) const
 
- Public Member Functions inherited from TrackPool
 TrackPool ()
 
 ~TrackPool ()
 
void AddTrack (const int number, const double mass, const RecMdcTrack *trk)
 
void AddTrack (const int number, const double mass, const RecMdcKalTrack *trk)
 
void AddTrack (const int number, const double mass, const RecEmcShower *trk)
 
void AddTrack (const int number, const WTrackParameter wtrk)
 
void AddMissTrack (const int number, const double mass)
 
void AddMissTrack (const int number, const double mass, const HepLorentzVector p4)
 
void AddMissTrack (const int number, const double mass, const RecEmcShower *trk)
 
void AddMissTrack (const int number, const RecEmcShower *trk)
 
void AddMissTrack (const int number, const HepLorentzVector p4)
 
void AddTrackVertex (const int number, const double mass, const RecEmcShower *trk)
 
std::vector< int > AddList (int n1)
 
std::vector< int > AddList (int n1, int n2)
 
std::vector< int > AddList (int n1, int n2, int n3)
 
std::vector< int > AddList (int n1, int n2, int n3, int n4)
 
std::vector< int > AddList (int n1, int n2, int n3, int n4, int n5)
 
std::vector< int > AddList (int n1, int n2, int n3, int n4, int n5, int n6)
 
std::vector< int > AddList (int n1, int n2, int n3, int n4, int n5, int n6, int n7)
 
std::vector< int > AddList (int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8)
 
std::vector< int > AddList (int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8, int n9)
 
std::vector< int > AddList (int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8, int n9, int n10)
 
std::vector< int > AddList (int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8, int n9, int n10, int n11)
 
std::vector< int > AddList (int n1, int n2, int n3, int n4, int n5, int n6, int n7, int n8, int n9, int n10, int n11, int n12)
 
std::vector< WTrackParameter > wTrackOrigin () const
 
std::vector< WTrackParameter > wTrackInfit () const
 
std::vector< int > wTrackList () const
 
WTrackParameter wTrackOrigin (int n) const
 
WTrackParameter wTrackInfit (int n) const
 
int wTrackList (int n) const
 
int numberWTrack () const
 
std::vector< GammaShape > GammaShapeValue () const
 
std::vector< int > GammaShapeList () const
 
GammaShape GammaShapeValue (int n) const
 
int GammaShapeList (int n) const
 
int numberGammaShape () const
 
void setWTrackOrigin (const int n, const WTrackParameter wtrk)
 
void setWTrackInfit (const int n, const WTrackParameter wtrk)
 
void setWTrackOrigin (const WTrackParameter wtrk)
 
void setWTrackInfit (const WTrackParameter wtrk)
 
void setWTrackList (const int n)
 
void clearWTrackOrigin ()
 
void clearWTrackInfit ()
 
void clearWTrackList ()
 
void clearone ()
 
void cleartwo ()
 
int numberone () const
 
int numbertwo () const
 
vector< int > mappositionA () const
 
vector< int > mappositionB () const
 
vector< int > mapkinematic () const
 
void clearMapkinematic ()
 
void clearMappositionA ()
 
void clearMappositionB ()
 
void setMapkinematic (const int n)
 
void setMappositionA (const int n)
 
void setMappositionB (const int n)
 
void setGammaShape (const int n, const GammaShape gammashape)
 
void setGammaShape (const GammaShape gammashape)
 
void setGammaShapeList (const int n)
 
void clearGammaShape ()
 
void clearGammaShapeList ()
 
void setBeamPosition (const HepPoint3D BeamPosition)
 
void setVBeamPosition (const HepSymMatrix VBeamPosition)
 
HepPoint3D getBeamPosition () const
 
HepSymMatrix getVBeamPosition () const
 

Static Public Member Functions

static VertexFit * instance ()
 

Friends

class VertexFitRefine
 

Detailed Description

Definition at line 18 of file VertexFit.h.

Constructor & Destructor Documentation

◆ ~VertexFit()

VertexFit::~VertexFit ( )

Definition at line 22 of file VertexFit.cxx.

23{
24 //if (m_pointer) delete m_pointer;
25}

Member Function Documentation

◆ AddBeamFit()

void VertexFit::AddBeamFit ( int number,
VertexParameter vpar,
int n )

Definition at line 74 of file VertexFit.cxx.

75{
76 std::vector<int> tlis = AddList(n);
77 VertexConstraints vc;
78 vc.FixedVertexConstraints(tlis);
79 m_vc.push_back(vc);
80 m_vpar_origin.push_back(vpar);
81 m_vpar_infit.push_back(vpar);
82 if ((unsigned int)number != m_vc.size() - 1)
83 std::cout << "wrong kinematic constraints index" << std::endl;
84}
n
const Int_t n
Definition DataBase/tau_mode.c:65
TrackPool::AddList
std::vector< int > AddList(int n1)
Definition TrackPool.cxx:483
VertexConstraints::FixedVertexConstraints
void FixedVertexConstraints(std::vector< int > tlis)
Definition VertexConstraints.cxx:29
VertexFit::vpar
VertexParameter vpar(int n) const
Definition VertexFit.h:89

Referenced by VertexFitRefine::AddBeamFit().

◆ AddVertex() [1/12]

void VertexFit::AddVertex ( int number,
VertexParameter vpar,
int n1,
int n2 )

Definition at line 108 of file VertexFit.cxx.

109{
110 std::vector<int> tlis = AddList(n1, n2);
111 AddVertex(number, vpar, tlis);
112}
n2
int n2
Definition SD0Tag.cxx:55
n1
int n1
Definition SD0Tag.cxx:54
VertexFit::AddVertex
void AddVertex(int number, VertexParameter vpar, std::vector< int > lis)
Definition VertexFit.cxx:89

◆ AddVertex() [2/12]

void VertexFit::AddVertex ( int number,
VertexParameter vpar,
int n1,
int n2,
int n3 )

Definition at line 115 of file VertexFit.cxx.

116{
117 std::vector<int> tlis = AddList(n1, n2, n3);
118 AddVertex(number, vpar, tlis);
119}

◆ AddVertex() [3/12]

void VertexFit::AddVertex ( int number,
VertexParameter vpar,
int n1,
int n2,
int n3,
int n4 )

Definition at line 121 of file VertexFit.cxx.

122{
123 std::vector<int> tlis = AddList(n1, n2, n3, n4);
124 AddVertex(number, vpar, tlis);
125}

◆ AddVertex() [4/12]

void VertexFit::AddVertex ( int number,
VertexParameter vpar,
int n1,
int n2,
int n3,
int n4,
int n5 )

Definition at line 127 of file VertexFit.cxx.

128{
129 std::vector<int> tlis = AddList(n1, n2, n3, n4, n5);
130 AddVertex(number, vpar, tlis);
131}

◆ AddVertex() [5/12]

void VertexFit::AddVertex ( int number,
VertexParameter vpar,
int n1,
int n2,
int n3,
int n4,
int n5,
int n6 )

Definition at line 133 of file VertexFit.cxx.

134{
135 std::vector<int> tlis = AddList(n1, n2, n3, n4, n5, n6);
136 AddVertex(number, vpar, tlis);
137}

◆ AddVertex() [6/12]

void VertexFit::AddVertex ( int number,
VertexParameter vpar,
int n1,
int n2,
int n3,
int n4,
int n5,
int n6,
int n7 )

Definition at line 139 of file VertexFit.cxx.

140{
141 std::vector<int> tlis = AddList(n1, n2, n3, n4, n5, n6, n7);
142 AddVertex(number, vpar, tlis);
143}

◆ AddVertex() [7/12]

void VertexFit::AddVertex ( int number,
VertexParameter vpar,
int n1,
int n2,
int n3,
int n4,
int n5,
int n6,
int n7,
int n8 )

Definition at line 145 of file VertexFit.cxx.

146{
147 std::vector<int> tlis = AddList(n1, n2, n3, n4, n5, n6, n7, n8);
148 AddVertex(number, vpar, tlis);
149}

◆ AddVertex() [8/12]

void VertexFit::AddVertex ( int number,
VertexParameter vpar,
int n1,
int n2,
int n3,
int n4,
int n5,
int n6,
int n7,
int n8,
int n9 )

Definition at line 151 of file VertexFit.cxx.

152{
153 std::vector<int> tlis = AddList(n1, n2, n3, n4, n5, n6, n7, n8, n9);
154 AddVertex(number, vpar, tlis);
155}

◆ AddVertex() [9/12]

void VertexFit::AddVertex ( int number,
VertexParameter vpar,
int n1,
int n2,
int n3,
int n4,
int n5,
int n6,
int n7,
int n8,
int n9,
int n10 )

Definition at line 157 of file VertexFit.cxx.

158{
159 std::vector<int> tlis = AddList(n1, n2, n3, n4, n5, n6, n7, n8, n9, n10);
160 AddVertex(number, vpar, tlis);
161}

◆ AddVertex() [10/12]

void VertexFit::AddVertex ( int number,
VertexParameter vpar,
int n1,
int n2,
int n3,
int n4,
int n5,
int n6,
int n7,
int n8,
int n9,
int n10,
int n11 )

Definition at line 163 of file VertexFit.cxx.

164{
165 std::vector<int> tlis = AddList(n1, n2, n3, n4, n5, n6, n7, n8, n9, n10, n11);
166 AddVertex(number, vpar, tlis);
167}

◆ AddVertex() [11/12]

void VertexFit::AddVertex ( int number,
VertexParameter vpar,
int n1,
int n2,
int n3,
int n4,
int n5,
int n6,
int n7,
int n8,
int n9,
int n10,
int n11,
int n12 )

Definition at line 169 of file VertexFit.cxx.

170{
171 std::vector<int> tlis = AddList(n1, n2, n3, n4, n5, n6, n7, n8, n9, n10, n11, n12);
172 AddVertex(number, vpar, tlis);
173}

◆ AddVertex() [12/12]

void VertexFit::AddVertex ( int number,
VertexParameter vpar,
std::vector< int > lis )

Definition at line 89 of file VertexFit.cxx.

90{
91 VertexConstraints vc;
92 vc.CommonVertexConstraints(tlis);
93 m_vc.push_back(vc);
94 HepVector vx(3, 0);
95 for (unsigned int i = 0; i < tlis.size(); i++)
96 vx += wTrackOrigin(tlis[i]).X();
97 vx = vx/tlis.size();
98 VertexParameter n_vpar = vpar;
99 n_vpar.setVx(vx);
100 m_vpar_origin.push_back(n_vpar);
101 m_vpar_infit.push_back(n_vpar);
102 WTrackParameter vwtrk;
103 m_virtual_wtrk.push_back(vwtrk);
104 if ((unsigned int)number != m_vc.size() - 1)
105 std::cout << "wrong kinematic constraints index" << std::endl;
106}
TrackPool::wTrackOrigin
std::vector< WTrackParameter > wTrackOrigin() const
Definition TrackPool.h:72
VertexConstraints::CommonVertexConstraints
void CommonVertexConstraints(std::vector< int > tlis)
Definition VertexConstraints.cxx:23
VertexFit::vx
HepPoint3D vx(int n) const
Definition VertexFit.h:85
VertexParameter::setVx
void setVx(const HepPoint3D &vx)
Definition VertexParameter.h:54

Referenced by AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), VertexFitRefine::AddVertex(), VertexFitRefine::AddVertex(), VertexFitRefine::AddVertex(), VertexFitRefine::AddVertex(), VertexFitRefine::AddVertex(), VertexFitRefine::AddVertex(), VertexFitRefine::AddVertex(), VertexFitRefine::AddVertex(), VertexFitRefine::AddVertex(), VertexFitRefine::AddVertex(), VertexFitRefine::AddVertex(), VertexFitRefine::AddVertex(), DQAKsInfo::calculate(), BeamParams::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), DQARhopi::execute(), Gam4pikp::execute(), inclks::execute(), incllambda::execute(), JsiLL::execute(), KShortReconstruction::execute(), LambdaReconstruction::execute(), Ppjrhopi::execute(), PrimaryVertex::execute(), Rhopi::execute(), rhopi::execute(), Signal::execute(), TestV::execute(), VertexFitRefine::Fit(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), K3pi::MTotal(), K3pipi0::MTotal(), Kk::MTotal(), Kkpi0::MTotal(), Kkpipi::MTotal(), Kpi::MTotal(), Kpipi0::MTotal(), Kpipi0pi0::MTotal(), Pipi::MTotal(), Pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), utility::SecondaryVFit_Lambda(), utility::vfit(), and utility::vfit().

◆ BeamFit()

bool VertexFit::BeamFit ( int n)

Definition at line 259 of file VertexFit.cxx.

260{
261 bool okfit = false;
262 if (n < 0 || (unsigned int)n >= m_vc.size())
263 return okfit;
264 for (unsigned int i = 0; i < (m_vc[n].Ltrk()).size(); i++)
265 {
266 int itk = (m_vc[n].Ltrk())[i];
267 setWTrackInfit(itk, wTrackOrigin(itk));
268 }
269 m_vpar_infit[n] = m_vpar_origin[n];
270
271 // iteration loop
272 std::vector<double> chisq;
273 chisq.clear();
274 for (int it = 0; it < m_niter; it++)
275 {
276 std::vector<WTrackParameter> wlis;
277 wlis.clear();
278 for (unsigned int i = 0; i < (m_vc[n].Ltrk()).size(); i++)
279 {
280 int itk = (m_vc[n].Ltrk())[i];
281 wlis.push_back(wTrackInfit(itk));
282 }
283 VertexParameter vpar = m_vpar_infit[n];
284 m_vc[n].UpdateConstraints(vpar, wlis);
285
286 fitBeam(n);
287 chisq.push_back(m_chisq[n]);
288 if (it > 0)
289 {
290 double delchi = chisq[it] - chisq[it-1];
291 if(fabs(delchi) < m_chiter)
292 break;
293 }
294 }
295 if(m_chisq[n]>=m_chicut) return okfit;
296 swimBeam(n);
297 okfit = true;
298 return okfit;
299}
TrackPool::wTrackInfit
std::vector< WTrackParameter > wTrackInfit() const
Definition TrackPool.h:73
TrackPool::setWTrackInfit
void setWTrackInfit(const int n, const WTrackParameter wtrk)
Definition TrackPool.h:106
VertexFit::chisq
double chisq() const
Definition VertexFit.h:66

◆ BuildVirtualParticle()

void VertexFit::BuildVirtualParticle ( int number)

Definition at line 619 of file VertexFit.cxx.

620{
621
622 vertexCovMatrix(n);
623 TStopwatch timer;
624 timer.Start();
625
626 HepMatrix A(NTRKPAR, NTRKPAR * m_nvtrk, 0);
627 HepMatrix B(NTRKPAR, NVTXPAR, 0);
628 VertexConstraints vc = m_vc[n];
629 unsigned int NSIZE = vc.Ltrk().size();
630 int charge = 0;
631
632 HepMatrix Ai(6, 6, 0);
633 Ai[0][0] = 1.0;
634 Ai[1][1] = 1.0;
635 Ai[2][2] = 1.0;
636 HepMatrix Bi(6, 3, 0);
637 Bi[3][0] = 1.0;
638 Bi[4][1] = 1.0;
639 Bi[5][2] = 1.0;
640 HepVector w1(6, 0);
641 HepVector w2(7, 0);
642 HepSymMatrix Ew(7, 0);
643 HepMatrix ew1(6, 6, 0);
644 HepMatrix ew2(7, 7, 0);
645 double totalE = 0;
646
647 for(unsigned int i = 0; i < NSIZE; i++)
648 {
649 int itk = vc.Ltrk()[i];
650 charge += wTrackInfit(itk).charge();
651// double afield = VertexFitBField::instance()->getCBz(m_xInfit, pInfit(itk).sub(4, 6));
652 double afield = m_factor * VertexFitBField::instance()->getCBz(m_xInfit, pInfit(itk).sub(4, 6));
653 double a = afield * wTrackInfit(itk).charge();
654
655 // totalE += wTrackOrigin(itk).w()[3];
656 totalE += wTrackInfit(itk).w()[3]; // update total energy after fit, commit by sunhk, thanks to Chengzhi He.
657 Ai[0][4] = a;
658 Ai[1][3] = -a;
659 Bi[0][1] = -a;
660 Bi[1][0] = a;
661 A.sub(1, NTRKPAR*itk + 1, Ai);
662 B += Bi;
663 }
664 B[3][0] = 1.0;
665 B[4][1] = 1.0;
666 B[5][2] = 1.0;
667
668 w1 = A * m_pInfit + B * m_xInfit;
669 ew1 = m_pcovInfit.similarity(A) + m_xcovInfit.similarity(B) + A*m_E*B.T()+B*(m_E.T())*A.T();
670
671 //convert w1(6x1) to w2(7x1)
672 w2[0] = w1[0];
673 w2[1] = w1[1];
674 w2[2] = w1[2];
675 w2[3] = totalE;
676 w2[4] = w1[3];
677 w2[5] = w1[4];
678 w2[6] = w1[5];
679 //convert ew1(6x6) to ew2(7x7)
680 m_TRB[3][0] = w1[0] / totalE;
681 m_TRB[3][1] = w1[1] / totalE;
682 m_TRB[3][2] = w1[2] / totalE;
683 ew2 = m_TRB * ew1 * m_TRB.T();
684 Ew.assign(ew2);
685 WTrackParameter vwtrk;
686 vwtrk.setCharge(charge);
687 vwtrk.setW(w2);
688 vwtrk.setEw(Ew);
689
690 m_virtual_wtrk[n] = vwtrk;
691 timer.Stop();
692 m_cpu[5] += timer.CpuTime();
693}
VertexConstraints::Ltrk
std::vector< int > Ltrk() const
Definition VertexConstraints.h:55
VertexFitBField::getCBz
double getCBz(const HepVector &vtx, const HepVector &trackPosition)
Definition Analysis/VertexFit/VertexFit-00-03-05/src/BField.cxx:63
VertexFit::Ew
HepSymMatrix Ew(int n) const
Definition VertexFit.h:78
WTrackParameter::setEw
void setEw(const HepSymMatrix &Ew)
Definition WTrackParameter.h:53
WTrackParameter::setCharge
void setCharge(const int charge)
Definition WTrackParameter.h:54
WTrackParameter::setW
void setW(const HepVector &w)
Definition WTrackParameter.h:51
charge
float charge
Definition prepare_data.cxx:28

Referenced by DQAKsInfo::calculate(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), inclks::execute(), incllambda::execute(), JsiLL::execute(), KShortReconstruction::execute(), LambdaReconstruction::execute(), Signal::execute(), VertexFitRefine::Fit(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), and utility::SecondaryVFit_Lambda().

◆ chisq() [1/2]

double VertexFit::chisq ( ) const
inline

Definition at line 66 of file VertexFit.h.

66{return m_chi;}

Referenced by BeamFit(), DQAKsInfo::calculate(), VertexFitRefine::chisq(), VertexFitRefine::chisq(), BeamParams::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), KShortReconstruction::execute(), LambdaReconstruction::execute(), PrimaryVertex::execute(), Signal::execute(), and Fit().

◆ chisq() [2/2]

double VertexFit::chisq ( int n) const
inline

Definition at line 67 of file VertexFit.h.

67{return m_chisq[n];}

◆ cpu()

HepVector VertexFit::cpu ( ) const
inline

Definition at line 82 of file VertexFit.h.

82{return m_cpu;}

◆ errorVx()

double VertexFit::errorVx ( int n,
int i ) const
inline

Definition at line 88 of file VertexFit.h.

88{ return sqrt((m_vpar_infit[n].Evx())[i][i]);}
VertexFit::Evx
HepSymMatrix Evx(int n) const
Definition VertexFit.h:87

Referenced by VertexFitRefine::errorVx(), and VertexFitRefine::Fit().

◆ Evx()

HepSymMatrix VertexFit::Evx ( int n) const
inline

Definition at line 87 of file VertexFit.h.

87{ return m_vpar_infit[n].Evx();}

Referenced by errorVx(), VertexFitRefine::Evx(), Gam4pikp::execute(), PrimaryVertex::execute(), and VertexFitRefine::Fit().

◆ Ew()

HepSymMatrix VertexFit::Ew ( int n) const
inline

Definition at line 78 of file VertexFit.h.

78{return wTrackInfit(n).Ew();}

Referenced by BuildVirtualParticle(), VertexFitRefine::Ew(), and Fit().

◆ factor()

double VertexFit::factor ( ) const
inline

Definition at line 51 of file VertexFit.h.

51{ return m_factor; }

Referenced by setMagCorrFactor().

◆ Fit() [1/2]

bool VertexFit::Fit ( )

Definition at line 301 of file VertexFit.cxx.

302{
303 bool okfit = false;
304 double mychi = 0;
305 for (unsigned int n = 0; n<(int)(m_vc.size()); n++)
306 {
307 Fit(n);
308 if (m_chisq[n] >= m_chicut) return okfit;
309 swimVertex(n);
310 mychi = mychi + m_chisq[n];
311 }
312 m_chi = mychi;
313 okfit = true;
314 return okfit;
315}
VertexFit::Fit
bool Fit()
Definition VertexFit.cxx:301

Referenced by DQAKsInfo::calculate(), BeamParams::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), DQARhopi::execute(), Gam4pikp::execute(), inclks::execute(), incllambda::execute(), JsiLL::execute(), KShortReconstruction::execute(), LambdaReconstruction::execute(), Ppjrhopi::execute(), PrimaryVertex::execute(), Rhopi::execute(), rhopi::execute(), Signal::execute(), TestV::execute(), Fit(), VertexFitRefine::Fit(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), K3pi::MTotal(), K3pipi0::MTotal(), Kk::MTotal(), Kkpi0::MTotal(), Kkpipi::MTotal(), Kpi::MTotal(), Kpipi0::MTotal(), Kpipi0pi0::MTotal(), Pipi::MTotal(), Pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), utility::SecondaryVFit_Lambda(), utility::vfit(), and utility::vfit().

◆ Fit() [2/2]

bool VertexFit::Fit ( int n)

Definition at line 175 of file VertexFit.cxx.

176{
177 bool okfit = false;
178 TStopwatch timer;
179 m_cpu = HepVector(10, 0);
180 if (n < 0 || (unsigned int)n >= m_vc.size()) return okfit;
181
182 timer.Start();
183 int ifail;
184 m_nvtrk = numberWTrack();
185
186 m_pOrigin = HepVector(m_nvtrk * NTRKPAR, 0);
187 m_pInfit = HepVector(m_nvtrk * NTRKPAR, 0);
188 m_pcovOrigin = HepSymMatrix(m_nvtrk * NTRKPAR, 0);
189 m_pcovInfit = HepSymMatrix(m_nvtrk * NTRKPAR, 0);
190
191 int ntrk = numberWTrack();
192 for(unsigned int itk = 0; itk < ntrk; itk++)
193 {
194 setWTrackInfit(itk, wTrackOrigin(itk));
195 setPOrigin(itk, Convert76(wTrackOrigin(itk).w()));
196 setPCovOrigin(itk, wTrackOrigin(itk).Ew().similarity(m_TRA));
197 }
198 m_pInfit = m_pOrigin;
199
200 m_xOrigin = HepVector(NVTXPAR, 0);
201 m_xInfit = HepVector(NVTXPAR, 0);
202 m_xcovOrigin = HepSymMatrix(NVTXPAR, 0);
203 m_xcovInfit = HepSymMatrix(NVTXPAR, 0);
204 m_xcovOriginInversed = HepSymMatrix(NVTXPAR, 0);
205 m_xcovInfitInversed = HepSymMatrix(NVTXPAR, 0);
206
207 m_xOrigin = m_vpar_origin[n].Vx();
208 m_xcovOrigin = m_vpar_origin[n].Evx();
209 m_xcovOriginInversed = m_xcovOrigin.inverse(ifail);
210 m_xInfit = m_xOrigin;
211
212 m_vpar_infit[n] = m_vpar_origin[n];
213
214 m_B = HepMatrix(NCONSTR*m_nvtrk, NVTXPAR, 0);
215 m_A = HepMatrix(NCONSTR*m_nvtrk, NTRKPAR*m_nvtrk, 0);
216 m_BT = HepMatrix(NVTXPAR,NCONSTR*m_nvtrk, 0);
217 m_AT = HepMatrix(NTRKPAR*m_nvtrk, NCONSTR*m_nvtrk, 0);
218 m_G = HepVector(NCONSTR*m_nvtrk, 0);
219 m_W = HepSymMatrix(NCONSTR*m_nvtrk, 0);
220 m_E = HepMatrix(NTRKPAR*m_nvtrk, NVTXPAR, 0);
221
222 timer.Stop();
223 m_cpu[0] += timer.CpuTime();
224
225 // iteration loop
226 std::vector<double> chisq;
227 chisq.clear();
228 for (int it = 0; it < m_niter; it++)
229 {
230 timer.Start();
231 UpdateConstraints(m_vc[n]);
232 timer.Stop();
233 m_cpu[1] += timer.CpuTime();
234 timer.Start();
235 fitVertex(n);
236 timer.Stop();
237 m_cpu[2] += timer.CpuTime();
238 chisq.push_back(m_chisq[n]);
239 if (it > 0)
240 {
241 double delchi = chisq[it] - chisq[it-1];
242 if (fabs(delchi) < m_chiter) break;
243 }
244 }
245
246 /*REVISED
247 if(m_chisq[n] >= m_chicut || m_chisq[n] < 0) return okfit;
248 REVISED*/
249 if (m_chisq[n] >= m_chicut) return okfit;
250
251 // update vertex and its covariance
252 m_vpar_infit[n].setVx(m_xInfit);
253 m_vpar_infit[n].setEvx(m_xcovInfit);
254
255 okfit = true;
256 return okfit;
257}
w
double w
Definition EvtPhokharaDef.hh:29
TrackPool::numberWTrack
int numberWTrack() const
Definition TrackPool.h:79

◆ init()

void VertexFit::init ( )

Definition at line 29 of file VertexFit.cxx.

30{
31 //derived from TrackPool
32 clearWTrackOrigin();
33 clearWTrackInfit();
34 clearWTrackList();
35 clearGammaShape();
36 clearGammaShapeList();
37 clearMapkinematic();
38 clearMappositionA();
39 clearMappositionB();
40 clearone();
41 cleartwo();
42
43 m_vpar_origin.clear();
44 m_vpar_infit.clear();
45 m_vc.clear();
46 m_chisq.clear();
47 m_chi = 9999.;
48 m_virtual_wtrk.clear();
49 m_niter = 10;
50 m_chiter = 1.0e-3;
51 m_chicut = 1000;
52
53 m_TRA = HepMatrix(6, 7, 0);
54 m_TRA[0][0] = 1.0;
55 m_TRA[1][1] = 1.0;
56 m_TRA[2][2] = 1.0;
57 m_TRA[3][4] = 1.0;
58 m_TRA[4][5] = 1.0;
59 m_TRA[5][6] = 1.0;
60 m_TRB = HepMatrix(7, 6, 0);
61 m_TRB[0][0] = 1.0;
62 m_TRB[1][1] = 1.0;
63 m_TRB[2][2] = 1.0;
64 m_TRB[4][3] = 1.0;
65 m_TRB[5][4] = 1.0;
66 m_TRB[6][5] = 1.0;
67
68 m_factor = 1.000;
69}
TrackPool::clearMapkinematic
void clearMapkinematic()
Definition TrackPool.h:124
TrackPool::clearGammaShapeList
void clearGammaShapeList()
Definition TrackPool.h:140
TrackPool::clearWTrackOrigin
void clearWTrackOrigin()
Definition TrackPool.h:111
TrackPool::clearone
void clearone()
Definition TrackPool.h:115
TrackPool::clearMappositionB
void clearMappositionB()
Definition TrackPool.h:126
TrackPool::clearWTrackList
void clearWTrackList()
Definition TrackPool.h:113
TrackPool::clearMappositionA
void clearMappositionA()
Definition TrackPool.h:125
TrackPool::cleartwo
void cleartwo()
Definition TrackPool.h:116
TrackPool::clearGammaShape
void clearGammaShape()
Definition TrackPool.h:139
TrackPool::clearWTrackInfit
void clearWTrackInfit()
Definition TrackPool.h:112

Referenced by DQAKsInfo::calculate(), BeamParams::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), DQARhopi::execute(), Gam4pikp::execute(), inclks::execute(), incllambda::execute(), JsiLL::execute(), KShortReconstruction::execute(), LambdaReconstruction::execute(), Ppjrhopi::execute(), PrimaryVertex::execute(), Rhopi::execute(), rhopi::execute(), Signal::execute(), TestV::execute(), VertexFitRefine::Fit(), DQAKsInfo::init(), VertexFitRefine::init(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), K3pi::MTotal(), K3pipi0::MTotal(), Kk::MTotal(), Kkpi0::MTotal(), Kkpipi::MTotal(), Kpi::MTotal(), Kpipi0::MTotal(), Kpipi0pi0::MTotal(), Pipi::MTotal(), Pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), utility::SecondaryVFit_Lambda(), utility::vfit(), utility::vfit(), and VertexFitRefine::~VertexFitRefine().

◆ instance()

VertexFit * VertexFit::instance ( )
static

Definition at line 15 of file VertexFit.cxx.

16{
17 if (m_pointer) return m_pointer;
18 m_pointer = new VertexFit();
19 return m_pointer;
20}

Referenced by DQADtag::DQADtag(), DQAKsInfo::DQAKsInfo(), DQAKsInfo::DQAKsInfo(), BeamParams::execute(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), DQARhopi::execute(), Gam4pikp::execute(), inclks::execute(), incllambda::execute(), JsiLL::execute(), KShortReconstruction::execute(), LambdaReconstruction::execute(), Ppjrhopi::execute(), PrimaryVertex::execute(), Rhopi::execute(), rhopi::execute(), Signal::execute(), TestV::execute(), VertexFitRefine::init(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), K3pi::MTotal(), K3pipi0::MTotal(), Kk::MTotal(), Kkpi0::MTotal(), Kkpipi::MTotal(), Kpi::MTotal(), Kpipi0::MTotal(), Kpipi0pi0::MTotal(), Pipi::MTotal(), Pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), utility::SecondaryVFit_Lambda(), utility::vfit(), and utility::vfit().

◆ pfit()

HepLorentzVector VertexFit::pfit ( int n) const
inline

Definition at line 75 of file VertexFit.h.

75{return wTrackInfit(n).p();}

Referenced by VertexFitRefine::Fit(), and VertexFitRefine::pfit().

◆ pull()

bool VertexFit::pull ( int n,
int itk,
HepVector & p )

Definition at line 457 of file VertexFit.cxx.

458{
459 assert(p.num_row() == 5);
460 vertexCovMatrix(n);
461
462 WTrackParameter wtrk0, wtrk1;
463 HepVector w1(6, 0);
464 HepVector w2(7, 0);
465 HepSymMatrix ew1(6, 0);
466 HepSymMatrix ew2(7, 0);
467 wtrk0 = wTrackOrigin(itk);
468 w1 = pInfit(itk);
469 ew1 = pcovInfit(itk);
470 w2 = Convert67(wtrk0.mass(), w1);
471 m_TRB[3][0] = w2[0] / w2[3];
472 m_TRB[3][1] = w2[1] / w2[3];
473 m_TRB[3][2] = w2[2] / w2[3];
474 ew2 = ew1.similarity(m_TRB);
475 wtrk1.setW(w2);
476 wtrk1.setEw(ew2);
477 wtrk1.setCharge(wtrk0.charge());
478
479 HTrackParameter htrk0(wtrk0);
480 HTrackParameter htrk1(wtrk1);
481 for (int i = 0; i < 5; i++)
482 {
483 double del = htrk0.eHel()[i][i] - htrk1.eHel()[i][i];
484 if (del == 0.0)
485 {
486 return false;
487 }
488 p[i] = (htrk0.helix()[i] - htrk1.helix()[i]) / sqrt(abs(del));
489 }
490 return true;
491}
WTrackParameter::charge
int charge() const
Definition WTrackParameter.h:62
WTrackParameter::mass
double mass() const
Definition WTrackParameter.h:63

Referenced by BeamParams::execute().

◆ setChisqCut()

void VertexFit::setChisqCut ( const double chicut = 1000,
const double chiter = 1.0e-3 )
inline

Definition at line 49 of file VertexFit.h.

49{ m_chicut = chicut;m_chiter = chiter; }

Referenced by DQAKsInfo::calculate(), and VertexFitRefine::setChisqCut().

◆ setIterNumber()

void VertexFit::setIterNumber ( const int niter = 10)
inline

Definition at line 48 of file VertexFit.h.

48{ m_niter = niter; }

Referenced by VertexFitRefine::setIterNumber().

◆ setMagCorrFactor()

void VertexFit::setMagCorrFactor ( const double factor = 1.000)
inline

Definition at line 50 of file VertexFit.h.

50{ m_factor = factor; }
VertexFit::factor
double factor() const
Definition VertexFit.h:51

Referenced by VertexFitRefine::setMagCorrFactor().

◆ Swim()

void VertexFit::Swim ( int n)
inline

Definition at line 59 of file VertexFit.h.

59{ vertexCovMatrix(n); swimVertex(n);}

Referenced by DQAKsInfo::calculate(), DQARhopi::execute(), Gam4pikp::execute(), JsiLL::execute(), Ppjrhopi::execute(), Rhopi::execute(), TestV::execute(), VertexFitRefine::Fit(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), K3pi::MTotal(), K3pipi0::MTotal(), Kk::MTotal(), Kkpi0::MTotal(), Kkpipi::MTotal(), Kpi::MTotal(), Kpipi0::MTotal(), Kpipi0pi0::MTotal(), Pipi::MTotal(), Pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), utility::SecondaryVFit_Lambda(), utility::vfit(), and utility::vfit().

◆ vpar()

VertexParameter VertexFit::vpar ( int n) const
inline

Definition at line 89 of file VertexFit.h.

89{return m_vpar_infit[n];}

Referenced by AddBeamFit(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), AddVertex(), BeamFit(), DQAKsInfo::calculate(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), inclks::execute(), incllambda::execute(), JsiLL::execute(), Signal::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), utility::SecondaryVFit_Lambda(), and VertexFitRefine::vpar().

◆ Vx()

HepVector VertexFit::Vx ( int n) const
inline

Definition at line 86 of file VertexFit.h.

86{ return m_vpar_infit[n].Vx();}

Referenced by BeamParams::execute(), PrimaryVertex::execute(), VertexFitRefine::Fit(), and VertexFitRefine::Vx().

◆ vx()

HepPoint3D VertexFit::vx ( int n) const
inline

Definition at line 85 of file VertexFit.h.

85{ return m_vpar_infit[n].vx();}

Referenced by AddVertex(), Gam4pikp::execute(), VertexFitRefine::Fit(), and VertexFitRefine::vx().

◆ w()

HepVector VertexFit::w ( int n) const
inline

Definition at line 77 of file VertexFit.h.

77{return wTrackInfit(n).w();}

Referenced by VertexFitRefine::w().

◆ wtrk()

WTrackParameter VertexFit::wtrk ( int n) const
inline

Definition at line 79 of file VertexFit.h.

79{return wTrackInfit(n);}

Referenced by DQARhopi::execute(), Gam4pikp::execute(), JsiLL::execute(), Ppjrhopi::execute(), Rhopi::execute(), rhopi::execute(), TestV::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), K3pi::MTotal(), K3pipi0::MTotal(), Kk::MTotal(), Kkpi0::MTotal(), Kkpipi::MTotal(), Kpi::MTotal(), Kpipi0::MTotal(), Kpipi0pi0::MTotal(), Pipi::MTotal(), Pipipi0::MTotal(), utility::vfit(), utility::vfit(), and VertexFitRefine::wtrk().

◆ wVirtualTrack()

WTrackParameter VertexFit::wVirtualTrack ( int n) const
inline

Definition at line 92 of file VertexFit.h.

92{return m_virtual_wtrk[n];}

Referenced by DQAKsInfo::calculate(), DQAKsKpi::execute(), DQAKsKpiDEDX::execute(), inclks::execute(), incllambda::execute(), JsiLL::execute(), KShortReconstruction::execute(), LambdaReconstruction::execute(), Signal::execute(), K0kk::MTotal(), K0kpi::MTotal(), K0pi0::MTotal(), K0pipi::MTotal(), K0pipipi0::MTotal(), LocalKsSelector::operator()(), LocalLambdaSelector::operator()(), utility::SecondaryVFit(), utility::SecondaryVFit_Lambda(), and VertexFitRefine::wVirtualTrack().

◆ xfit()

HepPoint3D VertexFit::xfit ( int n) const
inline

Definition at line 76 of file VertexFit.h.

76{return wTrackInfit(n).x();}

Referenced by VertexFitRefine::Fit(), and VertexFitRefine::xfit().

Friends And Related Symbol Documentation

◆ VertexFitRefine

friend class VertexFitRefine
friend

Definition at line 20 of file VertexFit.h.

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