#include <GsfConstraintAtVertex.h>

Public Member Functions
TrajectoryStateOnSurface	constrainAtBeamSpot (const reco::GsfTrack &, const reco::BeamSpot &) const
	(multi)TSOS after including the beamspot More...

TrajectoryStateOnSurface	constrainAtPoint (const reco::GsfTrack &, const GlobalPoint &globalPosition, const GlobalError &globalError) const
	(multi)TSOS after including a point with covariance matrix More...

TrajectoryStateOnSurface	constrainAtVertex (const reco::GsfTrack &, const reco::Vertex &) const
	(multi)TSOS after include a vertex More...

	GsfConstraintAtVertex (const edm::EventSetup &)

	~GsfConstraintAtVertex ()

Private Attributes
const TrackerGeometry *	geometry_

GsfPropagatorAdapter *	gsfPropagator_

GsfMultiStateUpdator	gsfUpdator_

const MagneticField *	magField_

MultiTrajectoryStateTransform	multiStateTransformer_

TransverseImpactPointExtrapolator *	tipExtrapolator_

Detailed Description

Definition at line 21 of file GsfConstraintAtVertex.h.

Constructor & Destructor Documentation

◆ GsfConstraintAtVertex()

GsfConstraintAtVertex::GsfConstraintAtVertex ( const edm::EventSetup & setup )

explicit

Definition at line 20 of file GsfConstraintAtVertex.cc.

                                                                        {
   edm::ESHandle<TrackerGeometry> geometryHandle;
   setup.get<TrackerDigiGeometryRecord>().get(geometryHandle);
   geometry_ = geometryHandle.product();
  
   edm::ESHandle<MagneticField> magFieldHandle;
   setup.get<IdealMagneticFieldRecord>().get(magFieldHandle);
   magField_ = magFieldHandle.product();
  
   //   edm::ESHandle<Propagator> propagatorHandle;
   //   setup.get<TrackingComponentsRecord>().get(propagatorName_,propagatorHandle);
   //   propagator_ = propagatorHandle.product();
   gsfPropagator_ = new GsfPropagatorAdapter(AnalyticalPropagator(magField_, anyDirection));
   tipExtrapolator_ = new TransverseImpactPointExtrapolator(*gsfPropagator_);
 }

References AnalyticalPropagator_cfi::AnalyticalPropagator, anyDirection, geometry_, get, gsfPropagator_, magField_, edm::ESHandle< T >::product(), singleTopDQM_cfi::setup, and tipExtrapolator_.

◆ ~GsfConstraintAtVertex()

GsfConstraintAtVertex::~GsfConstraintAtVertex ( )

Definition at line 36 of file GsfConstraintAtVertex.cc.

                                               {
   delete tipExtrapolator_;
   delete gsfPropagator_;
 }

References gsfPropagator_, and tipExtrapolator_.

Member Function Documentation

◆ constrainAtBeamSpot()

TrajectoryStateOnSurface GsfConstraintAtVertex::constrainAtBeamSpot	(	const reco::GsfTrack &	track,
		const reco::BeamSpot &	beamSpot
	)		const

(multi)TSOS after including the beamspot

Definition at line 41 of file GsfConstraintAtVertex.cc.

                                                                                                         {
   //
   // Beamspot (global co-ordinates)
   //
   GlobalPoint bsPosGlobal(beamSpot.x0(), beamSpot.y0(), beamSpot.z0());
   GlobalError bsCovGlobal(beamSpot.rotatedCovariance3D());
   //
   return constrainAtPoint(track, bsPosGlobal, bsCovGlobal);
 }

References pwdgSkimBPark_cfi::beamSpot, constrainAtPoint(), and HLT_2018_cff::track.

Referenced by GsfElectronAlgo::ElectronData::calculateTSOS().

◆ constrainAtPoint()

TrajectoryStateOnSurface GsfConstraintAtVertex::constrainAtPoint	(	const reco::GsfTrack &	track,
		const GlobalPoint &	globalPosition,
		const GlobalError &	globalError
	)		const

(multi)TSOS after including a point with covariance matrix

Definition at line 63 of file GsfConstraintAtVertex.cc.

                                                                                                        {
   //
   // Track on TIP plane
   //
   TrajectoryStateOnSurface innerState = multiStateTransformer_.innerStateOnSurface(track, *geometry_, magField_);
   if (!innerState.isValid())
     return TrajectoryStateOnSurface();
   TrajectoryStateOnSurface tipState = tipExtrapolator_->extrapolate(innerState, globalPosition);
   if (!tipState.isValid())
     return TrajectoryStateOnSurface();
   //
   // RecHit from beam spot
   //
   LocalError bsCovLocal = ErrorFrameTransformer().transform(globalError, tipState.surface());
   auto bsHit =
       TRecHit2DPosConstraint::build(tipState.surface().toLocal(globalPosition), bsCovLocal, &tipState.surface());
   //
   // update with constraint
   //
   TrajectoryStateOnSurface updatedState = gsfUpdator_.update(tipState, *bsHit);
   if (!updatedState.isValid()) {
     edm::LogWarning("GsfConstraintAtVertex") << " GSF update with vertex constraint failed";
     return TrajectoryStateOnSurface();
   }
  
   return updatedState;
 }

References TRecHit2DPosConstraint::build(), TransverseImpactPointExtrapolator::extrapolate(), geometry_, gsfUpdator_, MultiTrajectoryStateTransform::innerStateOnSurface(), TrajectoryStateOnSurface::isValid(), magField_, multiStateTransformer_, TrajectoryStateOnSurface::surface(), tipExtrapolator_, GloballyPositioned< T >::toLocal(), HLT_2018_cff::track, ErrorFrameTransformer::transform(), GsfMultiStateUpdator::update(), and TrackInfoProducer_cfi::updatedState.

Referenced by constrainAtBeamSpot(), and constrainAtVertex().

◆ constrainAtVertex()

TrajectoryStateOnSurface GsfConstraintAtVertex::constrainAtVertex	(	const reco::GsfTrack &	track,
		const reco::Vertex &	vertex
	)		const

(multi)TSOS after include a vertex

Definition at line 52 of file GsfConstraintAtVertex.cc.

                                                                                                   {
   //
   // Beamspot (global co-ordinates)
   //
   GlobalPoint vtxPosGlobal(vertex.position().x(), vertex.position().y(), vertex.position().z());
   GlobalError vtxCovGlobal(vertex.covariance());
   //
   return constrainAtPoint(track, vtxPosGlobal, vtxCovGlobal);
 }