TrigVertexFitConstraint Class Reference

#include <TrigL2Vertex.h>

Inheritance diagram for TrigVertexFitConstraint:

Collaboration diagram for TrigVertexFitConstraint:

Public Member Functions
	TrigVertexFitConstraint (double, const TrigVertexFitInputTrack *, const TrigVertexFitInputTrack *)
	two-track mass constraint More...
	TrigVertexFitConstraint (double, const TrigVertexFitInputTrack *, const TrigVertexFitInputTrack *, const TrigVertexFitInputTrack *)
	three-track mass constraint More...
	~TrigVertexFitConstraint ()
virtual double	getChi2Distance (class TrigL2Vertex *)
	implementation of abstract method from the base class More...
virtual void	updateVertex (class TrigL2Vertex *)
	implementation of abstract method from the base class More...
virtual MsgStream &	report (MsgStream &) const
double	getValue ()
	returns a mass of the constraint More...

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

Private Member Functions
double	calculateInvariantMass (TrigL2Vertex *pV)

Private Attributes
std::list< const TrigVertexFitInputTrack * >	m_trackList
double	m_value

Detailed Description

Definition at line 104 of file TrigL2Vertex.h.

Constructor & Destructor Documentation

TrigVertexFitConstraint() [1/2]

TrigVertexFitConstraint::TrigVertexFitConstraint	(	double	m,
		const TrigVertexFitInputTrack *	pT1,
		const TrigVertexFitInputTrack *	pT2
	)

two-track mass constraint

Definition at line 342 of file TrigL2Vertex.cxx.

                                                             : m_value(m)
{
  m_trackList.clear();
  m_trackList.push_back(pT1);
  m_trackList.push_back(pT2);
}

TrigVertexFitConstraint() [2/2]

TrigVertexFitConstraint::TrigVertexFitConstraint	(	double	m,
		const TrigVertexFitInputTrack *	pT1,
		const TrigVertexFitInputTrack *	pT2,
		const TrigVertexFitInputTrack *	pT3
	)

three-track mass constraint

Definition at line 351 of file TrigL2Vertex.cxx.

                                                             : m_value(m)
{
  m_trackList.clear();
  m_trackList.push_back(pT1);m_trackList.push_back(pT2);m_trackList.push_back(pT3);
}

~TrigVertexFitConstraint()

TrigVertexFitConstraint::~TrigVertexFitConstraint

(

)

Definition at line 361 of file TrigL2Vertex.cxx.

{
  m_trackList.clear();
}

Member Function Documentation

calculateInvariantMass()

double TrigVertexFitConstraint::calculateInvariantMass ( TrigL2Vertex *  pV )

private

Definition at line 366 of file TrigL2Vertex.cxx.

{
  const double C=0.029997;
  const double B=20.84;
 
  double invMass=0.0,alpha=C*B/1000.0;
  double P[3];double E=0.0;
  
  double* Rk = pV->getParametersVector();
 
  int offset=0;
  for(int i=0;i<3;i++) P[i]=0.0;
 
  for(std::list<const TrigVertexFitInputTrack*>::iterator it=m_trackList.begin();it!=m_trackList.end();++it)
    {
      offset=3+3*(*it)->getIndex();
      double mass=(*it)->getMass()/1000.0;
      double pT=fabs(Rk[offset+2]);
      double p=pT/sin(Rk[offset+1]);
 
      double psi=-asin(alpha*(Rk[0]*cos(Rk[offset])+Rk[1]*sin(Rk[offset]))/Rk[offset+2]);
      double phiV=Rk[offset]+psi;
      P[0]+=pT*cos(phiV);
      P[1]+=pT*sin(phiV);
      P[2]+=pT*cos(Rk[offset+1])/sin(Rk[offset+1]);
      E+=sqrt(mass*mass+p*p);
    }
  invMass=sqrt(E*E-P[0]*P[0]-P[1]*P[1]-P[2]*P[2]);
 
  return invMass;
}

getChi2Distance()

double TrigVertexFitConstraint::getChi2Distance ( class TrigL2Vertex *  pV )

virtual

implementation of abstract method from the base class

Implements TrigVertexFittingNode.

Definition at line 399 of file TrigL2Vertex.cxx.

{
  const double C=0.029997;
  const double B=20.84;
 
  double invMass=0.0,alpha=C*B/1000.0;
  double P[3];
  double E=0.0;
  int i{0},j{0};
  bool linFailed=false;
  double* Rk = pV->getParametersVector();
 
  int offset=0;
  for(i=0;i<3;i++) P[i]=0.0;
 
  int nSize=3+3*pV->getTracks()->size();
 
  for(std::list<const TrigVertexFitInputTrack*>::iterator it=m_trackList.begin();it!=m_trackList.end();++it)
    {
      offset=3+3*(*it)->getIndex();
      double mass=(*it)->getMass()/1000.0;
      double pT=fabs(Rk[offset+2]);
      double p=pT/sin(Rk[offset+1]);
 
      double psi=-asin(alpha*(Rk[0]*cos(Rk[offset])+Rk[1]*sin(Rk[offset]))/Rk[offset+2]);
      double phiV=Rk[offset]+psi;
      P[0]+=pT*cos(phiV);
      P[1]+=pT*sin(phiV);
      P[2]+=pT*cos(Rk[offset+1])/sin(Rk[offset+1]);
      E+=sqrt(mass*mass+p*p);
    }
  invMass=sqrt(E*E-P[0]*P[0]-P[1]*P[1]-P[2]*P[2]);
 
  m_resid[0]=m_value/1000.0-invMass;
  m_D[0][0]=0.0;m_D[0][1]=0.0;m_D[0][2]=0.0;
 
  for(std::list<const TrigVertexFitInputTrack*>::iterator it=m_trackList.begin();it!=m_trackList.end();++it)
    {
      offset=3+3*(*it)->getIndex();
 
      double mass=(*it)->getMass()/1000.0;
      double Ck=(Rk[offset+2]<0.0)?-1.0:1.0;
      double sinT=sin(Rk[offset+1]);
      double cosT=cos(Rk[offset+1]);
      double pT=fabs(Rk[offset+2]);
      double p=pT/sinT;
      double e=sqrt(p*p+mass*mass);
      double sinF=sin(Rk[offset]);
      double cosF=cos(Rk[offset]);
 
      double sinPsi=-alpha*(Rk[0]*cosF+Rk[1]*sinF)/Rk[offset+2];
      if(fabs(sinPsi)>1.0)
    {
      linFailed=true;
      break;
    }
      double psi=asin(sinPsi);
      double cosPsi=sqrt(1.0-sinPsi*sinPsi);
      double phiV=Rk[offset]+psi;
      
      double aCos=alpha*Ck/cosPsi;
      double dP=P[1]*cos(phiV)-P[0]*sin(phiV);
      double eE=E*Rk[offset+2]/(e*sinT);
 
      m_D[0][0]+=dP*cosF*aCos;
      m_D[0][1]+=dP*sinF*aCos;
      m_D[0][offset]=-dP*Ck*(Rk[offset+2]-Ck*aCos*(Rk[1]*cosF-Rk[0]*sinF));
      m_D[0][offset+1]=(Rk[offset+2]/(sinT*sinT))*(P[2]*Ck-eE*cosT);
      m_D[0][offset+2]=eE/sinT-Ck*(P[0]*cos(phiV)+P[1]*sin(phiV)+P[2]*cosT/sinT)+dP*Ck*sinPsi/cosPsi;     
    }
  for(i=0;i<nSize;i++) m_D[0][i]/=invMass;
  if(linFailed) return -999.9;
  /*
  for(i=0;i<nSize;i++) printf("%E ",m_D[0][i]);
  printf("\n");
  printf("----- numerical Jacobian -----\n");
 
  double oldPar, newPar, Delta, newMass, der;
  for(i=0;i<nSize;i++)
    {
      oldPar=Rk[i];
      Delta=0.00001*oldPar;
      newPar=oldPar+Delta;
      Rk[i]=newPar;
      newMass=calculateInvariantMass(pV);
      der=(newMass-invMass)/Delta;
      m_D[1][i]=der;
      Rk[i]=oldPar;
    }
  for(i=0;i<nSize;i++) printf("%E ",m_D[1][i]);
  printf("\n");
  printf("------------------------------\n");
  */
 
  double covM=1e-12;
 
  for(i=0;i<nSize;i++)
    {
      for(j=i;j<nSize;j++) 
    {
      double dCov=pV->m_Gk[i][j]*m_D[0][i]*m_D[0][j];
      if(i!=j) dCov*=2.0;
      covM+=dCov;
    }
    }
  m_V[0][0]=1.0/covM;
  return ((m_resid[0]*m_resid[0])*m_V[0][0]);
}

getValue()

double TrigVertexFitConstraint::getValue

(

)

returns a mass of the constraint

Definition at line 539 of file TrigL2Vertex.cxx.

{
  return m_value;
}

report()

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

virtual

Implements TrigVertexFittingNode.

Definition at line 533 of file TrigL2Vertex.cxx.

{
  out<<"Mass constraint with "<<m_trackList.size()<<" : mass="<<m_value<<endmsg;
  return out;
}

updateVertex()

void TrigVertexFitConstraint::updateVertex ( class TrigL2Vertex *  pV )

virtual

implementation of abstract method from the base class

Implements TrigVertexFittingNode.

Definition at line 508 of file TrigL2Vertex.cxx.

{
  double K[MAX_SIZE_VERT_COVM];
  int i{0},j{0};
  const int nSize=3+3*pV->getTracks()->size();
  double gain{0};
 
  for(i=0;i<nSize;i++)
    {
      gain=0.0;
      for(j=0;j<nSize;j++) gain+=pV->m_Gk[i][j]*m_D[0][j];
      m_D[1][i]=gain;
      K[i]=gain*m_V[0][0];
    }
  for(i=0;i<nSize;i++)
    {
      pV->getParametersVector()[i]+=K[i]*m_resid[0];
      for(j=i;j<nSize;j++)
    {
      pV->m_Gk[i][j]-=K[i]*m_D[1][j];
      pV->m_Gk[j][i]=pV->m_Gk[i][j];
    }
    }
}

Member Data Documentation

m_D

double TrigVertexFittingNode::m_D[2][MAX_SIZE_VERT_COVM]

protectedinherited

Definition at line 48 of file TrigL2Vertex.h.

m_resid

double TrigVertexFittingNode::m_resid[2]

protectedinherited

Definition at line 46 of file TrigL2Vertex.h.

m_trackList

std::list<const TrigVertexFitInputTrack*> TrigVertexFitConstraint::m_trackList

private

Definition at line 117 of file TrigL2Vertex.h.

m_V

double TrigVertexFittingNode::m_V[2][2]

protectedinherited

Definition at line 47 of file TrigL2Vertex.h.

m_value

double TrigVertexFitConstraint::m_value

private

Definition at line 118 of file TrigL2Vertex.h.

---

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

• TrigL2Vertex.h
• TrigL2Vertex.cxx