BasicTrajectoryState Class Reference

#include <BasicTrajectoryState.h>

Inheritance diagram for BasicTrajectoryState:

Public Types
typedef BasicTrajectoryState	BTSOS
typedef Proxy::pointer	pointer
typedef ProxyBase11< BTSOS >	Proxy
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 LocalTrajectoryParameters &par, const LocalTrajectoryError &err, const SurfaceType &aSurface, const MagneticField *field, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
	BasicTrajectoryState (const LocalTrajectoryParameters &par, const SurfaceType &aSurface, const MagneticField *field, 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)
	BasicTrajectoryState (const GlobalTrajectoryParameters &par, const SurfaceType &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
template<typename... Args>
	BasicTrajectoryState (double iweight, Args &&...args)
virtual bool	canUpdateLocalParameters () const
const CartesianTrajectoryError	cartesianError () const
TrackCharge	charge () const
virtual pointer	clone () const =0
virtual std::vector < TrajectoryStateOnSurface >	components () const =0
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)
CurvilinearTrajectoryError &	setCurvilinearError ()
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 SurfaceSide side) GCC11_FINAL
virtual void	update (double weight, const LocalTrajectoryParameters &p, const LocalTrajectoryError &err, const SurfaceType &aSurface, const MagneticField *field, const SurfaceSide side)
virtual void	update (const LocalTrajectoryParameters &p, const LocalTrajectoryError &err, const SurfaceSide side) GCC11_FINAL
double	weight () const
virtual	~BasicTrajectoryState ()

Static Public Member Functions
template<typename T , typename... Args>
static std::shared_ptr< BTSOS >	build (Args &&...args)
template<typename T , typename... Args>
static std::shared_ptr< BTSOS >	churn (Args &&...args)

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 =0.

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 70 of file BasicTrajectoryState.h.

Member Typedef Documentation

typedef BasicTrajectoryState BasicTrajectoryState::BTSOS

Definition at line 73 of file BasicTrajectoryState.h.

typedef Proxy::pointer BasicTrajectoryState::pointer

Definition at line 75 of file BasicTrajectoryState.h.

typedef ProxyBase11<BTSOS> BasicTrajectoryState::Proxy

Definition at line 74 of file BasicTrajectoryState.h.

typedef SurfaceSideDefinition::SurfaceSide BasicTrajectoryState::SurfaceSide

Definition at line 76 of file BasicTrajectoryState.h.

typedef Surface BasicTrajectoryState::SurfaceType

Definition at line 77 of file BasicTrajectoryState.h.

Constructor & Destructor Documentation

BasicTrajectoryState::BasicTrajectoryState ( )

inline

Definition at line 82 of file BasicTrajectoryState.h.

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

BasicTrajectoryState::BasicTrajectoryState ( const SurfaceType &  aSurface )

explicit

construct invalid trajectory state (without parameters)

Definition at line 42 of file BasicTrajectoryState.cc.

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

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 78 of file BasicTrajectoryState.cc.

                                      :
  theFreeState(fts),
  theLocalError(InvalidError()),
  theLocalParameters(),
  theLocalParametersValid(false),
  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
	)

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 110 of file BasicTrajectoryState.cc.

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

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

inline

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

Definition at line 124 of file BasicTrajectoryState.h.

                                                                             :
    BasicTrajectoryState(par, InvalidError(), aSurface, field, side){}

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 93 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  )

inline

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 145 of file BasicTrajectoryState.h.

                                                                             :
  theFreeState(par, err),
  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  )

inline

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

Definition at line 163 of file BasicTrajectoryState.h.

                                                                             :
    BasicTrajectoryState(par, InvalidError(), aSurface, side){}

template<typename... Args>

BasicTrajectoryState::BasicTrajectoryState ( double  iweight, Args &&...  args  )

inline

Definition at line 171 of file BasicTrajectoryState.h.

References theWeight.

                                                        : BasicTrajectoryState(std::forward<Args>(args)...) {
    theWeight = iweight;
  }

Member Function Documentation

template<typename T , typename... Args>

static std::shared_ptr<BTSOS> BasicTrajectoryState::build ( Args &&...  args )

inlinestatic

Definition at line 94 of file BasicTrajectoryState.h.

References ExtractAppInfoFromXML::args.

{ return std::make_shared<T>(std::forward<Args>(args)...);}

virtual bool BasicTrajectoryState::canUpdateLocalParameters ( ) const

inlinevirtual

Reimplemented in GCC11_FINAL< T, TOPO >.

Definition at line 305 of file BasicTrajectoryState.h.

{ return true; }

const CartesianTrajectoryError BasicTrajectoryState::cartesianError ( ) const

inline

Definition at line 230 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 220 of file BasicTrajectoryState.h.

References FreeTrajectoryState::charge(), and theFreeState.

                             {
    return theFreeState.charge();
  }

void BasicTrajectoryState::checkCurvilinError ( ) const

private

Definition at line 165 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,theFreeState);
  verifyCurvErr(cov,theFreeState); 
}

template<typename T , typename... Args>

static std::shared_ptr<BTSOS> BasicTrajectoryState::churn ( Args &&...  args )

inlinestatic

Definition at line 97 of file BasicTrajectoryState.h.

References ExtractAppInfoFromXML::args.

{ return std::allocate_shared<T>(churn_allocator<T>(),std::forward<Args>(args)...);}

virtual pointer 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 >, and GCC11_FINAL< T, TOPO >.

virtual std::vector<TrajectoryStateOnSurface> BasicTrajectoryState::components ( ) 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 >, and GCC11_FINAL< T, TOPO >.

void BasicTrajectoryState::createLocalError ( ) const

private

Definition at line 194 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 201 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();
  
  theLocalError = ROOT::Math::Similarity(jac, theFreeState.curvilinearError().matrix());

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

}

void BasicTrajectoryState::createLocalParameters ( ) const

private

Definition at line 183 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 237 of file BasicTrajectoryState.h.

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

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

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

inline

Definition at line 248 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 217 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 214 of file BasicTrajectoryState.h.

References FreeTrajectoryState::momentum(), and theFreeState.

                                      {
    return theFreeState.momentum();
  }

const GlobalTrajectoryParameters& BasicTrajectoryState::globalParameters ( void  ) const

inline

Definition at line 208 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 211 of file BasicTrajectoryState.h.

References FreeTrajectoryState::position(), and theFreeState.

                                     {
    return theFreeState.position();
  }

bool BasicTrajectoryState::hasError ( void  ) const

inline

Definition at line 300 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 204 of file BasicTrajectoryState.h.

References theValid.

Referenced by freeTrajectoryState(), and localParameters().

{ return theValid; }

LocalVector BasicTrajectoryState::localDirection ( ) const

inline

Definition at line 272 of file BasicTrajectoryState.h.

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

                                     {
    return localMomentum().unit();
  }

const LocalTrajectoryError& BasicTrajectoryState::localError ( ) const

inline

Definition at line 276 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 269 of file BasicTrajectoryState.h.

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

Referenced by localDirection().

                                    {
    return localParameters().momentum();
  }

const LocalTrajectoryParameters& BasicTrajectoryState::localParameters ( ) const

inline

Definition at line 260 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 266 of file BasicTrajectoryState.h.

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

                                   {
    return localParameters().position();
  }

const MagneticField* BasicTrajectoryState::magneticField ( ) const

inline

Definition at line 257 of file BasicTrajectoryState.h.

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

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

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

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

private

Definition at line 151 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 129 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 282 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);
  }
}

CurvilinearTrajectoryError& BasicTrajectoryState::setCurvilinearError ( )

inline

Definition at line 329 of file BasicTrajectoryState.h.

References FreeTrajectoryState::setCurvilinearError(), and theFreeState.

                                                    {
    return theFreeState.setCurvilinearError();
 }

double BasicTrajectoryState::signedInverseMomentum ( ) const

inline

Definition at line 223 of file BasicTrajectoryState.h.

References FreeTrajectoryState::signedInverseMomentum(), and theFreeState.

                                       {
    return theFreeState.signedInverseMomentum();
  }

const SurfaceType& BasicTrajectoryState::surface ( ) const

inline

Definition at line 285 of file BasicTrajectoryState.h.

References theSurfaceP.

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

                                     {
    return *theSurfaceP;
  }

SurfaceSide BasicTrajectoryState::surfaceSide ( ) const

inline

Position relative to material, defined relative to momentum vector.

Definition at line 296 of file BasicTrajectoryState.h.

References theSurfaceSide.

                                  {
    return theSurfaceSide;
  }

double BasicTrajectoryState::transverseCurvature ( ) const

inline

Definition at line 226 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 230 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 SurfaceSide  side  )

virtual

Definition at line 217 of file BasicTrajectoryState.cc.

References magneticField(), GlobalTrajectoryParameters::magneticFieldInTesla(), AlCaHLTBitMon_ParallelJobs::p, FreeTrajectoryState::parameters(), surface(), theFreeState, theLocalError, theLocalParameters, theLocalParametersValid, theSurfaceSide, and theValid.

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;
   theSurfaceSide = side;
   theLocalError = InvalidError();
   theFreeState=makeFTS(p,surface(),magneticField(), theFreeState.parameters().magneticFieldInTesla());
  
   theValid   = true;
   theLocalParametersValid  = true;

}

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

virtual

Reimplemented in GCC11_FINAL< T, TOPO >.

Definition at line 247 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;
}

void BasicTrajectoryState::update ( const LocalTrajectoryParameters &  p, const LocalTrajectoryError &  err, const SurfaceSide  side  )

virtual

Definition at line 267 of file BasicTrajectoryState.cc.

References magneticField(), GlobalTrajectoryParameters::magneticFieldInTesla(), AlCaHLTBitMon_ParallelJobs::p, FreeTrajectoryState::parameters(), surface(), theFreeState, theLocalError, theLocalParameters, theLocalParametersValid, theSurfaceSide, and theValid.

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;
    theSurfaceSide = side;
    theFreeState=   theFreeState=makeFTS(p,surface(),magneticField(), theFreeState.parameters().magneticFieldInTesla());

    theValid   = true;
    theLocalParametersValid  = true;
}

double BasicTrajectoryState::weight ( void  ) const

inline

Definition at line 289 of file BasicTrajectoryState.h.

References theWeight.

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

{return theWeight;} 

Member Data Documentation

FreeTrajectoryState BasicTrajectoryState::theFreeState

mutableprivate

Definition at line 354 of file BasicTrajectoryState.h.

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

LocalTrajectoryError BasicTrajectoryState::theLocalError

mutableprivate

Definition at line 356 of file BasicTrajectoryState.h.

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

LocalTrajectoryParameters BasicTrajectoryState::theLocalParameters

mutableprivate

Definition at line 357 of file BasicTrajectoryState.h.

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

bool BasicTrajectoryState::theLocalParametersValid

mutableprivate

Definition at line 359 of file BasicTrajectoryState.h.

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

ConstReferenceCountingPointer<SurfaceType> BasicTrajectoryState::theSurfaceP

private

Definition at line 364 of file BasicTrajectoryState.h.

Referenced by surface(), and update().

SurfaceSide BasicTrajectoryState::theSurfaceSide

private

Definition at line 363 of file BasicTrajectoryState.h.

Referenced by surfaceSide(), and update().

bool BasicTrajectoryState::theValid

mutableprivate

Definition at line 360 of file BasicTrajectoryState.h.

Referenced by isValid(), and update().

double BasicTrajectoryState::theWeight =0.

private

Definition at line 369 of file BasicTrajectoryState.h.

Referenced by BasicTrajectoryState(), update(), and weight().