A factory that produces instances of class ReferenceTrajectory from a given TrajTrackPairCollection. More...

Inheritance diagram for DualTrajectoryFactory:

Classes
struct	DualTrajectoryInput

Public Member Functions
virtual DualTrajectoryFactory *	clone () const override

	DualTrajectoryFactory (const edm::ParameterSet &config)

virtual const ReferenceTrajectoryCollection	trajectories (const edm::EventSetup &setup, const ConstTrajTrackPairCollection &tracks, const reco::BeamSpot &beamSpot) const override
	Produce the reference trajectories. More...

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

virtual	~DualTrajectoryFactory ()

Public Member Functions inherited from TrajectoryFactoryBase
const edm::ParameterSet &	configuration () const

MaterialEffects	materialEffects (void) const

PropagationDirection	propagationDirection (void) const

	TrajectoryFactoryBase (const edm::ParameterSet &config)

virtual	~TrajectoryFactoryBase (void)

Protected Member Functions
const TrajectoryStateOnSurface	propagateExternal (const TrajectoryStateOnSurface &external, const Surface &surface, const MagneticField *magField) const

const DualTrajectoryInput	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
double	theMass

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 instances of class ReferenceTrajectory from a given TrajTrackPairCollection.

Definition at line 18 of file DualTrajectoryFactory.cc.

Constructor & Destructor Documentation

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

Definition at line 57 of file DualTrajectoryFactory.cc.

References theMass.

Referenced by clone().

                                                                              :
   TrajectoryFactoryBase( config ),
   theMass(config.getParameter<double>("ParticleMass"))
 {
   edm::LogInfo("Alignment") << "@SUB=DualTrajectoryFactory"
                             << "mass: " << theMass;
 }

DualTrajectoryFactory::~DualTrajectoryFactory ( void )

virtual

Definition at line 66 of file DualTrajectoryFactory.cc.

66 {}

Member Function Documentation

virtual DualTrajectoryFactory* DualTrajectoryFactory::clone ( void ) const

inlineoverridevirtual

Implements TrajectoryFactoryBase.

Definition at line 34 of file DualTrajectoryFactory.cc.

References DualTrajectoryFactory().

34 { return new DualTrajectoryFactory(*this); }

DualTrajectoryFactory::DualTrajectoryFactory

DualTrajectoryFactory(const edm::ParameterSet &config)

Definition: DualTrajectoryFactory.cc:57

const TrajectoryStateOnSurface DualTrajectoryFactory::propagateExternal	(	const TrajectoryStateOnSurface &	external,
		const Surface &	surface,
		const MagneticField *	magField
	)		const

protected

Definition at line 236 of file DualTrajectoryFactory.cc.

References AnalyticalPropagator_cfi::AnalyticalPropagator, and anyDirection.

Referenced by trajectories().

 {
   AnalyticalPropagator propagator( magField, anyDirection );
   const std::pair< TrajectoryStateOnSurface, double > tsosWithPath =
     propagator.propagateWithPath( external, surface );
   return tsosWithPath.first;
 }

const DualTrajectoryFactory::DualTrajectoryInput DualTrajectoryFactory::referenceStateAndRecHits ( const ConstTrajTrackPair & track ) const

protected

Definition at line 189 of file DualTrajectoryFactory.cc.

References DualTrajectoryFactory::DualTrajectoryInput::bwdRecHits, DualTrajectoryFactory::DualTrajectoryInput::fwdRecHits, input, TrajectoryFactoryBase::orderedTrajectoryMeasurements(), DualTrajectoryFactory::DualTrajectoryInput::refTsos, and TrajectoryFactoryBase::useRecHit().

Referenced by trajectories().

 {
   DualTrajectoryInput input;
  
   // get the trajectory measurements in the correct order, i.e. reverse if needed
   Trajectory::DataContainer allTrajMeas = this->orderedTrajectoryMeasurements( *track.first );
   Trajectory::DataContainer usedTrajMeas;
   Trajectory::DataContainer::iterator itM;
   // get all relevant trajectory measurements
   for ( itM = allTrajMeas.begin(); itM != allTrajMeas.end(); itM++ )
   {
     if ( useRecHit( ( *itM ).recHit() ) ) usedTrajMeas.push_back( *itM );
   }
 
   unsigned int iMeas = 0;
   unsigned int nMeas = usedTrajMeas.size();
   unsigned int nRefStateMeas = nMeas/2;
   // get the valid RecHits
   for ( itM = usedTrajMeas.begin(); itM != usedTrajMeas.end(); itM++, iMeas++ )
   {
     TransientTrackingRecHit::ConstRecHitPointer aRecHit = ( *itM ).recHit();
 
     if ( iMeas < nRefStateMeas ) {
       input.bwdRecHits.push_back( aRecHit );
     } else if ( iMeas > nRefStateMeas ) {
       input.fwdRecHits.push_back( aRecHit );
     } else { // iMeas == nRefStateMeas
       if ( ( *itM ).updatedState().isValid() )
       {
     input.refTsos = ( *itM ).updatedState();
     input.bwdRecHits.push_back( aRecHit );
     input.fwdRecHits.push_back( aRecHit );
       } else {
     // if the tsos of the middle hit is not valid, try the next one ...
     nRefStateMeas++;
     input.bwdRecHits.push_back( aRecHit );
       }
     }
   }
 
   // bring input.fwdRecHits into correct order
   std::reverse( input.bwdRecHits.begin(), input.bwdRecHits.end() );
 
   return input;
 }

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

overridevirtual

Produce the reference trajectories.

Implements TrajectoryFactoryBase.

Definition at line 70 of file DualTrajectoryFactory.cc.

References SiPixelRawToDigiRegional_cfi::beamSpot, DualTrajectoryFactory::DualTrajectoryInput::bwdRecHits, DualTrajectoryFactory::DualTrajectoryInput::fwdRecHits, edm::EventSetup::get(), input, TrajectoryStateOnSurface::isValid(), TrajectoryFactoryBase::materialEffects(), edm::ESHandle< class >::product(), TrajectoryFactoryBase::propagationDirection(), referenceStateAndRecHits(), DualTrajectoryFactory::DualTrajectoryInput::refTsos, theMass, and TrajectoryFactoryBase::theUseBeamSpot.

Referenced by trajectories().

 {
   ReferenceTrajectoryCollection trajectories;
 
   edm::ESHandle< MagneticField > magneticField;
   setup.get< IdealMagneticFieldRecord >().get( magneticField );
   if (magneticField->inTesla(GlobalPoint(0.,0.,0.)).mag2() < 1.e-6) {
     edm::LogWarning("Alignment") << "@SUB=DualTrajectoryFactory::trajectories"
                                  << "B-field in z is " << magneticField->inTesla(GlobalPoint(0.,0.,0.)).z()
                                  << ": You should probably use the DualBzeroTrajectoryFactory\n" 
                                  << "or fix this code to switch automatically as the ReferenceTrajectoryFactory does.";
   }
 
   ConstTrajTrackPairCollection::const_iterator itTracks = tracks.begin();
 
   while ( itTracks != tracks.end() )
   { 
     const DualTrajectoryInput input = this->referenceStateAndRecHits( *itTracks );
     // Check input: If all hits were rejected, the TSOS is initialized as invalid.
     if ( input.refTsos.isValid() )
     {
       ReferenceTrajectoryPtr ptr( new DualReferenceTrajectory( input.refTsos,
                                    input.fwdRecHits,
                                    input.bwdRecHits,
                                    magneticField.product(),
                                    materialEffects(),
                                    propagationDirection(),
                                    theMass,
                                    theUseBeamSpot, beamSpot ) );
       trajectories.push_back( ptr );
     }
 
     ++itTracks;
   }
 
   return trajectories;
 }

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

overridevirtual

Implements TrajectoryFactoryBase.

Definition at line 111 of file DualTrajectoryFactory.cc.

References SiPixelRawToDigiRegional_cfi::beamSpot, DualTrajectoryFactory::DualTrajectoryInput::bwdRecHits, DualTrajectoryFactory::DualTrajectoryInput::fwdRecHits, edm::EventSetup::get(), input, TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localError(), TrajectoryFactoryBase::materialEffects(), LocalTrajectoryError::matrix(), edm::ESHandle< class >::product(), propagateExternal(), TrajectoryFactoryBase::propagationDirection(), referenceStateAndRecHits(), DualTrajectoryFactory::DualTrajectoryInput::refTsos, TrajectoryStateOnSurface::surface(), theMass, TrajectoryFactoryBase::theUseBeamSpot, and trajectories().

 {
   ReferenceTrajectoryCollection trajectories;
 
   if ( tracks.size() != external.size() )
   {
     edm::LogInfo("ReferenceTrajectories") << "@SUB=DualTrajectoryFactory::trajectories"
                       << "Inconsistent input:\n"
                       << "\tnumber of tracks = " << tracks.size()
                       << "\tnumber of external predictions = " << external.size();
     return trajectories;
   }
 
   edm::ESHandle< MagneticField > magneticField;
   setup.get< IdealMagneticFieldRecord >().get( magneticField );
   if (magneticField->inTesla(GlobalPoint(0.,0.,0.)).mag2() < 1.e-6) {
     edm::LogWarning("Alignment") << "@SUB=DualTrajectoryFactory::trajectories"
                                  << "B-field in z is " << magneticField->inTesla(GlobalPoint(0.,0.,0.)).z()
                                  << ": You should probably use the DualBzeroTrajectoryFactory\n" 
                                  << "or fix this code to switch automatically as the ReferenceTrajectoryFactory does.";
   }
 
   ConstTrajTrackPairCollection::const_iterator itTracks = tracks.begin();
   ExternalPredictionCollection::const_iterator itExternal = external.begin();
 
   while ( itTracks != tracks.end() )
   {
     const DualTrajectoryInput input = referenceStateAndRecHits( *itTracks );
     // Check input: If all hits were rejected, the TSOS is initialized as invalid.
     if ( input.refTsos.isValid() )
     {
       if ( (*itExternal).isValid() )
       {
     TrajectoryStateOnSurface propExternal =
       propagateExternal( *itExternal, input.refTsos.surface(), magneticField.product() );
 
     if ( !propExternal.isValid() ) continue;
 
     // set the flag for reversing the RecHits to false, since they are already in the correct order.
     ReferenceTrajectoryPtr ptr( new DualReferenceTrajectory( propExternal,
                                  input.fwdRecHits,
                                  input.bwdRecHits,
                                  magneticField.product(),
                                  materialEffects(),
                                  propagationDirection(),
                                  theMass, 
                                  theUseBeamSpot, beamSpot ) );
 
     AlgebraicSymMatrix externalParamErrors( asHepMatrix<5>( propExternal.localError().matrix() ) );
     ptr->setParameterErrors( externalParamErrors );
     trajectories.push_back( ptr );
       }
       else
       {
     ReferenceTrajectoryPtr ptr( new DualReferenceTrajectory( input.refTsos,
                                  input.fwdRecHits,
                                  input.bwdRecHits,
                                  magneticField.product(),
                                  materialEffects(),
                                  propagationDirection(),
                                  theMass, 
                                  theUseBeamSpot, beamSpot ) );
     trajectories.push_back( ptr );
       }
     }
 
     ++itTracks;
     ++itExternal;
   }
 
   return trajectories;
 }

Member Data Documentation

double DualTrajectoryFactory::theMass

protected

Definition at line 50 of file DualTrajectoryFactory.cc.

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

Classes

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation