#include <TSGFromPropagation.h>

Inheritance diagram for TSGFromPropagation:

Classes
struct	increasingEstimate

struct	isInvalid

Public Member Functions
void	init (const MuonServiceProxy *) override
	initialize More...

void	setEvent (const edm::Event &) override
	set an event More...

void	trackerSeeds (const TrackCand &, const TrackingRegion &, const TrackerTopology *, std::vector< TrajectorySeed > &) override
	generate seed(s) for a track More...

	TSGFromPropagation (const edm::ParameterSet &par, edm::ConsumesCollector &iC, const MuonServiceProxy *)

	TSGFromPropagation (const edm::ParameterSet &pset, edm::ConsumesCollector &iC)
	constructor More...

	~TSGFromPropagation () override
	destructor More...

Public Member Functions inherited from TrackerSeedGenerator
const edm::Event *	getEvent () const

	TrackerSeedGenerator ()

virtual void	trackerSeeds (const TrackCand &, const TrackingRegion &, const TrackerTopology *, BTSeedCollection &)

virtual	~TrackerSeedGenerator ()
	destructor More...

Private Member Functions
void	adjust (FreeTrajectoryState &) const
	adjust the error matrix of the FTS More...

void	adjust (TrajectoryStateOnSurface &) const
	adjust the error matrix of the TSOS More...

TrajectorySeed	createSeed (const TrajectoryStateOnSurface &, const DetId &) const
	create a hitless seed from a trajectory state More...

TrajectorySeed	createSeed (const TrajectoryStateOnSurface &tsos, const edm::OwnVector< TrackingRecHit > &container, const DetId &id) const
	create a seed from a trajectory state More...

double	dxyDis (const TrajectoryStateOnSurface &tsos) const

const Chi2MeasurementEstimator *	estimator () const

std::vector< TrajectoryMeasurement >	findMeasurements (const DetLayer *, const TrajectoryStateOnSurface &) const
	look for measurements on the first compatible layer More...

void	getRescalingFactor (const TrackCand &staMuon)

TrajectoryStateOnSurface	innerState (const TrackCand &) const

TrajectoryStateOnSurface	outerTkState (const TrackCand &) const

bool	passSelection (const TrajectoryStateOnSurface &) const
	check some quantity and beam-spot compatibility and decide to continue More...

edm::ESHandle< Propagator >	propagator () const

const TrajectoryStateUpdator *	updator () const

void	validMeasurements (std::vector< TrajectoryMeasurement > &) const
	select valid measurements More...

double	zDis (const TrajectoryStateOnSurface &tsos) const

Private Attributes
edm::Handle< reco::BeamSpot >	beamSpot

bool	theAdjustAtIp

edm::InputTag	theBeamSpotInputTag

edm::EDGetTokenT< reco::BeamSpot >	theBeamSpotToken

unsigned long long	theCacheId_MT

unsigned long long	theCacheId_TG

std::string	theCategory

edm::ParameterSet	theConfig

MuonErrorMatrix *	theErrorMatrixAdjuster

const Chi2MeasurementEstimator *	theEstimator

double	theFixedErrorRescaling

double	theFlexErrorRescaling

double	theMaxChi2

edm::ESHandle< MeasurementTracker >	theMeasTracker

edm::Handle< MeasurementTrackerEvent >	theMeasTrackerEvent

std::string	theMeasTrackerName

edm::InputTag	theMeasurementTrackerEventTag

edm::EDGetTokenT< MeasurementTrackerEvent >	theMeasurementTrackerEventToken

const DirectTrackerNavigation *	theNavigation

std::string	thePropagatorName

std::string	theResetMethod

bool	theSelectStateFlag

const MuonServiceProxy *	theService

double	theSigmaZ

edm::ESHandle< GeometricSearchTracker >	theTracker

TrajectoryStateTransform *	theTSTransformer

bool	theUpdateStateFlag

const TrajectoryStateUpdator *	theUpdator

bool	theUseVertexStateFlag

Additional Inherited Members
Public Types inherited from TrackerSeedGenerator
typedef std::vector< TrajectorySeed >	BTSeedCollection

typedef std::pair< const Trajectory *, reco::TrackRef >	TrackCand

Protected Attributes inherited from TrackerSeedGenerator
const edm::Event *	theEvent

const MuonServiceProxy *	theProxyService

Detailed Description

Author: Chang Liu - Purdue University

Tracker Seed Generator by propagating and updating a standAlone muon to the first 2 (or 1) rechits it meets in tracker system

Author: Chang Liu - Purdue University

Definition at line 31 of file TSGFromPropagation.h.

Constructor & Destructor Documentation

◆ TSGFromPropagation() [1/2]

TSGFromPropagation::TSGFromPropagation	(	const edm::ParameterSet &	pset,
		edm::ConsumesCollector &	iC
	)

constructor

Definition at line 29 of file TSGFromPropagation.cc.

30 : TSGFromPropagation(iConfig, iC, nullptr) {}

◆ TSGFromPropagation() [2/2]

TSGFromPropagation::TSGFromPropagation	(	const edm::ParameterSet &	par,
		edm::ConsumesCollector &	iC,
		const MuonServiceProxy *	service
	)

Definition at line 32 of file TSGFromPropagation.cc.

     : theTracker(nullptr),
       theMeasTracker(nullptr),
       theNavigation(nullptr),
       theService(service),
       theUpdator(nullptr),
       theEstimator(nullptr),
       theTSTransformer(nullptr),
       theSigmaZ(0),
       theConfig(iConfig) {
   theCategory = "Muon|RecoMuon|TSGFromPropagation";
   theMeasTrackerName = iConfig.getParameter<std::string>("MeasurementTrackerName");
   theMeasurementTrackerEventTag = iConfig.getParameter<edm::InputTag>("MeasurementTrackerEvent");
   theBeamSpotInputTag = theConfig.getParameter<edm::InputTag>("beamSpot");
   theBeamSpotToken = iC.consumes<reco::BeamSpot>(theBeamSpotInputTag);
   theMeasurementTrackerEventToken = iC.consumes<MeasurementTrackerEvent>(theMeasurementTrackerEventTag);
 }

References edm::ConsumesCollector::consumes(), edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, theBeamSpotInputTag, theBeamSpotToken, theCategory, theConfig, theMeasTrackerName, theMeasurementTrackerEventTag, and theMeasurementTrackerEventToken.

◆ ~TSGFromPropagation()

TSGFromPropagation::~TSGFromPropagation ( )

override

destructor

Definition at line 52 of file TSGFromPropagation.cc.

                                         {
   LogTrace(theCategory) << " TSGFromPropagation dtor called ";
   if (theNavigation)
     delete theNavigation;
   if (theUpdator)
     delete theUpdator;
   if (theEstimator)
     delete theEstimator;
   if (theErrorMatrixAdjuster)
     delete theErrorMatrixAdjuster;
 }

References LogTrace, theCategory, theErrorMatrixAdjuster, theEstimator, theNavigation, and theUpdator.

Member Function Documentation

◆ adjust() [1/2]

void TSGFromPropagation::adjust ( FreeTrajectoryState & state ) const

private

adjust the error matrix of the FTS

Definition at line 356 of file TSGFromPropagation.cc.

                                                                 {
   //rescale the error
   if (theResetMethod == "discreate") {
     state.rescaleError(theFlexErrorRescaling);
     return;
   }
  
   //rescale the error
   if (theResetMethod == "fixed" || !theErrorMatrixAdjuster) {
     state.rescaleError(theFixedErrorRescaling);
     return;
   }
  
   CurvilinearTrajectoryError oMat = state.curvilinearError();
   CurvilinearTrajectoryError sfMat = theErrorMatrixAdjuster->get(state.momentum());  //FIXME with position
   MuonErrorMatrix::multiply(oMat, sfMat);
  
   state = FreeTrajectoryState(state.parameters(), oMat);
 }

References FreeTrajectoryState::curvilinearError(), MuonErrorMatrix::get(), FreeTrajectoryState::momentum(), MuonErrorMatrix::multiply(), FreeTrajectoryState::parameters(), FreeTrajectoryState::rescaleError(), theErrorMatrixAdjuster, theFixedErrorRescaling, theFlexErrorRescaling, and theResetMethod.

Referenced by innerState(), and outerTkState().

◆ adjust() [2/2]

void TSGFromPropagation::adjust ( TrajectoryStateOnSurface & state ) const

private

adjust the error matrix of the TSOS

Definition at line 376 of file TSGFromPropagation.cc.

                                                                      {
   //rescale the error
   if (theResetMethod == "discreate") {
     state.rescaleError(theFlexErrorRescaling);
     return;
   }
  
   if (theResetMethod == "fixed" || !theErrorMatrixAdjuster) {
     state.rescaleError(theFixedErrorRescaling);
     return;
   }
  
   CurvilinearTrajectoryError oMat = state.curvilinearError();
   CurvilinearTrajectoryError sfMat = theErrorMatrixAdjuster->get(state.globalMomentum());  //FIXME with position
   MuonErrorMatrix::multiply(oMat, sfMat);
  
   state =
       TrajectoryStateOnSurface(state.weight(), state.globalParameters(), oMat, state.surface(), state.surfaceSide());
 }

References TrajectoryStateOnSurface::curvilinearError(), MuonErrorMatrix::get(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalParameters(), MuonErrorMatrix::multiply(), TrajectoryStateOnSurface::rescaleError(), TrajectoryStateOnSurface::surface(), TrajectoryStateOnSurface::surfaceSide(), theErrorMatrixAdjuster, theFixedErrorRescaling, theFlexErrorRescaling, theResetMethod, and TrajectoryStateOnSurface::weight().

◆ createSeed() [1/2]

TrajectorySeed TSGFromPropagation::createSeed	(	const TrajectoryStateOnSurface &	tsos,
		const DetId &	id
	)		const

private

create a hitless seed from a trajectory state

Definition at line 274 of file TSGFromPropagation.cc.

                                                                                                          {
   edm::OwnVector<TrackingRecHit> container;
   return createSeed(tsos, container, id);
 }

Referenced by trackerSeeds().

◆ createSeed() [2/2]

TrajectorySeed TSGFromPropagation::createSeed	(	const TrajectoryStateOnSurface &	tsos,
		const edm::OwnVector< TrackingRecHit > &	container,
		const DetId &	id
	)		const

private

create a seed from a trajectory state

Definition at line 279 of file TSGFromPropagation.cc.

                                                                      {
   PTrajectoryStateOnDet const& seedTSOS = trajectoryStateTransform::persistentState(tsos, id.rawId());
   return TrajectorySeed(seedTSOS, container, oppositeToMomentum);
 }

References oppositeToMomentum, and trajectoryStateTransform::persistentState().

◆ dxyDis()

double TSGFromPropagation::dxyDis ( const TrajectoryStateOnSurface & tsos ) const

private

Definition at line 334 of file TSGFromPropagation.cc.

                                                                             {
   return fabs(
       (-tsos.globalPosition().x() * tsos.globalMomentum().y() + tsos.globalPosition().y() * tsos.globalMomentum().x()) /
       tsos.globalMomentum().perp());
 }

References TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

◆ estimator()

const Chi2MeasurementEstimator* TSGFromPropagation::estimator ( void ) const

inlineprivate

Definition at line 60 of file TSGFromPropagation.h.

60 { return theEstimator; }

References theEstimator.

Referenced by findMeasurements(), and trackerSeeds().

◆ findMeasurements()

std::vector< TrajectoryMeasurement > TSGFromPropagation::findMeasurements	(	const DetLayer *	nl,
		const TrajectoryStateOnSurface &	staState
	)		const

private

look for measurements on the first compatible layer

Definition at line 292 of file TSGFromPropagation.cc.

                                                                         {
   std::vector<TrajectoryMeasurement> result;
  
   std::vector<DetLayer::DetWithState> compatDets = nl->compatibleDets(staState, *propagator(), *estimator());
   if (compatDets.empty())
     return result;
  
   for (std::vector<DetLayer::DetWithState>::const_iterator idws = compatDets.begin(); idws != compatDets.end();
        ++idws) {
     if (idws->second.isValid() && (idws->first)) {
       std::vector<TrajectoryMeasurement> tmptm =
           theMeasTrackerEvent->idToDet(idws->first->geographicalId())
               .fastMeasurements(idws->second, idws->second, *propagator(), *estimator());
       validMeasurements(tmptm);
       //         if ( tmptm.size() > 2 ) {
       //            std::stable_sort(tmptm.begin(),tmptm.end(),increasingEstimate());
       //            result.insert(result.end(),tmptm.begin(), tmptm.begin()+2);
       //         } else {
       result.insert(result.end(), tmptm.begin(), tmptm.end());
       //         }
     }
   }
  
   return result;
 }

References GeometricSearchDet::compatibleDets(), estimator(), MeasurementTrackerEvent::idToDet(), propagator(), mps_fire::result, theMeasTrackerEvent, and validMeasurements().

Referenced by trackerSeeds().

◆ getRescalingFactor()

void TSGFromPropagation::getRescalingFactor ( const TrackCand & staMuon )

private

Definition at line 345 of file TSGFromPropagation.cc.

                                                                     {
   float pt = (staMuon.second)->pt();
   if (pt < 13.0)
     theFlexErrorRescaling = 3;
   else if (pt < 30.0)
     theFlexErrorRescaling = 5;
   else
     theFlexErrorRescaling = 10;
   return;
 }

References DiDispStaMuonMonitor_cfi::pt, and theFlexErrorRescaling.

Referenced by trackerSeeds().

◆ init()

void TSGFromPropagation::init ( const MuonServiceProxy * service )

overridevirtual

initialize

Reimplemented from TrackerSeedGenerator.

Definition at line 151 of file TSGFromPropagation.cc.

                                                              {
   theMaxChi2 = theConfig.getParameter<double>("MaxChi2");
  
   theFixedErrorRescaling = theConfig.getParameter<double>("ErrorRescaling");
  
   theFlexErrorRescaling = 1.0;
  
   theResetMethod = theConfig.getParameter<std::string>("ResetMethod");
  
   if (theResetMethod != "discrete" && theResetMethod != "fixed" && theResetMethod != "matrix") {
     edm::LogError("TSGFromPropagation") << "Wrong error rescaling method: " << theResetMethod << "\n"
                                         << "Possible choices are: discrete, fixed, matrix.\n"
                                         << "Use discrete method" << std::endl;
     theResetMethod = "discrete";
   }
  
   theEstimator = new Chi2MeasurementEstimator(theMaxChi2);
  
   theCacheId_MT = 0;
  
   theCacheId_TG = 0;
  
   thePropagatorName = theConfig.getParameter<std::string>("Propagator");
  
   theService = service;
  
   theUseVertexStateFlag = theConfig.getParameter<bool>("UseVertexState");
  
   theUpdateStateFlag = theConfig.getParameter<bool>("UpdateState");
  
   theSelectStateFlag = theConfig.getParameter<bool>("SelectState");
  
   theUpdator = new KFUpdator();
  
   theSigmaZ = theConfig.getParameter<double>("SigmaZ");
  
   //theBeamSpotInputTag = theConfig.getParameter<edm::InputTag>("beamSpot");
  
   edm::ParameterSet errorMatrixPset = theConfig.getParameter<edm::ParameterSet>("errorMatrixPset");
   if (theResetMethod == "matrix" && !errorMatrixPset.empty()) {
     theAdjustAtIp = errorMatrixPset.getParameter<bool>("atIP");
     theErrorMatrixAdjuster = new MuonErrorMatrix(errorMatrixPset);
   } else {
     theAdjustAtIp = false;
     theErrorMatrixAdjuster = nullptr;
   }
  
   theService->eventSetup().get<TrackerRecoGeometryRecord>().get(theTracker);
   theNavigation = new DirectTrackerNavigation(theTracker);
 }

References Chi2MeasurementEstimator_cfi::Chi2MeasurementEstimator, HLT_2018_cff::errorMatrixPset, MuonServiceProxy::eventSetup(), edm::EventSetup::get(), get, edm::ParameterSet::getParameter(), AlCaHLTBitMon_QueryRunRegistry::string, theAdjustAtIp, theCacheId_MT, theCacheId_TG, theConfig, theErrorMatrixAdjuster, theEstimator, theFixedErrorRescaling, theFlexErrorRescaling, theMaxChi2, theNavigation, thePropagatorName, theResetMethod, theSelectStateFlag, theService, theSigmaZ, theTracker, theUpdateStateFlag, theUpdator, and theUseVertexStateFlag.

◆ innerState()

TrajectoryStateOnSurface TSGFromPropagation::innerState ( const TrackCand & staMuon ) const

private

Definition at line 235 of file TSGFromPropagation.cc.

                                                                                       {
   TrajectoryStateOnSurface innerTS;
  
   if (staMuon.first && staMuon.first->isValid()) {
     if (staMuon.first->direction() == alongMomentum) {
       innerTS = staMuon.first->firstMeasurement().updatedState();
     } else if (staMuon.first->direction() == oppositeToMomentum) {
       innerTS = staMuon.first->lastMeasurement().updatedState();
     }
   } else {
     innerTS = trajectoryStateTransform::innerStateOnSurface(
         *(staMuon.second), *theService->trackingGeometry(), &*theService->magneticField());
   }
   //rescale the error
   adjust(innerTS);
  
   return innerTS;
  
   //    return trajectoryStateTransform::innerStateOnSurface(*(staMuon.second),*theService->trackingGeometry(), &*theService->magneticField());
 }

References adjust(), alongMomentum, trajectoryStateTransform::innerStateOnSurface(), MuonServiceProxy::magneticField(), oppositeToMomentum, theService, and MuonServiceProxy::trackingGeometry().

Referenced by outerTkState().

◆ outerTkState()

TrajectoryStateOnSurface TSGFromPropagation::outerTkState ( const TrackCand & staMuon ) const

private

Definition at line 256 of file TSGFromPropagation.cc.

                                                                                         {
   TrajectoryStateOnSurface result;
  
   if (theUseVertexStateFlag && staMuon.second->pt() > 1.0) {
     FreeTrajectoryState iniState =
         trajectoryStateTransform::initialFreeState(*(staMuon.second), &*theService->magneticField());
     //rescale the error at IP
     adjust(iniState);
  
     StateOnTrackerBound fromInside(&*(theService->propagator("PropagatorWithMaterial")));
     result = fromInside(iniState);
   } else {
     StateOnTrackerBound fromOutside(&*propagator());
     result = fromOutside(innerState(staMuon));
   }
   return result;
 }

References adjust(), trajectoryStateTransform::initialFreeState(), innerState(), MuonServiceProxy::magneticField(), propagator(), MuonServiceProxy::propagator(), mps_fire::result, theService, and theUseVertexStateFlag.

Referenced by trackerSeeds().

◆ passSelection()

bool TSGFromPropagation::passSelection ( const TrajectoryStateOnSurface & tsos ) const

private

check some quantity and beam-spot compatibility and decide to continue

Definition at line 319 of file TSGFromPropagation.cc.

                                                                                  {
   if (!theSelectStateFlag)
     return true;
   else {
     if (beamSpot.isValid()) {
       return ((fabs(zDis(tsos) - beamSpot->z0()) < theSigmaZ));
  
     } else {
       return ((fabs(zDis(tsos)) < theSigmaZ));
       //      double theDxyCut = 100;
       //      return ( (zDis(tsos) < theSigmaZ) && (dxyDis(tsos) < theDxyCut) );
     }
   }
 }

References beamSpot, edm::HandleBase::isValid(), theSelectStateFlag, theSigmaZ, reco::BeamSpot::z0(), and zDis().

Referenced by trackerSeeds().

◆ propagator()

edm::ESHandle<Propagator> TSGFromPropagation::propagator ( ) const

inlineprivate

Definition at line 62 of file TSGFromPropagation.h.

62 { return theService->propagator(thePropagatorName); }

References MuonServiceProxy::propagator(), thePropagatorName, and theService.

Referenced by findMeasurements(), outerTkState(), and trackerSeeds().

◆ setEvent()

void TSGFromPropagation::setEvent ( const edm::Event & iEvent )

overridevirtual

set an event

Reimplemented from TrackerSeedGenerator.

Definition at line 202 of file TSGFromPropagation.cc.

                                                         {
   iEvent.getByToken(theBeamSpotToken, beamSpot);
  
   unsigned long long newCacheId_MT = theService->eventSetup().get<CkfComponentsRecord>().cacheIdentifier();
  
   if (theUpdateStateFlag && newCacheId_MT != theCacheId_MT) {
     LogTrace(theCategory) << "Measurment Tracker Geometry changed!";
     theCacheId_MT = newCacheId_MT;
     theService->eventSetup().get<CkfComponentsRecord>().get(theMeasTrackerName, theMeasTracker);
   }
  
   if (theUpdateStateFlag) {
     iEvent.getByToken(theMeasurementTrackerEventToken, theMeasTrackerEvent);
   }
  
   bool trackerGeomChanged = false;
  
   unsigned long long newCacheId_TG = theService->eventSetup().get<TrackerRecoGeometryRecord>().cacheIdentifier();
  
   if (newCacheId_TG != theCacheId_TG) {
     LogTrace(theCategory) << "Tracker Reco Geometry changed!";
     theCacheId_TG = newCacheId_TG;
     theService->eventSetup().get<TrackerRecoGeometryRecord>().get(theTracker);
     trackerGeomChanged = true;
   }
  
   if (trackerGeomChanged && (theTracker.product() != nullptr)) {
     if (theNavigation)
       delete theNavigation;
     theNavigation = new DirectTrackerNavigation(theTracker);
   }
 }

References beamSpot, MuonServiceProxy::eventSetup(), edm::EventSetup::get(), get, iEvent, LogTrace, edm::ESHandle< T >::product(), theBeamSpotToken, theCacheId_MT, theCacheId_TG, theCategory, theMeasTracker, theMeasTrackerEvent, theMeasTrackerName, theMeasurementTrackerEventToken, theNavigation, theService, theTracker, and theUpdateStateFlag.

◆ trackerSeeds()

void TSGFromPropagation::trackerSeeds	(	const TrackCand &	staMuon,
		const TrackingRegion &	region,
		const TrackerTopology *	tTopo,
		std::vector< TrajectorySeed > &	result
	)

override

generate seed(s) for a track

Definition at line 64 of file TSGFromPropagation.cc.

                                                                          {
   if (theResetMethod == "discrete")
     getRescalingFactor(staMuon);
  
   TrajectoryStateOnSurface staState = outerTkState(staMuon);
  
   if (!staState.isValid()) {
     LogTrace(theCategory) << "Error: initial state from L2 muon is invalid.";
     return;
   }
  
   LogTrace(theCategory) << "begin of trackerSeed:\n staState pos: " << staState.globalPosition()
                         << " mom: " << staState.globalMomentum() << "pos eta: " << staState.globalPosition().eta()
                         << "mom eta: " << staState.globalMomentum().eta();
  
   std::vector<const DetLayer*> nls = theNavigation->compatibleLayers(*(staState.freeState()), oppositeToMomentum);
  
   LogTrace(theCategory) << " compatible layers: " << nls.size();
  
   if (nls.empty())
     return;
  
   int ndesLayer = 0;
  
   bool usePredictedState = false;
  
   if (theUpdateStateFlag) {  //use updated states
     std::vector<TrajectoryMeasurement> alltm;
  
     for (std::vector<const DetLayer*>::const_iterator inl = nls.begin(); inl != nls.end(); inl++, ndesLayer++) {
       if ((*inl == nullptr))
         break;
       //         if ( (inl != nls.end()-1 ) && ( (*inl)->subDetector() == GeomDetEnumerators::TEC ) && ( (*(inl+1))->subDetector() == GeomDetEnumerators::TOB ) ) continue;
       alltm = findMeasurements(*inl, staState);
       if ((!alltm.empty())) {
         LogTrace(theCategory) << "final compatible layer: " << ndesLayer;
         break;
       }
     }
  
     if (alltm.empty()) {
       LogTrace(theCategory) << " NO Measurements Found: eta: " << staState.globalPosition().eta() << "pt "
                             << staState.globalMomentum().perp();
       usePredictedState = true;
     } else {
       LogTrace(theCategory) << " Measurements for seeds: " << alltm.size();
       std::stable_sort(alltm.begin(), alltm.end(), increasingEstimate());
       if (alltm.size() > 5)
         alltm.erase(alltm.begin() + 5, alltm.end());
  
       int i = 0;
       for (std::vector<TrajectoryMeasurement>::const_iterator itm = alltm.begin(); itm != alltm.end(); itm++, i++) {
         TrajectoryStateOnSurface updatedTSOS = updator()->update(itm->predictedState(), *(itm->recHit()));
         if (updatedTSOS.isValid() && passSelection(updatedTSOS)) {
           edm::OwnVector<TrackingRecHit> container;
           container.push_back(itm->recHit()->hit()->clone());
           TrajectorySeed ts = createSeed(updatedTSOS, container, itm->recHit()->geographicalId());
           result.push_back(ts);
         }
       }
       LogTrace(theCategory) << "result: " << result.size();
       return;
     }
   }
  
   if (!theUpdateStateFlag || usePredictedState) {  //use predicted states
     LogTrace(theCategory) << "use predicted state: ";
     for (std::vector<const DetLayer*>::const_iterator inl = nls.begin(); inl != nls.end(); inl++) {
       if (!result.empty() || *inl == nullptr) {
         break;
       }
       std::vector<DetLayer::DetWithState> compatDets = (*inl)->compatibleDets(staState, *propagator(), *estimator());
       LogTrace(theCategory) << " compatDets " << compatDets.size();
       if (compatDets.empty())
         continue;
       TrajectorySeed ts = createSeed(compatDets.front().second, compatDets.front().first->geographicalId());
       result.push_back(ts);
     }
     LogTrace(theCategory) << "result: " << result.size();
     return;
   }
   return;
 }

References DirectTrackerNavigation::compatibleLayers(), createSeed(), estimator(), PV3DBase< T, PVType, FrameType >::eta(), findMeasurements(), TrajectoryStateOnSurface::freeState(), getRescalingFactor(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), mps_fire::i, TrajectoryStateOnSurface::isValid(), LogTrace, oppositeToMomentum, outerTkState(), passSelection(), PV3DBase< T, PVType, FrameType >::perp(), propagator(), edm::OwnVector< T, P >::push_back(), mps_fire::result, theCategory, theNavigation, theResetMethod, theUpdateStateFlag, TrajectoryStateUpdator::update(), and updator().

◆ updator()

const TrajectoryStateUpdator* TSGFromPropagation::updator ( ) const

inlineprivate

Definition at line 58 of file TSGFromPropagation.h.

58 { return theUpdator; }

References theUpdator.

Referenced by trackerSeeds().

◆ validMeasurements()

void TSGFromPropagation::validMeasurements ( std::vector< TrajectoryMeasurement > & tms ) const

private

select valid measurements

Definition at line 286 of file TSGFromPropagation.cc.

                                                                                       {
   std::vector<TrajectoryMeasurement>::iterator tmsend = std::remove_if(tms.begin(), tms.end(), isInvalid());
   tms.erase(tmsend, tms.end());
   return;
 }

Referenced by findMeasurements().

◆ zDis()

double TSGFromPropagation::zDis ( const TrajectoryStateOnSurface & tsos ) const

private

Definition at line 340 of file TSGFromPropagation.cc.

                                                                           {
   return tsos.globalPosition().z() -
          tsos.globalPosition().perp() * tsos.globalMomentum().z() / tsos.globalMomentum().perp();
 }

References TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), PV3DBase< T, PVType, FrameType >::perp(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by passSelection().

Member Data Documentation

◆ beamSpot

edm::Handle<reco::BeamSpot> TSGFromPropagation::beamSpot

private

Definition at line 152 of file TSGFromPropagation.h.

Referenced by passSelection(), and setEvent().

◆ theAdjustAtIp

bool TSGFromPropagation::theAdjustAtIp

private

Definition at line 146 of file TSGFromPropagation.h.

Referenced by init().

◆ theBeamSpotInputTag

edm::InputTag TSGFromPropagation::theBeamSpotInputTag

private

Definition at line 153 of file TSGFromPropagation.h.

Referenced by TSGFromPropagation().

◆ theBeamSpotToken

edm::EDGetTokenT<reco::BeamSpot> TSGFromPropagation::theBeamSpotToken

private

Definition at line 154 of file TSGFromPropagation.h.

Referenced by setEvent(), and TSGFromPropagation().

◆ theCacheId_MT

unsigned long long TSGFromPropagation::theCacheId_MT

private

Definition at line 106 of file TSGFromPropagation.h.

Referenced by init(), and setEvent().

◆ theCacheId_TG

unsigned long long TSGFromPropagation::theCacheId_TG

private

Definition at line 107 of file TSGFromPropagation.h.

Referenced by init(), and setEvent().

◆ theCategory

std::string TSGFromPropagation::theCategory

private

Definition at line 109 of file TSGFromPropagation.h.

Referenced by setEvent(), trackerSeeds(), TSGFromPropagation(), and ~TSGFromPropagation().

◆ theConfig

edm::ParameterSet TSGFromPropagation::theConfig

private

Definition at line 150 of file TSGFromPropagation.h.

Referenced by init(), and TSGFromPropagation().

◆ theErrorMatrixAdjuster

MuonErrorMatrix* TSGFromPropagation::theErrorMatrixAdjuster

private

Definition at line 144 of file TSGFromPropagation.h.

Referenced by adjust(), init(), and ~TSGFromPropagation().

◆ theEstimator

const Chi2MeasurementEstimator* TSGFromPropagation::theEstimator

private

Definition at line 124 of file TSGFromPropagation.h.

Referenced by estimator(), init(), and ~TSGFromPropagation().

◆ theFixedErrorRescaling

double TSGFromPropagation::theFixedErrorRescaling

private

Definition at line 132 of file TSGFromPropagation.h.

Referenced by adjust(), and init().

◆ theFlexErrorRescaling

double TSGFromPropagation::theFlexErrorRescaling

private

Definition at line 130 of file TSGFromPropagation.h.

Referenced by adjust(), getRescalingFactor(), and init().

◆ theMaxChi2

double TSGFromPropagation::theMaxChi2

private

Definition at line 128 of file TSGFromPropagation.h.

Referenced by init().

◆ theMeasTracker

edm::ESHandle<MeasurementTracker> TSGFromPropagation::theMeasTracker

private

Definition at line 114 of file TSGFromPropagation.h.

Referenced by setEvent().

◆ theMeasTrackerEvent

edm::Handle<MeasurementTrackerEvent> TSGFromPropagation::theMeasTrackerEvent

private

Definition at line 116 of file TSGFromPropagation.h.

Referenced by findMeasurements(), and setEvent().

◆ theMeasTrackerName

std::string TSGFromPropagation::theMeasTrackerName

private

Definition at line 113 of file TSGFromPropagation.h.

Referenced by setEvent(), and TSGFromPropagation().

◆ theMeasurementTrackerEventTag

edm::InputTag TSGFromPropagation::theMeasurementTrackerEventTag

private

Definition at line 115 of file TSGFromPropagation.h.

Referenced by TSGFromPropagation().

◆ theMeasurementTrackerEventToken

edm::EDGetTokenT<MeasurementTrackerEvent> TSGFromPropagation::theMeasurementTrackerEventToken

private

Definition at line 155 of file TSGFromPropagation.h.

Referenced by setEvent(), and TSGFromPropagation().

◆ theNavigation

const DirectTrackerNavigation* TSGFromPropagation::theNavigation

private

Definition at line 118 of file TSGFromPropagation.h.

Referenced by init(), setEvent(), trackerSeeds(), and ~TSGFromPropagation().

◆ thePropagatorName

std::string TSGFromPropagation::thePropagatorName

private

Definition at line 142 of file TSGFromPropagation.h.

Referenced by init(), and propagator().

◆ theResetMethod

std::string TSGFromPropagation::theResetMethod

private

Definition at line 138 of file TSGFromPropagation.h.

Referenced by adjust(), init(), and trackerSeeds().

◆ theSelectStateFlag

bool TSGFromPropagation::theSelectStateFlag

private

Definition at line 140 of file TSGFromPropagation.h.

Referenced by init(), and passSelection().

◆ theService

const MuonServiceProxy* TSGFromPropagation::theService

private

Definition at line 120 of file TSGFromPropagation.h.

Referenced by init(), innerState(), outerTkState(), propagator(), and setEvent().

◆ theSigmaZ

double TSGFromPropagation::theSigmaZ

private

Definition at line 148 of file TSGFromPropagation.h.

Referenced by init(), and passSelection().

◆ theTracker

edm::ESHandle<GeometricSearchTracker> TSGFromPropagation::theTracker

private

Definition at line 111 of file TSGFromPropagation.h.

Referenced by init(), and setEvent().

◆ theTSTransformer

TrajectoryStateTransform* TSGFromPropagation::theTSTransformer

private

Definition at line 126 of file TSGFromPropagation.h.

◆ theUpdateStateFlag

bool TSGFromPropagation::theUpdateStateFlag

private

Definition at line 136 of file TSGFromPropagation.h.

Referenced by init(), setEvent(), and trackerSeeds().

◆ theUpdator

const TrajectoryStateUpdator* TSGFromPropagation::theUpdator

private

Definition at line 122 of file TSGFromPropagation.h.

Referenced by init(), updator(), and ~TSGFromPropagation().

◆ theUseVertexStateFlag

bool TSGFromPropagation::theUseVertexStateFlag

private

Definition at line 134 of file TSGFromPropagation.h.

Referenced by init(), and outerTkState().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TSGFromPropagation() [1/2]

◆ TSGFromPropagation() [2/2]

◆ ~TSGFromPropagation()

Member Function Documentation

◆ adjust() [1/2]

◆ adjust() [2/2]

◆ createSeed() [1/2]

◆ createSeed() [2/2]

◆ dxyDis()

◆ estimator()

◆ findMeasurements()

◆ getRescalingFactor()

◆ init()

◆ innerState()

◆ outerTkState()

◆ passSelection()

◆ propagator()

◆ setEvent()

◆ trackerSeeds()

◆ updator()

◆ validMeasurements()

◆ zDis()

Member Data Documentation

◆ beamSpot

◆ theAdjustAtIp

◆ theBeamSpotInputTag

◆ theBeamSpotToken

◆ theCacheId_MT

◆ theCacheId_TG

◆ theCategory

◆ theConfig

◆ theErrorMatrixAdjuster

◆ theEstimator

◆ theFixedErrorRescaling

◆ theFlexErrorRescaling

◆ theMaxChi2

◆ theMeasTracker

◆ theMeasTrackerEvent

◆ theMeasTrackerName

◆ theMeasurementTrackerEventTag

◆ theMeasurementTrackerEventToken

◆ theNavigation

◆ thePropagatorName

◆ theResetMethod

◆ theSelectStateFlag

◆ theService

◆ theSigmaZ

◆ theTracker

◆ theTSTransformer

◆ theUpdateStateFlag

◆ theUpdator

◆ theUseVertexStateFlag