#include <TrigInDetTrackFollowingTool.h>

Collaboration diagram for TrigFTF_ExtendedTrackState:

Public Member Functions
	TrigFTF_ExtendedTrackState ()=delete
	~TrigFTF_ExtendedTrackState ()
TrigFTF_ExtendedTrackState &	operator= (const TrigFTF_ExtendedTrackState &)=delete
	TrigFTF_ExtendedTrackState (double const , const Trk::PlaneSurface )
void	AddHole ()
void	AddHit (const Trk::PrepRawData *, double, int)
void	SwapTheEnds ()
void	correctAngles ()
void	report () const

Public Attributes
double	m_Xk [10]
double	m_Gk [10][10]
double	m_chi2
double	m_ndof
const Trk::PlaneSurface *	m_pS
const Trk::PlaneSurface *	m_pO
std::list< TrigFTF_HitAssignment >	m_track
int	m_nClusters
int	m_nHoles
bool	m_isSwapped

Detailed Description

Definition at line 52 of file TrigInDetTrackFollowingTool.h.

Constructor & Destructor Documentation

◆ TrigFTF_ExtendedTrackState() [1/2]

TrigFTF_ExtendedTrackState::TrigFTF_ExtendedTrackState ( )

delete

◆ ~TrigFTF_ExtendedTrackState()

TrigFTF_ExtendedTrackState::~TrigFTF_ExtendedTrackState ( )

inline

Definition at line 55 of file TrigInDetTrackFollowingTool.h.

55{};

◆ TrigFTF_ExtendedTrackState() [2/2]

TrigFTF_ExtendedTrackState::TrigFTF_ExtendedTrackState	(	double const *	P,
		const Trk::PlaneSurface *	pS )

Definition at line 40 of file TrigInDetTrackFollowingTool.cxx.

                                                                                                 : m_chi2(0), m_ndof(-5), m_pS(pS), m_pO(pS), m_nClusters(0), m_nHoles(0), m_isSwapped(false) {
 
  for(int i=0;i<5;i++) m_Xk[i] = m_Xk[i+5] = P[i];//duplicate the lower half of the state vector
 
  memset(&m_Gk[0][0],0,sizeof(m_Gk));
 
  //covariance block C
 
  m_Gk[0][0] = 1.0;
  m_Gk[1][1] = 1.0;
  m_Gk[2][2] = 0.0001;
  m_Gk[3][3] = 0.0001;
  m_Gk[4][4] = P[5];
 
  //triplicate the lower block
  for(int i=0;i<5;i++) {
    m_Gk[i+5][i+5] = m_Gk[i][i+5] = m_Gk[i+5][i] = m_Gk[i][i];
  }
 
  m_track.clear();
}

Member Function Documentation

◆ AddHit()

void TrigFTF_ExtendedTrackState::AddHit	(	const Trk::PrepRawData *	pPRD,
		double	dchi2,
		int	ndof )

Definition at line 84 of file TrigInDetTrackFollowingTool.cxx.

                                                                                          {
  
  correctAngles();
 
  double cov[15];
  int idx=0;
  for(int i=0;i<5;i++) 
    for(int j=0;j<=i;j++) 
      cov[idx++] = m_Gk[i][j];
 
  if(m_isSwapped) {
    m_track.emplace_front(pPRD, m_Xk, &cov[0], dchi2, ndof);
  }
  else {
    m_track.emplace_back(pPRD, m_Xk, &cov[0], dchi2, ndof);
  }
  
  m_chi2 += dchi2;
  m_ndof += ndof;
  m_nClusters++;
  m_nHoles = 0;//reset the counter
}

◆ AddHole()

void TrigFTF_ExtendedTrackState::AddHole ( )

Definition at line 62 of file TrigInDetTrackFollowingTool.cxx.

                                         {
  if(m_isSwapped) {
    m_track.emplace_front(nullptr, nullptr, nullptr, 0, -1);
  }
  else {
    m_track.emplace_back(nullptr, nullptr, nullptr, 0, -1);
  }
  m_nHoles++;
}

◆ correctAngles()

void TrigFTF_ExtendedTrackState::correctAngles ( )

Definition at line 72 of file TrigInDetTrackFollowingTool.cxx.

                                               {
  if (m_Xk[2] >  M_PI) m_Xk[2] -= 2 * M_PI;
  if (m_Xk[2] < -M_PI) m_Xk[2] += 2 * M_PI;
  if (m_Xk[3] < 0.0)   m_Xk[3] += M_PI;
  if (m_Xk[3] > M_PI)  m_Xk[3] -= M_PI;
  if (m_Xk[7] >  M_PI) m_Xk[7] -= 2 * M_PI;
  if (m_Xk[7] < -M_PI) m_Xk[7] += 2 * M_PI;
  if (m_Xk[8] < 0.0)   m_Xk[8] += M_PI;
  if (m_Xk[8] > M_PI)  m_Xk[8] -= M_PI;
 
}

◆ operator=()

TrigFTF_ExtendedTrackState & TrigFTF_ExtendedTrackState::operator= ( const TrigFTF_ExtendedTrackState & )

delete

◆ report()

void TrigFTF_ExtendedTrackState::report ( ) const

Definition at line 136 of file TrigInDetTrackFollowingTool.cxx.

                                              {
  std::cout<<"L: ";
  for(int i=0;i<4;i++) std::cout<<m_Xk[i]<<" ";
  std::cout<<1/m_Xk[4]<<" "<<std::sin(m_Xk[3])/m_Xk[4]<<std::endl;
  
 
  std::cout<<"Covariance at last point:"<<std::endl;
 
  for(int i=0;i<5;i++) {
    for(int j=0;j<5;j++) std::cout<<m_Gk[i][j]<<" ";
    std::cout<<std::endl;
  }
 
  std::cout<<"F: ";
  for(int i=0;i<4;i++) std::cout<<m_Xk[i+5]<<" ";
  std::cout<<1/m_Xk[4+5]<<" "<<std::sin(m_Xk[3+5])/m_Xk[4+5]<<std::endl;
  std::cout<<std::endl;
 
  std::cout<<"Covariance at the first point:"<<std::endl;
 
  for(int i=0;i<5;i++) {
    for(int j=0;j<5;j++) std::cout<<m_Gk[i+5][j+5]<<" ";
    std::cout<<std::endl;
  }
 
  std::cout<<"chi2="<<m_chi2<<" ndof="<<m_ndof<<" nClusters="<<m_nClusters<<" nHoles="<<m_nHoles<<std::endl;
 
}

◆ SwapTheEnds()

void TrigFTF_ExtendedTrackState::SwapTheEnds ( )

Definition at line 107 of file TrigInDetTrackFollowingTool.cxx.

                                             {
 
  m_isSwapped = !m_isSwapped;
 
  std::swap(m_pS, m_pO);
  
  double tmpX[10];
  memcpy(&tmpX[0], &m_Xk[0], sizeof(tmpX));
 
  for(int i=0;i<5;i++) {
    m_Xk[i]   = tmpX[i+5];
    m_Xk[i+5] = tmpX[i];
  }
 
  double tmpG[10][10];
  memcpy(&tmpG[0][0], &m_Gk[0][0], sizeof(tmpG));
 
  //covariance blocks are moved around to reflect the state vectors swap done above
 
  for(int i=0;i<5;i++) {
    for(int j=0;j<5;j++) {
      m_Gk[i+5][j+5] = tmpG[i][j];    //C -> E
      m_Gk[i][j]     = tmpG[i+5][j+5];//E -> C
      m_Gk[i][j+5]   = tmpG[i+5][j];  //D -> B
      m_Gk[i+5][j]   = tmpG[i][j+5];  //B -> D
    }
  }
}