Track Struct Reference

#include <Track.h>

Inheritance diagram for Track:

Collaboration diagram for Track:

Public Member Functions
void	addGrain (Grain gr)
bool	evaluateEstimators ()
unsigned int	getAreaSum () const
double	getBottomZ () const
unsigned short	getCount () const
double	getErrX () const
double	getErrY () const
double	getErrZ () const
std::vector< Grain >	getGrains () const
double	getSigma () const
double	getSX () const
double	getSY () const
double	getSZ () const
double	getTopZ () const
double	getX () const
double	getY () const
double	getZ () const
void	setAreaSum (unsigned int areaSum)
void	setBottomZ (double z)
void	setCount (unsigned short count)
void	setErrX (double errx)
void	setErrY (double erry)
void	setErrZ (double errz)
void	setSigma (double sigma)
void	setSX (double sx)
void	setSY (double sy)
void	setSZ (double sz)
void	setTopZ (double z)
void	setX (double x)
void	setY (double y)
void	setZ (double z)
	Track ()
	~Track ()

Public Attributes
double	dz
int	Field
float	FirstZ
TVector	Intercept
TVector	InterceptErrors
float	LastZ
double	meanDeltaTheta3D
double	meanDeltaThetaXZ
double	meanDeltaThetaYZ
double	meanDistanceClustersTo3DLine
double	meanGapClusterToCluster
double	mx
double	mxSigma
double	my
double	mySigma
BYTE *	pCorrection
char *	pCorrection
unsigned	PointsN
TVector *	pPoints
double	qx
double	qy
double	rx
double	ry
float	Sigma
double	sigmaDeltaTheta3D
double	sigmaDeltaThetaXZ
double	sigmaDeltaThetaYZ
TVector	Slope
TVector	SlopeErrors
boolean	Valid
bool	Valid
double	xChi2
double	yChi2

Private Member Functions
double	pointTo3DlineDistance (double x0, double y0, double z0)

Private Attributes
unsigned int	_areaSum
double	_bottomZ
unsigned short	_count
double	_errx
double	_erry
double	_errz
std::vector< Grain >	_grains
double	_sigma
double	_sx
double	_sy
double	_sz
double	_topZ
double	_x
double	_y
double	_z

Constructor & Destructor Documentation

Track()

Track::Track

(

)

{
    meanDeltaTheta3D = sigmaDeltaTheta3D =  meanDeltaThetaXZ = sigmaDeltaThetaXZ = meanDeltaThetaYZ = sigmaDeltaThetaYZ = 0;
    meanDistanceClustersTo3DLine = meanGapClusterToCluster = 0;
    mx =  qx = rx = xChi2 = 0;
    my =  qy = ry = yChi2 = 0;
    dz = 0;
    // to be set as input...
    _errx = _erry = 0.3;
    _errz = 3.;
}

~Track()

Track::~Track ( )

inline

{};

Member Function Documentation

addGrain()

void Track::addGrain ( Grain  gr )

inline

{_grains.push_back(gr);}

evaluateEstimators()

bool Track::evaluateEstimators

(

)

{
    //std::ofstream out("test.txt",std::ios::app); // used now for debugging --> later to be removed
    
    if (_grains.size() == 0)  
        return false;
    
    unsigned short nGrains = _count;
    //double s = 0, sx = 0, sy = 0, sz = 0, sxx = 0, syy = 0, szz = 0, szx = 0, szy = 0;
    //double d = _errz/_errx;
 
    Fitter fitXZ; fitXZ.setErrXY(3.,0.33);  //error on z coordinate is ~3 um, in x or y is 0.33 (pixel dimension)
    Fitter fitYZ; fitYZ.setErrXY(3.,0.33);
    for (unsigned short i = 0; i < nGrains; i++)
    {
        double x = _grains.at(i).x;
        double y = _grains.at(i).y;
        double z = _grains.at(i).z;
 
        fitXZ.addPoint(z,x);
        fitYZ.addPoint(z,y);
        
        meanDistanceClustersTo3DLine+= pointTo3DlineDistance(x,y,z);
 
        if (i < nGrains -1)
        {
            double dx = _grains.at(i+1).x-x;
            double dy = _grains.at(i+1).y-y;
            double dz = _grains.at(i+1).z-z;
            meanGapClusterToCluster += sqrt(dx*dx+dy*dy+dz*dz);
        }
 
        if (i < nGrains-2)
        {
            double dx1 = _grains.at(i+1).x-x;
            double dy1 = _grains.at(i+1).y-y;
            double dz1 = _grains.at(i+1).z-z;
            double dx2 = _grains.at(i+2).x-_grains.at(i+1).x;
            double dy2 = _grains.at(i+2).y-_grains.at(i+1).y;
            double dz2 = _grains.at(i+2).z-_grains.at(i+1).z;
 
            double sx1  = dx1/dz1; // tan_x
            double sy1  = dy1/dz1; // tan_y
 
            double sx2  = dx2/dz2; // tan_x
            double sy2  = dy2/dz2; // tan_y
            
            double valueXZ = acos((dx1*dx2+dz1*dz2)/(sqrt(dx1*dx1+dz1*dz1)*sqrt(dx2*dx2+dz2*dz2)));
            double valueYZ = acos((dy1*dy2+dz1*dz2)/(sqrt(dy1*dy1+dz1*dz1)*sqrt(dy2*dy2+dz2*dz2)));
            double value3D = acos((dx1*dx2+dy1*dy2+dz1*dz2)/(sqrt(dx1*dx1+dy1*dy1+dz1*dz1)*sqrt(dx2*dx2+dy2*dy2+dz2*dz2)));
            //double value3D = acos((sx1*sx2+sy1*sy2+1)/(sqrt(sx1*sx1+sy1*sy1+1)*sqrt(sx2*sx2+sy2*sy2+1)));
 
            double signXZ = 1.;
            double signYZ = 1.;
 
            if (sx2<sx1)
                signXZ*=-1;
            if (sy2<sy1)
                signYZ*=-1;
 
            meanDeltaThetaXZ+= valueXZ * signXZ;
            meanDeltaThetaYZ+= valueYZ * signYZ;
            meanDeltaTheta3D+= value3D;
 
            sigmaDeltaThetaXZ+=valueXZ*valueXZ;
            sigmaDeltaThetaYZ+=valueYZ*valueYZ;
            sigmaDeltaTheta3D+=value3D*value3D;
            
        }
    }
    meanDeltaThetaXZ/=(nGrains-2.);
    meanDeltaThetaYZ/=(nGrains-2.);
    meanDeltaTheta3D/=(nGrains-2.);
 
    meanGapClusterToCluster/=(nGrains-1.);
    
    sigmaDeltaThetaXZ/=(nGrains-2.);
    sigmaDeltaThetaYZ/=(nGrains-2.);
    sigmaDeltaTheta3D/=(nGrains-2.);
 
    sigmaDeltaThetaXZ-=pow(meanDeltaThetaXZ,2);
    sigmaDeltaThetaYZ-=pow(meanDeltaThetaYZ,2);
    sigmaDeltaTheta3D-=pow(meanDeltaTheta3D,2);
 
    sigmaDeltaThetaXZ = sqrt(sigmaDeltaThetaXZ);
    sigmaDeltaThetaYZ = sqrt(sigmaDeltaThetaYZ);
    sigmaDeltaTheta3D = sqrt(sigmaDeltaTheta3D);
 
    
    
    // fit line grains on plane x-z
    fitXZ.evaluateLinearFit();
    mx = fitXZ.getSlope();
    qx = fitXZ.getIntercept();
    rx = fitXZ.getPearsonCoefficent();
    
    
    // fit line grains on plane y-z
    fitYZ.evaluateLinearFit();
    my = fitYZ.getSlope();
    qy = fitYZ.getIntercept();
    ry = fitYZ.getPearsonCoefficent();
    
    xChi2 = fitXZ.getChi2();
    yChi2 = fitYZ.getChi2();
 
    // lenght of the segment (in z).
    dz = _grains.at(0).z-_grains.at(nGrains-1).z;
    meanDistanceClustersTo3DLine/=nGrains;
 
    return true;
}

getAreaSum()

unsigned int Track::getAreaSum ( ) const

inline

{return _areaSum;}

getBottomZ()

double Track::getBottomZ ( ) const

inline

{return _bottomZ;}

getCount()

unsigned short Track::getCount ( ) const

inline

{return _count;}

getErrX()

double Track::getErrX ( ) const

inline

{return _errx;}

getErrY()

double Track::getErrY ( ) const

inline

{return _erry;}

getErrZ()

double Track::getErrZ ( ) const

inline

{return _errz;}

getGrains()

std::vector< Grain > Track::getGrains ( ) const

inline

{return _grains;}

getSigma()

double Track::getSigma ( ) const

inline

{return _sigma;}

getSX()

double Track::getSX ( ) const

inline

{return _sx;}

getSY()

double Track::getSY ( ) const

inline

{return _sy;}

getSZ()

double Track::getSZ ( ) const

inline

{return _sz;}

getTopZ()

double Track::getTopZ ( ) const

inline

{return _topZ;}

getX()

double Track::getX ( ) const

inline

{return _x;}

getY()

double Track::getY ( ) const

inline

{return _y;}

getZ()

double Track::getZ ( ) const

inline

{return _z;}

pointTo3DlineDistance()

double Track::pointTo3DlineDistance ( double  x0, double  y0, double  z0  )

private

{
    double dist = 0;
    double qp_x = _x-x0;
    double qp_y = _y-y0;
    double qp_z = _z-z0;
    double vx = _sx;
    double vy = _sy;
    double vz = 1;
    double qr = (qp_x*vx+qp_y*vy+qp_z*vz)/sqrt(vx*vx+vy*vy+vz*vz);
    double qp =sqrt(qp_x*qp_x+qp_y*qp_y+qp_z*qp_z);
    dist = sqrt(qp*qp-qr*qr);
    //std::cout << qr << "\t" << qp << "\t" << dist << std::endl;
    return dist;
}

setAreaSum()

void Track::setAreaSum ( unsigned int  areaSum )

inline

{_areaSum = areaSum;}

setBottomZ()

void Track::setBottomZ ( double  z )

inline

{_bottomZ = z;}

setCount()

void Track::setCount ( unsigned short  count )

inline

{_count = count;}

setErrX()

void Track::setErrX ( double  errx )

inline

{_errx = errx;}

setErrY()

void Track::setErrY ( double  erry )

inline

{_erry = erry;}

setErrZ()

void Track::setErrZ ( double  errz )

inline

{_errz = errz;}

setSigma()

void Track::setSigma ( double  sigma )

inline

{_sigma = sigma;}

setSX()

void Track::setSX ( double  sx )

inline

{_sx = sx;}

setSY()

void Track::setSY ( double  sy )

inline

{_sy = sy;}

setSZ()

void Track::setSZ ( double  sz )

inline

{_sz = sz;}

setTopZ()

void Track::setTopZ ( double  z )

inline

{_topZ = z;}

setX()

void Track::setX ( double  x )

inline

{_x = x;}

setY()

void Track::setY ( double  y )

inline

{_y = y;}

setZ()

void Track::setZ ( double  z )

inline

{_z = z;}

Member Data Documentation

_areaSum

unsigned int Track::_areaSum

private

_bottomZ

double Track::_bottomZ

private

_count

unsigned short Track::_count

private

_errx

double Track::_errx

private

_erry

double Track::_erry

private

_errz

double Track::_errz

private

_grains

std::vector<Grain> Track::_grains

private

_sigma

double Track::_sigma

private

_sx

double Track::_sx

private

_sy

double Track::_sy

private

_sz

double Track::_sz

private

_topZ

double Track::_topZ

private

_x

double Track::_x

private

_y

double Track::_y

private

_z

double Track::_z

private

dz

double Track::dz

Field

int Track::Field

FirstZ

float Track::FirstZ

Intercept

TVector Track::Intercept

InterceptErrors

TVector Track::InterceptErrors

LastZ

float Track::LastZ

meanDeltaTheta3D

double Track::meanDeltaTheta3D

meanDeltaThetaXZ

double Track::meanDeltaThetaXZ

meanDeltaThetaYZ

double Track::meanDeltaThetaYZ

meanDistanceClustersTo3DLine

double Track::meanDistanceClustersTo3DLine

meanGapClusterToCluster

double Track::meanGapClusterToCluster

mx

double Track::mx

mxSigma

double Track::mxSigma

my

double Track::my

mySigma

double Track::mySigma

pCorrection [1/2]

BYTE* Track::pCorrection

pCorrection [2/2]

char* Track::pCorrection

PointsN

unsigned Track::PointsN

pPoints

TVector * Track::pPoints

qx

double Track::qx

qy

double Track::qy

rx

double Track::rx

ry

double Track::ry

Sigma

float Track::Sigma

sigmaDeltaTheta3D

double Track::sigmaDeltaTheta3D

sigmaDeltaThetaXZ

double Track::sigmaDeltaThetaXZ

sigmaDeltaThetaYZ

double Track::sigmaDeltaThetaYZ

Slope

TVector Track::Slope

SlopeErrors

TVector Track::SlopeErrors

Valid [1/2]

boolean Track::Valid

Valid [2/2]

bool Track::Valid

xChi2

double Track::xChi2

yChi2

double Track::yChi2

---

The documentation for this struct was generated from the following files:

• /home/antonio/fedra_doxygen/src/appl/rwcToEdb/Track.h
• /home/antonio/fedra_doxygen/src/libDataConversion/dataIO/Track.h
• /home/antonio/fedra_doxygen/src/libDataConversion/dataIO/Track.linux.h
• /home/antonio/fedra_doxygen/src/appl/rwcToEdb/Track.cpp