BasicSingleTrajectoryState Class Reference

vvv DEBUG More...

#include <BasicSingleTrajectoryState.h>

Inheritance diagram for BasicSingleTrajectoryState:

Public Member Functions
	BasicSingleTrajectoryState (const FreeTrajectoryState &fts, const Surface &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
	BasicSingleTrajectoryState (const GlobalTrajectoryParameters &par, const Surface &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
	BasicSingleTrajectoryState (const GlobalTrajectoryParameters &par, const CartesianTrajectoryError &err, const Surface &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
	BasicSingleTrajectoryState (const GlobalTrajectoryParameters &par, const CurvilinearTrajectoryError &err, const Surface &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface, double weight=1.)
	BasicSingleTrajectoryState (const GlobalTrajectoryParameters &par, const CurvilinearTrajectoryError &err, const Surface &aSurface, double weight)
	BasicSingleTrajectoryState (const LocalTrajectoryParameters &par, const Surface &aSurface, const MagneticField *field, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
	BasicSingleTrajectoryState (const LocalTrajectoryParameters &par, const LocalTrajectoryError &err, const Surface &aSurface, const MagneticField *field, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface, double weight=1.)
	BasicSingleTrajectoryState (const LocalTrajectoryParameters &par, const LocalTrajectoryError &err, const Surface &aSurface, const MagneticField *field, double weight)
	BasicSingleTrajectoryState (const Surface &aSurface)
	construct invalid trajectory state (without parameters) More...
virtual bool	canUpdateLocalParameters () const
const CartesianTrajectoryError &	cartesianError () const
TrackCharge	charge () const
BasicSingleTrajectoryState *	clone () 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	missingError (char const *where) const
void	rescaleError (double factor)
double	signedInverseMomentum () const
const Surface &	surface () const
virtual SurfaceSide	surfaceSide () const
	Position relative to material, defined relative to momentum vector. More...
double	transverseCurvature () const
virtual void	update (const LocalTrajectoryParameters &p, const Surface &aSurface, const MagneticField *field, const SurfaceSide side)
virtual void	update (const LocalTrajectoryParameters &p, const LocalTrajectoryError &err, const Surface &aSurface, const MagneticField *field, const SurfaceSide side, double weight)
virtual double	weight () const
virtual	~BasicSingleTrajectoryState ()
Public Member Functions inherited from BasicTrajectoryState
virtual std::vector < TrajectoryStateOnSurface >	components () const
virtual	~BasicTrajectoryState ()
Public Member Functions inherited from ReferenceCountedPoolAllocated
void	addReference () const
const ReferenceCountedPoolAllocated &	operator= (const ReferenceCountedPoolAllocated &)
	ReferenceCountedPoolAllocated ()
	ReferenceCountedPoolAllocated (const ReferenceCountedPoolAllocated &iRHS)
unsigned int	references () const
void	removeReference () const
virtual	~ReferenceCountedPoolAllocated ()
Public Member Functions inherited from BlockWipedPoolAllocated
virtual	~BlockWipedPoolAllocated ()

Static Public Member Functions
static void	notValid ()
Static Public Member Functions inherited from BlockWipedPoolAllocated
static BlockWipedAllocator &	allocator (size_t s)
static void	operator delete (void *p, size_t s)
static void *	operator new (size_t s, void *p)
static void *	operator new (size_t s)
static BlockWipedAllocator::Stat	stat (size_t s)
static void	usePool ()

Private Member Functions
void	checkCartesianError () const
void	checkCurvilinError () const
void	checkGlobalParameters () const
void	createLocalError () const
void	createLocalErrorFromCartesianError () const
void	createLocalErrorFromCurvilinearError () const
void	createLocalParameters () const

Private Attributes
bool	theCartesianErrorUp2Date
bool	theCurvilinErrorUp2Date
const MagneticField *	theField
DeepCopyPointer < FreeTrajectoryState >	theFreeState
bool	theGlobalParamsUp2Date
LocalTrajectoryError	theLocalError
bool	theLocalErrorValid
LocalTrajectoryParameters	theLocalParameters
bool	theLocalParametersValid
ConstReferenceCountingPointer < Surface >	theSurfaceP
SurfaceSide	theSurfaceSide
double	theWeight

Additional Inherited Members
Public Types inherited from BasicTrajectoryState
typedef BasicTrajectoryState	BTSOS
typedef ProxyBase< BTSOS, CopyUsingClone< BTSOS > >	Proxy
typedef ReferenceCountingPointer < BasicTrajectoryState >	RCPtr
typedef SurfaceSideDefinition::SurfaceSide	SurfaceSide
Static Public Attributes inherited from ReferenceCountedPoolAllocated
static int	s_alive =0
static int	s_referenced =0
Static Public Attributes inherited from BlockWipedPoolAllocated
static int	s_alive =0

Detailed Description

vvv DEBUG

Concrete implementation for the state of one trajectory on a surface.

Definition at line 24 of file BasicSingleTrajectoryState.h.

Constructor & Destructor Documentation

BasicSingleTrajectoryState::BasicSingleTrajectoryState	(	const FreeTrajectoryState &	fts,
		const Surface &	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 11 of file BasicSingleTrajectoryState.cc.

Referenced by clone().

                                        :
  theFreeState( new FreeTrajectoryState(fts)),
  theLocalError(),
  theLocalParameters(),
  theLocalParametersValid(false),
  theLocalErrorValid(false),
  theGlobalParamsUp2Date(true),
  theCartesianErrorUp2Date(fts.hasCartesianError()),
  theCurvilinErrorUp2Date(fts.hasCurvilinearError()),
  theSurfaceSide(side), 
  theSurfaceP( &aSurface), 
  theWeight(1.),
  theField( &fts.parameters().magneticField())
{}    

BasicSingleTrajectoryState::BasicSingleTrajectoryState	(	const GlobalTrajectoryParameters &	par,
		const Surface &	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 29 of file BasicSingleTrajectoryState.cc.

                                        :
  theFreeState( new FreeTrajectoryState(par)),
  theLocalError(),
  theLocalParameters(),
  theLocalParametersValid(false),
  theLocalErrorValid(false),
  theGlobalParamsUp2Date(true),
  theCartesianErrorUp2Date(false),
  theCurvilinErrorUp2Date(false),
  theSurfaceSide(side), 
  theSurfaceP( &aSurface), 
  theWeight(1.),
  theField( &par.magneticField())
{}

BasicSingleTrajectoryState::BasicSingleTrajectoryState	(	const GlobalTrajectoryParameters &	par,
		const CartesianTrajectoryError &	err,
		const Surface &	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 47 of file BasicSingleTrajectoryState.cc.

                                        :
  theFreeState( new FreeTrajectoryState(par, err)),
  theLocalError(),
  theLocalParameters(),
  theLocalParametersValid(false),
  theLocalErrorValid(false),
  theGlobalParamsUp2Date(true),
  theCartesianErrorUp2Date(true),
  theCurvilinErrorUp2Date(false),
  theSurfaceSide(side), 
  theSurfaceP( &aSurface), 
  theWeight(1.),
  theField( &par.magneticField())
{}

BasicSingleTrajectoryState::BasicSingleTrajectoryState	(	const GlobalTrajectoryParameters &	par,
		const CurvilinearTrajectoryError &	err,
		const Surface &	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 66 of file BasicSingleTrajectoryState.cc.

                               :
  theFreeState( new FreeTrajectoryState(par, err)),
  theLocalError(),
  theLocalParameters(),
  theLocalParametersValid(false),
  theLocalErrorValid(false),
  theGlobalParamsUp2Date(true),
  theCartesianErrorUp2Date(false),
  theCurvilinErrorUp2Date(true),
  theSurfaceSide(side), 
  theSurfaceP( &aSurface), 
  theWeight(weight),
  theField( &par.magneticField())
{}

BasicSingleTrajectoryState::BasicSingleTrajectoryState	(	const GlobalTrajectoryParameters &	par,
		const CurvilinearTrajectoryError &	err,
		const Surface &	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 86 of file BasicSingleTrajectoryState.cc.

                               :
  theFreeState( new FreeTrajectoryState(par, err)),
  theLocalError(),
  theLocalParameters(),
  theLocalParametersValid(false),
  theLocalErrorValid(false),
  theGlobalParamsUp2Date(true),
  theCartesianErrorUp2Date(false),
  theCurvilinErrorUp2Date(true),
  theSurfaceSide(SurfaceSideDefinition::atCenterOfSurface), 
  theSurfaceP( &aSurface), 
  theWeight(weight),
  theField( &par.magneticField())
{}

BasicSingleTrajectoryState::BasicSingleTrajectoryState	(	const LocalTrajectoryParameters &	par,
		const Surface &	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 105 of file BasicSingleTrajectoryState.cc.

                                        :
  theFreeState(0),
  theLocalError(),
  theLocalParameters(par),
  theLocalParametersValid(true),
  theLocalErrorValid(false),
  theGlobalParamsUp2Date(false),
  theCartesianErrorUp2Date(false),
  theCurvilinErrorUp2Date(false),
  theSurfaceSide(side),
  theSurfaceP( &aSurface), 
  theWeight(1.),
  theField(field) 
{}

BasicSingleTrajectoryState::BasicSingleTrajectoryState	(	const LocalTrajectoryParameters &	par,
		const LocalTrajectoryError &	err,
		const Surface &	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 124 of file BasicSingleTrajectoryState.cc.

                               :
  theFreeState(0),
  theLocalError(err),
  theLocalParameters(par),
  theLocalParametersValid(true),
  theLocalErrorValid(true),
  theGlobalParamsUp2Date(false),
  theCartesianErrorUp2Date(false),
  theCurvilinErrorUp2Date(false),
  theSurfaceSide(side), 
  theSurfaceP( &aSurface),
  theWeight(weight),
  theField(field)
{}

BasicSingleTrajectoryState::BasicSingleTrajectoryState	(	const LocalTrajectoryParameters &	par,
		const LocalTrajectoryError &	err,
		const Surface &	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 145 of file BasicSingleTrajectoryState.cc.

                               :
  theFreeState(0),
  theLocalError(err),
  theLocalParameters(par),
  theLocalParametersValid(true),
  theLocalErrorValid(true),
  theGlobalParamsUp2Date(false),
  theCartesianErrorUp2Date(false),
  theCurvilinErrorUp2Date(false),
  theSurfaceSide(SurfaceSideDefinition::atCenterOfSurface),
  theSurfaceP( &aSurface), 
  theWeight(weight),
  theField(field)
{}

BasicSingleTrajectoryState::BasicSingleTrajectoryState

(

const Surface &

aSurface

)

construct invalid trajectory state (without parameters)

Definition at line 165 of file BasicSingleTrajectoryState.cc.

                                                    :
  theFreeState(0),
  theLocalError(),
  theLocalParameters(),
  theLocalParametersValid(false),
  theLocalErrorValid(false),
  theGlobalParamsUp2Date(false),
  theCartesianErrorUp2Date(false),
  theCurvilinErrorUp2Date(false),
  theSurfaceSide(SurfaceSideDefinition::atCenterOfSurface), 
  theSurfaceP( &aSurface), 
  theWeight(0.),
  theField(0)
{}

BasicSingleTrajectoryState::~BasicSingleTrajectoryState ( )

virtual

Definition at line 181 of file BasicSingleTrajectoryState.cc.

{}

Member Function Documentation

virtual bool BasicSingleTrajectoryState::canUpdateLocalParameters ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 182 of file BasicSingleTrajectoryState.h.

{ return true; }

const CartesianTrajectoryError& BasicSingleTrajectoryState::cartesianError ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 126 of file BasicSingleTrajectoryState.h.

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

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

TrackCharge BasicSingleTrajectoryState::charge ( void  ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 117 of file BasicSingleTrajectoryState.h.

References FreeTrajectoryState::charge(), and freeTrajectoryState().

Referenced by PixelTrackBuilder::build().

                             {
    return freeTrajectoryState(false)->charge();
  }

void BasicSingleTrajectoryState::checkCartesianError ( ) const

private

Definition at line 239 of file BasicSingleTrajectoryState.cc.

References createLocalError(), createLocalParameters(), JacobianLocalToCartesian::jacobian(), localParameters(), LocalTrajectoryError::matrix(), FreeTrajectoryState::setCartesianError(), surface(), theCartesianErrorUp2Date, theFreeState, theLocalError, theLocalErrorValid, and theLocalParametersValid.

Referenced by freeTrajectoryState().

                                                           {
  if(!theCartesianErrorUp2Date){
    if(!theLocalParametersValid) createLocalParameters();
    if(!theLocalErrorValid) createLocalError();
    theCartesianErrorUp2Date = true;
   
    JacobianLocalToCartesian loc2Cart(surface(), localParameters());
    const AlgebraicMatrix65& jac = loc2Cart.jacobian();

    const AlgebraicSymMatrix66 &cov = ROOT::Math::Similarity(jac, theLocalError.matrix());
    
    //theFreeState->setCartesianError( CartesianTrajectoryError(cov) );
    theFreeState->setCartesianError( cov );
  }
}

void BasicSingleTrajectoryState::checkCurvilinError ( ) const

private

Definition at line 221 of file BasicSingleTrajectoryState.cc.

References createLocalError(), createLocalParameters(), JacobianLocalToCurvilinear::jacobian(), localParameters(), LocalTrajectoryError::matrix(), FreeTrajectoryState::setCurvilinearError(), surface(), theCurvilinErrorUp2Date, theField, theFreeState, theLocalError, theLocalErrorValid, and theLocalParametersValid.

Referenced by freeTrajectoryState().

                                                          {
  if(!theCurvilinErrorUp2Date){
    if(!theLocalParametersValid) createLocalParameters();
    // after createLocalParameters we can be sure theFreeState is not null
    if(!theLocalErrorValid) createLocalError();
    //    cout<<"!theCurvilinErrorUp2Date: create curviError from local"<<endl;
    theCurvilinErrorUp2Date = true;
    
    JacobianLocalToCurvilinear loc2Curv(surface(), localParameters(), *theField);
    const AlgebraicMatrix55& jac = loc2Curv.jacobian();
    
    const AlgebraicSymMatrix55 &cov = ROOT::Math::Similarity(jac, theLocalError.matrix());

    //theFreeState->setCurvilinearError( CurvilinearTrajectoryError(cov) );
    theFreeState->setCurvilinearError( cov );
  }
}

void BasicSingleTrajectoryState::checkGlobalParameters ( ) const

private

Definition at line 199 of file BasicSingleTrajectoryState.cc.

References LocalTrajectoryParameters::charge(), LocalTrajectoryParameters::momentum(), L1TEmulatorMonitor_cff::p, LocalTrajectoryParameters::position(), DeepCopyPointer< T >::replaceWith(), surface(), theCartesianErrorUp2Date, theCurvilinErrorUp2Date, theField, theFreeState, theGlobalParamsUp2Date, theLocalParameters, Surface::toGlobal(), and x.

Referenced by freeTrajectoryState().

                                                             {
  if(!theGlobalParamsUp2Date){
    //    cout<<"!theGlobalParamsUp2Date"<<endl;
    theGlobalParamsUp2Date = true;
    theCurvilinErrorUp2Date = false;
    theCartesianErrorUp2Date = false;
    // calculate global parameters from local
    GlobalPoint  x = surface().toGlobal(theLocalParameters.position());
    GlobalVector p = surface().toGlobal(theLocalParameters.momentum());
    // replace in place
    FreeTrajectoryState * fts = &(*theFreeState);
    if (fts) { 
        fts->~FreeTrajectoryState();
        new(fts) FreeTrajectoryState(x, p, theLocalParameters.charge(), theField);
    }else {
      theFreeState.replaceWith(new FreeTrajectoryState(x, p,
                                                       theLocalParameters.charge(),
                                                       theField));
    } 
  }
}

BasicSingleTrajectoryState* BasicSingleTrajectoryState::clone ( void  ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 173 of file BasicSingleTrajectoryState.h.

References BasicSingleTrajectoryState().

                                            {
    return new BasicSingleTrajectoryState(*this);
  }

void BasicSingleTrajectoryState::createLocalError ( ) const

private

Definition at line 267 of file BasicSingleTrajectoryState.cc.

References createLocalErrorFromCartesianError(), createLocalErrorFromCurvilinearError(), FreeTrajectoryState::hasCartesianError(), FreeTrajectoryState::hasCurvilinearError(), theFreeState, and theLocalErrorValid.

Referenced by checkCartesianError(), checkCurvilinError(), and localError().

                                                        {
  if(theFreeState->hasCartesianError())
    createLocalErrorFromCartesianError();
  else if(theFreeState->hasCurvilinearError())
    createLocalErrorFromCurvilinearError();
  else
    theLocalErrorValid = false;
}

void BasicSingleTrajectoryState::createLocalErrorFromCartesianError ( ) const

private

Definition at line 290 of file BasicSingleTrajectoryState.cc.

References funct::C, FreeTrajectoryState::cartesianError(), JacobianCartesianToLocal::jacobian(), localParameters(), CartesianTrajectoryError::matrix(), surface(), theFreeState, theLocalError, and theLocalErrorValid.

Referenced by createLocalError().

                                                                     {

  JacobianCartesianToLocal cart2Loc(surface(), localParameters());
  const AlgebraicMatrix56& jac = cart2Loc.jacobian();
    

  const AlgebraicSymMatrix55 &C = 
    ROOT::Math::Similarity(jac, theFreeState->cartesianError().matrix());
  theLocalError = LocalTrajectoryError(C);
  theLocalErrorValid = true;
}

void BasicSingleTrajectoryState::createLocalErrorFromCurvilinearError ( ) const

private

Definition at line 277 of file BasicSingleTrajectoryState.cc.

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

Referenced by createLocalError().

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

void BasicSingleTrajectoryState::createLocalParameters ( ) const

private

Definition at line 256 of file BasicSingleTrajectoryState.cc.

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

Referenced by checkCartesianError(), 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& BasicSingleTrajectoryState::curvilinearError ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 130 of file BasicSingleTrajectoryState.h.

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

Referenced by PixelTrackBuilder::build().

                                                             {
    if(!hasError()) missingError(" accesing curvilinearerror.");
    return freeTrajectoryState()->curvilinearError();
  }

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

virtual

Implements BasicTrajectoryState.

Definition at line 366 of file BasicSingleTrajectoryState.cc.

References checkCartesianError(), checkCurvilinError(), checkGlobalParameters(), hasError(), isValid(), and notValid().

Referenced by cartesianError(), charge(), curvilinearError(), globalDirection(), globalMomentum(), globalParameters(), globalPosition(), signedInverseMomentum(), and transverseCurvature().

                                                                     {
  if(!isValid()) notValid();
  checkGlobalParameters();
  //if(hasError()) { // let's start like this to see if we alloc less
  if(withErrors && hasError()) { // this is the right thing
    checkCartesianError();
    checkCurvilinError();
  }
  return &(*theFreeState);
}

GlobalVector BasicSingleTrajectoryState::globalDirection ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 114 of file BasicSingleTrajectoryState.h.

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

                                       {
    return freeTrajectoryState(false)->momentum().unit();
  }  

GlobalVector BasicSingleTrajectoryState::globalMomentum ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 111 of file BasicSingleTrajectoryState.h.

References freeTrajectoryState(), and FreeTrajectoryState::momentum().

Referenced by PixelTrackBuilder::build().

                                      {
    return freeTrajectoryState(false)->momentum();
  }

const GlobalTrajectoryParameters& BasicSingleTrajectoryState::globalParameters ( void  ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 105 of file BasicSingleTrajectoryState.h.

References freeTrajectoryState(), and FreeTrajectoryState::parameters().

                                                             {
    return freeTrajectoryState(false)->parameters();
  }

GlobalPoint BasicSingleTrajectoryState::globalPosition ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 108 of file BasicSingleTrajectoryState.h.

References freeTrajectoryState(), and FreeTrajectoryState::position().

Referenced by PixelTrackBuilder::build().

                                     {
    return freeTrajectoryState(false)->position();
  }

bool BasicSingleTrajectoryState::hasError ( void  ) const

virtual

Implements BasicTrajectoryState.

Definition at line 378 of file BasicSingleTrajectoryState.cc.

References FreeTrajectoryState::hasError(), theFreeState, and theLocalErrorValid.

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

                                           {
  return (theFreeState && theFreeState->hasError()) || theLocalErrorValid;
}

bool BasicSingleTrajectoryState::isValid ( void  ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 95 of file BasicSingleTrajectoryState.h.

References theFreeState, and theLocalParametersValid.

Referenced by freeTrajectoryState(), and localParameters().

                       {
    return theFreeState || theLocalParametersValid;
  }

LocalVector BasicSingleTrajectoryState::localDirection ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 153 of file BasicSingleTrajectoryState.h.

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

                                     {
    return localMomentum().unit();
  }

const LocalTrajectoryError& BasicSingleTrajectoryState::localError ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 157 of file BasicSingleTrajectoryState.h.

References createLocalError(), hasError(), missingError(), theLocalError, and theLocalErrorValid.

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

LocalVector BasicSingleTrajectoryState::localMomentum ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 150 of file BasicSingleTrajectoryState.h.

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

Referenced by localDirection().

                                    {
    return localParameters().momentum();
  }

const LocalTrajectoryParameters& BasicSingleTrajectoryState::localParameters ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 141 of file BasicSingleTrajectoryState.h.

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

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

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

LocalPoint BasicSingleTrajectoryState::localPosition ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 147 of file BasicSingleTrajectoryState.h.

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

                                   {
    return localParameters().position();
  }

const MagneticField* BasicSingleTrajectoryState::magneticField ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 138 of file BasicSingleTrajectoryState.h.

References theField.

{ return theField; }

void BasicSingleTrajectoryState::missingError

(

char const *

where

)

const

Definition at line 188 of file BasicSingleTrajectoryState.cc.

References theFreeState, and theLocalErrorValid.

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

                                                                     {
  std::stringstream form;
  form<<"TrajectoryStateOnSurface: attempt to access errors when none available "
      <<where<<".\nfreestate pointer: "
      <<theFreeState<<"\nlocal error valid :"<<theLocalErrorValid ;
  throw TrajectoryStateException(form.str());
}

void BasicSingleTrajectoryState::notValid ( )

static

Definition at line 183 of file BasicSingleTrajectoryState.cc.

Referenced by freeTrajectoryState(), and localParameters().

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

void BasicSingleTrajectoryState::rescaleError ( double  factor )

virtual

Implements BasicTrajectoryState.

Definition at line 343 of file BasicSingleTrajectoryState.cc.

References benchmark_cfg::errors, hasError(), i, MagneticField::inInverseGeV(), j, PV3DBase< T, PVType, FrameType >::mag2(), LocalTrajectoryError::matrix(), missingError(), FreeTrajectoryState::rescaleError(), theField, theFreeState, theLocalError, and theLocalErrorValid.

                                                      {
  if (!hasError()) missingError(" trying to rescale");    
  if (theFreeState)
    theFreeState->rescaleError(factor);
  
  if (theLocalErrorValid){
    //do it by hand if the free state is not around.
    bool zeroField =theField->inInverseGeV(GlobalPoint(0,0,0)).mag2()==0;
    if (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 BasicSingleTrajectoryState::signedInverseMomentum ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 120 of file BasicSingleTrajectoryState.h.

References freeTrajectoryState(), and FreeTrajectoryState::signedInverseMomentum().

                                       {
    return freeTrajectoryState(false)->signedInverseMomentum();
  }

const Surface& BasicSingleTrajectoryState::surface ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 164 of file BasicSingleTrajectoryState.h.

References theSurfaceP.

Referenced by checkCartesianError(), checkCurvilinError(), checkGlobalParameters(), createLocalErrorFromCartesianError(), createLocalErrorFromCurvilinearError(), and createLocalParameters().

                                 {
    return *theSurfaceP;
  }

virtual SurfaceSide BasicSingleTrajectoryState::surfaceSide ( ) const

inlinevirtual

Position relative to material, defined relative to momentum vector.

Implements BasicTrajectoryState.

Definition at line 178 of file BasicSingleTrajectoryState.h.

References theSurfaceSide.

                                          {
    return theSurfaceSide;
  }

double BasicSingleTrajectoryState::transverseCurvature ( ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 123 of file BasicSingleTrajectoryState.h.

References freeTrajectoryState(), and FreeTrajectoryState::transverseCurvature().

                                     {
    return freeTrajectoryState(false)->transverseCurvature();
  }

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

virtual

Implements BasicTrajectoryState.

Definition at line 303 of file BasicSingleTrajectoryState.cc.

References L1TEmulatorMonitor_cff::p, theCartesianErrorUp2Date, theCurvilinErrorUp2Date, theGlobalParamsUp2Date, theLocalErrorValid, theLocalParameters, theLocalParametersValid, theSurfaceP, theSurfaceSide, and theWeight.

Referenced by python.Vispa.Gui.VispaWidget.VispaWidget::autosize(), python.Vispa.Views.LineDecayView.LineDecayContainer::createObject(), python.Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), python.Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), python.Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), python.Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), python.Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), python.Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), python.Vispa.Gui.FindDialog.FindDialog::reset(), python.Vispa.Gui.PortConnection.PointToPointConnection::select(), python.Vispa.Gui.VispaWidget.VispaWidget::select(), python.Vispa.Views.LineDecayView.LineDecayContainer::select(), python.Vispa.Gui.VispaWidget.VispaWidget::setText(), python.Vispa.Gui.VispaWidget.VispaWidget::setTitle(), python.Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), python.Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and python.Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

{
    theLocalParameters = p;
    if (&aSurface != &*theSurfaceP) theSurfaceP.reset(&aSurface);
    theSurfaceSide = side;
    theWeight      = 1.0; 

    theGlobalParamsUp2Date   = false;
    theCartesianErrorUp2Date = false;
    theCurvilinErrorUp2Date  = false; 
    theLocalErrorValid       = false;
    theLocalParametersValid  = true;
}

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

virtual

Implements BasicTrajectoryState.

Definition at line 321 of file BasicSingleTrajectoryState.cc.

References L1TEmulatorMonitor_cff::p, theCartesianErrorUp2Date, theCurvilinErrorUp2Date, theGlobalParamsUp2Date, theLocalError, theLocalErrorValid, theLocalParameters, theLocalParametersValid, theSurfaceP, theSurfaceSide, theWeight, and weight().

Referenced by python.Vispa.Gui.VispaWidget.VispaWidget::autosize(), python.Vispa.Views.LineDecayView.LineDecayContainer::createObject(), python.Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), python.Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), python.Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), python.Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), python.Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), python.Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), python.Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), python.Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), python.Vispa.Gui.FindDialog.FindDialog::reset(), python.Vispa.Gui.PortConnection.PointToPointConnection::select(), python.Vispa.Gui.VispaWidget.VispaWidget::select(), python.Vispa.Views.LineDecayView.LineDecayContainer::select(), python.Vispa.Gui.VispaWidget.VispaWidget::setText(), python.Vispa.Gui.VispaWidget.VispaWidget::setTitle(), python.Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), python.Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and python.Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

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

    theGlobalParamsUp2Date   = false;
    theCartesianErrorUp2Date = false;
    theCurvilinErrorUp2Date  = false; 
    theLocalErrorValid       = true;
    theLocalParametersValid  = true;

}

virtual double BasicSingleTrajectoryState::weight ( ) const

inlinevirtual

Reimplemented from BasicTrajectoryState.

Definition at line 168 of file BasicSingleTrajectoryState.h.

References theWeight.

Referenced by update().

{return theWeight;} 

Member Data Documentation

bool BasicSingleTrajectoryState::theCartesianErrorUp2Date

mutableprivate

Definition at line 217 of file BasicSingleTrajectoryState.h.

Referenced by checkCartesianError(), checkGlobalParameters(), and update().

bool BasicSingleTrajectoryState::theCurvilinErrorUp2Date

mutableprivate

Definition at line 218 of file BasicSingleTrajectoryState.h.

Referenced by checkCurvilinError(), checkGlobalParameters(), and update().

const MagneticField* BasicSingleTrajectoryState::theField

private

Definition at line 225 of file BasicSingleTrajectoryState.h.

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

DeepCopyPointer<FreeTrajectoryState> BasicSingleTrajectoryState::theFreeState

mutableprivate

Definition at line 209 of file BasicSingleTrajectoryState.h.

Referenced by checkCartesianError(), checkCurvilinError(), checkGlobalParameters(), createLocalError(), createLocalErrorFromCartesianError(), createLocalErrorFromCurvilinearError(), createLocalParameters(), hasError(), isValid(), missingError(), and rescaleError().

bool BasicSingleTrajectoryState::theGlobalParamsUp2Date

mutableprivate

Definition at line 216 of file BasicSingleTrajectoryState.h.

Referenced by checkGlobalParameters(), and update().

LocalTrajectoryError BasicSingleTrajectoryState::theLocalError

mutableprivate

Definition at line 211 of file BasicSingleTrajectoryState.h.

Referenced by checkCartesianError(), checkCurvilinError(), createLocalErrorFromCartesianError(), createLocalErrorFromCurvilinearError(), localError(), rescaleError(), and update().

bool BasicSingleTrajectoryState::theLocalErrorValid

mutableprivate

Definition at line 214 of file BasicSingleTrajectoryState.h.

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

LocalTrajectoryParameters BasicSingleTrajectoryState::theLocalParameters

mutableprivate

Definition at line 212 of file BasicSingleTrajectoryState.h.

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

bool BasicSingleTrajectoryState::theLocalParametersValid

mutableprivate

Definition at line 213 of file BasicSingleTrajectoryState.h.

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

ConstReferenceCountingPointer<Surface> BasicSingleTrajectoryState::theSurfaceP

private

Definition at line 222 of file BasicSingleTrajectoryState.h.

Referenced by surface(), and update().

SurfaceSide BasicSingleTrajectoryState::theSurfaceSide

private

Definition at line 221 of file BasicSingleTrajectoryState.h.

Referenced by surfaceSide(), and update().

double BasicSingleTrajectoryState::theWeight

private

Definition at line 224 of file BasicSingleTrajectoryState.h.

Referenced by update(), and weight().