TrigPrimaryVertexTrack Class Reference

#include <TrigPrimaryVertexTrack.h>

Inheritance diagram for TrigPrimaryVertexTrack:

Collaboration diagram for TrigPrimaryVertexTrack:

Public Member Functions
	TrigPrimaryVertexTrack (const Trk::Track *)
	constructor for EF (offline) tracks
	~TrigPrimaryVertexTrack ()
const Trk::Track *	getTrkTrack ()
	getter for EF (offline) tracks
virtual double	getChi2Distance (TrigL2Vertex *)
	implementation of abstract method from the base class
virtual void	updateVertex (TrigL2Vertex *)
	implementation of abstract method from the base class
virtual MsgStream &	report (MsgStream &) const
void	setIndex (int)
	to be used by TrigVertexingTool
int	getIndex () const
	to be used by TrigVertexingTool
int	getTrackType ()
	0: L2 track, 1: EF(offline) track
bool	isActive ()
	if true this track will be used in the vertex fit otherwise it will be masked
void	activate ()
	sets m_isActive to true
void	mask ()
	sets m_isActive to false
const double *	Perigee () const
	track parameters at the perigee
double	PerigeeCovariance (int, int) const
	covariance of track parameters at the perigee
double	getChi2Contribution ()
	chi2-contribution to the vertex fit

Protected Attributes
double	m_resid [2] {}
double	m_V [2][2] {}
double	m_D [2][MAX_SIZE_VERT_COVM] {}

Private Attributes
const Trk::Track *	m_pTrkTrack
int	m_nTrackType
int	m_index = 0
double	m_Vqq [3][3] {}
double	m_Vuq [2][3] {}
double	m_Vuu [2][2] {}
double	m_u [2] {}
double	m_q [3] {}
double	m_Perigee [5] {}
double	m_PerigeeCovariance [5][5] {}
double	m_A [2][3] {}
double	m_B [2][3] {}
double	m_dChi2
bool	m_active

Detailed Description

Definition at line 13 of file TrigPrimaryVertexTrack.h.

Constructor & Destructor Documentation

TrigPrimaryVertexTrack()

TrigPrimaryVertexTrack::TrigPrimaryVertexTrack

(

const Trk::Track *

pT

)

constructor for EF (offline) tracks

Definition at line 11 of file TrigPrimaryVertexTrack.cxx.

{
  m_nTrackType=2;m_active=true;
  m_pTrkTrack=pT;m_dChi2=-100.0;
 
  const Trk::Perigee* pP= pT->perigeeParameters();
 
  if(pP!=NULL) {
    m_u[0]=pP->parameters()[Trk::d0];
    m_u[1]=pP->parameters()[Trk::z0];
    m_q[0]=pP->parameters()[Trk::phi0];
    m_q[1]=pP->parameters()[Trk::theta];
 
    const double ptC=1000.0*pP->parameters()[Trk::qOverP];
    const double inv_ptC = 1. / ptC;
 
    m_q[2]=sin(pP->parameters()[Trk::theta])*inv_ptC;
 
    const AmgSymMatrix(5)& TC= (*pP->covariance());
 
    const double a=cos(pP->parameters()[Trk::theta])*inv_ptC;
    const double b=-sin(pP->parameters()[Trk::theta])/(pP->parameters()[Trk::qOverP]*ptC);
 
         
    const double Ck[5][5] = {{TC(0,0), TC(1,0), TC(2,0), TC(3,0), a*TC(3,0) + b*TC(4,0)},
                             {TC(1,0), TC(1,1), TC(2,1), TC(3,1), a*TC(3,1) + b*TC(4,1)},
                             {TC(2,0), TC(2,1), TC(2,2), TC(3,2), a*TC(3,2) + b*TC(4,2)},
                             {TC(3,0), TC(3,1), TC(3,2), TC(3,3), a*TC(3,3) + b*TC(4,3)},
                             {a*TC(3,0) + b*TC(4,0), a*TC(3,1) + b*TC(4,1), a*TC(3,2) + b*TC(4,2), a*TC(3,3) + b*TC(4,3), a*a*TC(3,3) + 2*a*b*TC(4,3) + b*b*TC(4,4)}};
 
    for(int i=0;i<2;i++) for(int j=0;j<2;j++) m_Vuu[i][j]=Ck[i][j];
    for(int i=0;i<3;i++) for(int j=0;j<3;j++) m_Vqq[i][j]=Ck[i+2][j+2];
    for(int i=0;i<2;i++) for(int j=0;j<3;j++) m_Vuq[i][j]=Ck[i][j+2];
 
    m_Perigee[0]=m_u[0];m_Perigee[1]=m_u[1];m_Perigee[2]=m_q[0];m_Perigee[3]=m_q[1];m_Perigee[4]=m_q[2];
  }
}

~TrigPrimaryVertexTrack()

TrigPrimaryVertexTrack::~TrigPrimaryVertexTrack

(

)

Definition at line 49 of file TrigPrimaryVertexTrack.cxx.

{
 
}

Member Function Documentation

activate()

void TrigPrimaryVertexTrack::activate

(

)

sets m_isActive to true

Definition at line 59 of file TrigPrimaryVertexTrack.cxx.

{
  m_active=true;
}

getChi2Contribution()

double TrigPrimaryVertexTrack::getChi2Contribution

(

)

chi2-contribution to the vertex fit

Definition at line 194 of file TrigPrimaryVertexTrack.cxx.

{
  return m_dChi2;
}

getChi2Distance()

double TrigPrimaryVertexTrack::getChi2Distance ( TrigL2Vertex * pV )

virtual

implementation of abstract method from the base class

Implements TrigVertexFittingNode.

Definition at line 100 of file TrigPrimaryVertexTrack.cxx.

{
  const double C=0.02997;
  const double B=20.84;
    
  double Sk[2][2],detr,chi2;
  double AC[2][3],BV[2][3],h[2];
  int i,j,k;
  double psi,sinPsi,xv,yv,zv,
    cosPhi0,sinPhi0,ctt,sint,phi0,theta0,P0;
  const double alpha = C*B*1e-3;
  const double inv_alpha = 1. / alpha;
 
  xv=pV->getParametersVector()[0];
  yv=pV->getParametersVector()[1];
  zv=pV->getParametersVector()[2];
 
  phi0=m_q[0];
  theta0=m_q[1];
  P0=m_q[2];
 
  cosPhi0=cos(phi0);sinPhi0=sin(phi0);
  sinPsi=-alpha*(xv*cosPhi0+yv*sinPhi0)/P0;
  if(fabs(sinPsi)>1.0) return -999.9;
  const double cosPsi=sqrt(1.0-sinPsi*sinPsi);
  const double inv_cosPsi = 1. / cosPsi;
  psi=asin(sinPsi);
  sint=sin(theta0);
  ctt=cos(theta0)/sint;
 
  m_A[0][0]=-sin(phi0+psi)*inv_cosPsi;
  m_A[0][1]= cos(phi0+psi)*inv_cosPsi;
  m_A[0][2]=0.0;
 
  m_A[1][0]=-ctt*cosPhi0*inv_cosPsi;
  m_A[1][1]=-ctt*sinPhi0*inv_cosPsi;
  m_A[1][2]=1.0;
 
  m_B[0][0]=-xv*m_A[0][1]+yv*m_A[0][0];
  m_B[0][1]=0.0;
  m_B[0][2]=(1.0-inv_cosPsi)*inv_alpha;
 
  m_B[1][0]=-xv*m_A[1][1]+yv*m_A[1][0];
  m_B[1][1]=-P0*psi/(alpha*sint*sint);
  m_B[1][2]=ctt*(psi-sinPsi*inv_cosPsi)*inv_alpha;
 
  h[0]=yv*cosPhi0-xv*sinPhi0+P0*(1-cosPsi)*inv_alpha;
  h[1]=zv+P0*ctt*psi*inv_alpha;
 
  m_resid[0]=m_u[0]-h[0];
  m_resid[1]=m_u[1]-h[1];
 
  for(i=0;i<2;i++) for(j=0;j<2;j++) Sk[i][j]=m_Vuu[i][j];
  for(i=0;i<2;i++) for(j=0;j<3;j++)
    {
      AC[i][j]=0.0;
      for(k=0;k<3;k++) AC[i][j]+=m_A[i][k]*pV->m_Gk[k][j];
    }
  for(i=0;i<2;i++) for(j=0;j<2;j++)
    {
      for(k=0;k<3;k++) Sk[i][j]+=AC[i][k]*m_A[j][k];
    }
  for(i=0;i<2;i++)
    for(j=0;j<3;j++)
      {
    BV[i][j]=0.0;
    for(k=0;k<3;k++) BV[i][j]+=m_B[i][k]*m_Vqq[k][j];
      }
  for(i=0;i<2;i++)
    for(j=0;j<2;j++)
      {
    for(k=0;k<3;k++) Sk[i][j]+=BV[i][k]*m_B[j][k];
      }
  Sk[0][0]-=2.0*(m_Vuq[0][0]*m_B[0][0]+m_Vuq[0][1]*m_B[0][1]+m_Vuq[0][2]*m_B[0][2]);
  Sk[1][1]-=2.0*(m_Vuq[1][0]*m_B[1][0]+m_Vuq[1][1]*m_B[1][1]+m_Vuq[1][2]*m_B[1][2]);
  Sk[0][1]-=m_Vuq[1][0]*m_B[0][0]+m_Vuq[1][1]*m_B[0][1]+m_Vuq[1][2]*m_B[0][2]+
    m_Vuq[0][0]*m_B[1][0]+m_Vuq[0][1]*m_B[1][1]+m_Vuq[0][2]*m_B[1][2];
  Sk[1][0]=Sk[0][1];
 
  detr=1.0/(Sk[0][0]*Sk[1][1]-Sk[0][1]*Sk[1][0]);
  m_V[0][0]=Sk[1][1]*detr;
  m_V[1][1]=Sk[0][0]*detr;
  m_V[0][1]=m_V[1][0]=-Sk[0][1]*detr;
  
  chi2=m_V[0][0]*m_resid[0]*m_resid[0]+m_V[1][1]*m_resid[1]*m_resid[1]+
    2.0*m_V[0][1]*m_resid[1]*m_resid[0];
 
  for(i=0;i<2;i++) 
    for(j=0;j<3;j++) m_D[i][j]=AC[i][j];
    
  m_dChi2=chi2;
  return chi2;
}

getIndex()

int TrigPrimaryVertexTrack::getIndex

(

)

const

to be used by TrigVertexingTool

Definition at line 85 of file TrigPrimaryVertexTrack.cxx.

{
  return m_index;
}

getTrackType()

int TrigPrimaryVertexTrack::getTrackType

(

)

0: L2 track, 1: EF(offline) track

Definition at line 69 of file TrigPrimaryVertexTrack.cxx.

{
  return m_nTrackType;
}

getTrkTrack()

const Trk::Track * TrigPrimaryVertexTrack::getTrkTrack

(

)

getter for EF (offline) tracks

Definition at line 75 of file TrigPrimaryVertexTrack.cxx.

{
    return m_pTrkTrack;
}

isActive()

bool TrigPrimaryVertexTrack::isActive

(

)

if true this track will be used in the vertex fit otherwise it will be masked

Definition at line 64 of file TrigPrimaryVertexTrack.cxx.

{
  return m_active;
}

mask()

void TrigPrimaryVertexTrack::mask

(

)

sets m_isActive to false

Definition at line 54 of file TrigPrimaryVertexTrack.cxx.

{
  m_active=false;
}

Perigee()

const double * TrigPrimaryVertexTrack::Perigee

(

)

const

track parameters at the perigee

Definition at line 90 of file TrigPrimaryVertexTrack.cxx.

{
  return &m_Perigee[0];
}

PerigeeCovariance()

double TrigPrimaryVertexTrack::PerigeeCovariance	(	int	i,
		int	j ) const

covariance of track parameters at the perigee

Definition at line 95 of file TrigPrimaryVertexTrack.cxx.

{
  return m_PerigeeCovariance[i][j];
}

report()

MsgStream & TrigPrimaryVertexTrack::report ( MsgStream & out ) const

virtual

Implements TrigVertexFittingNode.

Definition at line 224 of file TrigPrimaryVertexTrack.cxx.

{
  int i;
 
  out<<"Primary track "<<m_index<<endmsg;
  for(i=0;i<2;i++)
    {
      out<<"  u"<<i<<" = "<<m_u[i]<<"      "<<m_Vuu[i][0]<<"   "<<m_Vuu[i][1]<<endmsg;
    }
  for(i=0;i<3;i++)
    {
      out<<"  q"<<i<<" = "<<m_q[i]<<"      "<<m_Vqq[i][0]<<"   "<<m_Vqq[i][1]<<"   "<<m_Vqq[i][2]<<endmsg;
    }
  return out;
}

setIndex()

void TrigPrimaryVertexTrack::setIndex

(

int

i

)

to be used by TrigVertexingTool

Definition at line 80 of file TrigPrimaryVertexTrack.cxx.

{
  m_index=i;
}

updateVertex()

void TrigPrimaryVertexTrack::updateVertex ( TrigL2Vertex * pV )

virtual

implementation of abstract method from the base class

Implements TrigVertexFittingNode.

Definition at line 199 of file TrigPrimaryVertexTrack.cxx.

{
  int i,j,k;
 
  double K[2][3];
 
  for(i=0;i<2;i++)
    {
      for(j=0;j<3;j++)
    {
      K[i][j]=0.0;
      for(k=0;k<2;k++) K[i][j]+=m_D[k][j]*m_V[k][i];
    }
    }
  for(i=0;i<3;i++)
    {
      pV->getParametersVector()[i]+=K[0][i]*m_resid[0]+K[1][i]*m_resid[1];
      for(j=i;j<3;j++)
    {
      pV->m_Gk[i][j]-=K[0][i]*m_D[0][j]+K[1][i]*m_D[1][j];
      pV->m_Gk[j][i]=pV->m_Gk[i][j];
    }
    }
}

Member Data Documentation

m_A

double TrigPrimaryVertexTrack::m_A[2][3] {}

private

Definition at line 43 of file TrigPrimaryVertexTrack.h.

{};

m_active

bool TrigPrimaryVertexTrack::m_active

private

Definition at line 46 of file TrigPrimaryVertexTrack.h.

m_B

double TrigPrimaryVertexTrack::m_B[2][3] {}

private

Definition at line 44 of file TrigPrimaryVertexTrack.h.

{};

m_D

double TrigVertexFittingNode::m_D[2][MAX_SIZE_VERT_COVM] {}

protectedinherited

Definition at line 48 of file TrigL2Vertex.h.

{};

m_dChi2

double TrigPrimaryVertexTrack::m_dChi2

private

Definition at line 45 of file TrigPrimaryVertexTrack.h.

m_index

int TrigPrimaryVertexTrack::m_index = 0

private

Definition at line 35 of file TrigPrimaryVertexTrack.h.

m_nTrackType

int TrigPrimaryVertexTrack::m_nTrackType

private

Definition at line 34 of file TrigPrimaryVertexTrack.h.

m_Perigee

double TrigPrimaryVertexTrack::m_Perigee[5] {}

private

Definition at line 41 of file TrigPrimaryVertexTrack.h.

{};

m_PerigeeCovariance

double TrigPrimaryVertexTrack::m_PerigeeCovariance[5][5] {}

private

Definition at line 42 of file TrigPrimaryVertexTrack.h.

{};

m_pTrkTrack

const Trk::Track* TrigPrimaryVertexTrack::m_pTrkTrack

private

Definition at line 33 of file TrigPrimaryVertexTrack.h.

m_q

double TrigPrimaryVertexTrack::m_q[3] {}

private

Definition at line 40 of file TrigPrimaryVertexTrack.h.

{};

m_resid

double TrigVertexFittingNode::m_resid[2] {}

protectedinherited

Definition at line 46 of file TrigL2Vertex.h.

{};

m_u

double TrigPrimaryVertexTrack::m_u[2] {}

private

Definition at line 39 of file TrigPrimaryVertexTrack.h.

{};

m_V

double TrigVertexFittingNode::m_V[2][2] {}

protectedinherited

Definition at line 47 of file TrigL2Vertex.h.

{};

m_Vqq

double TrigPrimaryVertexTrack::m_Vqq[3][3] {}

private

Definition at line 36 of file TrigPrimaryVertexTrack.h.

{};

m_Vuq

double TrigPrimaryVertexTrack::m_Vuq[2][3] {}

private

Definition at line 37 of file TrigPrimaryVertexTrack.h.

{};

m_Vuu

double TrigPrimaryVertexTrack::m_Vuu[2][2] {}

private

Definition at line 38 of file TrigPrimaryVertexTrack.h.

{};

---

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

• TrigPrimaryVertexTrack.h
• TrigPrimaryVertexTrack.cxx