#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
Public Member Functions
	BasicTrajectoryState ()
	BasicTrajectoryState (const Surface &aSurface)
	construct invalid trajectory state (without parameters)
	BasicTrajectoryState (const FreeTrajectoryState &fts, const Surface &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
	BasicTrajectoryState (const GlobalTrajectoryParameters &par, const CurvilinearTrajectoryError &err, const Surface &aSurface, double weight)
	BasicTrajectoryState (const LocalTrajectoryParameters &par, const Surface &aSurface, const MagneticField *field, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
	BasicTrajectoryState (const GlobalTrajectoryParameters &par, const Surface &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
	BasicTrajectoryState (const LocalTrajectoryParameters &par, const LocalTrajectoryError &err, const Surface &aSurface, const MagneticField *field, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface, double weight=1.)
	BasicTrajectoryState (const LocalTrajectoryParameters &par, const LocalTrajectoryError &err, const Surface &aSurface, const MagneticField *field, double weight)
	BasicTrajectoryState (const GlobalTrajectoryParameters &par, const CartesianTrajectoryError &err, const Surface &aSurface, const SurfaceSide side=SurfaceSideDefinition::atCenterOfSurface)
	BasicTrajectoryState (const GlobalTrajectoryParameters &par, const CurvilinearTrajectoryError &err, const Surface &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 Surface &	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 Surface &aSurface, const MagneticField *field, const SurfaceSide side, double weight)
virtual void	update (const LocalTrajectoryParameters &p, const Surface &aSurface, const MagneticField *field, const SurfaceSide side)
double	weight () const
virtual	~BasicTrajectoryState ()
Private Member Functions
void	checkCurvilinError () const dso_internal
void	checkGlobalParameters () const dso_internal
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
const MagneticField *	theField
DeepCopyPointer < FreeTrajectoryState >	theFreeState
bool	theGlobalParamsUp2Date
LocalTrajectoryError	theLocalError
LocalTrajectoryParameters	theLocalParameters
bool	theLocalParametersValid
ConstReferenceCountingPointer < Surface >	theSurfaceP
SurfaceSide	theSurfaceSide
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 47 of file BasicTrajectoryState.h.

Member Typedef Documentation

typedef BasicTrajectoryState BasicTrajectoryState::BTSOS

Definition at line 50 of file BasicTrajectoryState.h.

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

Definition at line 51 of file BasicTrajectoryState.h.

typedef ReferenceCountingPointer<BasicTrajectoryState> BasicTrajectoryState::RCPtr

Definition at line 52 of file BasicTrajectoryState.h.

typedef SurfaceSideDefinition::SurfaceSide BasicTrajectoryState::SurfaceSide

Definition at line 53 of file BasicTrajectoryState.h.

Constructor & Destructor Documentation

BasicTrajectoryState::BasicTrajectoryState ( ) [inline]

Definition at line 61 of file BasicTrajectoryState.h.

{}

BasicTrajectoryState::BasicTrajectoryState ( const Surface & aSurface ) [explicit]

construct invalid trajectory state (without parameters)

Definition at line 144 of file BasicTrajectoryState.cc.

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

BasicTrajectoryState::~BasicTrajectoryState ( ) [virtual]

Definition at line 157 of file BasicTrajectoryState.cc.

{}

BasicTrajectoryState::BasicTrajectoryState	(	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 14 of file BasicTrajectoryState.cc.

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

BasicTrajectoryState::BasicTrajectoryState	(	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 BasicTrajectoryState.cc.

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

BasicTrajectoryState::BasicTrajectoryState	(	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 44 of file BasicTrajectoryState.cc.

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

BasicTrajectoryState::BasicTrajectoryState	(	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 60 of file BasicTrajectoryState.cc.

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

BasicTrajectoryState::BasicTrajectoryState	(	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 77 of file BasicTrajectoryState.cc.

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

BasicTrajectoryState::BasicTrajectoryState	(	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 93 of file BasicTrajectoryState.cc.

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

BasicTrajectoryState::BasicTrajectoryState	(	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 109 of file BasicTrajectoryState.cc.

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

BasicTrajectoryState::BasicTrajectoryState	(	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 127 of file BasicTrajectoryState.cc.

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

Member Function Documentation

virtual bool BasicTrajectoryState::canUpdateLocalParameters ( ) const [inline, virtual]

Reimplemented in GCC11_FINAL< T >.

Definition at line 230 of file BasicTrajectoryState.h.

{ return true; }

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

Definition at line 161 of file BasicTrajectoryState.h.

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

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

TrackCharge BasicTrajectoryState::charge ( void ) const [inline]

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

Definition at line 151 of file BasicTrajectoryState.h.

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

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

void BasicTrajectoryState::checkCurvilinError ( ) const [private]

Definition at line 214 of file BasicTrajectoryState.cc.

References createLocalParameters(), JacobianLocalToCurvilinear::jacobian(), likely, localParameters(), LocalTrajectoryError::matrix(), surface(), theField, theFreeState, theLocalError, theLocalParametersValid, and unlikely.

Referenced by freeTrajectoryState().

                                                    {
  if likely(theFreeState->hasCurvilinearError()) return;

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

   theFreeState->setCurvilinearError( cov );
 
   verifyLocalErr(theLocalError);
   verifyCurvErr(cov); 
}

void BasicTrajectoryState::checkGlobalParameters ( ) const [private]

Definition at line 192 of file BasicTrajectoryState.cc.

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

Referenced by freeTrajectoryState().

                                                       {
  if likely(theGlobalParamsUp2Date) return;
 
  //    cout<<"!theGlobalParamsUp2Date"<<endl;
  theGlobalParamsUp2Date = true;
  // 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));
  } 
}

virtual BasicTrajectoryState* BasicTrajectoryState::clone ( ) const [pure virtual]

std::vector< TrajectoryStateOnSurface > BasicTrajectoryState::components ( ) const [virtual]

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

Definition at line 343 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 245 of file BasicTrajectoryState.cc.

References createLocalErrorFromCurvilinearError(), likely, theFreeState, and theLocalError.

Referenced by localError().

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

void BasicTrajectoryState::createLocalErrorFromCurvilinearError ( ) const [private]

Definition at line 252 of file BasicTrajectoryState.cc.

References JacobianCurvilinearToLocal::jacobian(), localParameters(), surface(), theField, theFreeState, and theLocalError.

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

  verifyCurvErr(theFreeState->curvilinearError());
  verifyLocalErr(theLocalError);

}

void BasicTrajectoryState::createLocalParameters ( ) const [private]

Definition at line 234 of file BasicTrajectoryState.cc.

References PV3DBase< T, PVType, FrameType >::mag(), AlCaHLTBitMon_ParallelJobs::p, surface(), theFreeState, theLocalParameters, theLocalParametersValid, GloballyPositioned< T >::toLocal(), PV3DBase< T, PVType, FrameType >::x(), 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 168 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()->curvilinearError();
  }

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

Definition at line 330 of file BasicTrajectoryState.cc.

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

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

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

GlobalVector BasicTrajectoryState::globalDirection ( ) const [inline]

Definition at line 148 of file BasicTrajectoryState.h.

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

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

GlobalVector BasicTrajectoryState::globalMomentum ( ) const [inline]

Definition at line 145 of file BasicTrajectoryState.h.

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

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

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

Definition at line 139 of file BasicTrajectoryState.h.

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

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

GlobalPoint BasicTrajectoryState::globalPosition ( ) const [inline]

Definition at line 142 of file BasicTrajectoryState.h.

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

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

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

Reimplemented in GCC11_FINAL< T >.

Definition at line 223 of file BasicTrajectoryState.h.

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

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

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

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

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

Definition at line 133 of file BasicTrajectoryState.h.

References theFreeState, and theLocalParametersValid.

Referenced by freeTrajectoryState(), and localParameters().

                       {
    return theFreeState || theLocalParametersValid;
  }

LocalVector BasicTrajectoryState::localDirection ( ) const [inline]

Reimplemented in GCC11_FINAL< T >.

Definition at line 195 of file BasicTrajectoryState.h.

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

                                     {
    return localMomentum().unit();
  }

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

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

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

Referenced by localDirection().

                                    {
    return localParameters().momentum();
  }

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

Definition at line 183 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 >, GCC11_FINAL< T >, GCC11_FINAL< T >, GCC11_FINAL< T >, GCC11_FINAL< T >, GCC11_FINAL< T >, GCC11_FINAL< T >, GCC11_FINAL< T >, GCC11_FINAL< T >, and GCC11_FINAL< T >.

Definition at line 189 of file BasicTrajectoryState.h.

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

                                   {
    return localParameters().position();
  }

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

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

Definition at line 180 of file BasicTrajectoryState.h.

References theField.

{ return theField; }

void BasicTrajectoryState::missingError ( char const * where ) const [private]

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

Definition at line 307 of file BasicTrajectoryState.cc.

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

                                                {
  if unlikely(!hasError()) missingError(" trying to rescale");    
  if (theFreeState)
    theFreeState->rescaleError(factor);
  
  if (theLocalError.valid()){
    //do it by hand if the free state is not around.
    bool zeroField = (theField->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 154 of file BasicTrajectoryState.h.

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

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

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

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

Definition at line 208 of file BasicTrajectoryState.h.

References theSurfaceP.

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

                                 {
    return *theSurfaceP;
  }

SurfaceSide BasicTrajectoryState::surfaceSide ( ) const [inline]

Position relative to material, defined relative to momentum vector.

Definition at line 219 of file BasicTrajectoryState.h.

References theSurfaceSide.

                                  {
    return theSurfaceSide;
  }

double BasicTrajectoryState::transverseCurvature ( ) const [inline]

Definition at line 157 of file BasicTrajectoryState.h.

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

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

void BasicTrajectoryState::update	(	const LocalTrajectoryParameters &	p,
		const Surface &	aSurface,
		const MagneticField *	field,
		const SurfaceSide	side
	)		`[virtual]`

Reimplemented in GCC11_FINAL< T >.

Definition at line 270 of file BasicTrajectoryState.cc.

References AlCaHLTBitMon_ParallelJobs::p, theField, theGlobalParamsUp2Date, theLocalError, theLocalParameters, theLocalParametersValid, theSurfaceP, theSurfaceSide, and theWeight.

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

    theGlobalParamsUp2Date   = false;
    theLocalParametersValid  = true;
}

void BasicTrajectoryState::update	(	const LocalTrajectoryParameters &	p,
		const LocalTrajectoryError &	err,
		const Surface &	aSurface,
		const MagneticField *	field,
		const SurfaceSide	side,
		double	weight
	)		`[virtual]`

Reimplemented in GCC11_FINAL< T >.

Definition at line 287 of file BasicTrajectoryState.cc.

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

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

    theGlobalParamsUp2Date   = false;
    theLocalParametersValid  = true;

}