#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)
	BasicTrajectoryState (const FreeTrajectoryState &fts, const SurfaceType &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
	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 GlobalTrajectoryParameters &par, 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, double weight=1.)
	BasicTrajectoryState (const LocalTrajectoryParameters &par, const LocalTrajectoryError &err, const SurfaceType &aSurface, const MagneticField *field, double weight)
	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.)
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 *	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.
double	transverseCurvature () const
virtual void	update (const LocalTrajectoryParameters &p, const LocalTrajectoryError &err, const SurfaceType &aSurface, const MagneticField *field, const SurfaceSide side, double weight)
virtual void	update (const LocalTrajectoryParameters &p, const SurfaceType &aSurface, const MagneticField *field, const SurfaceSide side)
double	weight () const
virtual	~BasicTrajectoryState ()
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 79 of file BasicTrajectoryState.h.

Member Typedef Documentation

typedef BasicTrajectoryState BasicTrajectoryState::BTSOS

Definition at line 82 of file BasicTrajectoryState.h.

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

Definition at line 83 of file BasicTrajectoryState.h.

typedef ReferenceCountingPointer<BasicTrajectoryState> BasicTrajectoryState::RCPtr

Definition at line 84 of file BasicTrajectoryState.h.

typedef SurfaceSideDefinition::SurfaceSide BasicTrajectoryState::SurfaceSide

Definition at line 85 of file BasicTrajectoryState.h.

typedef Surface BasicTrajectoryState::SurfaceType

Definition at line 86 of file BasicTrajectoryState.h.

Constructor & Destructor Documentation

BasicTrajectoryState::BasicTrajectoryState ( ) [inline]

Definition at line 94 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 [inline, virtual]

Reimplemented in GCC11_FINAL< T, TOPO >.

Definition at line 268 of file BasicTrajectoryState.h.

{ return true; }

const CartesianTrajectoryError BasicTrajectoryState::cartesianError ( ) const [inline]

Definition at line 191 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]

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

Definition at line 181 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]

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 >, 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 198 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* BasicTrajectoryState::freeTrajectoryState ( bool withErrors = true ) const [inline]

Definition at line 209 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 178 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 175 of file BasicTrajectoryState.h.

References FreeTrajectoryState::momentum(), and theFreeState.

                                      {
    return theFreeState.momentum();
  }

const GlobalTrajectoryParameters& BasicTrajectoryState::globalParameters ( void ) const [inline]

Definition at line 169 of file BasicTrajectoryState.h.

References FreeTrajectoryState::parameters(), and theFreeState.

Referenced by checkCurvilinError(), and createLocalErrorFromCurvilinearError().

                                                             {
    return theFreeState.parameters();
  }

GlobalPoint BasicTrajectoryState::globalPosition ( ) const [inline]

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

Definition at line 172 of file BasicTrajectoryState.h.

References FreeTrajectoryState::position(), and theFreeState.

                                     {
    return theFreeState.position();
  }

bool BasicTrajectoryState::hasError ( void ) const [inline]

Reimplemented in GCC11_FINAL< T, TOPO >.

Definition at line 261 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]

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

Definition at line 165 of file BasicTrajectoryState.h.

References theValid.

Referenced by freeTrajectoryState(), and localParameters().

{ return theValid; }

LocalVector BasicTrajectoryState::localDirection ( ) const [inline]

Reimplemented in GCC11_FINAL< T, TOPO >.

Definition at line 233 of file BasicTrajectoryState.h.

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

                                     {
    return localMomentum().unit();
  }

const LocalTrajectoryError& BasicTrajectoryState::localError ( ) const [inline]

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

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

Referenced by localDirection().

                                    {
    return localParameters().momentum();
  }

const LocalTrajectoryParameters& BasicTrajectoryState::localParameters ( ) const [inline]

Definition at line 221 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]

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

Definition at line 227 of file BasicTrajectoryState.h.

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

                                   {
    return localParameters().position();
  }

const MagneticField* BasicTrajectoryState::magneticField ( ) const [inline]

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

Definition at line 218 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 ( ) [static, private]

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 )

Reimplemented in GCC11_FINAL< T, TOPO >.

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

References FreeTrajectoryState::signedInverseMomentum(), and theFreeState.

                                       {
    return theFreeState.signedInverseMomentum();
  }

const SurfaceType& BasicTrajectoryState::surface ( ) const [inline]

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

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

References theSurfaceSide.

                                  {
    return theSurfaceSide;
  }

double BasicTrajectoryState::transverseCurvature ( ) const [inline]

Definition at line 187 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.

{
    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().

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

    theValid   = true;
    theLocalParametersValid  = true;
}