#include <TrigPrimaryVertexTrack.h>

Inheritance diagram for TrigPrimaryVertexTrack:

Collaboration diagram for TrigPrimaryVertexTrack:

Public Member Functions
	TrigPrimaryVertexTrack (const TrigInDetTrack *)
	constructor for L2 tracks
	TrigPrimaryVertexTrack (const Trk::Track *)
	constructor for EF (offline) tracks
	~TrigPrimaryVertexTrack ()
const TrigInDetTrack *	getTrigTrack ()
	getter for L2 tracks
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 TrigInDetTrack *	m_pTrigTrack
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 14 of file TrigPrimaryVertexTrack.h.

Constructor & Destructor Documentation

◆ TrigPrimaryVertexTrack() [1/2]

TrigPrimaryVertexTrack::TrigPrimaryVertexTrack ( const TrigInDetTrack * pT )

constructor for L2 tracks

Definition at line 11 of file TrigPrimaryVertexTrack.cxx.

{
  m_nTrackType=1;m_active=true;
  m_pTrigTrack=pT;m_pTrkTrack=NULL;m_dChi2=-100.0;
  double Ck[5][5];
  int i,j;
 
  const TrigInDetTrackFitPar* p=pT->param();
 
  double a=-2.0*exp(-p->eta())/(1.0+exp(-2.0*p->eta()));
  m_u[0]=p->a0();
  m_u[1]=p->z0();
  m_q[0]=p->phi0();
  m_q[1]=2.0*atan(exp(-p->eta()));
  m_q[2]=p->pT()*1e-3;
 
 
  Ck[0][0]=(*p->cov())[0];Ck[0][1]=Ck[1][0]=(*p->cov())[2];
  Ck[0][2]=Ck[2][0]=(*p->cov())[1];Ck[0][3]=Ck[3][0]=(*p->cov())[3];
  Ck[0][4]=Ck[4][0]=(*p->cov())[4];Ck[1][1]=(*p->cov())[9];
  Ck[1][2]=Ck[2][1]=(*p->cov())[6];Ck[1][3]=Ck[3][1]=(*p->cov())[10];
  Ck[1][4]=Ck[4][1]=(*p->cov())[11];Ck[2][2]=(*p->cov())[5];
  Ck[2][3]=Ck[3][2]=(*p->cov())[7];Ck[2][4]=Ck[4][2]=(*p->cov())[8];
  Ck[3][3]=(*p->cov())[12];Ck[3][4]=Ck[4][3]=(*p->cov())[13];
  Ck[4][4]=(*p->cov())[14];
  for(i=0;i<5;i++)
    {
      Ck[3][i]=a*Ck[3][i];Ck[i][3]=Ck[3][i];
    }
  Ck[3][3]*=a;
  for(i=0;i<5;i++)                    
    {
      Ck[4][i]=Ck[4][i]*1e-3;   Ck[i][4]=Ck[4][i];
    }
  Ck[4][4]*=1e-3;
 
 
  for(i=0;i<2;i++) for(j=0;j<2;j++) m_Vuu[i][j]=Ck[i][j];
  for(i=0;i<2;i++) for(j=0;j<3;j++) m_Vuq[i][j]=Ck[i][j+2];
  for(i=0;i<3;i++) for(j=0;j<3;j++) m_Vqq[i][j]=Ck[i+2][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() [2/2]

TrigPrimaryVertexTrack::TrigPrimaryVertexTrack ( const Trk::Track * pT )

constructor for EF (offline) tracks

Definition at line 55 of file TrigPrimaryVertexTrack.cxx.

{
  m_nTrackType=2;m_active=true;
  m_pTrkTrack=pT;m_pTrigTrack=NULL;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 93 of file TrigPrimaryVertexTrack.cxx.

{
 
}

Member Function Documentation

◆ activate()

void TrigPrimaryVertexTrack::activate ( )

sets m_isActive to true

Definition at line 103 of file TrigPrimaryVertexTrack.cxx.

{
  m_active=true;
}

◆ getChi2Contribution()

double TrigPrimaryVertexTrack::getChi2Contribution ( )

chi2-contribution to the vertex fit

Definition at line 243 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 149 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 134 of file TrigPrimaryVertexTrack.cxx.

{
  return m_index;
}

◆ getTrackType()

int TrigPrimaryVertexTrack::getTrackType ( )

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

Definition at line 113 of file TrigPrimaryVertexTrack.cxx.

{
  return m_nTrackType;
}

◆ getTrigTrack()

const TrigInDetTrack * TrigPrimaryVertexTrack::getTrigTrack ( )

getter for L2 tracks

Definition at line 119 of file TrigPrimaryVertexTrack.cxx.

{
  return m_pTrigTrack;
}

◆ getTrkTrack()

const Trk::Track * TrigPrimaryVertexTrack::getTrkTrack ( )

getter for EF (offline) tracks

Definition at line 124 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 108 of file TrigPrimaryVertexTrack.cxx.

{
  return m_active;
}

◆ mask()

void TrigPrimaryVertexTrack::mask ( )

sets m_isActive to false

Definition at line 98 of file TrigPrimaryVertexTrack.cxx.

{
  m_active=false;
}

◆ Perigee()

const double * TrigPrimaryVertexTrack::Perigee ( ) const

track parameters at the perigee

Definition at line 139 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 144 of file TrigPrimaryVertexTrack.cxx.

{
  return m_PerigeeCovariance[i][j];
}

◆ report()

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

virtual

Implements TrigVertexFittingNode.

Definition at line 273 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 129 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 248 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 47 of file TrigPrimaryVertexTrack.h.

47{};

◆ m_active

bool TrigPrimaryVertexTrack::m_active

private

Definition at line 50 of file TrigPrimaryVertexTrack.h.

◆ m_B

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

private

Definition at line 48 of file TrigPrimaryVertexTrack.h.

48{};

◆ m_D

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

protectedinherited

Definition at line 48 of file TrigL2Vertex.h.

48{};

◆ m_dChi2

double TrigPrimaryVertexTrack::m_dChi2

private

Definition at line 49 of file TrigPrimaryVertexTrack.h.

◆ m_index

int TrigPrimaryVertexTrack::m_index = 0

private

Definition at line 39 of file TrigPrimaryVertexTrack.h.

◆ m_nTrackType

int TrigPrimaryVertexTrack::m_nTrackType

private

Definition at line 38 of file TrigPrimaryVertexTrack.h.

◆ m_Perigee

double TrigPrimaryVertexTrack::m_Perigee[5] {}

private

Definition at line 45 of file TrigPrimaryVertexTrack.h.

45{};

◆ m_PerigeeCovariance

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

private

Definition at line 46 of file TrigPrimaryVertexTrack.h.

46{};

◆ m_pTrigTrack

const TrigInDetTrack* TrigPrimaryVertexTrack::m_pTrigTrack

private

Definition at line 36 of file TrigPrimaryVertexTrack.h.

◆ m_pTrkTrack

const Trk::Track* TrigPrimaryVertexTrack::m_pTrkTrack

private

Definition at line 37 of file TrigPrimaryVertexTrack.h.

◆ m_q

double TrigPrimaryVertexTrack::m_q[3] {}

private

Definition at line 44 of file TrigPrimaryVertexTrack.h.

44{};

◆ m_resid

double TrigVertexFittingNode::m_resid[2] {}

protectedinherited

Definition at line 46 of file TrigL2Vertex.h.

46{};

◆ m_u

double TrigPrimaryVertexTrack::m_u[2] {}

private

Definition at line 43 of file TrigPrimaryVertexTrack.h.

43{};

◆ m_V

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

protectedinherited

Definition at line 47 of file TrigL2Vertex.h.

47{};

◆ m_Vqq

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

private

Definition at line 40 of file TrigPrimaryVertexTrack.h.

40{};

◆ m_Vuq

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

private

Definition at line 41 of file TrigPrimaryVertexTrack.h.

41{};

◆ m_Vuu

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

private

Definition at line 42 of file TrigPrimaryVertexTrack.h.

42{};

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

Public Member Functions

Protected Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TrigPrimaryVertexTrack() [1/2]

◆ TrigPrimaryVertexTrack() [2/2]

◆ ~TrigPrimaryVertexTrack()

Member Function Documentation

◆ activate()

◆ getChi2Contribution()

◆ getChi2Distance()

◆ getIndex()

◆ getTrackType()

◆ getTrigTrack()

◆ getTrkTrack()

◆ isActive()

◆ mask()

◆ Perigee()

◆ PerigeeCovariance()

◆ report()

◆ setIndex()

◆ updateVertex()

Member Data Documentation

◆ m_A

◆ m_active

◆ m_B

◆ m_D

◆ m_dChi2

◆ m_index

◆ m_nTrackType

◆ m_Perigee

◆ m_PerigeeCovariance

◆ m_pTrigTrack

◆ m_pTrkTrack

◆ m_q

◆ m_resid

◆ m_u

◆ m_V

◆ m_Vqq

◆ m_Vuq

◆ m_Vuu