DualKalmanFactory Class Reference

Inheritance diagram for DualKalmanFactory:

Classes
struct	DualKalmanInput

Public Member Functions
virtual DualKalmanFactory *	clone () const
	DualKalmanFactory (const edm::ParameterSet &config)
virtual const ReferenceTrajectoryCollection	trajectories (const edm::EventSetup &setup, const ConstTrajTrackPairCollection &tracks, const reco::BeamSpot &beamSpot) const
	Produce the reference trajectories. More...
virtual const ReferenceTrajectoryCollection	trajectories (const edm::EventSetup &setup, const ConstTrajTrackPairCollection &tracks, const ExternalPredictionCollection &external, const reco::BeamSpot &beamSpot) const
virtual	~DualKalmanFactory ()
Public Member Functions inherited from TrajectoryFactoryBase
MaterialEffects	materialEffects (void) const
PropagationDirection	propagationDirection (void) const
	TrajectoryFactoryBase (const edm::ParameterSet &config)
virtual	~TrajectoryFactoryBase (void)

Protected Member Functions
const DualKalmanInput	referenceStateAndRecHits (const ConstTrajTrackPair &track) const
Protected Member Functions inherited from TrajectoryFactoryBase
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
const double	theMass
const int	theResidMethod
Protected Attributes inherited from TrajectoryFactoryBase
bool	theUseBeamSpot

Additional Inherited Members
Public Types inherited from TrajectoryFactoryBase
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

Detailed Description

A factory that produces reference trajectory instances of class DualKalmanTrajectory from a given TrajTrackPairCollection.

Currently two methods to set residual and error can be choosen via cfg: 1: the unbiased residal approach 2: the pull approach

Author

: Gero Flucke date : October 2008

Revision:

1.4

Date:

2011/12/19 14:40:01

(last update by

Author:

mussgill

)

Definition at line 32 of file DualKalmanFactory.cc.

Constructor & Destructor Documentation

DualKalmanFactory::DualKalmanFactory

(

const edm::ParameterSet &

config

)

Definition at line 75 of file DualKalmanFactory.cc.

Referenced by clone().

  : TrajectoryFactoryBase(config), theMass(config.getParameter<double>("ParticleMass")),
    theResidMethod(config.getParameter<int>("ResidualMethod"))
{
  // Since theResidMethod is passed to DualKalmanTrajectory, valid values are checked there.
  edm::LogInfo("Alignment") << "@SUB=DualKalmanFactory" << "Factory created.";
}

DualKalmanFactory::~DualKalmanFactory ( )

virtual

Definition at line 85 of file DualKalmanFactory.cc.

{}

Member Function Documentation

virtual DualKalmanFactory* DualKalmanFactory::clone ( void  ) const

inlinevirtual

Implements TrajectoryFactoryBase.

Definition at line 49 of file DualKalmanFactory.cc.

References DualKalmanFactory().

{ return new DualKalmanFactory(*this); }

const DualKalmanFactory::DualKalmanInput DualKalmanFactory::referenceStateAndRecHits ( const ConstTrajTrackPair &  track ) const

protected

Definition at line 139 of file DualKalmanFactory.cc.

References DualKalmanFactory::DualKalmanInput::bwdRecHitNums, DualKalmanFactory::DualKalmanInput::fwdRecHitNums, LaserDQM_cfg::input, TrajectoryFactoryBase::orderedTrajectoryMeasurements(), DualKalmanFactory::DualKalmanInput::refTsos, DualKalmanFactory::DualKalmanInput::trajMeasurements, and TrajectoryFactoryBase::useRecHit().

Referenced by trajectories().

{
  // Same idea as in DualTrajectoryFactory::referenceStateAndRecHits(..):
  // Split trajectory in the middle, take middle as reference and provide first
  // and second half of hits, each starting from this middle hit.
  // In contrast to DualTrajectoryFactory we deal here with indices and not rechits directly
  // to be able to get measurements and uncertainties later from the Trajectory that is
  // provided by the Kalman track fit.

  DualKalmanInput input;
 
  // get the trajectory measurements in the correct order, i.e. reverse if needed
  input.trajMeasurements = this->orderedTrajectoryMeasurements(*track.first);

  // get indices of relevant trajectory measurements to find middle of them
  std::vector<unsigned int> usedTrajMeasNums;
  for (unsigned int iM = 0; iM < input.trajMeasurements.size(); ++iM) {
    if (this->useRecHit(input.trajMeasurements[iM].recHit())) usedTrajMeasNums.push_back(iM);
  }
  unsigned int nRefStateMeas = usedTrajMeasNums.size()/2;

  // get the valid RecHits numbers
  for (unsigned int iMeas = 0; iMeas < usedTrajMeasNums.size(); ++iMeas) {
    if (iMeas < nRefStateMeas) {
      input.bwdRecHitNums.push_back(usedTrajMeasNums[iMeas]);
    } else if (iMeas > nRefStateMeas) {
      input.fwdRecHitNums.push_back(usedTrajMeasNums[iMeas]);
    } else { // iMeas == nRefStateMeas
      if (input.trajMeasurements[usedTrajMeasNums[iMeas]].updatedState().isValid()) {
    input.refTsos = input.trajMeasurements[usedTrajMeasNums[iMeas]].updatedState();
    input.bwdRecHitNums.push_back(usedTrajMeasNums[iMeas]);
    input.fwdRecHitNums.push_back(usedTrajMeasNums[iMeas]);
      } else {
    // if the hit/tsos of the middle hit is not valid, try the next one...
    ++nRefStateMeas; // but keep hit if only TSOS bad
    input.bwdRecHitNums.push_back(usedTrajMeasNums[iMeas]);
      }
    }
  }

  // bring input.fwdRecHits into correct order
  std::reverse(input.bwdRecHitNums.begin(), input.bwdRecHitNums.end());

  return input;
}

const DualKalmanFactory::ReferenceTrajectoryCollection DualKalmanFactory::trajectories ( const edm::EventSetup &  setup, const ConstTrajTrackPairCollection &  tracks, const reco::BeamSpot &  beamSpot  ) const

virtual

Produce the reference trajectories.

Implements TrajectoryFactoryBase.

Definition at line 90 of file DualKalmanFactory.cc.

References SiPixelRawToDigiRegional_cfi::beamSpot, DualKalmanFactory::DualKalmanInput::bwdRecHitNums, DualKalmanFactory::DualKalmanInput::fwdRecHitNums, edm::EventSetup::get(), LaserDQM_cfg::input, TrajectoryStateOnSurface::isValid(), TrajectoryFactoryBase::materialEffects(), edm::ESHandle< class >::product(), TrajectoryFactoryBase::propagationDirection(), referenceStateAndRecHits(), DualKalmanFactory::DualKalmanInput::refTsos, theMass, theResidMethod, TrajectoryFactoryBase::theUseBeamSpot, and DualKalmanFactory::DualKalmanInput::trajMeasurements.

Referenced by trajectories().

{
  ReferenceTrajectoryCollection trajectories;

  edm::ESHandle<MagneticField> magneticField;
  setup.get<IdealMagneticFieldRecord>().get(magneticField);

  ConstTrajTrackPairCollection::const_iterator itTracks = tracks.begin();

  while (itTracks != tracks.end()) { 
    const DualKalmanInput input = this->referenceStateAndRecHits(*itTracks);
    // Check input: If all hits were rejected, the TSOS is initialized as invalid.
    if (input.refTsos.isValid()) {
      ReferenceTrajectoryPtr ptr(new DualKalmanTrajectory(input.trajMeasurements,
                              input.refTsos,
                              input.fwdRecHitNums,
                              input.bwdRecHitNums,
                              magneticField.product(),
                              this->materialEffects(),
                              this->propagationDirection(),
                              theMass, theUseBeamSpot, beamSpot,
                              theResidMethod));
      trajectories.push_back(ptr);
    }
    ++itTracks;
  }

  return trajectories;
}

const DualKalmanFactory::ReferenceTrajectoryCollection DualKalmanFactory::trajectories ( const edm::EventSetup &  setup, const ConstTrajTrackPairCollection &  tracks, const ExternalPredictionCollection &  external, const reco::BeamSpot &  beamSpot  ) const

virtual

Implements TrajectoryFactoryBase.

Definition at line 124 of file DualKalmanFactory.cc.

References trajectories().

{
  ReferenceTrajectoryCollection trajectories;

  edm::LogError("Alignment") << "@SUB=DualKalmanFactory::trajectories" 
                 << "Not implemented with ExternalPrediction.";
  return trajectories;
}

Member Data Documentation

const double DualKalmanFactory::theMass

protected

Definition at line 67 of file DualKalmanFactory.cc.

Referenced by twikiExport.Constituent::__str__(), and trajectories().

const int DualKalmanFactory::theResidMethod

protected

Definition at line 68 of file DualKalmanFactory.cc.

Referenced by trajectories().