#include <MultiRefittedTS.h>

Inheritance diagram for MultiRefittedTS:

Public Types
typedef ReferenceCountingPointer< LinearizedTrackState< 5 > >	RefCountedLinearizedTrackState

typedef ReferenceCountingPointer< RefittedTrackState< 5 > >	RefCountedRefittedTrackState

Public Types inherited from RefittedTrackState< 5 >
typedef ROOT::Math::SMatrix< double, N, N, ROOT::Math::MatRepSym< double, N > >	AlgebraicSymMatrixNN

typedef ROOT::Math::SVector< double, N-2 >	AlgebraicVectorM

typedef ROOT::Math::SVector< double, N >	AlgebraicVectorN

Public Member Functions
std::vector< ReferenceCountingPointer< RefittedTrackState< 5 > > >	components () const override

AlgebraicSymMatrixNN	covariance () const override

FreeTrajectoryState	freeTrajectoryState () const override

AlgebraicVectorM	momentumVector () const override

	MultiRefittedTS (const std::vector< RefCountedRefittedTrackState > &prtsComp, const Surface &referenceSurface)

	MultiRefittedTS (const std::vector< RefCountedRefittedTrackState > &prtsComp, const GlobalPoint &referencePosition)

AlgebraicVectorN	parameters () const override

GlobalPoint	position () const override

ReferenceCountingPointer< RefittedTrackState< 5 > >	stateWithNewWeight (const double newWeight) const override

TrajectoryStateOnSurface	trajectoryStateOnSurface (const Surface &surface) const override

TrajectoryStateOnSurface	trajectoryStateOnSurface (const Surface &surface, const Propagator &propagator) const override

reco::TransientTrack	transientTrack () const override

double	weight () const override

	~MultiRefittedTS () override

Public Member Functions inherited from RefittedTrackState< 5 >
	~RefittedTrackState () override

Private Types
typedef std::vector< RefCountedRefittedTrackState >	RTSvector

Private Member Functions
void	computeFreeTrajectoryState () const

Private Attributes
FreeTrajectoryState	fts

bool	ftsAvailable

const GlobalPoint	refPosition

ConstReferenceCountingPointer< Surface >	refSurface

const bool	surf

RTSvector	theComponents

double	totalWeight

bool	totalWeightAvailable

Detailed Description

Definition at line 17 of file MultiRefittedTS.h.

Member Typedef Documentation

typedef ReferenceCountingPointer<LinearizedTrackState<5> > MultiRefittedTS::RefCountedLinearizedTrackState

Definition at line 22 of file MultiRefittedTS.h.

typedef ReferenceCountingPointer<RefittedTrackState<5> > MultiRefittedTS::RefCountedRefittedTrackState

Definition at line 21 of file MultiRefittedTS.h.

typedef std::vector<RefCountedRefittedTrackState > MultiRefittedTS::RTSvector

private

Definition at line 121 of file MultiRefittedTS.h.

Constructor & Destructor Documentation

MultiRefittedTS::MultiRefittedTS	(	const std::vector< RefCountedRefittedTrackState > &	prtsComp,
		const Surface &	referenceSurface
	)

Constructor with a reference surface, to be used to assemble the TSOS components on one identical surface.

Definition at line 12 of file MultiRefittedTS.cc.

Referenced by stateWithNewWeight().

                                       :
     theComponents(prtsComp), ftsAvailable(false), refSurface(&referenceSurface),
     surf(true) {}

MultiRefittedTS::MultiRefittedTS	(	const std::vector< RefCountedRefittedTrackState > &	prtsComp,
		const GlobalPoint &	referencePosition
	)

Constructor with a reference position. The surface which is going to be usedto assemble the TSOS components will be the surface perpendicular to the PCA of the state with the highest weight to the reference point.

Definition at line 18 of file MultiRefittedTS.cc.

                                            :
     theComponents(prtsComp), ftsAvailable(false), refPosition(referencePosition),
     surf(false) {}

MultiRefittedTS::~MultiRefittedTS ( )

inlineoverride

Definition at line 39 of file MultiRefittedTS.h.

References covariance(), freeTrajectoryState(), momentumVector(), parameters(), position(), PhotonConversionTrajectorySeedProducerFromQuadruplets_cfi::propagator, trajectoryStateOnSurface(), transientTrack(), and weight().

39 {}

Member Function Documentation

std::vector<ReferenceCountingPointer<RefittedTrackState<5> > > MultiRefittedTS::components ( ) const

inlineoverridevirtual

Implements RefittedTrackState< 5 >.

Definition at line 99 of file MultiRefittedTS.h.

References computeFreeTrajectoryState(), stateWithNewWeight(), and theComponents.

   {
     return theComponents;
   }

void MultiRefittedTS::computeFreeTrajectoryState ( ) const

private

Definition at line 33 of file MultiRefittedTS.cc.

References TransverseImpactPointExtrapolator::extrapolate(), freeTrajectoryState(), TrajectoryStateOnSurface::freeTrajectoryState(), fts, ftsAvailable, refPosition, surf, TrajectoryStateOnSurface::surface(), theComponents, and trajectoryStateOnSurface().

Referenced by components(), and freeTrajectoryState().

 {
   if (surf) {
     fts =  *(trajectoryStateOnSurface(*refSurface).freeTrajectoryState());
   } else {
     double maxWeight = -1.;
     RTSvector::const_iterator maxIt;
     for (RTSvector::const_iterator it = theComponents.begin(); 
       it != theComponents.end(); it++) {
       if ( (**it).weight() > maxWeight ) {
     maxWeight = (**it).weight();
     maxIt = it;
       }
     }
 
     TransverseImpactPointExtrapolator tipe(&((**maxIt).freeTrajectoryState().parameters().magneticField()));
     TrajectoryStateOnSurface initialTSOS = tipe.extrapolate((**maxIt).freeTrajectoryState(), refPosition);
 
     fts = *(trajectoryStateOnSurface(initialTSOS.surface()).freeTrajectoryState());
   }
   ftsAvailable = true;
 }

MultiRefittedTS::AlgebraicSymMatrixNN MultiRefittedTS::covariance ( void ) const

overridevirtual

The covariance matrix. Not possible yet for a multi-state, throws an exception.

The covariance matrix

Implements RefittedTrackState< 5 >.

Definition at line 72 of file MultiRefittedTS.cc.

Referenced by ~MultiRefittedTS().

 {
   throw VertexException
     ("MultiRefittedTS::freeTrajectoryState(): Don't know how to do that yet...");
 }

FreeTrajectoryState MultiRefittedTS::freeTrajectoryState ( ) const

overridevirtual

Returns a FreeTrajectoryState. It will be the FTS of the single, collapsed state.

Transformation into a FreeTrajectoryState

Implements RefittedTrackState< 5 >.

Definition at line 27 of file MultiRefittedTS.cc.

References computeFreeTrajectoryState(), fts, and ftsAvailable.

Referenced by computeFreeTrajectoryState(), transientTrack(), and ~MultiRefittedTS().

 {
   if (!ftsAvailable) computeFreeTrajectoryState();
   return fts;
 }

MultiRefittedTS::AlgebraicVectorM MultiRefittedTS::momentumVector ( ) const

overridevirtual

Vector containing the parameters describing the momentum as the vertex. These are (signed transverse curvature, theta, phi). Not possible yet for a multi-state, throws an exception.

Vector containing the parameters describing the momentum as the vertex. These are (signed transverse curvature, theta, phi)

Implements RefittedTrackState< 5 >.

Definition at line 93 of file MultiRefittedTS.cc.

Referenced by ~MultiRefittedTS().

 {
   throw VertexException
     ("MultiRefittedTS::freeTrajectoryState(): Don't know how to do that yet...");
 }

MultiRefittedTS::AlgebraicVectorN MultiRefittedTS::parameters ( void ) const

overridevirtual

Vector containing the refitted track parameters. Not possible yet for a multi-state, throws an exception.

Vector containing the refitted track parameters.
These are (signed transverse curvature, theta, phi, (signed) transverse , longitudinal impact parameter)

Implements RefittedTrackState< 5 >.

Definition at line 62 of file MultiRefittedTS.cc.

Referenced by ~MultiRefittedTS().

 {
   throw VertexException
     ("MultiRefittedTS::freeTrajectoryState(): Don't know how to do that yet...");
 }

GlobalPoint MultiRefittedTS::position ( ) const

overridevirtual

Position at which the momentum is defined. Not possible yet for a multi-state, throws an exception.

Position at which the momentum is defined.

Implements RefittedTrackState< 5 >.

Definition at line 82 of file MultiRefittedTS.cc.

Referenced by ~MultiRefittedTS().

 {
   throw VertexException
     ("MultiRefittedTS::freeTrajectoryState(): Don't know how to do that yet...");
 }

ReferenceCountingPointer< RefittedTrackState< 5 > > MultiRefittedTS::stateWithNewWeight ( const double newWeight ) const

overridevirtual

This method is meant to returns a new refitted state of the same type, but with another weight. As we can have several components, each component of the new multi-state will be reweighted so that the sum of all weights is equal to the specified weight. The current state is unchanged.

Implements RefittedTrackState< 5 >.

Definition at line 117 of file MultiRefittedTS.cc.

References MultiRefittedTS(), refPosition, refSurface, surf, theComponents, and weight().

Referenced by components().

 {
   if (weight() < DBL_MIN) {
   throw VertexException
     ("MultiRefittedTS::stateWithNewWeight(): Can not reweight multi-state with total weight < DBL_MIN");
   }
   double factor = newWeight/weight();
 
   RTSvector reWeightedRTSC;
   reWeightedRTSC.reserve(theComponents.size());
 
   for (RTSvector::const_iterator it = theComponents.begin(); 
     it != theComponents.end(); it++) {
     reWeightedRTSC.push_back((**it).stateWithNewWeight((**it).weight()*factor));
   }
   if (surf) {
     return RefCountedRefittedTrackState(new MultiRefittedTS(reWeightedRTSC, *refSurface));
   } else {
     return RefCountedRefittedTrackState(new MultiRefittedTS(reWeightedRTSC, refPosition));
   }
 }

TrajectoryStateOnSurface MultiRefittedTS::trajectoryStateOnSurface ( const Surface & surface ) const

overridevirtual

Returns a multi-state TSOS at a given surface

Transformation into a TSOS at a given surface

Implements RefittedTrackState< 5 >.

Definition at line 142 of file MultiRefittedTS.cc.

References theComponents.

Referenced by computeFreeTrajectoryState(), and ~MultiRefittedTS().

 {
   vector<TrajectoryStateOnSurface> tsosComponents;
   tsosComponents.reserve(theComponents.size());
   for (RTSvector::const_iterator it = theComponents.begin(); 
     it != theComponents.end(); it++) {
     tsosComponents.push_back((**it).trajectoryStateOnSurface(surface));
   }
   return TrajectoryStateOnSurface((BasicTrajectoryState*)new BasicMultiTrajectoryState(tsosComponents));
 }

TrajectoryStateOnSurface MultiRefittedTS::trajectoryStateOnSurface	(	const Surface &	surface,
		const Propagator &	propagator
	)		const

overridevirtual

Returns a multi-state TSOS at a given surface, with a given propagator

Implements RefittedTrackState< 5 >.

Definition at line 154 of file MultiRefittedTS.cc.

References theComponents.

 { //fixme... is the propagation done correctly? Is there a gsf propagator?
   vector<TrajectoryStateOnSurface> tsosComponents;
   tsosComponents.reserve(theComponents.size());
   for (RTSvector::const_iterator it = theComponents.begin(); 
     it != theComponents.end(); it++) {
     tsosComponents.push_back((**it).trajectoryStateOnSurface(surface, propagator));
   }
   return TrajectoryStateOnSurface((BasicTrajectoryState*)new BasicMultiTrajectoryState(tsosComponents));
 }

reco::TransientTrack MultiRefittedTS::transientTrack ( ) const

overridevirtual

Returns a new reco::Track, which can then be made persistent. The parameters are taken from FTS described above.

Implements RefittedTrackState< 5 >.

Definition at line 166 of file MultiRefittedTS.cc.

References TransientTrackFromFTSFactory::build(), and freeTrajectoryState().

Referenced by ~MultiRefittedTS().

 {
   TransientTrackFromFTSFactory factory;
   return factory.build(freeTrajectoryState());
 }

double MultiRefittedTS::weight ( ) const

overridevirtual

The weight of this component in a mixture

Implements RefittedTrackState< 5 >.

Definition at line 99 of file MultiRefittedTS.cc.

References gather_cfg::cout, theComponents, totalWeight, and totalWeightAvailable.

Referenced by stateWithNewWeight(), and ~MultiRefittedTS().

 {
   if (!totalWeightAvailable)
   {
     totalWeight = 0.;
     if (theComponents.empty()) {
       cout << "Asking for weight of empty MultiRefittedTS, returning zero!" << endl;
     }
     for (RTSvector::const_iterator it = theComponents.begin(); 
       it != theComponents.end(); it++) {
       totalWeight += (**it).weight();
     }
   }
   return totalWeight;
 }