#include <BasicMultiTrajectoryState.h>

Inheritance diagram for BasicMultiTrajectoryState:

Public Member Functions
	BasicMultiTrajectoryState (const std::vector< TSOS > &tsvec)

	BasicMultiTrajectoryState ()

virtual bool	canUpdateLocalParameters () const

virtual BasicMultiTrajectoryState *	clone () const

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

void	rescaleError (double factor)

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)

Public Member Functions inherited from BasicTrajectoryState
	BasicTrajectoryState ()

	BasicTrajectoryState (const Surface &aSurface)
	construct invalid trajectory state (without parameters) More...

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

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

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

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

const CartesianTrajectoryError	cartesianError () const

TrackCharge	charge () 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. More...

double	transverseCurvature () const

virtual double	weight () 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 ()

Private Types
typedef TrajectoryStateOnSurface	TSOS

Private Member Functions
void	combine ()

Private Attributes
std::vector< TSOS >	theStates

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

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

static bool	s_usePool =false

Detailed Description

Class which combines a set of components of a Gaussian mixture into a single component. Given all the components of a mixture, it calculates the mean and covariance matrix of the entire mixture. This combiner class can also be used in the process of transforming a Gaussian mixture into another Gaussian mixture with a smaller number of components. The relevant formulas can be found in R. Fruhwirth, Computer Physics Communications 100 (1997), 1.

Definition at line 17 of file BasicMultiTrajectoryState.h.

Member Typedef Documentation

typedef TrajectoryStateOnSurface BasicMultiTrajectoryState::TSOS

private

Definition at line 19 of file BasicMultiTrajectoryState.h.

Constructor & Destructor Documentation

BasicMultiTrajectoryState::BasicMultiTrajectoryState ( const std::vector< TSOS > & tsvec )

Definition at line 8 of file BasicMultiTrajectoryState.cc.

References combine(), edm::hlt::Exception, i, theStates, and unlikely.

                                                                                   :
   BasicTrajectoryState(tsvec.front().surface()) {
  
   // only accept planes!!
   const BoundPlane* bp = dynamic_cast<const BoundPlane*>(&tsvec.begin()->surface());
   if unlikely( bp==0 )
            throw cms::Exception("LogicError") << "MultiTrajectoryState constructed on cylinder";
    
   theStates.reserve(tsvec.size());
   for (std::vector<TSOS>::const_iterator i=tsvec.begin(); i!=tsvec.end(); i++) {
     if unlikely(!i->isValid()) {
       throw cms::Exception("LogicError") << "MultiTrajectoryState constructed with invalid state";
     }
     if unlikely(i->hasError() != tsvec.front().hasError()) {
       throw cms::Exception("LogicError") << "MultiTrajectoryState mixes states with and without errors";
     }
     if unlikely( &i->surface() != &tsvec.front().surface()) {
       throw cms::Exception("LogicError") << "MultiTrajectoryState mixes states with different surfaces";
     }
     if unlikely( i->surfaceSide() != tsvec.front().surfaceSide()) {
       throw cms::Exception("LogicError") 
     << "MultiTrajectoryState mixes states defined before and after material";
     }
     if unlikely( i->localParameters().pzSign()*tsvec.front().localParameters().pzSign()<0. ) {
       throw cms::Exception("LogicError") 
     << "MultiTrajectoryState mixes states with different signs of local p_z";
     }
 
     theStates.push_back( *i);
   }
   combine();
 }

BasicMultiTrajectoryState::BasicMultiTrajectoryState ( )

inline

Definition at line 25 of file BasicMultiTrajectoryState.h.

Referenced by clone().

25 {}

Member Function Documentation

virtual bool BasicMultiTrajectoryState::canUpdateLocalParameters ( ) const

inlinevirtual

Reimplemented from BasicTrajectoryState.

Definition at line 45 of file BasicMultiTrajectoryState.h.

45 { return false; }

virtual BasicMultiTrajectoryState* BasicMultiTrajectoryState::clone ( void ) const

inlinevirtual

Implements BasicTrajectoryState.

Definition at line 36 of file BasicMultiTrajectoryState.h.

References BasicMultiTrajectoryState().

                                                    {
     return new BasicMultiTrajectoryState(*this);
   }

void BasicMultiTrajectoryState::combine ( )

private

Definition at line 57 of file BasicMultiTrajectoryState.cc.

References diffTreeTool::diff, timingPdfMaker::mean, alignCSCRings::s, theStates, unlikely, BasicTrajectoryState::update(), and BasicTrajectoryState::weight().

Referenced by BasicMultiTrajectoryState(), and rescaleError().

                                     {
   const std::vector<TrajectoryStateOnSurface>& tsos = theStates;
 
   if unlikely(tsos.empty()) {
     edm::LogError("MultiTrajectoryStateCombiner") 
       << "Trying to collapse empty set of trajectory states!";
     return;
   }
 
   double pzSign = tsos.front().localParameters().pzSign();
   for (std::vector<TrajectoryStateOnSurface>::const_iterator it = tsos.begin(); 
        it != tsos.end(); it++) {
     if unlikely(it->localParameters().pzSign() != pzSign) {
       edm::LogError("MultiTrajectoryStateCombiner") 
     << "Trying to collapse trajectory states with different signs on p_z!";
       return;
     }
   }
   
   if unlikely(tsos.size() == 1) {
     BasicTrajectoryState::update(tsos.front().localParameters(), 
                  tsos.front().localError(), 
                  tsos.front().surface(), 
                  tsos.front().magneticField(),
                  tsos.front().surfaceSide(), 
                  tsos.front().weight()
                  );
     return;
   }
   
   double sumw = 0.;
   //int dim = tsos.front().localParameters().vector().num_row();
   AlgebraicVector5 mean;
   AlgebraicSymMatrix55 covarPart1, covarPart2;
   for (std::vector<TrajectoryStateOnSurface>::const_iterator it1 = tsos.begin(); 
        it1 != tsos.end(); it1++) {
     double weight = it1->weight();
     AlgebraicVector5 param = it1->localParameters().vector();
     sumw += weight;
     mean += weight * param;
     covarPart1 += weight * it1->localError().matrix();
     for (std::vector<TrajectoryStateOnSurface>::const_iterator it2 = it1 + 1; 
      it2 != tsos.end(); it2++) {
       AlgebraicVector5 diff = param - it2->localParameters().vector();
       AlgebraicSymMatrix11 s = AlgebraicMatrixID(); //stupid trick to make CLHEP work decently
       covarPart2 += (weight * it2->weight()) * 
                     ROOT::Math::Similarity(AlgebraicMatrix51(diff.Array(), 5), s);
                         //FIXME: we can surely write this thing in a better way
     }   
   }
   double sumwI = 1.0/sumw;
   mean *= sumwI;
   covarPart1 *= sumwI; covarPart2 *= (sumwI*sumwI);
   AlgebraicSymMatrix55 covar = covarPart1 + covarPart2;
 
   BasicTrajectoryState::update(LocalTrajectoryParameters(mean, pzSign), 
                    LocalTrajectoryError(covar), 
                    tsos.front().surface(), 
                    tsos.front().magneticField(),
                    tsos.front().surfaceSide(), 
                    sumw);
 }

virtual std::vector<TrajectoryStateOnSurface> BasicMultiTrajectoryState::components ( ) const

inlinevirtual

Reimplemented from BasicTrajectoryState.

Definition at line 40 of file BasicMultiTrajectoryState.h.

References theStates.

                                                                  {
     return theStates;
   }

void BasicMultiTrajectoryState::rescaleError ( double factor )

Rescaling the error of the mixture with a given factor. Please note that this rescaling is imposed on each of the components of the mixture and does therefore not exactly correspond to rescaling theCombinedState with the same factor.

Definition at line 43 of file BasicMultiTrajectoryState.cc.

References combine(), theStates, and unlikely.

                                                           {
 
   if unlikely(theStates.empty()) {
     edm::LogError("BasicMultiTrajectoryState") << "Trying to rescale errors of empty MultiTrajectoryState!";
     return;
   }
   
   for (std::vector<TSOS>::iterator it = theStates.begin(); it != theStates.end(); it++) {
     it->rescaleError(factor);
   }
   combine();
 }

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

virtual

Reimplemented from BasicTrajectoryState.

Definition at line 122 of file BasicMultiTrajectoryState.cc.

References edm::hlt::Exception.

 {
   throw cms::Exception("LogicError", 
                        "BasicMultiTrajectoryState::update(LocalTrajectoryParameters, Surface, ...) called even if canUpdateLocalParameters() is false");
 }

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

virtual

Reimplemented from BasicTrajectoryState.

Definition at line 133 of file BasicMultiTrajectoryState.cc.

References edm::hlt::Exception.

 {
   throw cms::Exception("LogicError", 
                        "BasicMultiTrajectoryState::update(LocalTrajectoryParameters, LocalTrajectoryError, ...) called even if canUpdateLocalParameters() is false");
 }

Member Data Documentation

std::vector<TSOS> BasicMultiTrajectoryState::theStates

private

Definition at line 58 of file BasicMultiTrajectoryState.h.

Referenced by BasicMultiTrajectoryState(), combine(), components(), and rescaleError().

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation