StraightLine Class Reference

#include <StraightLine.h>

Public Member Functions
	StraightLine (const HepPoint3D &pivot, const HepVector &a)
	StraightLine (double dr, double phi0, double dz, double tanl)
	StraightLine (HepPoint3D pos1, HepPoint3D pos2)
virtual	~StraightLine ()
	Destructor.
const HepPoint3D &	pivot (void) const
	returns pivot position.
const HepVector &	a (void) const
	returns helix parameters.
double	dr (void) const
	returns an element of parameters.
double	phi0 (void) const
double	dz (void) const
double	tanl (void) const
double	sinPhi0 (void) const
double	cosPhi0 (void) const
double	lambda () const
double	theta () const
double	cosTheta () const
HepPoint3D	x (double s=0.) const
	returns position after moving s in downwoards
double	sAtR (double R, int direction=1) const
HepPoint3D	xAtR (double R, int direction=1) const
double	sAtPhi (double phi) const
HepPoint3D	xAtPhi (double phi) const
void	set (const HepVector &a)
	returns position and convariance matrix(Ex) after moving s.
void	set (const HepPoint3D &pivot)
void	set (const HepPoint3D &pivot, const HepVector &a)
	sets pivot position, parameters.
void	print ()

Detailed Description

Definition at line 23 of file StraightLine.h.

Constructor & Destructor Documentation

StraightLine() [1/3]

StraightLine::StraightLine	(	const HepPoint3D &	pivot,
		const HepVector &	a
	)

Definition at line 8 of file StraightLine.cxx.

 :m_pivot(pivot),m_a(a)
{
}

StraightLine() [2/3]

StraightLine::StraightLine	(	double	dr,
		double	phi0,
		double	dz,
		double	tanl
	)

Definition at line 13 of file StraightLine.cxx.

 :m_a(HepVector(4,0))
{
     m_a[0]=dr;
     m_a[1]=phi0;
     m_a[2]=dz;
     m_a[3]=tanl;
 
     HepPoint3D pivot(0,0,0);
     m_pivot=pivot;
}

StraightLine() [3/3]

StraightLine::StraightLine	(	HepPoint3D	pos1,
		HepPoint3D	pos2
	)

Definition at line 25 of file StraightLine.cxx.

 :m_a(HepVector(4,0))
{
     // cout << "debug int StraightLine" << endl;
     // cout << "2points: " << setw(15) << pos1.x() << setw(15) << pos1.y() << setw(15) << pos1.z()
     //       << setw(15) << pos2.x() << setw(15) << pos2.y() << setw(15) << pos2.z() << endl;
     HepPoint3D pivot(0,0,0);
     m_pivot=pivot;
 
     if(pos1.y()<pos2.y()) {
      HepPoint3D pos_tmp=pos1;
      pos1=pos2;
      pos2=pos_tmp;
     }
     else if(pos1.y()==pos2.y()&&pos1.x()>pos2.x()) {
      HepPoint3D pos_tmp = pos1;
      pos1=pos2;
      pos2=pos_tmp;
     }
     //cout<<__FILE__<<": pos1, pos2 = "<<pos1<<", "<<pos2<<endl;
     HepPoint3D vec = pos2-pos1;
     double phi_vec = vec.phi();// (pi,2pi]
     //cout<<__FILE__<<": phi_vec="<<phi_vec<<endl;
     double phi0=phi_vec+M_PI/2.;// automatically in (-pi/2,pi/2)
     //cout<<__FILE__<<": phi0   ="<<phi0   <<endl;
     double dphi1=pos1.phi()-phi0;
     double dr = pos1.perp()*cos(dphi1);
 
     double s1(0.),s2(0.);
     if(pos1.y()!=pos2.y()) {
      s1=(pos1.y()-dr*sin(phi0))/sin(phi_vec);
      s2=(pos2.y()-dr*sin(phi0))/sin(phi_vec);
     }
     else {
      s1=(pos1.x()-dr*cos(phi0))/cos(phi_vec);
      s2=(pos2.x()-dr*cos(phi0))/cos(phi_vec);
     }
    
     double dz(0.),tanl(0.);
     if(s1>=s2) {
      cout<<"StraightLine: s1>=s2!!!"<<endl;
     }
     else {
      tanl=(pos2.z()-pos1.z())/(s2-s1);
      dz=-s1*tanl+pos1.z();
     }
     m_a[0]=dr;
     m_a[1]=phi0;
     m_a[2]=dz;
     m_a[3]=tanl;
     // cout << "line: " << setw(15) << dr << setw(15) << phi0 << setw(15) << dz << setw(15) << tanl << endl;
}

~StraightLine()

StraightLine::~StraightLine ( )

virtual

Destructor.

Definition at line 78 of file StraightLine.cxx.

{}

Member Function Documentation

a()

const HepVector & StraightLine::a ( void  ) const

inline

returns helix parameters.

Definition at line 43 of file StraightLine.h.

{return m_a;};

Referenced by set().

cosPhi0()

double StraightLine::cosPhi0 ( void  ) const

inline

Definition at line 51 of file StraightLine.h.

{return cos(phi0());};

cosTheta()

double StraightLine::cosTheta ( ) const

inline

Definition at line 54 of file StraightLine.h.

{return cos(theta());};

dr()

double StraightLine::dr ( void  ) const

inline

returns an element of parameters.

Definition at line 45 of file StraightLine.h.

{return m_a[0];};

Referenced by CgemLineFit::execute(), CgemLineFit::Fit(), CgemLineFit::Get4Par(), CgemLineFit::Get_OtherIniPar(), CgemMidDriftPlane::getPointAligned(), BesCgemSD::ProcessHits(), and StraightLine().

dz()

double StraightLine::dz ( void  ) const

inline

Definition at line 48 of file StraightLine.h.

{return m_a[2];};

Referenced by CgemLineFit::execute(), CgemLineFit::Fit(), CgemLineFit::Get4Par(), CgemLineFit::Get_OtherIniPar(), CgemMidDriftPlane::getPointAligned(), and StraightLine().

lambda()

double StraightLine::lambda ( ) const

inline

Definition at line 52 of file StraightLine.h.

{return atan(m_a[3]);}; 

Referenced by theta().

phi0()

double StraightLine::phi0 ( void  ) const

inline

Definition at line 46 of file StraightLine.h.

{return m_a[1];};

Referenced by cosPhi0(), CgemLineFit::execute(), CgemLineFit::Fit(), CgemLineFit::Get4Par(), CgemLineFit::Get_OtherIniPar(), CgemMidDriftPlane::getPointAligned(), sAtPhi(), sinPhi0(), and StraightLine().

pivot()

const HepPoint3D & StraightLine::pivot ( void  ) const

inline

returns pivot position.

Definition at line 40 of file StraightLine.h.

{return m_pivot;};

Referenced by set(), and StraightLine().

print()

void StraightLine::print

(

)

Definition at line 152 of file StraightLine.cxx.

{
     cout<<"line pivot: "<<m_pivot<<endl;
     cout<<"     par: "<<m_a<<endl;
}

sAtPhi()

double StraightLine::sAtPhi

(

double

phi

)

const

Definition at line 108 of file StraightLine.cxx.

{
     // if pivot at (0,0,0)
     double ss=999999.;
     //double       phi0_global  = m_a[1];//(-pi/2,pi/2]
     if(m_a[0]<0) {
      //phi0_global+ = M_PI;  //(pi/2,3/2*pi]
      phi+=M_PI;
     }
     double dphi=phi0()-phi;
     while(dphi> M_PI) dphi-=2*M_PI;
     while(dphi<-M_PI) dphi+=2*M_PI;
     if(fabs(dphi)<M_PI/2)
     {
      ss=m_a[0]*tan(dphi);
     }
     return ss;
}

Referenced by xAtPhi().

sAtR()

double StraightLine::sAtR	(	double	R,
		int	direction = 1
	)		const

Definition at line 89 of file StraightLine.cxx.

{
     // if pivot at (0,0,0)
     double ss=999999.;
     if(fabs(m_a[0])<=R) {
      ss=sqrt(R*R-m_a[0]*m_a[0]);
      ss*=direction;
     }
     return ss;
}

Referenced by CgemGeoAlign::StraightLineConversion(), and xAtR().

set() [1/3]

void StraightLine::set

(

const HepPoint3D &

pivot

)

Definition at line 140 of file StraightLine.cxx.

{
     m_pivot = pivot;
     // m_a to be updated // FIXME
}

set() [2/3]

void StraightLine::set	(	const HepPoint3D &	pivot,
		const HepVector &	a
	)

sets pivot position, parameters.

Definition at line 146 of file StraightLine.cxx.

{
     m_pivot = pivot;
     m_a=a;
}

set() [3/3]

void StraightLine::set

(

const HepVector &

a

)

returns position and convariance matrix(Ex) after moving s.

Definition at line 135 of file StraightLine.cxx.

{
     m_a=a;
}

sinPhi0()

double StraightLine::sinPhi0 ( void  ) const

inline

Definition at line 50 of file StraightLine.h.

{return sin(phi0());};

tanl()

double StraightLine::tanl ( void  ) const

inline

Definition at line 49 of file StraightLine.h.

{return m_a[3];};

Referenced by CgemLineFit::execute(), CgemLineFit::Fit(), CgemLineFit::Get4Par(), CgemLineFit::Get_OtherIniPar(), CgemMidDriftPlane::getPointAligned(), and StraightLine().

theta()

double StraightLine::theta ( ) const

inline

Definition at line 53 of file StraightLine.h.

{return M_PI/2.-lambda();};

Referenced by cosTheta().

x()

HepPoint3D StraightLine::x

(

double

s = 0.

)

const

returns position after moving s in downwoards

Definition at line 81 of file StraightLine.cxx.

{
     double x=m_a[0]*cos(m_a[1])+s*cos(m_a[1]-M_PI/2.);
     double y=m_a[0]*sin(m_a[1])+s*sin(m_a[1]-M_PI/2.);
     double z=m_a[2]+s*m_a[3];
     return HepPoint3D(x,y,z);
}

Referenced by CgemLineFit::IncidentAngle(), CgemGeoAlign::StraightLineConversion(), x(), xAtPhi(), and xAtR().

xAtPhi()

HepPoint3D StraightLine::xAtPhi

(

double

phi

)

const

Definition at line 127 of file StraightLine.cxx.

{
     HepPoint3D pos(999999,999999,999999);
     double ss=sAtPhi(phi);
     if(ss!=999999.) pos=x(ss);
     return pos;
}

xAtR()

HepPoint3D StraightLine::xAtR	(	double	R,
		int	direction = 1
	)		const

Definition at line 100 of file StraightLine.cxx.

{
     HepPoint3D pos(999999,999999,999999);
     double ss=sAtR(R, direction);
     if(ss!=999999.) pos=x(ss);
     return pos;
}

Referenced by CgemMidDriftPlane::getPointIdealGeom(), BesCgemSD::ProcessHits(), and CgemGeoAlign::StraightLineConversion().

---

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

• StraightLine.h
• StraightLine.cxx