MultiRefittedTS Class Reference

#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

RefCountedLinearizedTrackState

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

Definition at line 20 of file MultiRefittedTS.h.

RefCountedRefittedTrackState

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

Definition at line 19 of file MultiRefittedTS.h.

RTSvector

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

private

Definition at line 109 of file MultiRefittedTS.h.

Constructor & Destructor Documentation

MultiRefittedTS() [1/2]

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 11 of file MultiRefittedTS.cc.

Referenced by stateWithNewWeight().

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

MultiRefittedTS() [2/2]

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 15 of file MultiRefittedTS.cc.

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

~MultiRefittedTS()

MultiRefittedTS::~MultiRefittedTS ( )

inlineoverride

Definition at line 35 of file MultiRefittedTS.h.

{}

Member Function Documentation

components()

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

inlineoverridevirtual

Implements RefittedTrackState< 5 >.

Definition at line 94 of file MultiRefittedTS.h.

References theComponents.

{ return theComponents; }

computeFreeTrajectoryState()

void MultiRefittedTS::computeFreeTrajectoryState ( ) const

private

Definition at line 29 of file MultiRefittedTS.cc.

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

Referenced by 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;
}

covariance()

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 64 of file MultiRefittedTS.cc.

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

freeTrajectoryState()

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 23 of file MultiRefittedTS.cc.

References computeFreeTrajectoryState(), fts, and ftsAvailable.

Referenced by computeFreeTrajectoryState(), and transientTrack().

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

momentumVector()

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 81 of file MultiRefittedTS.cc.

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

parameters()

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 56 of file MultiRefittedTS.cc.

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

position()

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 72 of file MultiRefittedTS.cc.

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

stateWithNewWeight()

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 98 of file MultiRefittedTS.cc.

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

                                                                                                                 {
  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() [1/2]

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 121 of file MultiRefittedTS.cc.

References theComponents.

Referenced by computeFreeTrajectoryState().

                                                                                               {
  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() [2/2]

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 130 of file MultiRefittedTS.cc.

References TrackCandidateProducer_cfi::propagator, and 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));
}

transientTrack()

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 140 of file MultiRefittedTS.cc.

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

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

weight()

double MultiRefittedTS::weight ( ) const

overridevirtual

The weight of this component in a mixture

Implements RefittedTrackState< 5 >.

Definition at line 85 of file MultiRefittedTS.cc.

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

Referenced by stateWithNewWeight().

                                     {
  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;
}

Member Data Documentation

fts

FreeTrajectoryState MultiRefittedTS::fts

mutableprivate

Definition at line 114 of file MultiRefittedTS.h.

Referenced by computeFreeTrajectoryState(), and freeTrajectoryState().

ftsAvailable

bool MultiRefittedTS::ftsAvailable

mutableprivate

Definition at line 112 of file MultiRefittedTS.h.

Referenced by computeFreeTrajectoryState(), and freeTrajectoryState().

refPosition

const GlobalPoint MultiRefittedTS::refPosition

private

Definition at line 115 of file MultiRefittedTS.h.

Referenced by computeFreeTrajectoryState(), and stateWithNewWeight().

refSurface

ConstReferenceCountingPointer<Surface> MultiRefittedTS::refSurface

private

Definition at line 116 of file MultiRefittedTS.h.

Referenced by stateWithNewWeight().

surf

const bool MultiRefittedTS::surf

private

Definition at line 117 of file MultiRefittedTS.h.

Referenced by computeFreeTrajectoryState(), and stateWithNewWeight().

theComponents

RTSvector MultiRefittedTS::theComponents

mutableprivate

Definition at line 111 of file MultiRefittedTS.h.

Referenced by components(), computeFreeTrajectoryState(), stateWithNewWeight(), trajectoryStateOnSurface(), and weight().

totalWeight

double MultiRefittedTS::totalWeight

mutableprivate

Definition at line 113 of file MultiRefittedTS.h.

Referenced by weight().

totalWeightAvailable

bool MultiRefittedTS::totalWeightAvailable

mutableprivate

Definition at line 112 of file MultiRefittedTS.h.

Referenced by weight().