reco::TrackResiduals Class Reference

#include <TrackResiduals.h>

Public Types
enum	ResidualType { X_Y_RESIDUALS, X_Y_PULLS }

Public Member Functions
void	print (std::ostream &stream=std::cout) const
void	print (const HitPattern &, std::ostream &stream=std::cout) const
double	residualX (int i, const HitPattern &) const
double	residualX (int i) const
double	residualY (int i, const HitPattern &) const
double	residualY (int i) const
void	setPullXY (int idx, double pullX, double pullY)
void	setResidualType (enum ResidualType)
void	setResidualXY (int idx, double residualX, double residualY)
	TrackResiduals ()
	TrackResiduals (enum ResidualType)

Protected Types
enum	{ numResiduals = 0x40 }
	number of residuals stored More...

Static Protected Member Functions
static unsigned char	pack_pull (double)
static unsigned char	pack_residual (double)
static double	unpack_pull (unsigned char)
static double	unpack_residual (unsigned char)

Private Attributes
unsigned char	residuals_ [numResiduals]
	residuals, bitpacked two hits to a char More...
char	residualType

Friends
class	Trajectory

Detailed Description

Definition at line 15 of file TrackResiduals.h.

Member Enumeration Documentation

anonymous enum

protected

number of residuals stored

Enumerator
numResiduals

Definition at line 45 of file TrackResiduals.h.

{ numResiduals = 0x40 };

enum reco::TrackResiduals::ResidualType

Enumerator
X_Y_RESIDUALS
X_Y_PULLS

Definition at line 22 of file TrackResiduals.h.

                      {
        X_Y_RESIDUALS,
        X_Y_PULLS
    };

Constructor & Destructor Documentation

TrackResiduals::TrackResiduals

(

)

Definition at line 11 of file TrackResiduals.cc.

References residuals_.

                                :
    residualType(X_Y_RESIDUALS)
{
    memset(residuals_, 0, sizeof(residuals_));
}

TrackResiduals::TrackResiduals

(

enum ResidualType

type

)

Definition at line 17 of file TrackResiduals.cc.

References residuals_.

                                                      :
    residualType(type)
{
    memset(residuals_, 0, sizeof(residuals_));
}

Member Function Documentation

unsigned char TrackResiduals::pack_pull ( double  pull )

staticprotected

Definition at line 132 of file TrackResiduals.cc.

References mag(), pull_char_to_double, and FWPFMaths::sgn().

Referenced by setPullXY().

{
     unsigned char sgn = (pull < 0) * 0x08; // 1xxx is -abs(0xxx)
     int mag = -1;
     while (++mag < 8 && pull_char_to_double[mag][1] < fabs(pull));
     return sgn + mag;
}

unsigned char TrackResiduals::pack_residual ( double  pull )

staticprotected

Definition at line 148 of file TrackResiduals.cc.

References mag(), pull_char_to_double, and FWPFMaths::sgn().

Referenced by setResidualXY().

{
    unsigned char sgn = (pull < 0) * 0x08; // 1xxx is -abs(0xxx)
    int mag = -1;
    while (++mag < 8 && pull_char_to_double[mag][1] < fabs(pull)) {
        // just counting
        ;
    }
    return sgn + mag;
}

void TrackResiduals::print

(

std::ostream &

stream = std::cout

)

const

Definition at line 159 of file TrackResiduals.cc.

References flags, i, numResiduals, and residuals_.

{
    stream << "TrackResiduals" << std::endl;
    std::ios_base::fmtflags flags = stream.flags();
    stream.setf (std::ios_base::hex, std::ios_base::basefield);
    stream.setf (std::ios_base::showbase);

    for (int i = 0; i < numResiduals; i++) {
        unsigned char residual = residuals_[i];
        printf("0x%x\n", residual);
        //stream << residual << std::endl;
    }

    stream.flags(flags);
}

void TrackResiduals::print	(	const HitPattern &	h,
		std::ostream &	stream = std::cout
	)		const

Definition at line 175 of file TrackResiduals.cc.

References reco::HitPattern::getHitPattern(), i, reco::HitPattern::numberOfHits(), residualX(), residualY(), reco::HitPattern::TRACK_HITS, and reco::HitPattern::validHitFilter().

{
    stream << "TrackResiduals" << std::endl;
    for (int i = 0; i < h.numberOfHits(HitPattern::TRACK_HITS); i++) {
        stream << (h.validHitFilter(h.getHitPattern(HitPattern::TRACK_HITS, i)) ? "valid hit:   " : "invalid hit: ")
               << "( " << residualX(i, h) << " , " << residualY(i, h) << " )"
               << std::endl;
    }
}

double TrackResiduals::residualX	(	int	i,
		const HitPattern &	h
	)		const

get the residual of the ith hit (needs the hit pattern to figure out which hits are valid)

Definition at line 83 of file TrackResiduals.cc.

References customizeTrackingMonitorSeedNumber::idx, and index_to_hitpattern().

Referenced by FWTrackResidualDetailView::prepareData(), and print().

{
    int idx = index_to_hitpattern(i, h);
    if (idx == -999) {
        return -999;
    }

    return residualX(idx);
}

double TrackResiduals::residualX

(

int

i

)

const

get the residual of the ith valid hit, with no regard for alignment with the HitPattern

Definition at line 38 of file TrackResiduals.cc.

References assert(), residuals_, residualType, unpack_pull(), unpack_residual(), X_Y_PULLS, and X_Y_RESIDUALS.

{
    switch (residualType) {
    case X_Y_RESIDUALS:
        return unpack_residual(residuals_[i] >> 4);
    case X_Y_PULLS:
        return unpack_pull(residuals_[i] >> 4);
    default:
        assert(0);
    }
    return 0;
}

double TrackResiduals::residualY	(	int	i,
		const HitPattern &	h
	)		const

Definition at line 93 of file TrackResiduals.cc.

References customizeTrackingMonitorSeedNumber::idx, and index_to_hitpattern().

Referenced by FWTrackResidualDetailView::prepareData(), and print().

{
    int idx = index_to_hitpattern(i, h);
    if (idx == -999) {
        return -999;
    }

    return residualY(idx);
}

double TrackResiduals::residualY

(

int

i

)

const

Definition at line 51 of file TrackResiduals.cc.

References assert(), residuals_, residualType, unpack_pull(), unpack_residual(), X_Y_PULLS, and X_Y_RESIDUALS.

{
    switch (residualType) {
    case X_Y_RESIDUALS:
        return unpack_residual(residuals_[i] & 0x0f);
    case X_Y_PULLS:
        return unpack_pull(residuals_[i] & 0x0f);
    default:
        assert(0);
    }
    return 0;
}

void TrackResiduals::setPullXY	(	int	idx,
		double	pullX,
		double	pullY
	)

Definition at line 103 of file TrackResiduals.cc.

References assert(), customizeTrackingMonitorSeedNumber::idx, numResiduals, pack_pull(), residuals_, residualType, and X_Y_PULLS.

Referenced by trajectoryToResiduals().

{
    assert(residualType == X_Y_PULLS);
    if (idx >= numResiduals) {
        edm::LogWarning("TrackResiduals") << " setting pulls over the array size.";
        return;
    }

    residuals_[idx] = (pack_pull(pullX) << 4) | pack_pull(pullY);
}

void TrackResiduals::setResidualType

(

enum ResidualType

type

)

Definition at line 23 of file TrackResiduals.cc.

References residualType.

{
    residualType = type;
}

void TrackResiduals::setResidualXY	(	int	idx,
		double	residualX,
		double	residualY
	)

Definition at line 28 of file TrackResiduals.cc.

References assert(), customizeTrackingMonitorSeedNumber::idx, numResiduals, pack_residual(), residuals_, residualType, and X_Y_RESIDUALS.

Referenced by trajectoryToResiduals().

{
    assert(residualType == X_Y_RESIDUALS);
    if (idx >= numResiduals) {
        edm::LogWarning("TrackResiduals") << " setting residual over the array size.";
        return;
    }
    residuals_[idx] = (pack_residual(residualX) << 4) | pack_residual(residualY);
}

double TrackResiduals::unpack_pull ( unsigned char  pull )

staticprotected

Definition at line 125 of file TrackResiduals.cc.

References mag(), pull_char_to_double, and FWPFMaths::sgn().

Referenced by residualX(), and residualY().

{
    int sgn = 1 - 2 * ((pull & 0x08) >> 3);
    unsigned char mag = pull & 0x07;
    return sgn * (pull_char_to_double[mag][0] + pull_char_to_double[mag][1]) / 2;
}

double TrackResiduals::unpack_residual ( unsigned char  pull )

staticprotected

Definition at line 141 of file TrackResiduals.cc.

References mag(), pull_char_to_double, and FWPFMaths::sgn().

Referenced by residualX(), and residualY().

{
    int sgn = 1 - 2 * ((pull & 0x08) >> 3);
    unsigned char mag = pull & 0x07;
    return sgn * (pull_char_to_double[mag][0] + pull_char_to_double[mag][1]) / 2;
}

Friends And Related Function Documentation

friend class Trajectory

friend

Definition at line 18 of file TrackResiduals.h.

Member Data Documentation

unsigned char reco::TrackResiduals::residuals_[numResiduals]

private

residuals, bitpacked two hits to a char

Definition at line 53 of file TrackResiduals.h.

Referenced by print(), residualX(), residualY(), setPullXY(), setResidualXY(), and TrackResiduals().

char reco::TrackResiduals::residualType

private

Definition at line 54 of file TrackResiduals.h.

Referenced by residualX(), residualY(), setPullXY(), setResidualType(), and setResidualXY().