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

Inheritance diagram for DualTrajectoryFactory:

Classes
struct	DualTrajectoryInput

Public Member Functions
DualTrajectoryFactory *	clone () const override

	DualTrajectoryFactory (const edm::ParameterSet &config)

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

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

	~DualTrajectoryFactory () override

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

MaterialEffects	materialEffects (void) const

PropagationDirection	propagationDirection (void) const

unsigned int	tracksPerTrajectory () const

	TrajectoryFactoryBase (const edm::ParameterSet &config)

	TrajectoryFactoryBase (const edm::ParameterSet &config, unsigned int tracksPerTrajectory)

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
const bool	allowZeroMaterial_

const bool	includeAPEs_

const bool	useBeamSpot_

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 17 of file DualTrajectoryFactory.cc.

Constructor & Destructor Documentation

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

Definition at line 54 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 )

override

Definition at line 60 of file DualTrajectoryFactory.cc.

60 {}

Member Function Documentation

DualTrajectoryFactory* DualTrajectoryFactory::clone ( void ) const

inlineoverridevirtual

Implements TrajectoryFactoryBase.

Definition at line 32 of file DualTrajectoryFactory.cc.

References DualTrajectoryFactory().

32 { return new DualTrajectoryFactory(*this); }

DualTrajectoryFactory::DualTrajectoryFactory

DualTrajectoryFactory(const edm::ParameterSet &config)

Definition: DualTrajectoryFactory.cc:54

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

protected

Definition at line 206 of file DualTrajectoryFactory.cc.

References anyDirection, DEFINE_EDM_PLUGIN, AnalyticalPropagator::propagateWithPath(), and TrackCandidateProducer_cfi::propagator.

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 162 of file DualTrajectoryFactory.cc.

References DualTrajectoryFactory::DualTrajectoryInput::bwdRecHits, DualTrajectoryFactory::DualTrajectoryInput::fwdRecHits, input, TrajectoryFactoryBase::orderedTrajectoryMeasurements(), DualTrajectoryFactory::DualTrajectoryInput::refTsos, groupFilesInBlocks::reverse, 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 62 of file DualTrajectoryFactory.cc.

References ReferenceTrajectoryBase::Config::allowZeroMaterial, TrajectoryFactoryBase::allowZeroMaterial_, pwdgSkimBPark_cfi::beamSpot, DualTrajectoryFactory::DualTrajectoryInput::bwdRecHits, looper::config, DualTrajectoryFactory::DualTrajectoryInput::fwdRecHits, edm::EventSetup::get(), ReferenceTrajectoryBase::Config::includeAPEs, TrajectoryFactoryBase::includeAPEs_, input, MagneticField::inTesla(), TrajectoryStateOnSurface::isValid(), PV3DBase< T, PVType, FrameType >::mag2(), HLT_2018_cff::magneticField, TrajectoryFactoryBase::materialEffects(), edm::ESHandle< T >::product(), TrajectoryFactoryBase::propagationDirection(), referenceStateAndRecHits(), DualTrajectoryFactory::DualTrajectoryInput::refTsos, theMass, ReferenceTrajectoryBase::Config::useBeamSpot, TrajectoryFactoryBase::useBeamSpot_, and PV3DBase< T, PVType, FrameType >::z().

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()) {
       ReferenceTrajectoryBase::Config config(materialEffects(), propagationDirection(), theMass);
       config.useBeamSpot = useBeamSpot_;
       config.includeAPEs = includeAPEs_;
       config.allowZeroMaterial = allowZeroMaterial_;
       ReferenceTrajectoryPtr ptr(new DualReferenceTrajectory(
           input.refTsos, input.fwdRecHits, input.bwdRecHits, magneticField.product(), beamSpot, config));
       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 96 of file DualTrajectoryFactory.cc.

References ReferenceTrajectoryBase::Config::allowZeroMaterial, TrajectoryFactoryBase::allowZeroMaterial_, pwdgSkimBPark_cfi::beamSpot, DualTrajectoryFactory::DualTrajectoryInput::bwdRecHits, looper::config, DualTrajectoryFactory::DualTrajectoryInput::fwdRecHits, edm::EventSetup::get(), ReferenceTrajectoryBase::Config::includeAPEs, TrajectoryFactoryBase::includeAPEs_, input, MagneticField::inTesla(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localError(), PV3DBase< T, PVType, FrameType >::mag2(), HLT_2018_cff::magneticField, TrajectoryFactoryBase::materialEffects(), LocalTrajectoryError::matrix(), edm::ESHandle< T >::product(), propagateExternal(), TrajectoryFactoryBase::propagationDirection(), referenceStateAndRecHits(), DualTrajectoryFactory::DualTrajectoryInput::refTsos, TrajectoryStateOnSurface::surface(), theMass, trajectories(), ReferenceTrajectoryBase::Config::useBeamSpot, TrajectoryFactoryBase::useBeamSpot_, and PV3DBase< T, PVType, FrameType >::z().

                                         {
   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;
 
         ReferenceTrajectoryBase::Config config(materialEffects(), propagationDirection(), theMass);
         config.useBeamSpot = useBeamSpot_;
         config.includeAPEs = includeAPEs_;
         config.allowZeroMaterial = allowZeroMaterial_;
         ReferenceTrajectoryPtr ptr(new DualReferenceTrajectory(
             propExternal, input.fwdRecHits, input.bwdRecHits, magneticField.product(), beamSpot, config));
 
         AlgebraicSymMatrix externalParamErrors(asHepMatrix<5>(propExternal.localError().matrix()));
         ptr->setParameterErrors(externalParamErrors);
         trajectories.push_back(ptr);
       } else {
         ReferenceTrajectoryBase::Config config(materialEffects(), propagationDirection(), theMass);
         config.useBeamSpot = useBeamSpot_;
         config.includeAPEs = includeAPEs_;
         config.allowZeroMaterial = allowZeroMaterial_;
         ReferenceTrajectoryPtr ptr(new DualReferenceTrajectory(
             input.refTsos, input.fwdRecHits, input.bwdRecHits, magneticField.product(), beamSpot, config));
         trajectories.push_back(ptr);
       }
     }
 
     ++itTracks;
     ++itExternal;
   }
 
   return trajectories;
 }

Member Data Documentation

double DualTrajectoryFactory::theMass

protected

Definition at line 47 of file DualTrajectoryFactory.cc.

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