BOSS 6.6.4.p01
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | List of all members
VertexConstraints Class Reference

#include <VertexConstraints.h>

Public Member Functions

 VertexConstraints ()
 
 ~VertexConstraints ()
 
 VertexConstraints (const VertexConstraints &)
 
void CommonVertexConstraints (std::vector< int > tlis)
 
void FixedVertexConstraints (std::vector< int > tlis)
 
void UpdateConstraints (VertexParameter vpar, std::vector< WTrackParameter > wlis)
 
void setCovax (const int i, const HepMatrix cov)
 
void setLambda (const int i, const HepVector L)
 
void setLambda0 (const int i, const HepVector L0)
 
void setVD (const int i, const HepSymMatrix vd)
 
void setEVDE (const int i, const HepSymMatrix evde)
 
void setType (const int type)
 
std::vector< HepMatrix > Ec () const
 
std::vector< HepMatrix > Dc () const
 
std::vector< HepVector > dc () const
 
std::vector< HepSymMatrix > VD () const
 
std::vector< HepSymMatrix > EVDE () const
 
std::vector< HepMatrix > covax () const
 
std::vector< HepVector > lambda0 () const
 
std::vector< HepVector > lambda () const
 
std::vector< int > Ltrk () const
 
int type () const
 
int commonVertex () const
 
int fixedVertex () const
 
VertexConstraints & operator= (const VertexConstraints &)
 

Detailed Description

Definition at line 30 of file VertexConstraints.h.

Constructor & Destructor Documentation

◆ VertexConstraints() [1/2]

VertexConstraints::VertexConstraints ( )

Definition at line 9 of file VertexConstraints.cxx.

10{
11 m_Ec.clear();
12 m_Dc.clear();
13 m_dc.clear();
14 m_VD.clear();
15 m_EVDE.clear();
16 m_lambda0.clear();
17 m_lambda.clear();
18 m_covax.clear();
19 m_ltrk.clear();
20 m_type = 0;
21}

◆ ~VertexConstraints()

VertexConstraints::~VertexConstraints ( )
inline

Definition at line 35 of file VertexConstraints.h.

35{;}

◆ VertexConstraints() [2/2]

VertexConstraints::VertexConstraints ( const VertexConstraints &  vc)
inline

Definition at line 77 of file VertexConstraints.h.

78{
79 m_Ec = vc.m_Ec;
80 m_Dc = vc.m_Dc;
81 m_dc = vc.m_dc;
82 m_VD = vc.m_VD;
83 m_EVDE = vc.m_EVDE;
84 m_covax = vc.m_covax;
85 m_lambda0 = vc.m_lambda0;
86 m_lambda = vc.m_lambda;
87 m_ltrk = vc.m_ltrk;
88 m_type = vc.m_type;
89}

Member Function Documentation

◆ commonVertex()

int VertexConstraints::commonVertex ( ) const
inline

Definition at line 57 of file VertexConstraints.h.

57{return CommonVertex;}

Referenced by CommonVertexConstraints().

◆ CommonVertexConstraints()

void VertexConstraints::CommonVertexConstraints ( std::vector< int >  tlis)

Definition at line 23 of file VertexConstraints.cxx.

24{
25 m_ltrk = tlis;
26 setType(commonVertex());
27}
VertexConstraints::commonVertex
int commonVertex() const
Definition: VertexConstraints.h:57
VertexConstraints::setType
void setType(const int type)
Definition: VertexConstraints.h:45

Referenced by VertexFit::AddVertex().

◆ covax()

std::vector< HepMatrix > VertexConstraints::covax ( ) const
inline

Definition at line 52 of file VertexConstraints.h.

52{return m_covax;}

Referenced by UpdateConstraints().

◆ Dc()

std::vector< HepMatrix > VertexConstraints::Dc ( ) const
inline

Definition at line 48 of file VertexConstraints.h.

48{return m_Dc;}

Referenced by UpdateConstraints().

◆ dc()

std::vector< HepVector > VertexConstraints::dc ( ) const
inline

Definition at line 49 of file VertexConstraints.h.

49{return m_dc;}

Referenced by UpdateConstraints().

◆ Ec()

std::vector< HepMatrix > VertexConstraints::Ec ( ) const
inline

Definition at line 47 of file VertexConstraints.h.

47{return m_Ec;}

Referenced by UpdateConstraints().

◆ EVDE()

std::vector< HepSymMatrix > VertexConstraints::EVDE ( ) const
inline

Definition at line 51 of file VertexConstraints.h.

51{return m_EVDE;}

Referenced by UpdateConstraints().

◆ fixedVertex()

int VertexConstraints::fixedVertex ( ) const
inline

Definition at line 58 of file VertexConstraints.h.

58{return FixedVertex;}

Referenced by FixedVertexConstraints().

◆ FixedVertexConstraints()

void VertexConstraints::FixedVertexConstraints ( std::vector< int >  tlis)

Definition at line 29 of file VertexConstraints.cxx.

30{
31 m_ltrk = tlis;
32 setType(fixedVertex());
33}
VertexConstraints::fixedVertex
int fixedVertex() const
Definition: VertexConstraints.h:58

Referenced by VertexFit::AddBeamFit().

◆ lambda()

std::vector< HepVector > VertexConstraints::lambda ( ) const
inline

Definition at line 54 of file VertexConstraints.h.

54{return m_lambda;}

Referenced by UpdateConstraints().

◆ lambda0()

std::vector< HepVector > VertexConstraints::lambda0 ( ) const
inline

Definition at line 53 of file VertexConstraints.h.

53{return m_lambda0;}

Referenced by UpdateConstraints().

◆ Ltrk()

std::vector< int > VertexConstraints::Ltrk ( ) const
inline

Definition at line 55 of file VertexConstraints.h.

55{return m_ltrk;}

Referenced by VertexFit::BuildVirtualParticle().

◆ operator=()

VertexConstraints & VertexConstraints::operator= ( const VertexConstraints &  vc)
inline

Definition at line 91 of file VertexConstraints.h.

92{
93 m_Ec = vc.m_Ec;
94 m_Dc = vc.m_Dc;
95 m_dc = vc.m_dc;
96 m_VD = vc.m_VD;
97 m_EVDE = vc.m_EVDE;
98 m_covax = vc.m_covax;
99 m_lambda0 = vc.m_lambda0;
100 m_lambda = vc.m_lambda;
101 m_ltrk = vc.m_ltrk;
102 m_type = vc.m_type;
103 return (*this);
104}

◆ setCovax()

void VertexConstraints::setCovax ( const int  i,
const HepMatrix  cov 
)
inline

Definition at line 40 of file VertexConstraints.h.

40{m_covax[i] = cov;}

◆ setEVDE()

void VertexConstraints::setEVDE ( const int  i,
const HepSymMatrix  evde 
)
inline

Definition at line 44 of file VertexConstraints.h.

44{m_EVDE[i] = evde;}

◆ setLambda()

void VertexConstraints::setLambda ( const int  i,
const HepVector  L 
)
inline

Definition at line 41 of file VertexConstraints.h.

41{m_lambda[i] = L;}

◆ setLambda0()

void VertexConstraints::setLambda0 ( const int  i,
const HepVector  L0 
)
inline

Definition at line 42 of file VertexConstraints.h.

42{m_lambda0[i] = L0;}

◆ setType()

void VertexConstraints::setType ( const int  type)
inline

Definition at line 45 of file VertexConstraints.h.

45{m_type = type;}

Referenced by CommonVertexConstraints(), and FixedVertexConstraints().

◆ setVD()

void VertexConstraints::setVD ( const int  i,
const HepSymMatrix  vd 
)
inline

Definition at line 43 of file VertexConstraints.h.

43{m_VD[i] = vd;}

◆ type()

int VertexConstraints::type ( ) const
inline

Definition at line 56 of file VertexConstraints.h.

56{return m_type;}

Referenced by setType().

◆ UpdateConstraints()

void VertexConstraints::UpdateConstraints ( VertexParameter  vpar,
std::vector< WTrackParameter >  wlis 
)

Definition at line 35 of file VertexConstraints.cxx.

36{
37 m_Ec.clear();
38 m_Dc.clear();
39 m_dc.clear();
40 m_VD.clear();
41 m_EVDE.clear();
42 m_lambda0.clear();
43 m_lambda.clear();
44 m_covax.clear();
45
46 switch(m_type)
47 {
48 case FixedVertex:
49 {
50 WTrackParameter wtrk;
51 for (unsigned int i = 0; i < wlis.size(); i++)
52 {
53 wtrk = wlis[i];
54 HepLorentzVector p = wtrk.p();
55 HepPoint3D x = wtrk.x();
56 HepPoint3D delx = vpar.vx() - x;
57
58 double afield = VertexFitBField::instance()->getCBz(vpar.Vx(), wtrk.X());
59 double a = afield * (0.0+wtrk.charge());
60 double J = a * (delx.x()*p.px() + delx.y()*p.py())/p.perp2();
61
62 J = std::min(J, 1-1e-4);
63 J = std::max(J, -1+1e-4);
64 double Rx = delx.x() - 2*p.px() * (delx.x()*p.px() + delx.y()*p.py()) / p.perp2();
65 double Ry = delx.y() - 2*p.py() * (delx.x()*p.px() + delx.y()*p.py()) / p.perp2();
66 double S = 1.0/sqrt(1-J*J)/p.perp2();
67
68 // dc
69 HepVector dc(2, 0);
70 dc[0] = delx.y() * p.px() - delx.x() * p.py() - 0.5 * a * (delx.x() * delx.x() + delx.y() * delx.y());
71 dc[1] = delx.z() - p.pz()/a*asin(J);
72 m_dc.push_back(dc);
73 // Ec
74 HepMatrix Ec(2, 3, 0);
75 m_Ec.push_back(Ec);
76 // Dc
77 HepMatrix Dc(2, 6, 0);
78 Dc[0][0] = delx.y();
79 Dc[0][1] = -delx.x();
80 Dc[0][2] = 0;
81 Dc[0][3] = p.py() + a * delx.x();
82 Dc[0][4] = -p.px() + a * delx.y();
83 Dc[0][5] = 0;
84 Dc[1][0] = -p.pz() * S * Rx;
85 Dc[1][1] = -p.pz() * S * Ry;
86 Dc[1][2] = - asin(J) / a;
87 Dc[1][3] = p.px() * p.pz() * S;
88 Dc[1][4] = p.py() * p.pz() * S;
89 Dc[1][5] = -1.0;
90 m_Dc.push_back(Dc);
91 // VD
92 HepSymMatrix VD(2, 0);
93 // int ifail;
94 // VD = ((wtrk.Ew()).similarity(Dc)).inverse(ifail); // 2x2 matrix (Dc Ew Dc.T())
95 m_VD.push_back(VD);
96 // EVDE
97 HepSymMatrix EVDE(3, 0);
98 m_EVDE.push_back(EVDE);
99 // lambda0
100 HepVector lambda0(2, 0);
101 m_lambda0.push_back(lambda0);
102 // lambda
103 HepVector lambda(2, 0);
104 m_lambda.push_back(lambda);
105 // covax
106 HepMatrix covax(6, 3, 0);
107 m_covax.push_back(covax);
108 break;
109 }
110 }
111 case CommonVertex:
112 default:
113 {
114 WTrackParameter wtrk;
115 for (unsigned int i = 0; i < wlis.size(); i++)
116 {
117 wtrk = wlis[i];
118 HepLorentzVector p = wtrk.p();
119 HepPoint3D x = wtrk.x();
120 HepPoint3D delx = vpar.vx() - x;
121 if (wtrk.charge() == 0)
122 {
123 //dc
124 HepVector dc(2,0);
125 dc[0] = p.pz()*delx.x() - p.px()*delx.z();
126 dc[1] = p.pz()*delx.y() - p.py()*delx.z();
127 m_dc.push_back(dc);
128 //Ec
129 HepMatrix Ec(2,3,0);
130 Ec[0][0] = p.pz();
131 Ec[0][1] = 0;
132 Ec[0][2] = -p.px();
133 Ec[1][0] = 0;
134 Ec[1][1] = p.pz();
135 Ec[1][2] = -p.py();
136 m_Ec.push_back(Ec);
137 //Dc
138 HepMatrix Dc(2,6,0);
139 Dc[0][0] = -delx.z();
140 Dc[0][1] = 0;
141 Dc[0][2] = delx.x();
142 Dc[0][3] = -p.pz();
143 Dc[0][4] = 0;
144 Dc[0][5] = p.px();
145 Dc[1][0] = 0;
146 Dc[1][1] = -delx.z();
147 Dc[1][2] = delx.y();
148 Dc[1][3] = 0;
149 Dc[1][4] = -p.pz();
150 Dc[1][5] = p.py();
151 m_Dc.push_back(Dc);
152 }
153 else
154 {
155 double afield = VertexFitBField::instance()->getCBz(vpar.Vx(), wtrk.X());
156 double a = afield * (0.0+wtrk.charge());
157 double J = a * (delx.x()*p.px() + delx.y()*p.py())/p.perp2();
158 J = std::min(J, 1-1e-4);
159 J = std::max(J, -1+1e-4);
160 double Rx = delx.x() - 2*p.px() * (delx.x()*p.px() + delx.y()*p.py()) / p.perp2();
161 double Ry = delx.y() - 2*p.py() * (delx.x()*p.px() + delx.y()*p.py()) / p.perp2();
162 double S = 1.0/sqrt(1-J*J)/p.perp2();
163 // dc
164 HepVector dc(2, 0);
165 dc[0] = delx.y() * p.px() - delx.x() * p.py() - 0.5 * a * (delx.x() * delx.x() + delx.y() * delx.y());
166 dc[1] = delx.z() - p.pz()/a*asin(J);
167 m_dc.push_back(dc);
168 // Ec
169 HepMatrix Ec(2, 3, 0);
170 Ec[0][0] = -p.py()- a * delx.x();
171 Ec[0][1] = p.px() - a * delx.y();
172 Ec[0][2] = 0;
173 Ec[1][0] = -p.px() * p.pz() * S ;
174 Ec[1][1] = -p.py() * p.pz() * S ;
175 Ec[1][2] = 1.0;
176 m_Ec.push_back(Ec);
177 // Dc
178 HepMatrix Dc(2, 6, 0);
179 Dc[0][0] = delx.y();
180 Dc[0][1] = -delx.x();
181 Dc[0][2] = 0;
182 Dc[0][3] = p.py() + a * delx.x();
183 Dc[0][4] = -p.px() + a * delx.y();
184 Dc[0][5] = 0;
185 Dc[1][0] = -p.pz() * S * Rx;
186 Dc[1][1] = -p.pz() * S * Ry;
187 Dc[1][2] = - asin(J) / a;
188 Dc[1][3] = p.px() * p.pz() * S;
189 Dc[1][4] = p.py() * p.pz() * S;
190 Dc[1][5] = -1.0;
191 m_Dc.push_back(Dc);
192 }
193 // VD
194 HepSymMatrix VD(2, 0);
195 m_VD.push_back(VD);
196 // EVDE
197 HepSymMatrix EVDE(3, 0);
198 m_EVDE.push_back(EVDE);
199 // lambda0
200 HepVector lambda0(2, 0);
201 m_lambda0.push_back(lambda0);
202 // lambda
203 HepVector lambda(2, 0);
204 m_lambda.push_back(lambda);
205 // covax
206 HepMatrix covax(6, 3, 0);
207 m_covax.push_back(covax);
208 }
209 break;
210 }
211 }
212}
x
Double_t x[10]
Definition: DataBase/tau_mode.c:57
VertexConstraints::Dc
std::vector< HepMatrix > Dc() const
Definition: VertexConstraints.h:48
VertexConstraints::EVDE
std::vector< HepSymMatrix > EVDE() const
Definition: VertexConstraints.h:51
VertexConstraints::VD
std::vector< HepSymMatrix > VD() const
Definition: VertexConstraints.h:50
VertexConstraints::covax
std::vector< HepMatrix > covax() const
Definition: VertexConstraints.h:52
VertexConstraints::Ec
std::vector< HepMatrix > Ec() const
Definition: VertexConstraints.h:47
VertexConstraints::lambda0
std::vector< HepVector > lambda0() const
Definition: VertexConstraints.h:53
VertexConstraints::dc
std::vector< HepVector > dc() const
Definition: VertexConstraints.h:49
VertexConstraints::lambda
std::vector< HepVector > lambda() const
Definition: VertexConstraints.h:54
VertexFitBField::getCBz
double getCBz(const HepVector &vtx, const HepVector &trackPosition)
Definition: Analysis/VertexFit/VertexFit-00-02-69/src/BField.cxx:63
VertexParameter::vx
HepPoint3D vx() const
Definition: VertexParameter.h:30
VertexParameter::Vx
HepVector Vx() const
Definition: VertexParameter.h:31
WTrackParameter::charge
int charge() const
Definition: WTrackParameter.h:62
WTrackParameter::p
HepLorentzVector p() const
Definition: WTrackParameter.h:69
WTrackParameter::x
HepPoint3D x() const
Definition: WTrackParameter.h:71
WTrackParameter::X
HepVector X() const
Definition: WTrackParameter.h:72

◆ VD()

std::vector< HepSymMatrix > VertexConstraints::VD ( ) const
inline

Definition at line 50 of file VertexConstraints.h.

50{return m_VD;}

Referenced by UpdateConstraints().

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