#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 51 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 54 of file TrigInDetTrackFollowingTool.h.

54 {};

◆ TrigFTF_ExtendedTrackState() [2/2]

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

Definition at line 39 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 83 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 61 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 71 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 135 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 106 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
     }
   }
 }