TrigVertexFitInputTrack Class Reference

#include <TrigL2Vertex.h>

Inheritance diagram for TrigVertexFitInputTrack:

Collaboration diagram for TrigVertexFitInputTrack:

Public Member Functions
	TrigVertexFitInputTrack (const TrigInDetTrack *, double)
	constructor for L2 tracks
	TrigVertexFitInputTrack (const Trk::Track *, double)
	constructor for EF (offline) tracks
	~TrigVertexFitInputTrack ()
const TrigInDetTrack *	getTrigTrack ()
	getter for L2 tracks
const Trk::Track *	getTrkTrack ()
	getter for EF (offline) tracks
void	initializeVertex (class TrigL2Vertex *)
	resets and fills its part of the covariance and parameter vector
bool	linearize (class TrigL2Vertex *)
	re-calculates linearization of the measurement model in the vicinity of parameters provided by input TrigL2Vertex
virtual double	getChi2Distance (class TrigL2Vertex *)
	implementation of abstract method from the base class
virtual void	updateVertex (class 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 isActive to true
void	mask ()
	sets isActive to false
void	setMass (double)
	sets a mass of this particle
double	getMass () const
	gets a mass of this particle
const double *	Perigee () const
	track parameters at the perigee
double	PerigeeCovariance (int, int) const
	covariance of track parameters at the perigee

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 {}
double	m_mass {}
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_h [2] {}
bool	m_active {true}

Detailed Description

Definition at line 61 of file TrigL2Vertex.h.

Constructor & Destructor Documentation

TrigVertexFitInputTrack() [1/2]

TrigVertexFitInputTrack::TrigVertexFitInputTrack	(	const TrigInDetTrack *	pT,
		double	mass = 0.0 )

constructor for L2 tracks

Definition at line 9 of file TrigL2Vertex.cxx.

                                                                                          :
  m_pTrigTrack(pT),
  m_nTrackType(1),
  m_mass(mass)
{
  double Ck[5][5];
  int i{0},j{0};
 
  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()/1000.0;      
 
 
  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]/1000.0;Ck[i][4]=Ck[4][i];
    }
  Ck[4][4]/=1000.0; 
 
 
  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_Vuu[0][0]=(*p->cov())[0];
  m_Vuu[0][1]=m_Vuu[1][0]=(*p->cov())[2];
  m_Vuu[1][1]=(*p->cov())[9];
 
  m_Vuq[0][0]=(*p->cov())[1];m_Vuq[0][1]=(*p->cov())[3]*a;
  m_Vuq[0][2]=(*p->cov())[4]*1e-3;m_Vuq[1][0]=(*p->cov())[6];m_Vuq[1][1]=(*p->cov())[10]*a;
  m_Vuq[1][2]=(*p->cov())[11]*1e-3;
 
  m_Vqq[0][0]=(*p->cov())[5];m_Vqq[0][1]=m_Vqq[1][0]=(*p->cov())[7]*a;
  m_Vqq[0][2]=m_Vqq[2][0]=(*p->cov())[8]*1e-3;
  m_Vqq[1][1]=(*p->cov())[12]*a*a;m_Vqq[1][2]=m_Vqq[2][1]=(*p->cov())[13]*a*1e-3;
  m_Vqq[2][2]=(*p->cov())[14]*1e-6;
  */
  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];
}

TrigVertexFitInputTrack() [2/2]

TrigVertexFitInputTrack::TrigVertexFitInputTrack	(	const Trk::Track *	pT,
		double	mass = 0.0 )

constructor for EF (offline) tracks

Definition at line 68 of file TrigL2Vertex.cxx.

                                                                                   :
  m_pTrkTrack(pT),
  m_nTrackType(2),
  m_mass(mass)
{
  double Ck[5][5];
 
  const Trk::Perigee* pP= pT->perigeeParameters();
 
  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];
 
  double ptC=1000.0*pP->parameters()[Trk::qOverP];
 
  m_q[2]=sin(pP->parameters()[Trk::theta])/ptC;
 
    if( pP->covariance() ){
        // has covariance matrix
        const AmgSymMatrix(5)* TC = pP->covariance();
 
        Ck[0][0]=(*TC)(Trk::d0,Trk::d0);
        Ck[0][1]=Ck[1][0]=(*TC)(Trk::d0,Trk::z0);
        Ck[0][2]=Ck[2][0]=(*TC)(Trk::d0,Trk::phi0);
        Ck[0][3]=Ck[3][0]=(*TC)(Trk::d0,Trk::theta);
        Ck[0][4]=Ck[4][0]=(*TC)(Trk::d0,Trk::qOverP);
        Ck[1][1]=(*TC)(Trk::z0,Trk::z0);
        Ck[1][2]=Ck[2][1]=(*TC)(Trk::z0,Trk::phi0);
        Ck[1][3]=Ck[3][1]=(*TC)(Trk::z0,Trk::theta);
        Ck[1][4]=Ck[4][1]=(*TC)(Trk::z0,Trk::qOverP);
        Ck[2][2]=(*TC)(Trk::phi0,Trk::phi0);
        Ck[2][3]=Ck[3][2]=(*TC)(Trk::phi0,Trk::theta);
        Ck[2][4]=Ck[4][2]=(*TC)(Trk::phi0,Trk::qOverP);
        Ck[3][3]=(*TC)(Trk::theta,Trk::theta);
        Ck[3][4]=Ck[4][3]=(*TC)(Trk::theta,Trk::qOverP);
        Ck[4][4]=(*TC)(Trk::qOverP,Trk::qOverP);
 
        const double a = cos(pP->parameters()[Trk::theta])/ptC;
        const double b = -sin(pP->parameters()[Trk::theta])/(pP->parameters()[Trk::qOverP]*ptC);
 
        Ck[3][3]=Ck[3][3]+2.0*a*Ck[3][4]+a*a*Ck[4][4];
        Ck[3][4]=Ck[4][3]=b*Ck[3][4]+a*b*Ck[4][4];
        Ck[4][4]=b*b*Ck[4][4];
        Ck[0][3]=Ck[3][0]=Ck[0][3]+a*Ck[0][4];Ck[0][4]*=b;Ck[4][0]=Ck[0][4];
        Ck[1][3]=Ck[3][1]=Ck[1][3]+a*Ck[1][4];Ck[1][4]*=b;Ck[4][1]=Ck[1][4];
        Ck[2][3]=Ck[3][2]=Ck[2][3]+a*Ck[2][4];Ck[2][4]*=b;Ck[4][2]=Ck[2][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];
    }
    else{
        m_active=false;
        // no covariance
    }
 
}

~TrigVertexFitInputTrack()

TrigVertexFitInputTrack::~TrigVertexFitInputTrack

(

)

Definition at line 128 of file TrigL2Vertex.cxx.

{
 
}

Member Function Documentation

activate()

void TrigVertexFitInputTrack::activate

(

)

sets isActive to true

Definition at line 143 of file TrigL2Vertex.cxx.

{
  m_active=true;
}

getChi2Distance()

double TrigVertexFitInputTrack::getChi2Distance ( class TrigL2Vertex * pV )

virtual

implementation of abstract method from the base class

Implements TrigVertexFittingNode.

Definition at line 196 of file TrigL2Vertex.cxx.

{
  const double C=0.02997;
  const double B=20.84;
  const double alpha=C*B/1000.0;
    
  double Sk[2][2],detr,chi2;
  double AC[2][3],BV[2][3];
  int i{0},j{0},k{0};
 
  const int shift=m_index*3;
 
  k=3+shift;
  for(i=0;i<3;i++)
    for(j=0;j<3;j++)
      pV->m_Gk[i+k][j+k]=m_Vqq[i][j];  
 
  initializeVertex(pV);
  const double xv=pV->getParametersVector()[0];
  const double yv=pV->getParametersVector()[1];
  const double zv=pV->getParametersVector()[2];
  const double phi0=pV->getParametersVector()[3+shift];
  const double theta0=pV->getParametersVector()[4+shift];
  const double P0=pV->getParametersVector()[5+shift];
 
  const double cosPhi0=cos(phi0);
  const double sinPhi0=sin(phi0);
  const double 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 psi=asin(sinPsi);
  const double sint=sin(theta0);
  const double ctt=cos(theta0)/sint;
 
  m_A[0][0]=-sin(phi0+psi)/cosPsi;
  m_A[0][1]= cos(phi0+psi)/cosPsi;
  m_A[0][2]=0.0;
 
  m_A[1][0]=-ctt*cosPhi0/cosPsi;
  m_A[1][1]=-ctt*sinPhi0/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-1.0/cosPsi)/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/cosPsi)/alpha;
 
  m_h[0]=yv*cosPhi0-xv*sinPhi0+P0*(1-cosPsi)/alpha;
  m_h[1]=zv+P0*ctt*psi/alpha;
 
  m_resid[0]=m_u[0]-m_h[0];
  m_resid[1]=m_u[1]-m_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];
      for(j=3;j<3+shift;j++) 
    {
      m_D[i][j]=0.0;
      for(k=0;k<3;k++)
        m_D[i][j]+=m_A[i][k]*pV->m_Gk[k][j];
    }
      for(j=0;j<3;j++) m_D[i][j+3+shift]=BV[i][j]-m_Vuq[i][j];
    }
    
  return chi2;
}

getIndex()

int TrigVertexFitInputTrack::getIndex

(

)

const

to be used by TrigVertexingTool

Definition at line 174 of file TrigL2Vertex.cxx.

{
  return m_index;
}

getMass()

double TrigVertexFitInputTrack::getMass ( ) const

inline

gets a mass of this particle

Definition at line 81 of file TrigL2Vertex.h.

getTrackType()

int TrigVertexFitInputTrack::getTrackType

(

)

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

Definition at line 153 of file TrigL2Vertex.cxx.

{
  return m_nTrackType;
}

getTrigTrack()

const TrigInDetTrack * TrigVertexFitInputTrack::getTrigTrack

(

)

getter for L2 tracks

Definition at line 159 of file TrigL2Vertex.cxx.

{
  return m_pTrigTrack;
}

getTrkTrack()

const Trk::Track * TrigVertexFitInputTrack::getTrkTrack

(

)

getter for EF (offline) tracks

Definition at line 164 of file TrigL2Vertex.cxx.

{
    return m_pTrkTrack;
}

initializeVertex()

void TrigVertexFitInputTrack::initializeVertex

(

class TrigL2Vertex *

pV

)

resets and fills its part of the covariance and parameter vector

Definition at line 189 of file TrigL2Vertex.cxx.

{
  const int i = 3+m_index*3;
  for(int j=0;j<3;j++) 
    pV->getParametersVector()[i+j]=m_q[j];  
}

isActive()

bool TrigVertexFitInputTrack::isActive

(

)

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

Definition at line 148 of file TrigL2Vertex.cxx.

{
  return m_active;
}

linearize()

bool TrigVertexFitInputTrack::linearize

(

class TrigL2Vertex *

)

re-calculates linearization of the measurement model in the vicinity of parameters provided by input TrigL2Vertex

mask()

void TrigVertexFitInputTrack::mask

(

)

sets isActive to false

Definition at line 138 of file TrigL2Vertex.cxx.

{
  m_active=false;
}

Perigee()

const double * TrigVertexFitInputTrack::Perigee

(

)

const

track parameters at the perigee

Definition at line 179 of file TrigL2Vertex.cxx.

{
  return &m_Perigee[0];
}

PerigeeCovariance()

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

covariance of track parameters at the perigee

Definition at line 184 of file TrigL2Vertex.cxx.

{
  return m_PerigeeCovariance[i][j];
}

report()

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

virtual

Implements TrigVertexFittingNode.

Definition at line 329 of file TrigL2Vertex.cxx.

{
  out<<"Track "<<m_index<<" : mass="<<m_mass<<endmsg;
  for(int i=0;i<2;i++)
    {
      out<<"  u"<<i<<" = "<<m_u[i]<<"      "<<m_Vuu[i][0]<<"   "<<m_Vuu[i][1]<<endmsg;
    }
  for(int 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 TrigVertexFitInputTrack::setIndex

(

int

i

)

to be used by TrigVertexingTool

Definition at line 169 of file TrigL2Vertex.cxx.

{
  m_index=i;
}

setMass()

void TrigVertexFitInputTrack::setMass

(

double

m

)

sets a mass of this particle

Definition at line 133 of file TrigL2Vertex.cxx.

{
  m_mass=m;
}

updateVertex()

void TrigVertexFitInputTrack::updateVertex ( class TrigL2Vertex * pV )

virtual

implementation of abstract method from the base class

Implements TrigVertexFittingNode.

Definition at line 303 of file TrigL2Vertex.cxx.

{
  int i{0},j{0},k{0};
  const int nSize=6+3*m_index;
 
  double K[2][MAX_SIZE_VERT_COVM];
 
  for(i=0;i<2;i++)
    {
      for(j=0;j<nSize;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<nSize;i++)
    {
      pV->getParametersVector()[i]+=K[0][i]*m_resid[0]+K[1][i]*m_resid[1];
      for(j=i;j<nSize;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 TrigVertexFitInputTrack::m_A[2][3] {}

private

Definition at line 97 of file TrigL2Vertex.h.

{};

m_active

bool TrigVertexFitInputTrack::m_active {true}

private

Definition at line 101 of file TrigL2Vertex.h.

{true};

m_B

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

private

Definition at line 98 of file TrigL2Vertex.h.

{};

m_D

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

protectedinherited

Definition at line 48 of file TrigL2Vertex.h.

{};

m_h

double TrigVertexFitInputTrack::m_h[2] {}

private

Definition at line 99 of file TrigL2Vertex.h.

{};

m_index

int TrigVertexFitInputTrack::m_index {}

private

Definition at line 88 of file TrigL2Vertex.h.

{};

m_mass

double TrigVertexFitInputTrack::m_mass {}

private

Definition at line 89 of file TrigL2Vertex.h.

{};

m_nTrackType

int TrigVertexFitInputTrack::m_nTrackType {}

private

Definition at line 87 of file TrigL2Vertex.h.

{};

m_Perigee

double TrigVertexFitInputTrack::m_Perigee[5] {}

private

Definition at line 95 of file TrigL2Vertex.h.

{};

m_PerigeeCovariance

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

private

Definition at line 96 of file TrigL2Vertex.h.

{};

m_pTrigTrack

const TrigInDetTrack* TrigVertexFitInputTrack::m_pTrigTrack {}

private

Definition at line 85 of file TrigL2Vertex.h.

{};

m_pTrkTrack

const Trk::Track* TrigVertexFitInputTrack::m_pTrkTrack {}

private

Definition at line 86 of file TrigL2Vertex.h.

{};

m_q

double TrigVertexFitInputTrack::m_q[3] {}

private

Definition at line 94 of file TrigL2Vertex.h.

{};

m_resid

double TrigVertexFittingNode::m_resid[2] {}

protectedinherited

Definition at line 46 of file TrigL2Vertex.h.

{};

m_u

double TrigVertexFitInputTrack::m_u[2] {}

private

Definition at line 93 of file TrigL2Vertex.h.

{};

m_V

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

protectedinherited

Definition at line 47 of file TrigL2Vertex.h.

{};

m_Vqq

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

private

Definition at line 90 of file TrigL2Vertex.h.

{};

m_Vuq

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

private

Definition at line 91 of file TrigL2Vertex.h.

{};

m_Vuu

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

private

Definition at line 92 of file TrigL2Vertex.h.

{};

---

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

• TrigL2Vertex.h
• TrigL2Vertex.cxx