#include <SteppingHelixStateInfo.h>

Classes
struct	VolumeBounds

Public Types
typedef CLHEP::Hep3Vector	Point

enum	Result { OK = 0, FAULT, APPROX, RANGEOUT, INACC, NOT_IMPLEMENTED, WRONG_VOLUME, UNDEFINED, MAX_RESULT }

typedef CLHEP::Hep3Vector	Vector

Public Member Functions
int	charge () const

void	getFreeState (FreeTrajectoryState &fts) const
	convert internal structure into the fts More...

TrajectoryStateOnSurface	getStateOnSurface (const Surface &surf, bool returnTangentPlane=false) const

bool	hasErrorPropagated () const

bool	isValid () const

GlobalVector	momentum () const

double	path () const

GlobalPoint	position () const

double	radPath () const

Result	status () const

	SteppingHelixStateInfo ()

	SteppingHelixStateInfo (const FreeTrajectoryState &fts)

Static Public Attributes
static const std::string	ResultName [MAX_RESULT]

Protected Attributes
Vector	bf

Vector	bfGradLoc

AlgebraicSymMatrix55	covCurv

double	dEdx

double	dEdXPrime

double	dir

const MagneticField *	field

bool	hasErrorPropagated_

bool	isComplete

bool	isValid_

bool	isYokeVol

const MagVolume *	magVol

AlgebraicSymMatrix55	matDCovCurv

Vector	p3

double	path_

int	q

Point	r3

double	radPath_

double	radX0

VolumeBounds	rzLims

Result	status_

Friends
class	SteppingHelixPropagator

Detailed Description

Holder of SteppingHelixState information

Author: Vyacheslav Krutelyov (slava77)

Implementation part of the stepping helix propagator state data structure

Author: Vyacheslav Krutelyov (slava77)

Definition at line 27 of file SteppingHelixStateInfo.h.

Member Typedef Documentation

◆ Point

typedef CLHEP::Hep3Vector SteppingHelixStateInfo::Point

Definition at line 32 of file SteppingHelixStateInfo.h.

◆ Vector

typedef CLHEP::Hep3Vector SteppingHelixStateInfo::Vector

Definition at line 31 of file SteppingHelixStateInfo.h.

Member Enumeration Documentation

◆ Result

enum SteppingHelixStateInfo::Result

Enumerator
OK
FAULT
APPROX
RANGEOUT
INACC
NOT_IMPLEMENTED
WRONG_VOLUME
UNDEFINED
MAX_RESULT

Definition at line 34 of file SteppingHelixStateInfo.h.

34 { OK = 0, FAULT, APPROX, RANGEOUT, INACC, NOT_IMPLEMENTED, WRONG_VOLUME, UNDEFINED, MAX_RESULT };

Constructor & Destructor Documentation

◆ SteppingHelixStateInfo() [1/2]

SteppingHelixStateInfo::SteppingHelixStateInfo ( )

inline

Definition at line 38 of file SteppingHelixStateInfo.h.

       : path_(0),
         radPath_(0),
         dir(0),
         magVol(nullptr),
         isYokeVol(false),
         field(nullptr),
         dEdx(0),
         dEdXPrime(0),
         radX0(1e12),
         isComplete(false),
         isValid_(false),
         hasErrorPropagated_(false),
         status_(UNDEFINED) {}

◆ SteppingHelixStateInfo() [2/2]

SteppingHelixStateInfo::SteppingHelixStateInfo ( const FreeTrajectoryState & fts )

Definition at line 27 of file SteppingHelixStateInfo.cc.

     : path_(0),
       radPath_(0),
       dir(0),
       magVol(nullptr),
       isYokeVol(false),
       field(nullptr),
       dEdx(0),
       dEdXPrime(0),
       radX0(1e12),
       status_(UNDEFINED) {
   p3.set(fts.momentum().x(), fts.momentum().y(), fts.momentum().z());
   r3.set(fts.position().x(), fts.position().y(), fts.position().z());
   q = fts.charge();
  
   if (fts.hasError()) {
     covCurv = fts.curvilinearError().matrix();
     hasErrorPropagated_ = true;
   } else {
     covCurv = AlgebraicSymMatrix55();
     hasErrorPropagated_ = false;
   }
   static const std::string metname = "SteppingHelixPropagator";
   if (fts.hasError()) {
     LogTrace(metname) << "Created SHPStateInfo from FTS\n" << fts;
     //    LogTrace(metname)<<"and cartesian error of\n"<<fts.cartesianError().matrix();
   } else
     LogTrace(metname) << "Created SHPStateInfo from FTS without errors";
  
   isComplete = false;
   isValid_ = true;
 }

References FreeTrajectoryState::charge(), covCurv, FreeTrajectoryState::curvilinearError(), FreeTrajectoryState::hasError(), hasErrorPropagated_, isComplete, isValid_, LogTrace, CurvilinearTrajectoryError::matrix(), metname, FreeTrajectoryState::momentum(), p3, FreeTrajectoryState::position(), q, r3, AlCaHLTBitMon_QueryRunRegistry::string, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Member Function Documentation

◆ charge()

int SteppingHelixStateInfo::charge ( void ) const

inline

Definition at line 61 of file SteppingHelixStateInfo.h.

61 { return q; }

References q.

◆ getFreeState()

void SteppingHelixStateInfo::getFreeState ( FreeTrajectoryState & fts ) const

convert internal structure into the fts

Definition at line 88 of file SteppingHelixStateInfo.cc.

                                                                         {
   if (isValid()) {
     GlobalVector p3GV(p3.x(), p3.y(), p3.z());
     GlobalPoint r3GP(r3.x(), r3.y(), r3.z());
     GlobalTrajectoryParameters tPars(r3GP, p3GV, q, field);
     //    CartesianTrajectoryError tCov(cov);
     //    CurvilinearTrajectoryError tCCov(ROOT::Math::Similarity(JacobianCartesianToCurvilinear(tPars).jacobian(), cov));
     CurvilinearTrajectoryError tCCov(covCurv);
  
     fts = (hasErrorPropagated_) ? FreeTrajectoryState(tPars, tCCov) : FreeTrajectoryState(tPars);
     //      ? FreeTrajectoryState(tPars, tCov, tCCov) : FreeTrajectoryState(tPars);
     //    if (fts.hasError()) fts.curvilinearError(); //call it so it gets created
   }
 }

References covCurv, field, hasErrorPropagated_, isValid(), p3, q, and r3.

Referenced by CachedTrajectory::propagate(), and SteppingHelixPropagator::propagateWithPath().

◆ getStateOnSurface()

TrajectoryStateOnSurface SteppingHelixStateInfo::getStateOnSurface	(	const Surface &	surf,
		bool	returnTangentPlane = `false`
	)		const

Definition at line 60 of file SteppingHelixStateInfo.cc.

                                                                                                                      {
   static const std::string metname = "SteppingHelixPropagator";
   if (!isValid())
     LogTrace(metname) << "Return TSOS is invalid";
   else
     LogTrace(metname) << "Return TSOS is valid";
   if (!isValid())
     return TrajectoryStateOnSurface();
   GlobalVector p3GV(p3.x(), p3.y(), p3.z());
   GlobalPoint r3GP(r3.x(), r3.y(), r3.z());
   GlobalTrajectoryParameters tPars(r3GP, p3GV, q, field);
   //  CartesianTrajectoryError tCov(cov);
  
   //  CurvilinearTrajectoryError tCCov(ROOT::Math::Similarity(JacobianCartesianToCurvilinear(tPars).jacobian(), cov));
   CurvilinearTrajectoryError tCCov(covCurv);
  
   FreeTrajectoryState fts(tPars, tCCov);
   if (!hasErrorPropagated_)
     fts = FreeTrajectoryState(tPars);
  
   SurfaceSideDefinition::SurfaceSide side = SurfaceSideDefinition::atCenterOfSurface;
   if (dir > 0)
     side = SurfaceSideDefinition::beforeSurface;
   if (dir < 0)
     side = SurfaceSideDefinition::afterSurface;
   return TrajectoryStateOnSurface(fts, returnTangentPlane ? *surf.tangentPlane(fts.position()) : surf, side);
 }