#include "KFitTrack.h"

Inheritance diagram for trkf::KFitTrack:

Public Types
enum	FitStatus { INVALID, UNKNOWN, FORWARD, FORWARD_PREDICTED, BACKWARD, BACKWARD_PREDICTED, OPTIMAL, OPTIMAL_PREDICTED }
	Fit status enum. More...

Public Member Functions
	KFitTrack ()
	Default constructor. More...

	KFitTrack (const KETrack &tre, double s=0., double chisq=0., FitStatus stat=INVALID)
	Initializing constructor. More...

virtual	~KFitTrack ()
	Destructor. More...

double	getPath () const
	Propagation distance. More...

double	getChisq () const
	Fit chisquare. More...

FitStatus	getStat () const
	Fit status. More...

void	setPath (double path)
	Set propagation distance. More...

void	setChisq (double chisq)
	Set chisquare. More...

void	setStat (FitStatus stat)
	Set fit status. More...

bool	combineFit (const KFitTrack &trf)
	Combine two tracks. More...

virtual std::ostream &	Print (std::ostream &out, bool doTitle=true) const
	Printout. More...

const TrackError &	getError () const
	Track error matrix. More...

TrackError &	getError ()
	Modifiable error matrix. More...

double	PointingError () const
	Pointing error (radians). More...

void	setError (const TrackError &err)
	Set error matrix. More...

boost::optional< double >	combineTrack (const KETrack &tre)
	Combine two tracks. More...

const std::shared_ptr< const Surface > &	getSurface () const
	Surface. More...

const TrackVector &	getVector () const
	Track state vector. More...

TrackVector &	getVector ()
	Modifiable state vector. More...

Surface::TrackDirection	getDirection () const
	Track direction. More...

int	PdgCode () const
	Pdg code. More...

double	Mass () const
	Based on pdg code. More...

void	setSurface (const std::shared_ptr< const Surface > &psurf)
	Set surface. More...

void	setVector (const TrackVector &vec)
	Set state vector. More...

void	setDirection (Surface::TrackDirection dir)
	Set direction. More...

void	setPdgCode (int pdg)
	Set pdg code. More...

bool	isValid () const
	Test if track is valid. More...

void	getPosition (double xyz[3]) const
	Get position of track. More...

void	getMomentum (double mom[3]) const
	Get momentum vector of track. More...

double	XLatitude () const
	Get x-latitude. More...

double	XLongitude () const
	Get x-longitude. More...

Private Attributes
double	fPath
	Propagation distance. More...

double	fChisq
	Fit chisquare. More...

FitStatus	fStat
	Fit status. More...

Detailed Description

Definition at line 36 of file KFitTrack.h.

Member Enumeration Documentation

enum trkf::KFitTrack::FitStatus

Fit status enum.

Enumerator
INVALID
UNKNOWN
FORWARD
FORWARD_PREDICTED
BACKWARD
BACKWARD_PREDICTED
OPTIMAL
OPTIMAL_PREDICTED

Definition at line 41 of file KFitTrack.h.

                    {
       INVALID,             // No valid fit information.
       UNKNOWN,             // Unknown.
       FORWARD,             // Fit based on past measurements, including current surface.
       FORWARD_PREDICTED,   // Fit based on past measurements, not including current surface.
       BACKWARD,            // Fit based on future measurements, including current surface.
       BACKWARD_PREDICTED,  // Fit based on future measurements, not including current surface.
       OPTIMAL,             // Fit based on all measurements, including current surface.
       OPTIMAL_PREDICTED    // Fit based on all measurements, except current surface.
     };

Constructor & Destructor Documentation

trkf::KFitTrack::KFitTrack ( )

Default constructor.

Definition at line 17 of file KFitTrack.cxx.

                        :
     fPath(0.),
     fChisq(0.),
     fStat(INVALID)
   {}

trkf::KFitTrack::KFitTrack	(	const KETrack &	tre,
		double	s = `0.`,
		double	chisq = `0.`,
		FitStatus	stat = `INVALID`
	)

Initializing constructor.

Arguments:

tre - KETrack. s - Path distance. chisq - Fit chisquare. stat - Fit status.

Definition at line 32 of file KFitTrack.cxx.

                                        :
     KETrack(tre),
     fPath(s),
     fChisq(chisq),
     fStat(stat)
   {}

trkf::KFitTrack::~KFitTrack ( )

virtual

Destructor.

Definition at line 43 of file KFitTrack.cxx.

44 {}

Member Function Documentation

bool trkf::KFitTrack::combineFit ( const KFitTrack & trf )

Combine two tracks.

Arguments:

trf - Another track.

Returns: True if success.

This method updates the current track to be the weighted average of itself and another track. Track parameters and error matrix are updated by method KETrack::combineTrack. The updated chissquare is the sum of the original chisquare, the chisquare from the other track, and the chisquare of the combination. The path distance is not updated.

The updated status is derived from the input status according to the following table.

FORWARD + BACKWARD_PREDICTED = OPTIMAL FORWARD_PREDICTED + BACKWARD = OPTIMAL FORWARD_PREDICTED + BACKWARD_PREDICTED = OPTIMAL_PREDICTED

Any other combination of input statuses will cause this method to return failure. This method can also return failure because of singular error matrices. In case of failure, the original track is not updated.

The two tracks being combined should be on the same surface (throw exception if not).

Definition at line 76 of file KFitTrack.cxx.

References BACKWARD, BACKWARD_PREDICTED, trkf::KETrack::combineTrack(), fChisq, FORWARD, FORWARD_PREDICTED, fStat, getChisq(), getStat(), trkf::KTrack::getSurface(), OPTIMAL, OPTIMAL_PREDICTED, and UNKNOWN.

Referenced by setStat(), and trkf::KalmanFilterAlg::smoothTrack().

   {
     // Make sure that the two track surfaces are the same.
     // Throw an exception if they are not.
 
     if(!getSurface()->isEqual(*trf.getSurface()))
       throw cet::exception("KFitTrack") << "Track combination surfaces are not the same.\n";
 
     // Default result failure.
 
     bool result = false;
 
     // Make sure input statuses are compatible and figure out the
     // updated status.
 
     FitStatus stat1 = getStat();
     FitStatus stat2 = trf.getStat();
     FitStatus statu = UNKNOWN;
     if((stat1 == FORWARD && stat2 == BACKWARD_PREDICTED) ||
        (stat1 == FORWARD_PREDICTED && stat2 == BACKWARD) ||
        (stat1 == BACKWARD && stat2 == FORWARD_PREDICTED) ||
        (stat1 == BACKWARD_PREDICTED && stat2 == FORWARD))
       statu = OPTIMAL;
     else if(stat1 == FORWARD_PREDICTED && stat2 == BACKWARD_PREDICTED)
       statu = OPTIMAL_PREDICTED;
     else
       statu = UNKNOWN;
     if(statu != UNKNOWN) {
 
       // Update track parameters and error matrix.
 
       boost::optional<double> chisq = combineTrack(trf);
 
       // Update status and chisquare.
 
       if(!!chisq) {
         result = true;
         fStat = statu;
         fChisq = fChisq + trf.getChisq() + *chisq;
       }
     }
 
     // Done.
 
     return result;
   }

boost::optional< double > trkf::KETrack::combineTrack ( const KETrack & tre )

inherited

Combine two tracks.

Arguments:

tre - Another track.

Returns: Chisquare + success flag.

This method updates the current track to be the weighted average of itself and another track. The chisquare of the combination is returned as the result value. The combination can fail because the sum of the two error matrices is singular, in which case the success flag embedded in the return value is false.

Definition at line 95 of file KETrack.cxx.

References trkf::KETrack::getError(), trkf::KTrack::getSurface(), trkf::KTrack::getVector(), trkf::KTrack::isValid(), trkf::KETrack::setError(), trkf::KTrack::setVector(), and trkf::syminvert().

Referenced by combineFit(), and trkf::KETrack::setError().

   {
     // Make sure that the two track surfaces are the same.
     // Throw an exception if they are not.
 
     if(!getSurface()->isEqual(*tre.getSurface()))
       throw cet::exception("KETrack") << "Track combination surfaces are not the same.\n";
 
     // Default result is failure.
 
     boost::optional<double> result(false, 0.);
 
     // We will use asymmetric versions of the updating formulas, such
     // that the result is calculated as a perturbation on the
     // better-measured track.  We define the better measured track as
     // the one with the smaller error matrix trace.
 
     // Extract the two state vectors and error matrices as pointers.
 
     const TrackVector* vec1 = &getVector();
     const TrackError* err1 = &getError();
     const TrackVector* vec2 = &tre.getVector();
     const TrackError* err2 = &tre.getError();
 
     // Calculate the traces of the error matrices.
 
     double tr1 = 0;
     for(unsigned int i=0; i<err1->size1(); ++i)
       tr1 += (*err1)(i,i);
 
     double tr2 = 0;
     for(unsigned int i=0; i<err2->size1(); ++i)
       tr2 += (*err2)(i,i);
 
     // Define vec1, err1 as belong to the better measured track.
     // Swap if necessary.
 
     if(tr1 > tr2) {
       const TrackVector* tvec = vec1;
       vec1 = vec2;
       vec2 = tvec;
       const TrackError* terr = err1;
       err1 = err2;
       err2 = terr;
     }
 
     // Calculate the difference vector and difference error matrix.
 
     TrackVector dvec = *vec1 - *vec2;
     TrackError derr = *err1 + *err2;
 
     // Invert the difference error matrix.
     // This is the only place where a detectable failure can occur.
 
     bool ok = syminvert(derr);
     if(ok) {
 
       // Calculate updated state vector.
       // vec1 = vec1 - err1 * derr * dvec
 
       TrackVector tvec1 = prod(derr, dvec);
       TrackVector tvec2 = prod(*err1, tvec1);
       TrackVector tvec3 = *vec1 - tvec2;
       setVector(tvec3);
 
       // Calculate updated error matrix.
       // err1 = err1 - err1 * derr * err1
 
       TrackMatrix terr1 = prod(derr, *err1);
       TrackMatrix terr2 = prod(*err1, terr1);
       TrackError terr2s = ublas::symmetric_adaptor<TrackMatrix>(terr2);
       TrackError terr3 = *err1 - terr2s;
       setError(terr3);
 
       // Calculate chisquare.
       // chisq = dvec^T * derr * dvec
 
       TrackVector dvec1 = prod(derr, dvec);
       double chisq = inner_prod(dvec, dvec1);
       result = boost::optional<double>(true, chisq);
     }
 
     // Final validity check.
 
     if(!isValid())
       result = boost::optional<double>(false, 0.);
 
     // Done.
 
     return result;
   }

double trkf::KFitTrack::getChisq ( ) const

inline

Fit chisquare.

Definition at line 67 of file KFitTrack.h.

References fChisq.

Referenced by trkf::KalmanFilterAlg::buildTrack(), combineFit(), trkf::KalmanFilterAlg::extendTrack(), trkf::KGTrack::fillTrack(), trkf::Track3DKalmanHitAlg::smoothandextendTrack(), and trkf::KalmanFilterAlg::smoothTrack().

Surface::TrackDirection trkf::KTrack::getDirection ( ) const

inherited

Track direction.

Track direction accessor. Track direction implied by track parameters has precedence over track direction attribute. If the surface pointer is null, return UNKNOWN.

Definition at line 72 of file KTrack.cxx.

References trkf::KTrack::fDir, trkf::KTrack::fSurf, trkf::KTrack::fVec, and trkf::Surface::UNKNOWN.

Referenced by trkf::KTrack::getMomentum(), trkf::KTrack::getVector(), trkf::KTrack::isValid(), trkf::Propagator::lin_prop(), trkf::InteractPlane::noise(), trkf::PropYZPlane::origin_vec_prop(), trkf::PropXYZPlane::origin_vec_prop(), trkf::PropYZLine::origin_vec_prop(), trkf::PropYZPlane::short_vec_prop(), and trkf::PropXYZPlane::short_vec_prop().

   {
     Surface::TrackDirection result = Surface::UNKNOWN;
     if(fSurf.get() != 0)
       result = fSurf->getDirection(fVec, fDir);
     return result;
   }

const TrackError& trkf::KETrack::getError ( ) const

inlineinherited

Track error matrix.

Definition at line 54 of file KETrack.h.

References trkf::KETrack::fErr, and trkf::KETrack::PointingError().

Referenced by trkf::KETrack::combineTrack(), trkf::Propagator::err_prop(), trkf::KalmanFilterAlg::extendTrack(), trkf::KGTrack::fillTrack(), trkf::KalmanFilterAlg::fitMomentum(), trkf::KalmanFilterAlg::fitMomentumMS(), trkf::KalmanFilterAlg::fitMomentumRange(), trkf::Propagator::noise_prop(), trkf::KHitMulti::predict(), trkf::KHitWireLine::subpredict(), trkf::KHitWireX::subpredict(), trkf::KHitMulti::update(), trkf::KHit< N >::update(), and trkf::KalmanFilterAlg::updateMomentum().

TrackError& trkf::KETrack::getError ( )

inlineinherited

Modifiable error matrix.

Definition at line 59 of file KETrack.h.

References trkf::KETrack::fErr.

void trkf::KTrack::getMomentum ( double mom[3] ) const

inherited

Get momentum vector of track.

Get momentum vector of track. Throw an exception if track is not valid.

Arguments:

mom - Momentum vector of track.

Definition at line 217 of file KTrack.cxx.

References dir, trkf::KTrack::fSurf, trkf::KTrack::fVec, trkf::KTrack::getDirection(), and trkf::KTrack::isValid().

Referenced by trkf::KGTrack::fillTrack(), trkf::InteractGeneral::noise(), trkf::KTrack::Print(), trkf::KGTrack::Print(), trkf::Track3DKalmanHitAlg::qualityCutsOnSeedTrack(), trkf::KTrack::setPdgCode(), trkf::Propagator::vec_prop(), trkf::KTrack::XLatitude(), and trkf::KTrack::XLongitude().

   {
     if(!isValid())
       throw cet::exception("KTrack") << "Momentum vector requested for invalid track.\n";
     Surface::TrackDirection dir = getDirection();
     fSurf->getMomentum(fVec, mom, dir);
   }

double trkf::KFitTrack::getPath ( ) const

inline

Propagation distance.

Definition at line 66 of file KFitTrack.h.

References fPath.

Referenced by trkf::KGTrack::addTrack(), trkf::KalmanFilterAlg::buildTrack(), trkf::KalmanFilterAlg::extendTrack(), trkf::KGTrack::recalibrate(), and trkf::KalmanFilterAlg::smoothTrack().

void trkf::KTrack::getPosition ( double xyz[3] ) const

inherited

Get position of track.

Get position of track. Throw an exception if track is not valid.

Arguments:

xyz - Position vector.

Definition at line 170 of file KTrack.cxx.

References trkf::KTrack::fSurf, trkf::KTrack::fVec, and trkf::KTrack::isValid().

Referenced by trkf::KGTrack::fillTrack(), trkf::InteractGeneral::noise(), trkf::PropYZLine::origin_vec_prop(), trkf::PropYZPlane::origin_vec_prop(), trkf::PropXYZPlane::origin_vec_prop(), trkf::KTrack::Print(), trkf::KGTrack::Print(), trkf::KTrack::setPdgCode(), trkf::PropXYZPlane::short_vec_prop(), trkf::PropYZPlane::short_vec_prop(), trkf::PropYZLine::short_vec_prop(), and trkf::Propagator::vec_prop().

   {
     if(!isValid())
       throw cet::exception("KTrack") << "Position requested for invalid track.\n";
     fSurf->getPosition(fVec, xyz);
   }

FitStatus trkf::KFitTrack::getStat ( ) const

inline

Fit status.

Definition at line 68 of file KFitTrack.h.

References fStat.

Referenced by combineFit(), trkf::KalmanFilterAlg::extendTrack(), trkf::KalmanFilterAlg::fitMomentumMS(), trkf::KalmanFilterAlg::fitMomentumRange(), trkf::KalmanFilterAlg::smoothTrack(), and trkf::KalmanFilterAlg::updateMomentum().

const std::shared_ptr<const Surface>& trkf::KTrack::getSurface ( ) const

inlineinherited

const TrackVector& trkf::KTrack::getVector ( ) const

inlineinherited

Track state vector.

Definition at line 56 of file KTrack.h.

References trkf::KTrack::fVec, and trkf::KTrack::getDirection().

Referenced by trkf::KETrack::combineTrack(), trkf::KalmanFilterAlg::extendTrack(), trkf::KalmanFilterAlg::fitMomentum(), trkf::KalmanFilterAlg::fitMomentumMS(), trkf::KalmanFilterAlg::fitMomentumRange(), trkf::Propagator::lin_prop(), trkf::InteractPlane::noise(), trkf::PropYZPlane::origin_vec_prop(), trkf::PropXYZPlane::origin_vec_prop(), trkf::PropYZLine::origin_vec_prop(), trkf::KETrack::PointingError(), trkf::KHitMulti::predict(), trkf::PropYZLine::short_vec_prop(), trkf::PropXYZPlane::short_vec_prop(), trkf::PropYZPlane::short_vec_prop(), trkf::KHitWireX::subpredict(), trkf::KHitWireLine::subpredict(), trkf::KHitMulti::update(), trkf::KHit< N >::update(), trkf::KalmanFilterAlg::updateMomentum(), and trkf::Propagator::vec_prop().

TrackVector& trkf::KTrack::getVector ( )

inlineinherited

Modifiable state vector.

Definition at line 63 of file KTrack.h.

References trkf::KTrack::fVec.

bool trkf::KTrack::isValid ( ) const

inherited

Test if track is valid.

A default-constructed or partially-constructed track, is invalid by virtue of having an unknown propagation direction or a null surface pointer.

Tracks can become invaliddynamically for other reasons. This method also does the following checks: a) Check for invalid floating point values (inf and nan). b) Surface-dependent checks via virtual method Surface::isTrackValid.

Definition at line 90 of file KTrack.cxx.

References trkf::KTrack::fSurf, trkf::KTrack::fVec, trkf::KTrack::getDirection(), geo::vect::isfinite(), and trkf::Surface::UNKNOWN.

Referenced by trkf::KGTrack::addTrack(), trkf::KalmanFilterAlg::buildTrack(), trkf::KETrack::combineTrack(), trkf::KalmanFilterAlg::extendTrack(), trkf::KTrack::getMomentum(), trkf::KTrack::getPosition(), trkf::Propagator::lin_prop(), trkf::PropYZLine::origin_vec_prop(), trkf::PropXYZPlane::origin_vec_prop(), trkf::PropYZPlane::origin_vec_prop(), trkf::KETrack::PointingError(), trkf::KTrack::setPdgCode(), and trkf::KalmanFilterAlg::smoothTrack().

   {
     bool result = true;
 
     // Check for valid direction.
 
     if(getDirection() == Surface::UNKNOWN)
       result = false;
 
     // Check for non-null surface pointer (for safety, should be redundant
     // with previous check.
 
     if(result && fSurf.get() == 0)
       result = false;
 
     // Check for track parameters containing invalid floating point
     // values.
 
     if(result) {
       for(unsigned int i=0; i<fVec.size(); ++i) {
         if(!std::isfinite(fVec(i))) {
           result = false;
           break;
         }
       }
     }
 
     // Surface-dependent check on track validity.
 
     if(result && !fSurf->isTrackValid(fVec))
       result = false;
 
     // Done.
 
     return result;
   }

double trkf::KTrack::Mass ( ) const

inherited

Based on pdg code.

Particle mass based on pdg code.

Definition at line 128 of file KTrack.cxx.

References trkf::KTrack::fPdgCode.

Referenced by trkf::KalmanFilterAlg::fitMomentumMS(), trkf::InteractPlane::noise(), trkf::KTrack::PdgCode(), trkf::PropXYZPlane::short_vec_prop(), trkf::PropYZLine::short_vec_prop(), trkf::PropYZPlane::short_vec_prop(), and trkf::Propagator::vec_prop().

   {
     double mass = 0.;
     int apdg = std::abs(fPdgCode);
 
     // Muon
 
     if(apdg == 13)
       mass = mumass;
 
     // Charged pion
 
     else if(apdg == 211)
       mass = pimass;
 
     // Charged kaon
 
     else if(apdg == 321)
       mass = kmass;
 
     // (Anti)proton
 
     else if(apdg == 2212)
       mass = pmass;
 
     // Anything else throw exception
 
     else
       throw cet::exception("KTrack") << "Mass requested for invalid pdg id = " << fPdgCode << "\n";
 
     // Done.
 
     return mass;
   }

int trkf::KTrack::PdgCode ( ) const

inlineinherited

Pdg code.

Definition at line 58 of file KTrack.h.

References trkf::KTrack::fPdgCode, and trkf::KTrack::Mass().

double trkf::KETrack::PointingError ( ) const

inherited

Pointing error (radians).

Calculate track pointing error (sigma, in radians).

This method calculates a single pointing error (sigma, in radians) based on the track parameters and error matrix. The actual calculation is done by the similarly names method of the surface class, since this class doesn't know what the track parameters mean.

Definition at line 74 of file KETrack.cxx.

References trkf::KETrack::fErr, trkf::KTrack::getSurface(), trkf::KTrack::getVector(), and trkf::KTrack::isValid().

Referenced by trkf::KalmanFilterAlg::buildTrack(), trkf::KalmanFilterAlg::extendTrack(), trkf::KETrack::getError(), and trkf::KETrack::Print().

   {
     if(!isValid())
       throw cet::exception("KETrack") << "Pointing error requested for invalid track.\n";
     return getSurface()->PointingError(getVector(), fErr);
   }

std::ostream & trkf::KFitTrack::Print	(	std::ostream &	out,
		bool	doTitle = `true`
	)		const

virtual

Printout.

Reimplemented from trkf::KETrack.

Reimplemented in trkf::KHitsTrack, and trkf::KHitTrack.

Definition at line 124 of file KFitTrack.cxx.

References BACKWARD, BACKWARD_PREDICTED, fChisq, FORWARD, FORWARD_PREDICTED, fPath, fStat, INVALID, OPTIMAL, trkf::KETrack::Print(), and UNKNOWN.

Referenced by trkf::KHitTrack::Print(), trkf::KHitsTrack::Print(), and setStat().

   {
     if(doTitle)
       out << "KFitTrack:\n";
 
     // Print information specific to this class.
 
     out << "  Distance = " << fPath << "\n";
     out << "  Chisquare = " << fChisq << "\n";
     out << "  Status = " << (fStat == INVALID ? "INVALID" :
                              (fStat == UNKNOWN ? "UNKNOWN" :
                               (fStat == FORWARD ? "FORWARD" :
                                (fStat == FORWARD_PREDICTED ? "FORWARD_PREDICTED" :
                                 (fStat == BACKWARD ? "BACKWARD" :
                                  (fStat == BACKWARD_PREDICTED ? "BACKWARD_PREDICTED" :
                                   (fStat == OPTIMAL ? "OPTIMAL" : 
                                    "OPTIMAL_PREDICTED"))))))) << "\n";
 
     // Print base class.
 
     KETrack::Print(out, false);
     return out;
   }

void trkf::KFitTrack::setChisq ( double chisq )

inline

Set chisquare.

Definition at line 73 of file KFitTrack.h.

References fChisq.

Referenced by trkf::KalmanFilterAlg::buildTrack(), trkf::KalmanFilterAlg::extendTrack(), and trkf::KalmanFilterAlg::smoothTrack().

void trkf::KTrack::setDirection ( Surface::TrackDirection dir )

inlineinherited

Set direction.

Definition at line 68 of file KTrack.h.

References dir, and trkf::KTrack::fDir.

Referenced by trkf::Propagator::lin_prop(), trkf::PropXYZPlane::origin_vec_prop(), trkf::PropYZLine::origin_vec_prop(), and trkf::PropYZPlane::origin_vec_prop().

void trkf::KETrack::setError ( const TrackError & err )

inlineinherited

Set error matrix.

Definition at line 60 of file KETrack.h.

References trkf::KETrack::combineTrack(), trkf::KETrack::fErr, and trkf::KETrack::Print().

Referenced by trkf::KETrack::combineTrack(), trkf::Propagator::err_prop(), trkf::Propagator::noise_prop(), trkf::KHitMulti::update(), and trkf::KHit< N >::update().

void trkf::KFitTrack::setPath ( double path )

inline

Set propagation distance.

Definition at line 72 of file KFitTrack.h.

References fPath.

Referenced by trkf::KalmanFilterAlg::buildTrack(), trkf::KalmanFilterAlg::extendTrack(), trkf::KGTrack::recalibrate(), and trkf::KalmanFilterAlg::smoothTrack().

void trkf::KTrack::setPdgCode ( int pdg )

inlineinherited

Set pdg code.

Definition at line 69 of file KTrack.h.

References trkf::KTrack::fPdgCode, trkf::KTrack::getMomentum(), trkf::KTrack::getPosition(), trkf::KTrack::isValid(), trkf::KTrack::Print(), trkf::KTrack::XLatitude(), and trkf::KTrack::XLongitude().

void trkf::KFitTrack::setStat ( FitStatus stat )

inline

Set fit status.

Definition at line 74 of file KFitTrack.h.

References combineFit(), fStat, and Print().

Referenced by trkf::KalmanFilterAlg::buildTrack(), trkf::KalmanFilterAlg::extendTrack(), and trkf::KalmanFilterAlg::smoothTrack().

void trkf::KTrack::setSurface ( const std::shared_ptr< const Surface > & psurf )

inlineinherited

Set surface.

Definition at line 66 of file KTrack.h.

References trkf::KTrack::fSurf.

Referenced by trkf::Propagator::lin_prop(), trkf::PropYZPlane::origin_vec_prop(), trkf::PropXYZPlane::origin_vec_prop(), trkf::PropYZLine::origin_vec_prop(), trkf::PropYZPlane::short_vec_prop(), trkf::PropYZLine::short_vec_prop(), and trkf::PropXYZPlane::short_vec_prop().

66 {fSurf = psurf;}

trkf::KTrack::fSurf

std::shared_ptr< const Surface > fSurf

Track surface.

Definition: KTrack.h:93

void trkf::KTrack::setVector ( const TrackVector & vec )

inlineinherited

Set state vector.

Definition at line 67 of file KTrack.h.

References trkf::KTrack::fVec.

Referenced by trkf::KETrack::combineTrack(), trkf::Propagator::lin_prop(), trkf::PropYZPlane::origin_vec_prop(), trkf::PropXYZPlane::origin_vec_prop(), trkf::PropYZLine::origin_vec_prop(), trkf::PropXYZPlane::short_vec_prop(), trkf::PropYZLine::short_vec_prop(), trkf::PropYZPlane::short_vec_prop(), trkf::KHitMulti::update(), and trkf::KHit< N >::update().

double trkf::KTrack::XLatitude ( ) const

inherited

Get x-latitude.

The x-latitude is the latitude defined with respect to the x-axis. The x-latitude is zero of the track is traveling parallel to the wire planes.

Definition at line 183 of file KTrack.cxx.

References trkf::KTrack::getMomentum().

Referenced by trkf::KTrack::Print(), and trkf::KTrack::setPdgCode().

   {
     double mom[3];
     getMomentum(mom);
     double ptx = std::sqrt(mom[1]*mom[1] + mom[2]*mom[2]);
     double result = 0.;
     if(ptx > 0. || mom[0] > 0.)
       result = atan2(mom[0], ptx);
     return result;
   }

double trkf::KTrack::XLongitude ( ) const

inherited

Get x-longitude.

The x-longitude is the longitude defined with respect to the y- and z-axes. The x-longitude is zero of the track is parallel to the z-axis in the yz-plane.

Definition at line 200 of file KTrack.cxx.

References trkf::KTrack::getMomentum().

Referenced by trkf::KTrack::Print(), and trkf::KTrack::setPdgCode().

   {
     double mom[3];
     getMomentum(mom);
     double result = 0.;
     if(mom[1] != 0. || mom[2] != 0.)
       result = atan2(mom[1], mom[2]);
     return result;
   }

Member Data Documentation

double trkf::KFitTrack::fChisq

private

Fit chisquare.

Definition at line 87 of file KFitTrack.h.

Referenced by combineFit(), getChisq(), Print(), and setChisq().

double trkf::KFitTrack::fPath

private

Propagation distance.

Definition at line 86 of file KFitTrack.h.

Referenced by getPath(), Print(), and setPath().

FitStatus trkf::KFitTrack::fStat

private

Fit status.

Definition at line 88 of file KFitTrack.h.

Referenced by combineFit(), getStat(), Print(), and setStat().

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

/cvmfs/larsoft.opensciencegrid.org/products/lardata/v07_01_01/source/lardata/RecoObjects/KFitTrack.h
/cvmfs/larsoft.opensciencegrid.org/products/lardata/v07_01_01/source/lardata/RecoObjects/KFitTrack.cxx

Public Types

Public Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation