TTrack.h File Reference

#include <string>
#include "CLHEP/Vector/ThreeVector.h"
#include "TrackUtil/Helix.h"
#include "TrkReco/TMDCUtil.h"
#include "TrkReco/TTrackBase.h"
#include "TrkReco/TMLink.h"
#include "TrkReco/TBuilder0.h"
#include "TrkReco/TBuilderCosmic.h"
#include "TrkReco/TBuilderCurl.h"
#include "TrkReco/TPoint2D.h"
#include "TrkReco/T3DLine.h"
#include "TrkReco/TRunge.h"
#include "GaudiKernel/IInterface.h"
#include "GaudiKernel/Kernel.h"
#include "GaudiKernel/Service.h"
#include "MagneticField/IMagneticFieldSvc.h"
#include "MagneticField/MagneticFieldSvc.h"
#include "CLHEP/Matrix/Vector.h"
#include "CLHEP/Matrix/SymMatrix.h"
#include "CLHEP/Vector/LorentzVector.h"
#include "CLHEP/Geometry/Point3D.h"

Go to the source code of this file.

Classes
class	TTrack
	A class to represent a track in tracking. More...

Macros
#define	TrackOldConformalFinder   1
#define	TrackFastFinder   2
#define	TrackSlowFinder   4
#define	TrackCurlFinder   8
#define	TrackTrackManager   16
#define	TrackPMCurlFinder   32
#define	TrackSVDAssociator   64
#define	TrackFinderMask   255
#define	TrackTypeUndefined   0
#define	TrackTypeNormal   1
#define	TrackTypeCurl   2
#define	TrackTypeCircle   4
#define	TrackTypeCosmic   8
#define	TrackTypeIncomingCosmic   8
#define	TrackTypeOutgoingCosmic   16
#define	TrackTypeKink   32
#define	TrackTypeSVDOnly   256
#define	TrackQualityOutsideCurler   1
#define	TrackQualityAfterKink   2
#define	TrackQualityCosmic   4
#define	TrackQuality2D   8
#define	TrackQualityShift   8
#define	TrackQualityMask   255
#define	TrackFitGlobal   1
#define	TrackFitCosmic   2
#define	TrackFitCdcKalman   4
#define	TrackFitSvdCdcKalman   8
#define	TrackFitShift   16
#define	TrackFitMask   255
#define	TrackHasMother   1
#define	TrackHasDaughter   2
#define	TrackRelationShift   24
#define	TrackRelationMask   255
#define	HEP_SHORT_NAMES
#define	OLD_STEREO   1
	calculates arc length and z for a stereo hit. uses these functions for curl tracks.
#define	TTrack_INLINE_DEFINE_HERE

Typedefs
typedef HepGeom::Point3D< double >	HepPoint3D

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

Macro Definition Documentation

HEP_SHORT_NAMES

#define HEP_SHORT_NAMES

Definition at line 67 of file TTrack.h.

OLD_STEREO

#define OLD_STEREO   1

calculates arc length and z for a stereo hit. uses these functions for curl tracks.

Definition at line 250 of file TTrack.h.

TrackCurlFinder

#define TrackCurlFinder   8

Definition at line 26 of file TTrack.h.

TrackFastFinder

#define TrackFastFinder   2

Definition at line 24 of file TTrack.h.

TrackFinderMask

#define TrackFinderMask   255

Definition at line 30 of file TTrack.h.

TrackFitCdcKalman

#define TrackFitCdcKalman   4

Definition at line 54 of file TTrack.h.

TrackFitCosmic

#define TrackFitCosmic   2

Definition at line 53 of file TTrack.h.

TrackFitGlobal

#define TrackFitGlobal   1

Definition at line 52 of file TTrack.h.

TrackFitMask

#define TrackFitMask   255

Definition at line 57 of file TTrack.h.

TrackFitShift

#define TrackFitShift   16

Definition at line 56 of file TTrack.h.

TrackFitSvdCdcKalman

#define TrackFitSvdCdcKalman   8

Definition at line 55 of file TTrack.h.

TrackHasDaughter

#define TrackHasDaughter   2

Definition at line 61 of file TTrack.h.

TrackHasMother

#define TrackHasMother   1

Definition at line 60 of file TTrack.h.

TrackOldConformalFinder

#define TrackOldConformalFinder   1

Definition at line 23 of file TTrack.h.

TrackPMCurlFinder

#define TrackPMCurlFinder   32

Definition at line 28 of file TTrack.h.

TrackQuality2D

#define TrackQuality2D   8

Definition at line 47 of file TTrack.h.

TrackQualityAfterKink

#define TrackQualityAfterKink   2

Definition at line 45 of file TTrack.h.

TrackQualityCosmic

#define TrackQualityCosmic   4

Definition at line 46 of file TTrack.h.

TrackQualityMask

#define TrackQualityMask   255

Definition at line 49 of file TTrack.h.

TrackQualityOutsideCurler

#define TrackQualityOutsideCurler   1

Definition at line 44 of file TTrack.h.

TrackQualityShift

#define TrackQualityShift   8

Definition at line 48 of file TTrack.h.

TrackRelationMask

#define TrackRelationMask   255

Definition at line 63 of file TTrack.h.

TrackRelationShift

#define TrackRelationShift   24

Definition at line 62 of file TTrack.h.

TrackSlowFinder

#define TrackSlowFinder   4

Definition at line 25 of file TTrack.h.

TrackSVDAssociator

#define TrackSVDAssociator   64

Definition at line 29 of file TTrack.h.

TrackTrackManager

#define TrackTrackManager   16

Definition at line 27 of file TTrack.h.

TrackTypeCircle

#define TrackTypeCircle   4

Definition at line 36 of file TTrack.h.

TrackTypeCosmic

#define TrackTypeCosmic   8

Definition at line 37 of file TTrack.h.

TrackTypeCurl

#define TrackTypeCurl   2

Definition at line 35 of file TTrack.h.

TrackTypeIncomingCosmic

#define TrackTypeIncomingCosmic   8

Definition at line 38 of file TTrack.h.

TrackTypeKink

#define TrackTypeKink   32

Definition at line 40 of file TTrack.h.

TrackTypeNormal

#define TrackTypeNormal   1

Definition at line 34 of file TTrack.h.

TrackTypeOutgoingCosmic

#define TrackTypeOutgoingCosmic   16

Definition at line 39 of file TTrack.h.

TrackTypeSVDOnly

#define TrackTypeSVDOnly   256

Definition at line 41 of file TTrack.h.

TrackTypeUndefined

#define TrackTypeUndefined   0

Definition at line 33 of file TTrack.h.

TTrack_INLINE_DEFINE_HERE

#define TTrack_INLINE_DEFINE_HERE

Definition at line 470 of file TTrack.h.

Typedef Documentation

HepPoint3D

typedef HepGeom::Point3D<double> HepPoint3D

Definition at line 103 of file TTrack.h.

Function Documentation

HelixHasNan()

bool HelixHasNan

(

const Helix &

h

)

Helix parameter validity.

Definition at line 3934 of file TTrack.cxx.

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

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

PositiveDefinite()

bool PositiveDefinite

(

const Helix &

h

)

Error matrix validity.

Definition at line 3918 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().

SortByPt()

int SortByPt	(	const void *	a,
		const void *	b
	)

Utility functions.

Definition at line 2530 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 3948 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 3684 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 3656 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 3712 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);
}

TrackDump()

std::string TrackDump ( const TTrack &  t )

inline

to dump a track.

Definition at line 697 of file TTrack.h.

                            {
    return t.name() + " " + TrackStatus(t) + " " + TrackKinematics(t) + " " +
    TrackInformation(t);
}

Referenced by TTrackManager::dump(), and TTrackManager::saveTables().

TrackInformation() [1/3]

std::string TrackInformation

(

const MdcRec_trk &

r

)

Definition at line 3868 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 3856 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(), and TrackInformation().

TrackInformation() [3/3]

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

Definition at line 3880 of file TTrack.cxx.

                                                                    {
    std::string s;
 
    s += "a" + itostring(int(nA));
    s += " s" + itostring(int(nS));
    s += " n" + itostring(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 += itostring(y) + ".";
    float z = fabs(x);
    for (unsigned j = 0; j < 1; j++) {
    z *= 10.;
    y = (int(z) % 10);
    s += itostring(y);
    }
 
    return s;
}

TrackKinematics() [1/2]

std::string TrackKinematics

(

const Helix &

h

)

Definition at line 3740 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 += itostring(y) + ".";
    float z = fabs(x[i]);
    for (unsigned j = 0; j < 3; j++) {
        z *= 10.;
        y = (int(z) % 10);
        s += itostring(y);
    }
    }
 
    return s;
}

TrackKinematics() [2/2]

std::string TrackKinematics ( const TTrack &  t )

inline

Definition at line 691 of file TTrack.h.

                                  {
    return TrackKinematics(t.helix());
}

Referenced by TrackDump(), and TrackKinematics().

TrackLayerUsage()

std::string TrackLayerUsage

(

const TTrack &

t

)

Definition at line 3905 of file TTrack.cxx.

                                  {
    unsigned n[11];
    NHitsSuperLayer(t.links(), n);
    std::string nh;
    for (unsigned i = 0; i < 11; i++) {
    nh += itostring(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 3789 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);
}

TrackStatus() [2/4]

std::string TrackStatus

(

const MdcRec_trk_add &

a

)

Definition at line 3796 of file TTrack.cxx.

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

TrackStatus() [3/4]

std::string TrackStatus

(

const TTrack &

t

)

returns string of track status.

Definition at line 3779 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 3806 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 + itostring(m) + p + itostring(d);
}

TrackType() [1/2]

std::string TrackType ( const TTrack &  t )

inline

returns string of track type.

Definition at line 685 of file TTrack.h.

                            {
    return TrackType(t.type());
}

TrackType() [2/2]

std::string TrackType ( unsigned  type )

inline

Definition at line 3030 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  ");
}