BasicTrajectoryState Class Reference

#include <BasicTrajectoryState.h>

Inheritance diagram for BasicTrajectoryState:

Public Types
typedef BasicTrajectoryState	BTSOS
typedef ProxyBase< BTSOS, CopyUsingClone< BTSOS > >	Proxy
typedef ReferenceCountingPointer < BasicTrajectoryState >	RCPtr
typedef SurfaceSideDefinition::SurfaceSide	SurfaceSide
typedef Surface	SurfaceType

Public Member Functions
	BasicTrajectoryState ()
	BasicTrajectoryState (const SurfaceType &aSurface)
	construct invalid trajectory state (without parameters) More...
	BasicTrajectoryState (const FreeTrajectoryState &fts, const SurfaceType &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
	BasicTrajectoryState (const GlobalTrajectoryParameters &par, const SurfaceType &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
	BasicTrajectoryState (const GlobalTrajectoryParameters &par, const CartesianTrajectoryError &err, const SurfaceType &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
	BasicTrajectoryState (const GlobalTrajectoryParameters &par, const CurvilinearTrajectoryError &err, const SurfaceType &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface, double weight=1.)
	BasicTrajectoryState (const GlobalTrajectoryParameters &par, const CurvilinearTrajectoryError &err, const SurfaceType &aSurface, double weight)
	BasicTrajectoryState (const LocalTrajectoryParameters &par, const SurfaceType &aSurface, const MagneticField *field, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
	BasicTrajectoryState (const LocalTrajectoryParameters &par, const LocalTrajectoryError &err, const SurfaceType &aSurface, const MagneticField *field, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface, double weight=1.)
	BasicTrajectoryState (const LocalTrajectoryParameters &par, const LocalTrajectoryError &err, const SurfaceType &aSurface, const MagneticField *field, double weight)
virtual bool	canUpdateLocalParameters () const
const CartesianTrajectoryError	cartesianError () const
TrackCharge	charge () const
virtual BasicTrajectoryState *	clone () const =0
virtual std::vector < TrajectoryStateOnSurface >	components () const
const CurvilinearTrajectoryError &	curvilinearError () const
FreeTrajectoryState const *	freeTrajectoryState (bool withErrors=true) const
GlobalVector	globalDirection () const
GlobalVector	globalMomentum () const
const GlobalTrajectoryParameters &	globalParameters () const
GlobalPoint	globalPosition () const
bool	hasError () const
bool	isValid () const
LocalVector	localDirection () const
const LocalTrajectoryError &	localError () const
LocalVector	localMomentum () const
const LocalTrajectoryParameters &	localParameters () const
LocalPoint	localPosition () const
const MagneticField *	magneticField () const
void	rescaleError (double factor)
double	signedInverseMomentum () const
const SurfaceType &	surface () const
SurfaceSide	surfaceSide () const
	Position relative to material, defined relative to momentum vector. More...
double	transverseCurvature () const
virtual void	update (const LocalTrajectoryParameters &p, const SurfaceType &aSurface, const MagneticField *field, const SurfaceSide side)
virtual void	update (const LocalTrajectoryParameters &p, const LocalTrajectoryError &err, const SurfaceType &aSurface, const MagneticField *field, const SurfaceSide side, double weight)
double	weight () const
virtual	~BasicTrajectoryState ()
Public Member Functions inherited from BasicReferenceCounted
void	addReference () const
	BasicReferenceCounted ()
	BasicReferenceCounted (const BasicReferenceCounted &iRHS)
const BasicReferenceCounted &	operator= (const BasicReferenceCounted &)
unsigned int	references () const
void	removeReference () const
virtual	~BasicReferenceCounted ()

Private Member Functions
void	checkCurvilinError () const
void	createLocalError () const
void	createLocalErrorFromCurvilinearError () const dso_internal
void	createLocalParameters () const
void	missingError (char const *where) const

Static Private Member Functions
static void	notValid ()

Private Attributes
FreeTrajectoryState	theFreeState
LocalTrajectoryError	theLocalError
LocalTrajectoryParameters	theLocalParameters
bool	theLocalParametersValid
ConstReferenceCountingPointer < SurfaceType >	theSurfaceP
SurfaceSide	theSurfaceSide
bool	theValid
double	theWeight

Friends
class	ProxyBase< BTSOS, CopyUsingClone< BTSOS > >
class	ReferenceCountingPointer< BasicTrajectoryState >

Detailed Description

No so Abstract (anyore) base class for TrajectoryState. It is ReferenceCounted.

VI 8/12/2011 content of BasicSingleTrajectoryState moved here.... fully devirtualized

Definition at line 75 of file BasicTrajectoryState.h.

Member Typedef Documentation

typedef BasicTrajectoryState BasicTrajectoryState::BTSOS

Definition at line 78 of file BasicTrajectoryState.h.

typedef ProxyBase< BTSOS, CopyUsingClone<BTSOS> > BasicTrajectoryState::Proxy

Definition at line 79 of file BasicTrajectoryState.h.

typedef ReferenceCountingPointer<BasicTrajectoryState> BasicTrajectoryState::RCPtr

Definition at line 80 of file BasicTrajectoryState.h.

typedef SurfaceSideDefinition::SurfaceSide BasicTrajectoryState::SurfaceSide

Definition at line 81 of file BasicTrajectoryState.h.

typedef Surface BasicTrajectoryState::SurfaceType

Definition at line 82 of file BasicTrajectoryState.h.

Constructor & Destructor Documentation

BasicTrajectoryState::BasicTrajectoryState ( )

inline

Definition at line 90 of file BasicTrajectoryState.h.

: theValid(false), theWeight(0){}

BasicTrajectoryState::BasicTrajectoryState ( const SurfaceType &  aSurface )

explicit

construct invalid trajectory state (without parameters)

Definition at line 176 of file BasicTrajectoryState.cc.

                                                  :
  theLocalError(InvalidError()),
  theLocalParameters(),
  theLocalParametersValid(false),
  theValid(false),
  theSurfaceSide(SurfaceSideDefinition::atCenterOfSurface), 
  theSurfaceP( &aSurface), 
  theWeight(0)
{}

BasicTrajectoryState::~BasicTrajectoryState ( )

virtual

Definition at line 40 of file BasicTrajectoryState.cc.

{}

BasicTrajectoryState::BasicTrajectoryState	(	const FreeTrajectoryState &	fts,
		const SurfaceType &	aSurface,
		const SurfaceSide	side = SurfaceSideDefinition::atCenterOfSurface
	)

Constructor from FTS and surface. For surfaces with material the side of the surface should be specified explicitely.

Definition at line 55 of file BasicTrajectoryState.cc.

                                        :
  theFreeState(fts),
  theLocalError(InvalidError()),
  theLocalParameters(),
  theLocalParametersValid(false),
  theValid(true),
  theSurfaceSide(side), 
  theSurfaceP( &aSurface), 
  theWeight(1.)
{}    

BasicTrajectoryState::BasicTrajectoryState	(	const GlobalTrajectoryParameters &	par,
		const SurfaceType &	aSurface,
		const SurfaceSide	side = SurfaceSideDefinition::atCenterOfSurface
	)

Constructor from global parameters and surface. For surfaces with material the side of the surface should be specified explicitely.

Definition at line 69 of file BasicTrajectoryState.cc.

                                        :
  theFreeState(par),
  theLocalError(InvalidError()),
  theLocalParameters(),
  theLocalParametersValid(false),
  theValid(true),
  theSurfaceSide(side), 
  theSurfaceP( &aSurface), 
  theWeight(1.)
{}

BasicTrajectoryState::BasicTrajectoryState	(	const GlobalTrajectoryParameters &	par,
		const CartesianTrajectoryError &	err,
		const SurfaceType &	aSurface,
		const SurfaceSide	side = SurfaceSideDefinition::atCenterOfSurface
	)

Constructor from global parameters, errors and surface. For surfaces with material the side of the surface should be specified explicitely.

Definition at line 83 of file BasicTrajectoryState.cc.

                                        :
  theFreeState(par, err),
  theLocalError(InvalidError()),
  theLocalParameters(),
  theLocalParametersValid(false),
  theValid(true),
  theSurfaceSide(side), 
  theSurfaceP( &aSurface), 
  theWeight(1.)
{}

BasicTrajectoryState::BasicTrajectoryState	(	const GlobalTrajectoryParameters &	par,
		const CurvilinearTrajectoryError &	err,
		const SurfaceType &	aSurface,
		const SurfaceSide	side = SurfaceSideDefinition::atCenterOfSurface,
		double	weight = 1.
	)

Constructor from global parameters, errors and surface. For surfaces with material the side of the surface should be specified explicitely. For multi-states the weight should be specified explicitely.

Definition at line 98 of file BasicTrajectoryState.cc.

                               :
  theFreeState(par, err),
  theLocalError(InvalidError()),
  theLocalParameters(),
  theLocalParametersValid(false),
  theValid(true),
  theSurfaceSide(side), 
  theSurfaceP( &aSurface), 
  theWeight(weight)
{}

BasicTrajectoryState::BasicTrajectoryState	(	const GlobalTrajectoryParameters &	par,
		const CurvilinearTrajectoryError &	err,
		const SurfaceType &	aSurface,
		double	weight
	)

Constructor from global parameters, errors and surface. For multi-states the weight should be specified explicitely. For backward compatibility without specification of the side of the surface.

Definition at line 114 of file BasicTrajectoryState.cc.

                               :
  theFreeState(par, err),
  theLocalError(InvalidError()),
  theLocalParameters(),
  theLocalParametersValid(false),
  theValid(true),
  theSurfaceSide(SurfaceSideDefinition::atCenterOfSurface), 
  theSurfaceP( &aSurface), 
  theWeight(weight)
{}

BasicTrajectoryState::BasicTrajectoryState	(	const LocalTrajectoryParameters &	par,
		const SurfaceType &	aSurface,
		const MagneticField *	field,
		const SurfaceSide	side = SurfaceSideDefinition::atCenterOfSurface
	)

Constructor from local parameters, errors and surface. For surfaces with material the side of the surface should be specified explicitely.

Definition at line 129 of file BasicTrajectoryState.cc.

                                        :
  theFreeState(makeFTS(par,aSurface,field)),
  theLocalError(InvalidError()),
  theLocalParameters(par),
  theLocalParametersValid(true),
  theValid(true),
   theSurfaceSide(side),
  theSurfaceP( &aSurface), 
  theWeight(1.)
{}

BasicTrajectoryState::BasicTrajectoryState	(	const LocalTrajectoryParameters &	par,
		const LocalTrajectoryError &	err,
		const SurfaceType &	aSurface,
		const MagneticField *	field,
		const SurfaceSide	side = SurfaceSideDefinition::atCenterOfSurface,
		double	weight = 1.
	)

Constructor from local parameters, errors and surface. For surfaces with material the side of the surface should be specified explicitely. For multi-states the weight should be specified explicitely.

Definition at line 144 of file BasicTrajectoryState.cc.

                               :
  theFreeState(makeFTS(par,aSurface,field)),
  theLocalError(err),
  theLocalParameters(par),
  theLocalParametersValid(true),
  theValid(true),
  theSurfaceSide(side), 
  theSurfaceP( &aSurface),
  theWeight(weight)
{}

BasicTrajectoryState::BasicTrajectoryState	(	const LocalTrajectoryParameters &	par,
		const LocalTrajectoryError &	err,
		const SurfaceType &	aSurface,
		const MagneticField *	field,
		double	weight
	)

Constructor from local parameters, errors and surface. For multi-states the weight should be specified explicitely. For backward compatibility without specification of the side of the surface.

Definition at line 161 of file BasicTrajectoryState.cc.

                               :
  theFreeState(makeFTS(par,aSurface,field)),
  theLocalError(err),
  theLocalParameters(par),
  theLocalParametersValid(true),
  theValid(true),
  theSurfaceSide(SurfaceSideDefinition::atCenterOfSurface),
  theSurfaceP( &aSurface), 
  theWeight(weight){}

Member Function Documentation

virtual bool BasicTrajectoryState::canUpdateLocalParameters ( ) const

inlinevirtual

Reimplemented in GCC11_FINAL< T, TOPO >.

Definition at line 264 of file BasicTrajectoryState.h.

{ return true; }

const CartesianTrajectoryError BasicTrajectoryState::cartesianError ( ) const

inline

Definition at line 187 of file BasicTrajectoryState.h.

References FreeTrajectoryState::cartesianError(), freeTrajectoryState(), hasError(), missingError(), and unlikely.

                                                        {
    if unlikely(!hasError()) {
    missingError(" accesing cartesian error.");
    return CartesianTrajectoryError();
      }
    return freeTrajectoryState(true)->cartesianError();
  }

TrackCharge BasicTrajectoryState::charge ( void  ) const

inline

Definition at line 177 of file BasicTrajectoryState.h.

References FreeTrajectoryState::charge(), and theFreeState.

                             {
    return theFreeState.charge();
  }

void BasicTrajectoryState::checkCurvilinError ( ) const

private

Definition at line 220 of file BasicTrajectoryState.cc.

References createLocalParameters(), globalParameters(), FreeTrajectoryState::hasCurvilinearError(), JacobianLocalToCurvilinear::jacobian(), likely, localParameters(), magneticField(), LocalTrajectoryError::matrix(), FreeTrajectoryState::setCurvilinearError(), surface(), theFreeState, theLocalError, theLocalParametersValid, and unlikely.

Referenced by freeTrajectoryState().

                                                    {
  if likely(theFreeState.hasCurvilinearError()) return;

  if unlikely(!theLocalParametersValid) createLocalParameters();
  
  JacobianLocalToCurvilinear loc2Curv(surface(), localParameters(), globalParameters(), *magneticField());
  const AlgebraicMatrix55& jac = loc2Curv.jacobian();
  const AlgebraicSymMatrix55 &cov = ROOT::Math::Similarity(jac, theLocalError.matrix());

  theFreeState.setCurvilinearError( cov );
  
  verifyLocalErr(theLocalError);
  verifyCurvErr(cov); 
}

virtual BasicTrajectoryState* BasicTrajectoryState::clone ( ) const

pure virtual

Implemented in GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, and GCC11_FINAL< T, TOPO >.

std::vector< TrajectoryStateOnSurface > BasicTrajectoryState::components ( ) const

virtual

Reimplemented in GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, GCC11_FINAL< T, TOPO >, and GCC11_FINAL< T, TOPO >.

Definition at line 336 of file BasicTrajectoryState.cc.

References query::result.

                                       {
  std::vector<TrajectoryStateOnSurface> result; result.reserve(1);
  result.push_back( const_cast<BasicTrajectoryState*>(this));
  return result;
}

void BasicTrajectoryState::createLocalError ( ) const

private

Definition at line 249 of file BasicTrajectoryState.cc.

References createLocalErrorFromCurvilinearError(), FreeTrajectoryState::hasCurvilinearError(), likely, theFreeState, and theLocalError.

Referenced by localError().

                                                  {
  if likely(theFreeState.hasCurvilinearError())
    createLocalErrorFromCurvilinearError();
  else theLocalError = InvalidError();
}

void BasicTrajectoryState::createLocalErrorFromCurvilinearError ( ) const

private

Definition at line 256 of file BasicTrajectoryState.cc.

References FreeTrajectoryState::curvilinearError(), globalParameters(), JacobianCurvilinearToLocal::jacobian(), localParameters(), magneticField(), CurvilinearTrajectoryError::matrix(), surface(), theFreeState, and theLocalError.

Referenced by createLocalError().

                                                                 {
  
  JacobianCurvilinearToLocal curv2Loc(surface(), localParameters(), globalParameters(), *magneticField());
  const AlgebraicMatrix55& jac = curv2Loc.jacobian();
  
  const AlgebraicSymMatrix55 &cov = 
    ROOT::Math::Similarity(jac, theFreeState.curvilinearError().matrix());
  //    cout<<"Clocal via curvilinear error"<<endl;
  theLocalError = LocalTrajectoryError(cov);

  verifyCurvErr(theFreeState.curvilinearError());
  verifyLocalErr(theLocalError);

}

void BasicTrajectoryState::createLocalParameters ( ) const

private

Definition at line 238 of file BasicTrajectoryState.cc.

References FreeTrajectoryState::charge(), PV3DBase< T, PVType, FrameType >::mag(), FreeTrajectoryState::momentum(), FreeTrajectoryState::position(), FreeTrajectoryState::signedInverseMomentum(), surface(), theFreeState, theLocalParameters, theLocalParametersValid, GloballyPositioned< T >::toLocal(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by checkCurvilinError(), and localParameters().

                                                       {
  LocalPoint  x = surface().toLocal(theFreeState.position());
  LocalVector p = surface().toLocal(theFreeState.momentum());
// believe p.z() never exactly equals 0.
  bool isCharged = theFreeState.charge()!=0;
  theLocalParameters =
    LocalTrajectoryParameters(isCharged?theFreeState.signedInverseMomentum():1./p.mag(),
      p.x()/p.z(), p.y()/p.z(), x.x(), x.y(), p.z()>0. ? 1.:-1., isCharged);
  theLocalParametersValid = true;
}

const CurvilinearTrajectoryError& BasicTrajectoryState::curvilinearError ( ) const

inline

Definition at line 194 of file BasicTrajectoryState.h.

References FreeTrajectoryState::curvilinearError(), freeTrajectoryState(), hasError(), missingError(), and unlikely.

                                                             {
    if unlikely(!hasError()) {
    missingError(" accesing curvilinearerror.");
    static CurvilinearTrajectoryError crap;
    return crap;
      }
    return freeTrajectoryState(true)->curvilinearError();
  }

FreeTrajectoryState const* BasicTrajectoryState::freeTrajectoryState ( bool  withErrors = true ) const

inline

Definition at line 205 of file BasicTrajectoryState.h.

References checkCurvilinError(), hasError(), isValid(), notValid(), theFreeState, and unlikely.

Referenced by cartesianError(), and curvilinearError().

                                                                             {
    if unlikely(!isValid()) notValid();
    if(withErrors && hasError()) { // this is the right thing
      checkCurvilinError();
    }
    return &theFreeState;
  }

GlobalVector BasicTrajectoryState::globalDirection ( ) const

inline

Definition at line 174 of file BasicTrajectoryState.h.

References FreeTrajectoryState::momentum(), theFreeState, and Vector3DBase< T, FrameTag >::unit().

                                       {
    return theFreeState.momentum().unit();
  }  

GlobalVector BasicTrajectoryState::globalMomentum ( ) const

inline

Definition at line 171 of file BasicTrajectoryState.h.

References FreeTrajectoryState::momentum(), and theFreeState.

                                      {
    return theFreeState.momentum();
  }

const GlobalTrajectoryParameters& BasicTrajectoryState::globalParameters ( void  ) const

inline

Definition at line 165 of file BasicTrajectoryState.h.

References FreeTrajectoryState::parameters(), and theFreeState.

Referenced by checkCurvilinError(), and createLocalErrorFromCurvilinearError().

                                                             {
    return theFreeState.parameters();
  }

GlobalPoint BasicTrajectoryState::globalPosition ( ) const

inline

Definition at line 168 of file BasicTrajectoryState.h.

References FreeTrajectoryState::position(), and theFreeState.

                                     {
    return theFreeState.position();
  }

bool BasicTrajectoryState::hasError ( void  ) const

inline

Definition at line 257 of file BasicTrajectoryState.h.

References FreeTrajectoryState::hasError(), theFreeState, theLocalError, and LocalTrajectoryError::valid().

Referenced by cartesianError(), curvilinearError(), freeTrajectoryState(), localError(), and rescaleError().

                        {
    return theFreeState.hasError() || theLocalError.valid();
  }

bool BasicTrajectoryState::isValid ( void  ) const

inline

Definition at line 161 of file BasicTrajectoryState.h.

References theValid.

Referenced by freeTrajectoryState(), and localParameters().

{ return theValid; }

LocalVector BasicTrajectoryState::localDirection ( ) const

inline

Definition at line 229 of file BasicTrajectoryState.h.

References localMomentum(), and Vector3DBase< T, FrameTag >::unit().

                                     {
    return localMomentum().unit();
  }

const LocalTrajectoryError& BasicTrajectoryState::localError ( ) const

inline

Definition at line 233 of file BasicTrajectoryState.h.

References createLocalError(), hasError(), LocalTrajectoryError::invalid(), missingError(), theLocalError, and unlikely.

                                                 {
    if unlikely(!hasError()) {
    missingError(" accessing local error.");
    return theLocalError;
      }
    if unlikely(theLocalError.invalid()) createLocalError();
    return theLocalError;
  }

LocalVector BasicTrajectoryState::localMomentum ( ) const

inline

Definition at line 226 of file BasicTrajectoryState.h.

References localParameters(), and LocalTrajectoryParameters::momentum().

Referenced by localDirection().

                                    {
    return localParameters().momentum();
  }

const LocalTrajectoryParameters& BasicTrajectoryState::localParameters ( ) const

inline

Definition at line 217 of file BasicTrajectoryState.h.

References createLocalParameters(), isValid(), notValid(), theLocalParameters, theLocalParametersValid, and unlikely.

Referenced by checkCurvilinError(), createLocalErrorFromCurvilinearError(), localMomentum(), and localPosition().

                                                           {
    if unlikely(!isValid()) notValid();
    if unlikely(!theLocalParametersValid)
      createLocalParameters();
    return theLocalParameters;
  }

LocalPoint BasicTrajectoryState::localPosition ( ) const

inline

Definition at line 223 of file BasicTrajectoryState.h.

References localParameters(), and LocalTrajectoryParameters::position().

                                   {
    return localParameters().position();
  }

const MagneticField* BasicTrajectoryState::magneticField ( ) const

inline

Definition at line 214 of file BasicTrajectoryState.h.

References GlobalTrajectoryParameters::magneticField(), FreeTrajectoryState::parameters(), and theFreeState.

Referenced by checkCurvilinError(), createLocalErrorFromCurvilinearError(), and rescaleError().

{ return &theFreeState.parameters().magneticField(); }

void BasicTrajectoryState::missingError ( char const *  where ) const

private

Definition at line 206 of file BasicTrajectoryState.cc.

References LocalTrajectoryError::matrix(), theFreeState, theLocalError, and LocalTrajectoryError::valid().

Referenced by cartesianError(), curvilinearError(), localError(), and rescaleError().

                                                               {
  std::stringstream form;
  form<<"BasicTrajectoryState: attempt to access errors when none available "
      <<where<<".\nfreestate pointer: " <<theFreeState
      <<"\nlocal error valid/values :"<< theLocalError.valid() << "\n" 
      <<  theLocalError.matrix();

  edm::LogWarning("BasicTrajectoryState") << form.str();

  // throw TrajectoryStateException(form.str());
}

void BasicTrajectoryState::notValid ( )

staticprivate

Definition at line 188 of file BasicTrajectoryState.cc.

Referenced by freeTrajectoryState(), and localParameters().

                                    {
  throw TrajectoryStateException("TrajectoryStateOnSurface is invalid and cannot return any parameters");
}

void BasicTrajectoryState::rescaleError

(

double

factor

)

Definition at line 309 of file BasicTrajectoryState.cc.

References benchmark_cfg::errors, hasError(), i, j, magneticField(), LocalTrajectoryError::matrix(), missingError(), MagneticField::nominalValue(), FreeTrajectoryState::rescaleError(), theFreeState, theLocalError, unlikely, and LocalTrajectoryError::valid().

                                                {
  if unlikely(!hasError()) missingError(" trying to rescale");    
  theFreeState.rescaleError(factor);
  
  if (theLocalError.valid()){
    //do it by hand if the free state is not around.
    bool zeroField = (magneticField()->nominalValue()==0);
    if unlikely(zeroField){
      AlgebraicSymMatrix55 errors=theLocalError.matrix();
      //scale the 0 indexed covariance by the square root of the factor
      for (unsigned int i=1;i!=5;++i)      errors(i,0)*=factor;
      double factor_squared=factor*factor;
      //scale all others by the scaled factor
      for (unsigned int i=1;i!=5;++i)  for (unsigned int j=i;j!=5;++j) errors(i,j)*=factor_squared;
      //term 0,0 is not scaled at all
      theLocalError = LocalTrajectoryError(errors);
    }
    else theLocalError *= (factor*factor);
  }
}

double BasicTrajectoryState::signedInverseMomentum ( ) const

inline

Definition at line 180 of file BasicTrajectoryState.h.

References FreeTrajectoryState::signedInverseMomentum(), and theFreeState.

                                       {
    return theFreeState.signedInverseMomentum();
  }

const SurfaceType& BasicTrajectoryState::surface ( ) const

inline

Definition at line 242 of file BasicTrajectoryState.h.

References theSurfaceP.

Referenced by checkCurvilinError(), createLocalErrorFromCurvilinearError(), and createLocalParameters().

                                     {
    return *theSurfaceP;
  }

SurfaceSide BasicTrajectoryState::surfaceSide ( ) const

inline

Position relative to material, defined relative to momentum vector.

Definition at line 253 of file BasicTrajectoryState.h.

References theSurfaceSide.

                                  {
    return theSurfaceSide;
  }

double BasicTrajectoryState::transverseCurvature ( ) const

inline

Definition at line 183 of file BasicTrajectoryState.h.

References theFreeState, and FreeTrajectoryState::transverseCurvature().

                                     {
    return theFreeState.transverseCurvature();
  }

void BasicTrajectoryState::update ( const LocalTrajectoryParameters &  p, const SurfaceType &  aSurface, const MagneticField *  field, const SurfaceSide  side  )

virtual

Reimplemented in GCC11_FINAL< T, TOPO >.

Definition at line 273 of file BasicTrajectoryState.cc.

References AlCaHLTBitMon_ParallelJobs::p, theFreeState, theLocalError, theLocalParameters, theLocalParametersValid, theSurfaceP, theSurfaceSide, theValid, and theWeight.

Referenced by progressbar.ProgressBar::__next__(), relval_steps.Matrix::__setitem__(), relval_steps.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), relval_steps.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

{
    theLocalParameters = p;
    if (&aSurface != &*theSurfaceP) theSurfaceP.reset(&aSurface);
    theSurfaceSide = side;
    theWeight      = 1.0; 
    theLocalError = InvalidError();
    theFreeState=makeFTS(p,aSurface,field);

    theValid   = true;
    theLocalParametersValid  = true;
}

void BasicTrajectoryState::update ( const LocalTrajectoryParameters &  p, const LocalTrajectoryError &  err, const SurfaceType &  aSurface, const MagneticField *  field, const SurfaceSide  side, double  weight  )

virtual

Reimplemented in GCC11_FINAL< T, TOPO >.

Definition at line 290 of file BasicTrajectoryState.cc.

References AlCaHLTBitMon_ParallelJobs::p, theFreeState, theLocalError, theLocalParameters, theLocalParametersValid, theSurfaceP, theSurfaceSide, theValid, theWeight, and weight().

Referenced by progressbar.ProgressBar::__next__(), relval_steps.Matrix::__setitem__(), relval_steps.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), relval_steps.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

{
    theLocalParameters = p;
    theLocalError      = err;
    if (&aSurface != &*theSurfaceP) theSurfaceP.reset(&aSurface);
    theSurfaceSide = side;
    theWeight      = weight; 
    theFreeState=makeFTS(p,aSurface,field);

    theValid   = true;
    theLocalParametersValid  = true;
}

double BasicTrajectoryState::weight ( void  ) const

inline

Definition at line 246 of file BasicTrajectoryState.h.

References theWeight.

Referenced by cuy.ValElement::__init__(), cuy.additionElement::__init__(), cuy.superimposeElement::__init__(), cuy.graphElement::__init__(), and update().

{return theWeight;} 

Friends And Related Function Documentation

friend class ProxyBase< BTSOS, CopyUsingClone< BTSOS > >

friend

Definition at line 85 of file BasicTrajectoryState.h.

friend class ReferenceCountingPointer< BasicTrajectoryState >

friend

Definition at line 86 of file BasicTrajectoryState.h.

Member Data Documentation

FreeTrajectoryState BasicTrajectoryState::theFreeState

mutableprivate

Definition at line 299 of file BasicTrajectoryState.h.

Referenced by charge(), checkCurvilinError(), createLocalError(), createLocalErrorFromCurvilinearError(), createLocalParameters(), freeTrajectoryState(), globalDirection(), globalMomentum(), globalParameters(), globalPosition(), hasError(), magneticField(), missingError(), rescaleError(), signedInverseMomentum(), transverseCurvature(), and update().

LocalTrajectoryError BasicTrajectoryState::theLocalError

mutableprivate

Definition at line 301 of file BasicTrajectoryState.h.

Referenced by checkCurvilinError(), createLocalError(), createLocalErrorFromCurvilinearError(), hasError(), localError(), missingError(), rescaleError(), and update().

LocalTrajectoryParameters BasicTrajectoryState::theLocalParameters

mutableprivate

Definition at line 302 of file BasicTrajectoryState.h.

Referenced by createLocalParameters(), localParameters(), and update().

bool BasicTrajectoryState::theLocalParametersValid

mutableprivate

Definition at line 304 of file BasicTrajectoryState.h.

Referenced by checkCurvilinError(), createLocalParameters(), localParameters(), and update().

ConstReferenceCountingPointer<SurfaceType> BasicTrajectoryState::theSurfaceP

private

Definition at line 309 of file BasicTrajectoryState.h.

Referenced by surface(), and update().

SurfaceSide BasicTrajectoryState::theSurfaceSide

private

Definition at line 308 of file BasicTrajectoryState.h.

Referenced by surfaceSide(), and update().

bool BasicTrajectoryState::theValid

mutableprivate

Definition at line 305 of file BasicTrajectoryState.h.

Referenced by isValid(), and update().

double BasicTrajectoryState::theWeight

private

Definition at line 311 of file BasicTrajectoryState.h.

Referenced by update(), and weight().