TTrack.cxx File Reference

#include <cfloat>
#include "CLHEP/Matrix/DiagMatrix.h"
#include "CLHEP/Matrix/Matrix.h"
#include "CLHEP/Matrix/SymMatrix.h"
#include "CLHEP/Matrix/Vector.h"
#include "TrkReco/TCircle.h"
#include "TrkReco/TMDCUtil.h"
#include "TrkReco/TMDCWire.h"
#include "TrkReco/TMDCWireHit.h"
#include "TrkReco/TMDCWireHitMC.h"
#include "TrkReco/TMLine.h"
#include "TrkReco/TMLink.h"
#include "TrkReco/TSegment.h"
#include "TrkReco/TTrack.h"
#include "MdcTables/HepevtTables.h"
#include "MdcTables/MdcTables.h"
#include "MdcTables/MdstTables.h"
#include "MdcTables/TrkTables.h"
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/IInterface.h"
#include "GaudiKernel/IMessageSvc.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/Kernel.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/Service.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/StatusCode.h"

Go to the source code of this file.

Macros
#define	OPTNK
#define	NEW_FIT2D   1
#define	JJTEST   0
#define	JJJTEST   0
#define	t_dot(a, b)
#define	t_dot2(a, b)
#define	t_print(a, b)

Functions
int	SortByPt (const void *av, const void *bv)
	Utility functions.
std::string	TrackType (unsigned type)
Helix	Track2Helix (const MdcTrk_localz &t)
Helix	Track2Helix (const MdcRec_trk &t)
	returns helix.
Helix	Track2Helix (const Mdst_trk_fit &t)
std::string	TrackKinematics (const Helix &h)
std::string	TrackStatus (const TTrack &t)
	returns string of track status.
std::string	TrackStatus (const MdcRec_trk &c)
std::string	TrackStatus (const MdcRec_trk_add &a)
std::string	TrackStatus (unsigned md, unsigned mk, unsigned mq, unsigned ms, unsigned mm, unsigned ma)
std::string	TrackInformation (const TTrack &t)
std::string	TrackInformation (const MdcRec_trk &r)
std::string	TrackInformation (unsigned nA, unsigned nS, unsigned n, float chisq)
std::string	TrackLayerUsage (const TTrack &t)
bool	PositiveDefinite (const Helix &h)
	Error matrix validity.
bool	HelixHasNan (const Helix &h)
	Helix parameter validity.
Helix	Track2Helix (const Gen_hepevt &t)

Macro Definition Documentation

JJJTEST

#define JJJTEST   0

JJTEST

#define JJTEST   0

NEW_FIT2D

#define NEW_FIT2D   1

Definition at line 895 of file TTrack.cxx.

OPTNK

#define OPTNK

Definition at line 1 of file TTrack.cxx.

t_dot

#define t_dot	(		a,
			b )

Value:

( a[0] * b[0] + a[1] * b[1] + a[2] * b[2] )

Definition at line 3068 of file TTrack.cxx.

t_dot2

#define t_dot2	(		a,
			b )

Value:

( a[0] * b[0] + a[1] * b[1] )

Definition at line 3069 of file TTrack.cxx.

t_print

#define t_print	(		a,
			b )

Value:

std::cout << b << " = " << a[0] << " " << a[1] << " " << a[2] << std::endl;

Definition at line 3070 of file TTrack.cxx.

#  define t_print( a, b )                                                                     \
    std::cout << b << " = " << a[0] << " " << a[1] << " " << a[2] << std::endl;

Function Documentation

HelixHasNan()

bool HelixHasNan

(

const Helix &

h

)

Helix parameter validity.

Definition at line 3939 of file TTrack.cxx.

                                   {
  const Vector& a = h.a();
  for ( unsigned i = 0; i < 5; i++ )
    // maqm if (isnan(a[i]))
    if ( std::isnan( a[i] ) ) return true;
  const SymMatrix& Ea = h.Ea();
  for ( unsigned i = 0; i < 5; i++ )
    for ( unsigned j = 0; j <= i; j++ )
      // maqm if (isnan(Ea[i][j]))
      if ( std::isnan( Ea[i][j] ) ) return true;
  return false;
}

Referenced by TTrackManager::goodTrack(), TrackInformation(), and TrackInformation().

PositiveDefinite()

bool PositiveDefinite

(

const Helix &

h

)

Error matrix validity.

Definition at line 3924 of file TTrack.cxx.

                                        {
  const SymMatrix& e  = h.Ea();
  SymMatrix        e2 = e.sub( 1, 2 );
  SymMatrix        e3 = e.sub( 1, 3 );
  SymMatrix        e4 = e.sub( 1, 4 );
 
  bool positive = true;
  if ( e[0][0] <= 0. ) positive = false;
  else if ( e2.determinant() <= 0. ) positive = false;
  else if ( e3.determinant() <= 0. ) positive = false;
  else if ( e4.determinant() <= 0. ) positive = false;
  else if ( e.determinant() <= 0. ) positive = false;
  return positive;
}

Referenced by TrackInformation(), and TrackInformation().

SortByPt()

int SortByPt	(	const void *	av,
		const void *	bv )

Utility functions.

Definition at line 2514 of file TTrack.cxx.

                                                          {
  const TTrack** a( (const TTrack**)av );
  const TTrack** b( (const TTrack**)bv );
  if ( ( *a )->pt() < ( *b )->pt() ) return 1;
  else if ( ( *a )->pt() == ( *b )->pt() ) return 0;
  else return -1;
}

Referenced by TTrackManager::append2D().

Track2Helix() [1/4]

Helix Track2Helix

(

const Gen_hepevt &

t

)

Definition at line 3952 of file TTrack.cxx.

                                         {
  float charge = 1;
  if ( t.idhep == 11 ) charge = -1;
  else if ( t.idhep == -11 ) charge = 1;
  else if ( t.idhep == 13 ) charge = -1;
  else if ( t.idhep == -13 ) charge = 1;
  else if ( t.idhep == 211 ) charge = 1;
  else if ( t.idhep == -211 ) charge = -1;
  else if ( t.idhep == 321 ) charge = 1;
  else if ( t.idhep == -321 ) charge = -1;
  else if ( t.idhep == 2212 ) charge = 1;
  else if ( t.idhep == -2212 ) charge = -1;
  else
  {
    std::cout << "Track2Helix(gen_hepevt) !!! charge of id=";
    std::cout << t.idhep << " is unknown" << std::endl;
  }
 
  Hep3Vector mom( t.P[0], t.P[1], t.P[2] );
  Hep3Vector pos( t.V[0] / 10., t.V[1] / 10., t.V[2] / 10. );
  return Helix( pos, mom, charge );
}

Track2Helix() [2/4]

Helix Track2Helix

(

const MdcRec_trk &

t

)

returns helix.

Definition at line 3707 of file TTrack.cxx.

                                         {
  HepVector    a( 5 );
  Hep3Vector   p( t.pivot[0], t.pivot[1], t.pivot[2] );
  HepSymMatrix er( 5, 0 );
  a( 1 )     = t.helix[0];
  a( 2 )     = t.helix[1];
  a( 3 )     = t.helix[2];
  a( 4 )     = t.helix[3];
  a( 5 )     = t.helix[4];
  er( 1, 1 ) = t.error[0];
  er( 2, 1 ) = t.error[1];
  er( 2, 2 ) = t.error[2];
  er( 3, 1 ) = t.error[3];
  er( 3, 2 ) = t.error[4];
  er( 3, 3 ) = t.error[5];
  er( 4, 1 ) = t.error[6];
  er( 4, 2 ) = t.error[7];
  er( 4, 3 ) = t.error[8];
  er( 4, 4 ) = t.error[9];
  er( 5, 1 ) = t.error[10];
  er( 5, 2 ) = t.error[11];
  er( 5, 3 ) = t.error[12];
  er( 5, 4 ) = t.error[13];
  er( 5, 5 ) = t.error[14];
  return Helix( p, a, er );
}

Track2Helix() [3/4]

Helix Track2Helix

(

const MdcTrk_localz &

t

)

Definition at line 3680 of file TTrack.cxx.

                                            {
  HepVector    a( 5 );
  Hep3Vector   p( t.pivot[0], t.pivot[1], t.pivot[2] );
  HepSymMatrix er( 5, 0 );
  a( 1 )     = t.helix[0];
  a( 2 )     = t.helix[1];
  a( 3 )     = t.helix[2];
  a( 4 )     = t.helix[3];
  a( 5 )     = t.helix[4];
  er( 1, 1 ) = t.error[0];
  er( 2, 1 ) = t.error[1];
  er( 2, 2 ) = t.error[2];
  er( 3, 1 ) = t.error[3];
  er( 3, 2 ) = t.error[4];
  er( 3, 3 ) = t.error[5];
  er( 4, 1 ) = t.error[6];
  er( 4, 2 ) = t.error[7];
  er( 4, 3 ) = t.error[8];
  er( 4, 4 ) = t.error[9];
  er( 5, 1 ) = t.error[10];
  er( 5, 2 ) = t.error[11];
  er( 5, 3 ) = t.error[12];
  er( 5, 4 ) = t.error[13];
  er( 5, 5 ) = t.error[14];
  return Helix( p, a, er );
}

Referenced by TrackInformation().

Track2Helix() [4/4]

Helix Track2Helix

(

const Mdst_trk_fit &

t

)

Definition at line 3734 of file TTrack.cxx.

                                           {
  HepVector    a( 5 );
  Hep3Vector   p( t.pivot_x, t.pivot_y, t.pivot_z );
  HepSymMatrix er( 5, 0 );
  a( 1 )     = t.helix[0];
  a( 2 )     = t.helix[1];
  a( 3 )     = t.helix[2];
  a( 4 )     = t.helix[3];
  a( 5 )     = t.helix[4];
  er( 1, 1 ) = t.error[0];
  er( 2, 1 ) = t.error[1];
  er( 2, 2 ) = t.error[2];
  er( 3, 1 ) = t.error[3];
  er( 3, 2 ) = t.error[4];
  er( 3, 3 ) = t.error[5];
  er( 4, 1 ) = t.error[6];
  er( 4, 2 ) = t.error[7];
  er( 4, 3 ) = t.error[8];
  er( 4, 4 ) = t.error[9];
  er( 5, 1 ) = t.error[10];
  er( 5, 2 ) = t.error[11];
  er( 5, 3 ) = t.error[12];
  er( 5, 4 ) = t.error[13];
  er( 5, 5 ) = t.error[14];
  return Helix( p, a, er );
}

TrackInformation() [1/3]

std::string TrackInformation

(

const MdcRec_trk &

r

)

Definition at line 3878 of file TTrack.cxx.

                                                  {
  std::string p;
  if ( PositiveDefinite( Track2Helix( r ) ) ) p = " posi";
  else p = " nega";
  if ( HelixHasNan( Track2Helix( r ) ) ) p += " with NaN";
  return TrackInformation( r.nhits - r.nster, r.nster, r.nhits, r.chiSq ) + p;
}

TrackInformation() [2/3]

std::string TrackInformation

(

const TTrack &

t

)

Definition at line 3867 of file TTrack.cxx.

                                              {
  const AList<TMLink> cores = t.cores();
  unsigned            n     = cores.length();
  unsigned            nS    = NStereoHits( cores );
  unsigned            nA    = n - nS;
  std::string         p;
  if ( !PositiveDefinite( t.helix() ) ) p = " negative";
  if ( HelixHasNan( t.helix() ) ) p += " NaN";
  return TrackInformation( nA, nS, n, t.chi2() ) + p;
}

Referenced by TrackDump(), TrackInformation(), and TrackInformation().

TrackInformation() [3/3]

std::string TrackInformation	(	unsigned	nA,
		unsigned	nS,
		unsigned	n,
		float	chisq )

Definition at line 3886 of file TTrack.cxx.

                                                                                {
  std::string s;
 
  s += "a" + std::to_string( int( nA ) );
  s += " s" + std::to_string( int( nS ) );
  s += " n" + std::to_string( int( n ) );
  // s += " ndf" + std::string(int(r.m_ndf));
  float x = chisq;
 
  if ( x < 0. ) s += " -";
  else s += " ";
 
  int y = int( fabs( x ) );
  s += std::to_string( y ) + ".";
  float z = fabs( x );
  for ( unsigned j = 0; j < 1; j++ )
  {
    z *= 10.;
    y = ( int( z ) % 10 );
    s += std::to_string( y );
  }
 
  return s;
}

TrackKinematics()

std::string TrackKinematics

(

const Helix &

h

)

Definition at line 3761 of file TTrack.cxx.

                                            {
  static const HepPoint3D IP( 0., 0., 0. );
  Helix                   hIp = h;
  hIp.pivot( IP );
 
  float       chrg = hIp.a()[2] / fabs( hIp.a()[2] );
  std::string s;
  if ( chrg > 0. ) s = "+";
  else s = "-";
 
  float x[4];
  x[0] = fabs( hIp.a()[0] );
  x[1] = hIp.a()[3];
  x[2] = 1. / fabs( hIp.a()[2] );
  x[3] = ( 1. / fabs( hIp.a()[2] ) ) * hIp.a()[4];
 
  if ( ( x[0] < 2. ) && ( fabs( x[1] ) < 4. ) ) s += "i ";
  else s += "  ";
 
  for ( unsigned i = 0; i < 4; i++ )
  {
    if ( i ) s += " ";
 
    if ( x[i] < 0. ) s += "-";
    else s += " ";
 
    int y = int( fabs( x[i] ) );
    s += std::to_string( y ) + ".";
    float z = fabs( x[i] );
    for ( unsigned j = 0; j < 3; j++ )
    {
      z *= 10.;
      y = ( int( z ) % 10 );
      s += std::to_string( y );
    }
  }
 
  return s;
}

TrackLayerUsage()

std::string TrackLayerUsage

(

const TTrack &

t

)

Definition at line 3911 of file TTrack.cxx.

                                             {
  unsigned n[11];
  NHitsSuperLayer( t.links(), n );
  std::string nh;
  for ( unsigned i = 0; i < 11; i++ )
  {
    nh += std::to_string( n[i] );
    if ( i % 2 ) nh += "-";
    else if ( i < 10 ) nh += ",";
  }
  return nh;
}

TrackStatus() [1/4]

std::string TrackStatus

(

const MdcRec_trk &

c

)

Definition at line 3805 of file TTrack.cxx.

                                             {
  //    const reccdc_trk_add & a =
  //    * (reccdc_trk_add *) BsGetEnt(RECMDC_TRK_ADD, c.m_ID, BBS_No_Index);
  //    return TrackStatus(a);
  return std::string( "" );
}

TrackStatus() [2/4]

std::string TrackStatus

(

const MdcRec_trk_add &

a

)

Definition at line 3812 of file TTrack.cxx.

                                                 {
  //    return TrackStatus(a.decision,
  //               a.kind,
  //               a.quality,
  //               a.stat,
  //               a.mother,
  //               a.daughter);
  return std::string( "" );
}

TrackStatus() [3/4]

std::string TrackStatus

(

const TTrack &

t

)

returns string of track status.

Definition at line 3801 of file TTrack.cxx.

                                         {
  return TrackStatus( t.finder(), t.type(), t.quality(), t.fitting(), 0, 0 );
}

Referenced by TrackDump(), and TrackStatus().

TrackStatus() [4/4]

std::string TrackStatus	(	unsigned	md,
		unsigned	mk,
		unsigned	mq,
		unsigned	ms,
		unsigned	mm,
		unsigned	ma )

Definition at line 3822 of file TTrack.cxx.

                                       {
 
  std::string f;
  if ( md & TrackOldConformalFinder ) f += "o";
  if ( md & TrackFastFinder ) f += "f";
  if ( md & TrackSlowFinder ) f += "s";
  if ( md & TrackCurlFinder ) f += "c";
  if ( md & TrackTrackManager ) f += "t";
  if ( f == "" ) f = "?";
 
  std::string k;
  if ( mk & TrackTypeNormal ) k += "Norm";
  if ( mk & TrackTypeCurl ) k += "Curl";
  if ( mk & TrackTypeCircle ) k += "Circ";
  if ( mk & TrackTypeIncomingCosmic ) k += "Inco";
  if ( mk & TrackTypeOutgoingCosmic ) k += "Outc";
  if ( mk & TrackTypeKink ) k += "Kink";
  if ( mk & TrackTypeSVDOnly ) k += "Svd";
  if ( k == "" ) k = "?";
 
  std::string b;
  if ( mq & TrackQualityOutsideCurler ) b += "Curlback";
  if ( mq & TrackQualityAfterKink ) b += "Afterkink";
  if ( mq & TrackQualityCosmic ) b += "Cosmic";
  if ( mq & TrackQuality2D ) b += "2D";
  if ( b == "" ) b = "ok";
 
  std::string s;
  if ( ms & TrackFitGlobal ) s += "HFit";
  if ( ms & TrackFitCosmic ) s += "CFit";
  if ( ms & TrackFitCdcKalman ) s += "CKal";
  if ( ms & TrackFitSvdCdcKalman ) s += "SKal";
  if ( s == "" ) s = "?";
 
  int m = mm;
  if ( m ) --m;
 
  int d = ma;
  if ( d ) --d;
 
  std::string p = " ";
  return f + p + k + p + b + p + s + p + std::to_string( m ) + p + std::to_string( d );
}

TrackType()

std::string TrackType ( unsigned type )

inline

Definition at line 3055 of file TTrack.cxx.

                                     {
  switch ( type )
  {
  case TrackTypeUndefined: return std::string( "undefined" );
  case TrackTypeNormal: return std::string( "normal" );
  case TrackTypeCurl: return std::string( "curl  " );
  case TrackTypeCircle: return std::string( "circle" );
  case TrackTypeCosmic: return std::string( "cosmic" );
  }
  return std::string( "unknown  " );
}