#include <TrajectoryFactoryBase.h>

Inheritance diagram for TrajectoryFactoryBase:

Public Types
typedef AlignmentAlgorithmBase::ConstTrajTrackPair	ConstTrajTrackPair
typedef AlignmentAlgorithmBase::ConstTrajTrackPairCollection	ConstTrajTrackPairCollection
typedef std::vector < TrajectoryStateOnSurface >	ExternalPredictionCollection
typedef ReferenceTrajectoryBase::MaterialEffects	MaterialEffects
typedef std::vector < ReferenceTrajectoryPtr >	ReferenceTrajectoryCollection
typedef ReferenceTrajectoryBase::ReferenceTrajectoryPtr	ReferenceTrajectoryPtr
typedef std::pair < TrajectoryStateOnSurface, TransientTrackingRecHit::ConstRecHitContainer >	TrajectoryInput
Public Member Functions
virtual TrajectoryFactoryBase *	clone (void) const =0
MaterialEffects	materialEffects (void) const
PropagationDirection	propagationDirection (void) const
virtual const ReferenceTrajectoryCollection	trajectories (const edm::EventSetup &setup, const ConstTrajTrackPairCollection &tracks, const ExternalPredictionCollection &external, const reco::BeamSpot &beamSpot) const =0
virtual const ReferenceTrajectoryCollection	trajectories (const edm::EventSetup &setup, const ConstTrajTrackPairCollection &tracks, const reco::BeamSpot &beamSpot) const =0
	TrajectoryFactoryBase (const edm::ParameterSet &config)
virtual	~TrajectoryFactoryBase (void)
Protected Member Functions
virtual const TrajectoryInput	innermostStateAndRecHits (const ConstTrajTrackPair &track) const
virtual const Trajectory::DataContainer	orderedTrajectoryMeasurements (const Trajectory &trajectory) const
bool	sameSurface (const Surface &s1, const Surface &s2) const
bool	useRecHit (const TransientTrackingRecHit::ConstRecHitPointer &hitPtr) const
Protected Attributes
bool	theUseBeamSpot
Private Member Functions
MaterialEffects	materialEffects (const std::string &strME) const
PropagationDirection	propagationDirection (const std::string &strPD) const
Private Attributes
MaterialEffects	theMaterialEffects
PropagationDirection	thePropDir
bool	theUseInvalidHits
bool	theUseProjectedHits
bool	theUseWithoutDet

Detailed Description

Definition at line 17 of file TrajectoryFactoryBase.h.

Member Typedef Documentation

typedef AlignmentAlgorithmBase::ConstTrajTrackPair TrajectoryFactoryBase::ConstTrajTrackPair

Definition at line 24 of file TrajectoryFactoryBase.h.

typedef AlignmentAlgorithmBase::ConstTrajTrackPairCollection TrajectoryFactoryBase::ConstTrajTrackPairCollection

Definition at line 25 of file TrajectoryFactoryBase.h.

typedef std::vector< TrajectoryStateOnSurface > TrajectoryFactoryBase::ExternalPredictionCollection

Definition at line 28 of file TrajectoryFactoryBase.h.

typedef ReferenceTrajectoryBase::MaterialEffects TrajectoryFactoryBase::MaterialEffects

Definition at line 23 of file TrajectoryFactoryBase.h.

typedef std::vector< ReferenceTrajectoryPtr > TrajectoryFactoryBase::ReferenceTrajectoryCollection

Definition at line 26 of file TrajectoryFactoryBase.h.

typedef ReferenceTrajectoryBase::ReferenceTrajectoryPtr TrajectoryFactoryBase::ReferenceTrajectoryPtr

Definition at line 22 of file TrajectoryFactoryBase.h.

typedef std::pair< TrajectoryStateOnSurface, TransientTrackingRecHit::ConstRecHitContainer > TrajectoryFactoryBase::TrajectoryInput

Definition at line 27 of file TrajectoryFactoryBase.h.

Constructor & Destructor Documentation

TrajectoryFactoryBase::TrajectoryFactoryBase ( const edm::ParameterSet & config )

Definition at line 10 of file TrajectoryFactoryBase.cc.

References edm::ParameterSet::getParameter(), materialEffects(), propagationDirection(), theMaterialEffects, thePropDir, theUseBeamSpot, theUseInvalidHits, theUseProjectedHits, and theUseWithoutDet.

{
  const std::string strMaterialEffects = config.getParameter< std::string >( "MaterialEffects" );
  theMaterialEffects = this->materialEffects( strMaterialEffects );
  
  const std::string strPropagationDirection = config.getParameter< std::string >( "PropagationDirection" );
  thePropDir = this->propagationDirection( strPropagationDirection );

  theUseWithoutDet = config.getParameter< bool >( "UseHitWithoutDet" );
  theUseInvalidHits = config.getParameter< bool >( "UseInvalidHits" );
  theUseProjectedHits = config.getParameter< bool >( "UseProjectedHits" );
  theUseBeamSpot = config.getParameter< bool >( "UseBeamSpot" );
}

TrajectoryFactoryBase::~TrajectoryFactoryBase ( void ) [virtual]

Definition at line 25 of file TrajectoryFactoryBase.cc.

{}

Member Function Documentation

virtual TrajectoryFactoryBase* TrajectoryFactoryBase::clone ( void ) const [pure virtual]

Implemented in BzeroReferenceTrajectoryFactory, CombinedTrajectoryFactory, DualBzeroTrajectoryFactory, DualKalmanFactory, DualTrajectoryFactory, ReferenceTrajectoryFactory, and TwoBodyDecayTrajectoryFactory.

const TrajectoryFactoryBase::TrajectoryInput TrajectoryFactoryBase::innermostStateAndRecHits ( const ConstTrajTrackPair & track ) const [protected, virtual]

Definition at line 29 of file TrajectoryFactoryBase.cc.

References orderedTrajectoryMeasurements(), query::result, and useRecHit().

Referenced by TwoBodyDecayTrajectoryFactory::constructTrajectories(), BzeroReferenceTrajectoryFactory::trajectories(), and ReferenceTrajectoryFactory::trajectories().

{
  TrajectoryInput result;

  // get the trajectory measurements in the correct order, i.e. reverse if needed
  Trajectory::DataContainer trajectoryMeasurements 
    = this->orderedTrajectoryMeasurements( *track.first );
  Trajectory::DataContainer::iterator itM = trajectoryMeasurements.begin();

  // get the innermost valid trajectory state - the corresponding hit must be o.k. as well
  while ( itM != trajectoryMeasurements.end() )
  {
    if ( ( *itM ).updatedState().isValid() && useRecHit( ( *itM ).recHit() ) ) break;
    ++itM;
  }
  if ( itM != trajectoryMeasurements.end() ) result.first = ( *itM ).updatedState();

  // get the valid RecHits
  while ( itM != trajectoryMeasurements.end() )
  {
    TransientTrackingRecHit::ConstRecHitPointer aRecHit = ( *itM ).recHit();
    if ( useRecHit( aRecHit ) ) result.second.push_back( aRecHit );
    ++itM;
  }

  return result;
}

MaterialEffects TrajectoryFactoryBase::materialEffects ( void ) const [inline]

Definition at line 44 of file TrajectoryFactoryBase.h.

References theMaterialEffects.

Referenced by TwoBodyDecayTrajectoryFactory::constructTrajectories(), BzeroReferenceTrajectoryFactory::trajectories(), DualKalmanFactory::trajectories(), ReferenceTrajectoryFactory::trajectories(), DualBzeroTrajectoryFactory::trajectories(), DualTrajectoryFactory::trajectories(), and TrajectoryFactoryBase().

{ return theMaterialEffects; }

TrajectoryFactoryBase::MaterialEffects TrajectoryFactoryBase::materialEffects ( const std::string & strME ) const [private]

Definition at line 112 of file TrajectoryFactoryBase.cc.

References ReferenceTrajectoryBase::breakPoints, ReferenceTrajectoryBase::brokenLinesCoarse, ReferenceTrajectoryBase::brokenLinesFine, runregparse::combined, ReferenceTrajectoryBase::energyLoss, Exception, ReferenceTrajectoryBase::multipleScattering, and btau_dummy::none.

{
  if ( strME == "MultipleScattering" ) return ReferenceTrajectoryBase::multipleScattering;
  if ( strME == "EnergyLoss" ) return ReferenceTrajectoryBase::energyLoss;
  if ( strME == "Combined" ) return ReferenceTrajectoryBase::combined;
  if ( strME == "None" ) return ReferenceTrajectoryBase::none;
  if ( strME == "BreakPoints" ) return ReferenceTrajectoryBase::breakPoints;
  if ( strME == "BrokenLines" ) return ReferenceTrajectoryBase::brokenLinesCoarse;
  if ( strME == "BrokenLinesCoarse" ) return ReferenceTrajectoryBase::brokenLinesCoarse;
  if ( strME == "BrokenLinesFine" ) return ReferenceTrajectoryBase::brokenLinesFine;
          
  throw cms::Exception("BadConfig")
    << "[TrajectoryFactoryBase::materialEffects] Unknown parameter: " << strME;
}

const Trajectory::DataContainer TrajectoryFactoryBase::orderedTrajectoryMeasurements ( const Trajectory & trajectory ) const [protected, virtual]

Definition at line 59 of file TrajectoryFactoryBase.cc.

References anyDirection, dir, Trajectory::direction(), Trajectory::measurements(), definitions::original, and thePropDir.

Referenced by innermostStateAndRecHits(), DualBzeroTrajectoryFactory::referenceStateAndRecHits(), DualTrajectoryFactory::referenceStateAndRecHits(), and DualKalmanFactory::referenceStateAndRecHits().

{
  const PropagationDirection dir = trajectory.direction();
  const bool hitsAreReverse = ( ( dir == thePropDir || thePropDir == anyDirection ) ? false : true );

  const Trajectory::DataContainer & original = trajectory.measurements();

  if ( hitsAreReverse )
  {
    // Simply use this line instead of the copying by hand?
    // const Trajectory::DataContainer reordered(original.rbegin(), original.rend());
    Trajectory::DataContainer reordered;
    reordered.reserve( original.size() );

    Trajectory::DataContainer::const_reverse_iterator itM;
    for ( itM = original.rbegin(); itM != original.rend(); ++itM )
    {
      reordered.push_back( *itM );
    }
    return reordered;
  }

  return original;
}

PropagationDirection TrajectoryFactoryBase::propagationDirection ( void ) const [inline]

Definition at line 45 of file TrajectoryFactoryBase.h.

References thePropDir.

Referenced by TwoBodyDecayTrajectoryFactory::constructTrajectories(), DualKalmanFactory::trajectories(), BzeroReferenceTrajectoryFactory::trajectories(), ReferenceTrajectoryFactory::trajectories(), DualBzeroTrajectoryFactory::trajectories(), DualTrajectoryFactory::trajectories(), and TrajectoryFactoryBase().

{ return thePropDir; }

PropagationDirection TrajectoryFactoryBase::propagationDirection ( const std::string & strPD ) const [private]

Definition at line 129 of file TrajectoryFactoryBase.cc.

References alongMomentum, anyDirection, Exception, and oppositeToMomentum.

{
  if ( strPD == "oppositeToMomentum" ) return oppositeToMomentum;
  if ( strPD == "alongMomentum" ) return alongMomentum;
  if ( strPD == "anyDirection" ) return anyDirection;

  throw cms::Exception("BadConfig")
    << "[TrajectoryFactoryBase::propagationDirection] Unknown parameter: " << strPD;
}

bool TrajectoryFactoryBase::sameSurface	(	const Surface &	s1,
		const Surface &	s2
	)		const `[protected]`

Definition at line 85 of file TrajectoryFactoryBase.cc.

References GloballyPositioned< T >::eta(), GloballyPositioned< T >::phi(), and GloballyPositioned< T >::position().

Referenced by BzeroReferenceTrajectoryFactory::trajectories(), and ReferenceTrajectoryFactory::trajectories().

{
  // - Should use perp2() instead of perp()
  // - Should not rely on floating point equality, but make a minimal range, e.g. 1.e-6 ?
  return ( s1.eta() == s2.eta() ) && ( s1.phi() == s2.phi() ) && ( s1.position().perp() == s2.position().perp() );
}

virtual const ReferenceTrajectoryCollection TrajectoryFactoryBase::trajectories	(	const edm::EventSetup &	setup,
		const ConstTrajTrackPairCollection &	tracks,
		const ExternalPredictionCollection &	external,
		const reco::BeamSpot &	beamSpot
	)		const `[pure virtual]`

Implemented in BzeroReferenceTrajectoryFactory, CombinedTrajectoryFactory, DualBzeroTrajectoryFactory, DualKalmanFactory, DualTrajectoryFactory, ReferenceTrajectoryFactory, and TwoBodyDecayTrajectoryFactory.

virtual const ReferenceTrajectoryCollection TrajectoryFactoryBase::trajectories	(	const edm::EventSetup &	setup,
		const ConstTrajTrackPairCollection &	tracks,
		const reco::BeamSpot &	beamSpot
	)		const `[pure virtual]`

Implemented in BzeroReferenceTrajectoryFactory, CombinedTrajectoryFactory, DualBzeroTrajectoryFactory, DualKalmanFactory, DualTrajectoryFactory, ReferenceTrajectoryFactory, and TwoBodyDecayTrajectoryFactory.

Referenced by MillePedeAlignmentAlgorithm::run().

bool TrajectoryFactoryBase::useRecHit ( const TransientTrackingRecHit::ConstRecHitPointer & hitPtr ) const [protected]

Definition at line 94 of file TrajectoryFactoryBase.cc.

References theUseInvalidHits, theUseProjectedHits, and theUseWithoutDet.

Referenced by innermostStateAndRecHits(), DualBzeroTrajectoryFactory::referenceStateAndRecHits(), DualTrajectoryFactory::referenceStateAndRecHits(), and DualKalmanFactory::referenceStateAndRecHits().

{
  const GeomDet* det = hitPtr->det();
  if ( !det && !theUseWithoutDet ) return false;
  
  if ( !( theUseInvalidHits || hitPtr->isValid() ) ) return false;

  if ( !theUseProjectedHits )
  {
    const ProjectedRecHit2D* projectedHit = dynamic_cast< const ProjectedRecHit2D* >( hitPtr.get() );
    if ( projectedHit != 0 ) return false;
  }

  return true;
}