TTrack.cxx File Reference

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

Go to the source code of this file.

Macros
#define	OPTNK
#define	HEP_SHORT_NAMES
#define	NEW_FIT2D   1
#define	JJTEST   0
#define	JJJTEST   0
#define	t_dot(a, b)   (a[0]*b[0]+a[1]*b[1]+a[2]*b[2])
#define	t_dot2(a, b)   (a[0]*b[0]+a[1]*b[1])
#define	t_print(a, b)   std::cout << b << " = " << a[0] << " " << a[1] << " " << a[2] << std::endl;

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

HEP_SHORT_NAMES

#define HEP_SHORT_NAMES

Definition at line 33 of file TTrack.cxx.

JJJTEST

#define JJJTEST   0

JJTEST

#define JJTEST   0

NEW_FIT2D

#define NEW_FIT2D   1

Definition at line 940 of file TTrack.cxx.

OPTNK

#define OPTNK

Definition at line 1 of file TTrack.cxx.

t_dot

#define t_dot	(		a,
			b
	)		(a[0]*b[0]+a[1]*b[1]+a[2]*b[2])

Definition at line 3047 of file TTrack.cxx.

t_dot2

#define t_dot2	(		a,
			b
	)		(a[0]*b[0]+a[1]*b[1])

Definition at line 3048 of file TTrack.cxx.

t_print

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

Definition at line 3049 of file TTrack.cxx.

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 *	av,
		const void *	bv
	)

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);
}

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()

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;
}

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()

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  ");
}