CGEM BOSS 6.6.5.h
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | List of all members
TrkCompTrk Class Reference

#include <TrkCompTrk.h>

+ Inheritance diagram for TrkCompTrk:

Public Member Functions

int charge () const
 
double chisq () const
 
int nDof () const
 
const TrkDifTraj & traj () const
 
const BField & bField () const
 
HepPoint3D position (double fltL) const
 
Hep3Vector direction (double fltL) const
 
Hep3Vector momentum (double fltL=0.) const
 
double pt (double fltL=0.) const
 
BesPointErr positionErr (double fltL) const
 
BesVectorErr directionErr (double fltL) const
 
BesVectorErr momentumErr (double fltL) const
 
double startValidRange () const
 
double endValidRange () const
 
virtual HepMatrix posmomCov (double fltL) const
 
virtual void getAllCovs (double fltL, HepSymMatrix &xxCov, HepSymMatrix &ppCov, HepMatrix &xpCov) const
 
virtual void getAllWeights (double fltL, HepVector &pos, HepVector &mom, HepSymMatrix &xxWeight, HepSymMatrix &ppWeight, HepMatrix &xpWeight) const
 
virtual void getAllWeights (const HepPoint3D &pt, HepVector &pos, HepVector &mom, HepSymMatrix &xxWeight, HepSymMatrix &ppWeight, HepMatrix &xpWeight) const
 
virtual void print (std::ostream &) const
 
virtual void printAll (std::ostream &) const
 
 TrkCompTrk (const BesPointErr &pos, const BesVectorErr &mom, const HepMatrix &xpCov, int charge, double chisq, int nDoF, const BField *bf)
 
 TrkCompTrk (const TrkCompTrk &right)
 
virtual ~TrkCompTrk ()
 
const TrkCompTrk & operator= (const TrkCompTrk &right)
 
- Public Member Functions inherited from TrkAbsFit

Additional Inherited Members

- Protected Member Functions inherited from TrkAbsFit
 TrkAbsFit ()
 
virtual ~TrkAbsFit ()
 

Detailed Description

Definition at line 24 of file TrkCompTrk.h.

Constructor & Destructor Documentation

◆ TrkCompTrk() [1/2]

TrkCompTrk::TrkCompTrk ( const BesPointErr & pos,
const BesVectorErr & mom,
const HepMatrix & xpCov,
int charge,
double chisq,
int nDoF,
const BField * bf )

Definition at line 28 of file TrkCompTrk.cxx.

31 :
32 _chisq(chisq),_nDof(nDof) {
33//---------------------------------------------------------------------------
34 _bf = bf;
35 _charge=charge;
36 if(_charge!=0) {
37 TrkExchangePar par1 = TrkHelixUtils::helixFromMomErr(pos,mom,xpCov,charge,bField());
38 _traj.reset( new HelixTraj(par1.params(), par1.covariance()) );
39 }
40 else {
41 _traj.reset( new NeutTraj(TrkHelixUtils::lineFromMomErr(pos,mom,xpCov,1,bField())) );
42 }
43}
TrkCompTrk::chisq
double chisq() const
Definition TrkCompTrk.cxx:80
TrkCompTrk::bField
const BField & bField() const
Definition TrkCompTrk.h:36
TrkCompTrk::charge
int charge() const
Definition TrkCompTrk.cxx:89
TrkCompTrk::nDof
int nDof() const
Definition TrkCompTrk.cxx:73
TrkExchangePar::params
const HepVector & params() const
Definition TrkExchangePar.h:48
TrkExchangePar::covariance
const HepSymMatrix & covariance() const
Definition TrkExchangePar.h:50
TrkHelixUtils::lineFromMomErr
static NeutParams lineFromMomErr(const BesPointErr &vertex, const BesVectorErr &p, const HepMatrix &cxp, double sign, const BField &)
Definition TrkHelixUtils.cxx:238
TrkHelixUtils::helixFromMomErr
static TrkExchangePar helixFromMomErr(const BesPointErr &vertex, const BesVectorErr &p, const HepMatrix &cxp, double sign, const BField &)
Definition TrkHelixUtils.cxx:146

◆ TrkCompTrk() [2/2]

TrkCompTrk::TrkCompTrk ( const TrkCompTrk & right)

Definition at line 47 of file TrkCompTrk.cxx.

47 : _bf(rhs._bf),_chisq(rhs._chisq),_nDof(rhs._nDof) {
48//---------------------------------------------
49 _traj.reset((TrkSimpTraj*)rhs.traj().clone());
50 _charge=rhs.charge();
51}

◆ ~TrkCompTrk()

TrkCompTrk::~TrkCompTrk ( )
virtual

Definition at line 54 of file TrkCompTrk.cxx.

54 {
55//----------------------------------------------------------------------
56}

Member Function Documentation

◆ bField()

const BField & TrkCompTrk::bField ( ) const
inline

Definition at line 36 of file TrkCompTrk.h.

36{return *_bf;}

Referenced by getAllCovs(), getAllWeights(), momentum(), momentumErr(), posmomCov(), pt(), and TrkCompTrk().

◆ charge()

int TrkCompTrk::charge ( ) const
virtual

Implements TrkAbsFit.

Definition at line 89 of file TrkCompTrk.cxx.

89 {
90//------------------------------------------------------------------------
91 return _charge;
92}

Referenced by operator=(), TrkCompTrk(), and TrkCompTrk().

◆ chisq()

double TrkCompTrk::chisq ( ) const
virtual

Implements TrkAbsFit.

Definition at line 80 of file TrkCompTrk.cxx.

80 {
81//------------------------------------------------------------------------
82 return _chisq;
83}

Referenced by operator=().

◆ direction()

Hep3Vector TrkCompTrk::direction ( double fltL) const
virtual

Implements TrkAbsFit.

Definition at line 111 of file TrkCompTrk.cxx.

111 {
112//------------------------------------------------------------------------
113 return traj().direction(fltL);
114}
Trajectory::direction
virtual Hep3Vector direction(double) const =0
TrkCompTrk::traj
const TrkDifTraj & traj() const
Definition TrkCompTrk.cxx:97

◆ directionErr()

BesVectorErr TrkCompTrk::directionErr ( double fltL) const
virtual

Implements TrkAbsFit.

Definition at line 138 of file TrkCompTrk.cxx.

139{
140 // If getDFInfo is split up, use faster version - done JEA 020300
141 DifPoint posD;
142 DifVector dirD;
143 traj().getDFInfo2(fltL, posD, dirD);
144 HepMatrix err = dirD.errorMatrix( dirD.x.indepPar()->covariance() );
145 BesError symErr(3);
146 symErr.assign(err);
147 Hep3Vector dir(dirD.x.number(), dirD.y.number(), dirD.z.number());
148 return BesVectorErr(dir, symErr);
149}
DifIndepPar::covariance
HepSymMatrix & covariance()
Definition DifIndepPar.h:53
DifNumber::number
double number() const
Definition DifNumber.h:87
DifNumber::indepPar
const DifIndepPar * indepPar() const
Definition DifNumber.h:93
DifVector::y
DifNumber y
Definition DifVector.h:149
DifVector::x
DifNumber x
Definition DifVector.h:148
DifVector::z
DifNumber z
Definition DifVector.h:150
DifVector::errorMatrix
HepSymMatrix errorMatrix(const HepSymMatrix &e) const
Definition DifVector.cxx:54
TrkDifTraj::getDFInfo2
virtual void getDFInfo2(double fltLen, DifPoint &pos, DifVector &direction) const
Definition TrkDifTraj.cxx:25

◆ endValidRange()

double TrkCompTrk::endValidRange ( ) const
virtual

Implements TrkAbsFit.

Definition at line 174 of file TrkCompTrk.cxx.

175{
176 return traj().hiRange();
177}
Trajectory::hiRange
double hiRange() const
Definition Trajectory.h:92

◆ getAllCovs()

void TrkCompTrk::getAllCovs ( double fltL,
HepSymMatrix & xxCov,
HepSymMatrix & ppCov,
HepMatrix & xpCov ) const
virtual

Implements TrkAbsFit.

Definition at line 215 of file TrkCompTrk.cxx.

219{
220 const BField& theField = bField();
221 TrkMomCalculator::getAllCovs(*_traj, theField, fltL,
222 xxCov,ppCov,xpCov);
223}
TrkMomCalculator::getAllCovs
static void getAllCovs(const TrkSimpTraj &, const BField &, double fltlen, HepSymMatrix &xxCov, HepSymMatrix &ppCov, HepMatrix &xpCov)
Definition TrkMomCalculator.cxx:207

◆ getAllWeights() [1/2]

void TrkCompTrk::getAllWeights ( const HepPoint3D & pt,
HepVector & pos,
HepVector & mom,
HepSymMatrix & xxWeight,
HepSymMatrix & ppWeight,
HepMatrix & xpWeight ) const
virtual

Definition at line 239 of file TrkCompTrk.cxx.

245{
246 double fltL=0;
247 TrkPoca poca(traj(),fltL, pt);
248 fltL = poca.flt1();
249 getAllWeights(fltL,pos,mom,xxWeight,ppWeight,xpWeight);
250}
TrkCompTrk::pt
double pt(double fltL=0.) const
Definition TrkCompTrk.cxx:158
TrkCompTrk::getAllWeights
virtual void getAllWeights(double fltL, HepVector &pos, HepVector &mom, HepSymMatrix &xxWeight, HepSymMatrix &ppWeight, HepMatrix &xpWeight) const
Definition TrkCompTrk.cxx:225

◆ getAllWeights() [2/2]

void TrkCompTrk::getAllWeights ( double fltL,
HepVector & pos,
HepVector & mom,
HepSymMatrix & xxWeight,
HepSymMatrix & ppWeight,
HepMatrix & xpWeight ) const
virtual

Implements TrkAbsFit.

Definition at line 225 of file TrkCompTrk.cxx.

231{
232 const BField& theField = bField();
233 TrkMomCalculator::getAllWeights(*_traj, theField, fltL,
234 pos,mom,xxWeight,ppWeight,xpWeight);
235
236}
TrkMomCalculator::getAllWeights
static void getAllWeights(const TrkSimpTraj &, const BField &, double fltlen, HepVector &pos, HepVector &mom, HepSymMatrix &xxWeight, HepSymMatrix &ppWeight, HepMatrix &xpWeight)
Definition TrkMomCalculator.cxx:252

Referenced by getAllWeights().

◆ momentum()

Hep3Vector TrkCompTrk::momentum ( double fltL = 0.) const
virtual

Implements TrkAbsFit.

Definition at line 117 of file TrkCompTrk.cxx.

118{
119 return TrkMomCalculator::vecMom(*_traj, bField(), fltL);
120
121}
TrkMomCalculator::vecMom
static Hep3Vector vecMom(const TrkSimpTraj &, const BField &, double fltlen)
Definition TrkMomCalculator.cxx:82

◆ momentumErr()

BesVectorErr TrkCompTrk::momentumErr ( double fltL) const
virtual

Implements TrkAbsFit.

Definition at line 152 of file TrkCompTrk.cxx.

153{
154 return TrkMomCalculator::errMom(*_traj, bField(), fltL);
155}
TrkMomCalculator::errMom
static BesVectorErr errMom(const TrkSimpTraj &, const BField &, double fltlen)
Definition TrkMomCalculator.cxx:112

◆ nDof()

int TrkCompTrk::nDof ( ) const
virtual

Implements TrkAbsFit.

Definition at line 73 of file TrkCompTrk.cxx.

73 {
74//------------------------------------------------------------------------
75 return _nDof;
76}

Referenced by operator=().

◆ operator=()

const TrkCompTrk & TrkCompTrk::operator= ( const TrkCompTrk & right)

Definition at line 60 of file TrkCompTrk.cxx.

60 {
61//----------------------------------------------------------------------
62 if (&right == this) return *this;
63 _bf = right._bf;
64 _traj.reset((TrkSimpTraj*)right.traj().clone());
65 _chisq=right.chisq();
66 _nDof=right.nDof();
67 _charge=right.charge();
68 return *this;
69}
Trajectory::clone
virtual Trajectory * clone() const =0

◆ position()

HepPoint3D TrkCompTrk::position ( double fltL) const
virtual

Implements TrkAbsFit.

Definition at line 104 of file TrkCompTrk.cxx.

104 {
105//------------------------------------------------------------------------
106 return traj().position(fltL);
107}
Trajectory::position
virtual HepPoint3D position(double) const =0

◆ positionErr()

BesPointErr TrkCompTrk::positionErr ( double fltL) const
virtual

Implements TrkAbsFit.

Definition at line 124 of file TrkCompTrk.cxx.

125{
126 // If getDFInfo is split up, use faster version - done JEA 020300
127 DifPoint posD;
128 DifVector dirD;
129 traj().getDFInfo2(fltL, posD, dirD);
130 HepMatrix err = posD.errorMatrix( posD.x.indepPar()->covariance() );
131 BesError symErr(3);
132 symErr.assign(err);
133 HepPoint3D point(posD.x.number(), posD.y.number(), posD.z.number());
134 return BesPointErr(point, symErr);
135}

◆ posmomCov()

HepMatrix TrkCompTrk::posmomCov ( double fltL) const
virtual

Implements TrkAbsFit.

Definition at line 209 of file TrkCompTrk.cxx.

210{
211 const BField& theField = bField();
212 return TrkMomCalculator::posmomCov(*_traj, theField, fltL);
213}
TrkMomCalculator::posmomCov
static HepMatrix posmomCov(const TrkSimpTraj &, const BField &, double fltlen)
Definition TrkMomCalculator.cxx:180

◆ print()

void TrkCompTrk::print ( std::ostream & ) const
virtual

Implements TrkAbsFit.

Definition at line 180 of file TrkCompTrk.cxx.

181{
182 ostr << "Traj: ";
183 if(_charge==-1 || _charge==1) {
184 ostr << "Charged Particle -> Helix Trajectory" << endl;
185 HelixTraj& theTraj=(HelixTraj&)traj();
186 ostr << " phi0=" << BesAngle(theTraj.phi0()).rad() << endl;
187 ostr << " d0=" << theTraj.d0() << endl;
188 ostr << " z0=" << theTraj.z0() << endl;
189 ostr << " omega=" << theTraj.omega() << endl;
190 ostr << " tanDip=" << theTraj.tanDip() << endl;
191 }
192 else {
193 ostr << "Neutral Particle -> NeutTraj" << endl;
194 NeutTraj& theTraj=(NeutTraj&)traj();
195 ostr << " phi0=" << BesAngle(theTraj.params().phi0()).rad() << endl;
196 ostr << " d0=" << theTraj.params().d0() << endl;
197 ostr << " z0=" << theTraj.params().z0() << endl;
198 ostr << " p=" << theTraj.params().p() << endl;
199 ostr << " tanDip=" << theTraj.params().tanDip() << endl;
200 }
201}
BesAngle::rad
double rad() const
Definition BesAngle.h:118
HelixTraj::omega
double omega() const
Definition HelixTraj.h:81
HelixTraj::d0
double d0() const
Definition HelixTraj.h:79
HelixTraj::phi0
double phi0() const
Definition HelixTraj.cxx:510
HelixTraj::z0
double z0() const
Definition HelixTraj.h:82
HelixTraj::tanDip
double tanDip() const
Definition HelixTraj.h:83
NeutParams::phi0
double & phi0()
Definition NeutParams.h:38
NeutParams::z0
double & z0()
Definition NeutParams.h:40
NeutParams::tanDip
double & tanDip()
Definition NeutParams.h:41
NeutParams::p
double & p()
Definition NeutParams.h:39
NeutParams::d0
double & d0()
Definition NeutParams.h:37
NeutTraj::params
NeutParams & params()
Definition NeutTraj.h:54

Referenced by printAll().

◆ printAll()

void TrkCompTrk::printAll ( std::ostream & ) const
virtual

Implements TrkAbsFit.

Definition at line 204 of file TrkCompTrk.cxx.

205{
206 print(ostr);
207}
TrkCompTrk::print
virtual void print(std::ostream &) const
Definition TrkCompTrk.cxx:180

◆ pt()

double TrkCompTrk::pt ( double fltL = 0.) const
virtual

Implements TrkAbsFit.

Definition at line 158 of file TrkCompTrk.cxx.

159{
160 return TrkMomCalculator::ptMom(*_traj, bField(), fltL);
161 // This could probably be made more efficient.
162// Hep3Vector pVec = momentum(fltL);
163// return sqrt( pVec.x() * pVec.x() + pVec.y() * pVec.y() );
164}
TrkMomCalculator::ptMom
static double ptMom(const TrkSimpTraj &, const BField &, double fltlen)
Definition TrkMomCalculator.cxx:47

Referenced by getAllWeights().

◆ startValidRange()

double TrkCompTrk::startValidRange ( ) const
virtual

Implements TrkAbsFit.

Definition at line 168 of file TrkCompTrk.cxx.

169{
170 return traj().lowRange();
171}
Trajectory::lowRange
double lowRange() const
Definition Trajectory.h:91

◆ traj()

const TrkDifTraj & TrkCompTrk::traj ( ) const
virtual

Implements TrkAbsFit.

Definition at line 97 of file TrkCompTrk.cxx.

97 {
98//------------------------------------------------------------------------
99 return *_traj;
100}

Referenced by direction(), directionErr(), endValidRange(), getAllWeights(), operator=(), position(), positionErr(), print(), startValidRange(), and TrkCompTrk().

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