#include <FastTSGFromPropagation.h>

Inheritance diagram for FastTSGFromPropagation:

Classes
struct	increasingEstimate

struct	isInvalid

Public Member Functions
	FastTSGFromPropagation (const edm::ParameterSet &par, const MuonServiceProxy *, edm::ConsumesCollector &iC)

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

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 *tTopo, std::vector< TrajectorySeed > &) override
	generate seed(s) for a track More...

	~FastTSGFromPropagation () 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

void	stateOnDet (const TrajectoryStateOnSurface &ts, unsigned int detid, PTrajectoryStateOnDet &pts) const
	A mere copy (without memory leak) of an existing tracking method. More...

const TrajectoryStateUpdator *	updator () const

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

double	zDis (const TrajectoryStateOnSurface &tsos) const

Private Attributes
edm::EDGetTokenT< reco::BeamSpot >	beamSpot_

edm::Handle< FastTrackerRecHitCombinationCollection >	recHitCombinations

edm::EDGetTokenT< FastTrackerRecHitCombinationCollection >	recHitCombinationsToken_

bool	theAdjustAtIp

edm::Handle< reco::BeamSpot >	theBeamSpot

unsigned long long	theCacheId_MT

unsigned long long	theCacheId_TG

std::string	theCategory

const edm::ParameterSet	theConfig

std::unique_ptr< MuonErrorMatrix >	theErrorMatrixAdjuster

std::unique_ptr< const Chi2MeasurementEstimator >	theEstimator

double	theFixedErrorRescaling

double	theFlexErrorRescaling

const TrackerGeometry *	theGeometry

double	theMaxChi2

edm::ESHandle< MeasurementTracker >	theMeasTracker

edm::Handle< MeasurementTrackerEvent >	theMeasTrackerEvent

edm::EDGetTokenT< MeasurementTrackerEvent >	theMeasurementTrackerEventToken_

std::unique_ptr< const DirectTrackerNavigation >	theNavigation

std::string	thePropagatorName

std::string	theResetMethod

bool	theSelectStateFlag

const MuonServiceProxy *	theService

double	theSigmaZ

edm::EDGetTokenT< edm::SimTrackContainer >	theSimTrackCollectionToken_

edm::Handle< edm::SimTrackContainer >	theSimTracks

edm::ESHandle< GeometricSearchTracker >	theTracker

edm::ESHandle< TransientTrackingRecHitBuilder >	theTTRHBuilder

bool	theUpdateStateFlag

std::unique_ptr< 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

Emulate TSGFromPropagation in RecoMuon

Author: Hwidong Yoo - Purdue University

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

Emulate TSGFromPropagation in RecoMuon

Author: Hwidong Yoo - Purdue University

Definition at line 48 of file FastTSGFromPropagation.h.

Constructor & Destructor Documentation

◆ FastTSGFromPropagation() [1/2]

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

constructor

Definition at line 46 of file FastTSGFromPropagation.cc.

47 : FastTSGFromPropagation(iConfig, nullptr, iC) {}

◆ FastTSGFromPropagation() [2/2]

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

Definition at line 49 of file FastTSGFromPropagation.cc.

     : theCategory("FastSimulation|Muons|FastTSGFromPropagation"),
       theTracker(),
       theNavigation(),
       theService(service),
       theUpdator(),
       theEstimator(),
       theSigmaZ(0.0),
       theConfig(iConfig),
       theSimTrackCollectionToken_(
           iC.consumes<edm::SimTrackContainer>(theConfig.getParameter<edm::InputTag>("SimTrackCollectionLabel"))),
       recHitCombinationsToken_(
           iC.consumes<FastTrackerRecHitCombinationCollection>(theConfig.getParameter<edm::InputTag>("HitProducer"))),
       beamSpot_(iC.consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamSpot"))),
       theMeasurementTrackerEventToken_(
           iC.consumes<MeasurementTrackerEvent>(iConfig.getParameter<edm::InputTag>("MeasurementTrackerEvent"))) {}

◆ ~FastTSGFromPropagation()

FastTSGFromPropagation::~FastTSGFromPropagation ( )

override

destructor

Definition at line 68 of file FastTSGFromPropagation.cc.

68 { LogTrace(theCategory) << " FastTSGFromPropagation dtor called "; }

References LogTrace, and theCategory.

Member Function Documentation

◆ adjust() [1/2]

void FastTSGFromPropagation::adjust ( FreeTrajectoryState & state ) const

private

adjust the error matrix of the FTS

Definition at line 512 of file FastTSGFromPropagation.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(), FreeTrajectoryState::momentum(), MuonErrorMatrix::multiply(), FreeTrajectoryState::parameters(), FreeTrajectoryState::rescaleError(), theErrorMatrixAdjuster, theFixedErrorRescaling, theFlexErrorRescaling, and theResetMethod.

Referenced by innerState(), and outerTkState().

◆ adjust() [2/2]

void FastTSGFromPropagation::adjust ( TrajectoryStateOnSurface & state ) const

private

adjust the error matrix of the TSOS

Definition at line 532 of file FastTSGFromPropagation.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(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalParameters(), MuonErrorMatrix::multiply(), TrajectoryStateOnSurface::rescaleError(), TrajectoryStateOnSurface::surface(), TrajectoryStateOnSurface::surfaceSide(), theErrorMatrixAdjuster, theFixedErrorRescaling, theFlexErrorRescaling, theResetMethod, and TrajectoryStateOnSurface::weight().

◆ createSeed() [1/2]

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

private

create a hitless seed from a trajectory state

Definition at line 430 of file FastTSGFromPropagation.cc.

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

Referenced by trackerSeeds().

◆ createSeed() [2/2]

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

private

create a seed from a trajectory state

Definition at line 435 of file FastTSGFromPropagation.cc.

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

References oppositeToMomentum, and trajectoryStateTransform::persistentState().

◆ dxyDis()

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

private

Definition at line 490 of file FastTSGFromPropagation.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* FastTSGFromPropagation::estimator ( void ) const

inlineprivate

Definition at line 80 of file FastTSGFromPropagation.h.

80 { return theEstimator.get(); }

References theEstimator.

Referenced by findMeasurements(), and trackerSeeds().

◆ findMeasurements()

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

private

look for measurements on the first compatible layer

Definition at line 448 of file FastTSGFromPropagation.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(), propagator(), mps_fire::result, and theMeasTrackerEvent.

Referenced by trackerSeeds().

◆ getRescalingFactor()

void FastTSGFromPropagation::getRescalingFactor ( const TrackCand & staMuon )

private

Definition at line 501 of file FastTSGFromPropagation.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 FastTSGFromPropagation::init ( const MuonServiceProxy * service )

overridevirtual

initialize

Reimplemented from TrackerSeedGenerator.

Definition at line 303 of file FastTSGFromPropagation.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("FastTSGFromPropagation") << "Wrong error rescaling method: " << theResetMethod << "\n"
                                             << "Possible choices are: discrete, fixed, matrix.\n"
                                             << "Use discrete method" << std::endl;
     theResetMethod = "discrete";
   }
  
   theEstimator.reset(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.reset(new KFUpdator());
  
   theSigmaZ = theConfig.getParameter<double>("SigmaZ");
  
   edm::ParameterSet errorMatrixPset = theConfig.getParameter<edm::ParameterSet>("errorMatrixPset");
   if (theResetMethod == "matrix" && !errorMatrixPset.empty()) {
     theAdjustAtIp = errorMatrixPset.getParameter<bool>("atIP");
     theErrorMatrixAdjuster.reset(new MuonErrorMatrix(errorMatrixPset));
   } else {
     theAdjustAtIp = false;
     theErrorMatrixAdjuster.reset();
   }
  
   theService->eventSetup().get<TrackerRecoGeometryRecord>().get(theTracker);
   theNavigation.reset(new DirectTrackerNavigation(theTracker));
  
   edm::ESHandle<TrackerGeometry> geometry;
   theService->eventSetup().get<TrackerDigiGeometryRecord>().get(geometry);
   theGeometry = &(*geometry);
  
   theService->eventSetup().get<TransientRecHitRecord>().get("WithTrackAngle", theTTRHBuilder);
 }

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

◆ innerState()

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

private

Definition at line 393 of file FastTSGFromPropagation.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;
 }

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

Referenced by outerTkState().

◆ outerTkState()

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

private

Definition at line 412 of file FastTSGFromPropagation.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(), propagator(), mps_fire::result, theService, and theUseVertexStateFlag.

Referenced by trackerSeeds().

◆ passSelection()

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

private

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

Definition at line 475 of file FastTSGFromPropagation.cc.

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

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

Referenced by trackerSeeds().

◆ propagator()

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

inlineprivate

Definition at line 82 of file FastTSGFromPropagation.h.

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

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

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

◆ setEvent()

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

overridevirtual

set an event

Reimplemented from TrackerSeedGenerator.

Definition at line 358 of file FastTSGFromPropagation.cc.

                                                             {
   iEvent.getByToken(beamSpot_, theBeamSpot);
  
   // retrieve the MC truth (SimTracks)
   iEvent.getByToken(theSimTrackCollectionToken_, theSimTracks);
   iEvent.getByToken(recHitCombinationsToken_, recHitCombinations);
  
   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(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()) {
     theNavigation.reset(new DirectTrackerNavigation(theTracker));
   }
 }

References beamSpot_, get, iEvent, LogTrace, edm::ESHandle< T >::product(), recHitCombinations, recHitCombinationsToken_, theBeamSpot, theCacheId_MT, theCacheId_TG, theCategory, theMeasTracker, theMeasTrackerEvent, theMeasurementTrackerEventToken_, theNavigation, theService, theSimTrackCollectionToken_, theSimTracks, theTracker, and theUpdateStateFlag.

◆ stateOnDet()

void FastTSGFromPropagation::stateOnDet	(	const TrajectoryStateOnSurface &	ts,
		unsigned int	detid,
		PTrajectoryStateOnDet &	pts
	)		const

private

A mere copy (without memory leak) of an existing tracking method.

should check if corresponds to m

Definition at line 552 of file FastTSGFromPropagation.cc.

                                                                           {
   const AlgebraicSymMatrix55& m = ts.localError().matrix();
   int dim = 5;  
   float localErrors[15];
   int k = 0;
   for (int i = 0; i < dim; ++i) {
     for (int j = 0; j <= i; ++j) {
       localErrors[k++] = m(i, j);
     }
   }
   int surfaceSide = static_cast<int>(ts.surfaceSide());
   pts = PTrajectoryStateOnDet(ts.localParameters(), ts.globalMomentum().perp(), localErrors, detid, surfaceSide);
 }

References TrajectoryStateOnSurface::globalMomentum(), mps_fire::i, dqmiolumiharvest::j, dqmdumpme::k, TrajectoryStateOnSurface::localError(), TrajectoryStateOnSurface::localParameters(), visualization-live-secondInstance_cfg::m, LocalTrajectoryError::matrix(), PV3DBase< T, PVType, FrameType >::perp(), RPCpg::pts, and TrajectoryStateOnSurface::surfaceSide().

◆ trackerSeeds()

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

override

generate seed(s) for a track

Definition at line 70 of file FastTSGFromPropagation.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());
  
       const edm::SimTrackContainer* simTracks = &(*theSimTracks);
       TrajectorySeedHitCandidate theSeedHits;
       std::vector<TrajectorySeedHitCandidate> outerHits;
  
       //std::vector<TrajectorySeed>  tmpTS;
       bool isMatch = false;
       for (std::vector<TrajectoryMeasurement>::const_iterator itm = alltm.begin(); itm != alltm.end(); itm++) {
         const TrajectoryStateOnSurface seedState = itm->predictedState();
         double preY = seedState.globalPosition().y();
  
         // Check SimTrack
         FreeTrajectoryState simtrack_trackerstate;
         for (unsigned icomb = 0; icomb < recHitCombinations->size(); ++icomb) {
           const auto& recHitCombination = (*recHitCombinations)[icomb];
           if (recHitCombination.empty())
             continue;
           int32_t simTrackId = recHitCombination.back()->simTrackId(0);
           const SimTrack& simtrack = (*simTracks)[simTrackId];
  
           GlobalPoint position(simtrack.trackerSurfacePosition().x(),
                                simtrack.trackerSurfacePosition().y(),
                                simtrack.trackerSurfacePosition().z());
           GlobalVector momentum(simtrack.trackerSurfaceMomentum().x(),
                                 simtrack.trackerSurfaceMomentum().y(),
                                 simtrack.trackerSurfaceMomentum().z());
           int charge = (int)simtrack.charge();
           GlobalTrajectoryParameters glb_parameters(
               position, momentum, charge, &*theService->magneticField().product());
           simtrack_trackerstate = FreeTrajectoryState(glb_parameters);
  
           unsigned int outerId = 0;
           for (const auto& recHitRef : recHitCombination) {
             theSeedHits = TrajectorySeedHitCandidate(recHitRef.get(), tTopo);
             unsigned int id = theSeedHits.hit()->geographicalId().rawId();
             if (preY < 0) {
               if (id > outerId)
                 outerId = id;
             } else {
               if (id > outerId)
                 outerId = id;
             }
           }
           for (const auto& recHitRef : recHitCombination) {
             theSeedHits = TrajectorySeedHitCandidate(recHitRef.get(), tTopo);
             if (itm->recHit()->hit()->geographicalId().rawId() == theSeedHits.hit()->geographicalId().rawId()) {
               auto aTrackingRecHit = std::unique_ptr<TrackingRecHit>(theSeedHits.hit()->clone());
               TransientTrackingRecHit::ConstRecHitPointer recHit = theTTRHBuilder->build(aTrackingRecHit.get());
               if (!recHit)
                 continue;
               TrajectoryStateOnSurface updatedTSOS = updator()->update(seedState, *(recHit));
               if (updatedTSOS.isValid() && passSelection(updatedTSOS)) {
                 edm::OwnVector<TrackingRecHit> container;
                 container.push_back(recHit->hit()->clone());
                 fastTrackingUtilities::setRecHitCombinationIndex(container, icomb);
                 TrajectorySeed ts = createSeed(updatedTSOS, container, recHit->geographicalId());
                 // check direction
                 const TrajectorySeed* aSeed = &ts;
                 PTrajectoryStateOnDet PTSOD = aSeed->startingState();
  
                 const GeomDet* g = theGeometry->idToDet(PTSOD.detId());
                 TrajectoryStateOnSurface tsos = trajectoryStateTransform::transientState(
                     PTSOD, &(g->surface()), &*theService->magneticField().product());
                 if (tsos.globalMomentum().basicVector() * seedState.globalMomentum().basicVector() < 0.)
                   continue;
                 result.push_back(ts);
                 isMatch = true;
               }
             }
           }
         }
       }
       if (!isMatch) {
         // if there is no hits w.r.t. TM, find outermost hit
         for (std::vector<TrajectoryMeasurement>::const_iterator itm = alltm.begin(); itm != alltm.end(); itm++) {
           const TrajectoryStateOnSurface seedState = itm->predictedState();
           double preY = seedState.globalPosition().y();
  
           // Check SimTrack
           TrackingRecHit* aTrackingRecHit;
           FreeTrajectoryState simtrack_trackerstate;
  
           for (unsigned icomb = 0; icomb < recHitCombinations->size(); ++icomb) {
             const auto& recHitCombination = (*recHitCombinations)[icomb];
             if (recHitCombination.empty())
               continue;
             int32_t simTrackId = recHitCombination.back()->simTrackId(0);
             const SimTrack& simtrack = (*simTracks)[simTrackId];
  
             GlobalPoint position(simtrack.trackerSurfacePosition().x(),
                                  simtrack.trackerSurfacePosition().y(),
                                  simtrack.trackerSurfacePosition().z());
             GlobalVector momentum(simtrack.trackerSurfaceMomentum().x(),
                                   simtrack.trackerSurfaceMomentum().y(),
                                   simtrack.trackerSurfaceMomentum().z());
             int charge = (int)simtrack.charge();
             GlobalTrajectoryParameters glb_parameters(
                 position, momentum, charge, &*theService->magneticField().product());
             simtrack_trackerstate = FreeTrajectoryState(glb_parameters);
  
             unsigned int outerId = 0;
             for (const auto& recHitRef : recHitCombination) {
               theSeedHits = TrajectorySeedHitCandidate(recHitRef.get(), tTopo);
               unsigned int id = theSeedHits.hit()->geographicalId().rawId();
               if (preY < 0) {
                 if (id > outerId)
                   outerId = id;
               } else {
                 if (id > outerId)
                   outerId = id;
               }
             }
             for (const auto& recHitRef : recHitCombination) {
               theSeedHits = TrajectorySeedHitCandidate(recHitRef.get(), tTopo);
               if (outerId == theSeedHits.hit()->geographicalId().rawId()) {
                 aTrackingRecHit = theSeedHits.hit()->clone();
                 TransientTrackingRecHit::ConstRecHitPointer recHit = theTTRHBuilder->build(aTrackingRecHit);
                 if (!recHit)
                   continue;
                 TrajectoryStateOnSurface updatedTSOS = updator()->update(seedState, *(recHit));
                 if (updatedTSOS.isValid() && passSelection(updatedTSOS)) {
                   edm::OwnVector<TrackingRecHit> container;
                   container.push_back(recHit->hit()->clone());
                   fastTrackingUtilities::setRecHitCombinationIndex(container, icomb);
                   TrajectorySeed ts = createSeed(updatedTSOS, container, recHit->geographicalId());
                   // check direction
                   const TrajectorySeed* aSeed = &ts;
                   PTrajectoryStateOnDet PTSOD = aSeed->startingState();
  
                   const GeomDet* g = theGeometry->idToDet(PTSOD.detId());
                   TrajectoryStateOnSurface tsos = trajectoryStateTransform::transientState(
                       PTSOD, &(g->surface()), &*theService->magneticField().product());
                   if (tsos.globalMomentum().basicVector() * seedState.globalMomentum().basicVector() < 0.)
                     continue;
                   result.push_back(ts);
                 }
               }
             }
           }
         }
       }
  
       /*
        for( unsigned ir = 0; ir < tmpTS.size(); ir++ ) {
      const BasicTrajectorySeed* aSeed = &((tmpTS)[ir]);
      PTrajectoryStateOnDet PTSOD = aSeed->startingState();
      
      DetId seedDetId(PTSOD.detId());
      const GeomDet * g = theGeometry->idToDet(seedDetId);
      TrajectoryStateOnSurface tsos = trajectoryStateTransform::transientState(PTSOD, &(g->surface()),  &*theService->magneticField().product());
          cout << "tsos3 = " << tsos.globalMomentum() << endl;
      if( _index == ir ) {
          cout << "tsos4 = " << tsos.globalMomentum() << endl;
          result.push_back(tmpTS[ir]);
      }
        }
        */
       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 PV3DBase< T, PVType, FrameType >::basicVector(), TransientTrackingRecHitBuilder::build(), CoreSimTrack::charge(), ALCARECOTkAlJpsiMuMu_cff::charge, FastTrackerRecHit::clone(), createSeed(), PTrajectoryStateOnDet::detId(), estimator(), PV3DBase< T, PVType, FrameType >::eta(), findMeasurements(), TrajectoryStateOnSurface::freeState(), g, TrackingRecHit::geographicalId(), getRescalingFactor(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), TrajectorySeedHitCandidate::hit(), triggerObjects_cff::id, TrackerGeometry::idToDet(), createfilelist::int, cms::dd::isMatch(), TrajectoryStateOnSurface::isValid(), LogTrace, oppositeToMomentum, outerTkState(), passSelection(), PV3DBase< T, PVType, FrameType >::perp(), position, propagator(), edm::OwnVector< T, P >::push_back(), DetId::rawId(), rpcPointValidation_cfi::recHit, recHitCombinations, mps_fire::result, fastTrackingUtilities::setRecHitCombinationIndex(), TrackCandidateProducer_cfi::simTracks, theCategory, theGeometry, theNavigation, theResetMethod, theService, theTTRHBuilder, theUpdateStateFlag, SimTrack::trackerSurfaceMomentum(), SimTrack::trackerSurfacePosition(), trajectoryStateTransform::transientState(), TrajectoryStateUpdator::update(), updator(), and PV3DBase< T, PVType, FrameType >::y().

◆ updator()

const TrajectoryStateUpdator* FastTSGFromPropagation::updator ( ) const

inlineprivate

Definition at line 78 of file FastTSGFromPropagation.h.

78 { return theUpdator.get(); }

References theUpdator.

Referenced by trackerSeeds().

◆ validMeasurements()

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

private

select valid measurements

Definition at line 442 of file FastTSGFromPropagation.cc.

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

◆ zDis()

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

private

Definition at line 496 of file FastTSGFromPropagation.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::EDGetTokenT<reco::BeamSpot> FastTSGFromPropagation::beamSpot_

private

Definition at line 170 of file FastTSGFromPropagation.h.

Referenced by setEvent().

◆ recHitCombinations

edm::Handle<FastTrackerRecHitCombinationCollection> FastTSGFromPropagation::recHitCombinations

private

Definition at line 175 of file FastTSGFromPropagation.h.

Referenced by setEvent(), and trackerSeeds().

◆ recHitCombinationsToken_

edm::EDGetTokenT<FastTrackerRecHitCombinationCollection> FastTSGFromPropagation::recHitCombinationsToken_

private

Definition at line 169 of file FastTSGFromPropagation.h.

Referenced by setEvent().

◆ theAdjustAtIp

bool FastTSGFromPropagation::theAdjustAtIp

private

Definition at line 163 of file FastTSGFromPropagation.h.

Referenced by init().

◆ theBeamSpot

edm::Handle<reco::BeamSpot> FastTSGFromPropagation::theBeamSpot

private

Definition at line 173 of file FastTSGFromPropagation.h.

Referenced by passSelection(), and setEvent().

◆ theCacheId_MT

unsigned long long FastTSGFromPropagation::theCacheId_MT

private

Definition at line 126 of file FastTSGFromPropagation.h.

Referenced by init(), and setEvent().

◆ theCacheId_TG

unsigned long long FastTSGFromPropagation::theCacheId_TG

private

Definition at line 127 of file FastTSGFromPropagation.h.

Referenced by init(), and setEvent().

◆ theCategory

std::string FastTSGFromPropagation::theCategory

private

Definition at line 129 of file FastTSGFromPropagation.h.

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

◆ theConfig

const edm::ParameterSet FastTSGFromPropagation::theConfig

private

Definition at line 167 of file FastTSGFromPropagation.h.

Referenced by init().

◆ theErrorMatrixAdjuster

std::unique_ptr<MuonErrorMatrix> FastTSGFromPropagation::theErrorMatrixAdjuster

private

Definition at line 161 of file FastTSGFromPropagation.h.

Referenced by adjust(), and init().

◆ theEstimator

std::unique_ptr<const Chi2MeasurementEstimator> FastTSGFromPropagation::theEstimator

private

Definition at line 143 of file FastTSGFromPropagation.h.

Referenced by estimator(), and init().

◆ theFixedErrorRescaling

double FastTSGFromPropagation::theFixedErrorRescaling

private

Definition at line 149 of file FastTSGFromPropagation.h.

Referenced by adjust(), and init().

◆ theFlexErrorRescaling

double FastTSGFromPropagation::theFlexErrorRescaling

private

Definition at line 147 of file FastTSGFromPropagation.h.

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

◆ theGeometry

const TrackerGeometry* FastTSGFromPropagation::theGeometry

private

Definition at line 137 of file FastTSGFromPropagation.h.

Referenced by init(), and trackerSeeds().

◆ theMaxChi2

double FastTSGFromPropagation::theMaxChi2

private

Definition at line 145 of file FastTSGFromPropagation.h.

Referenced by init().

◆ theMeasTracker

edm::ESHandle<MeasurementTracker> FastTSGFromPropagation::theMeasTracker

private

Definition at line 133 of file FastTSGFromPropagation.h.

Referenced by setEvent().

◆ theMeasTrackerEvent

edm::Handle<MeasurementTrackerEvent> FastTSGFromPropagation::theMeasTrackerEvent

private

Definition at line 176 of file FastTSGFromPropagation.h.

Referenced by findMeasurements(), and setEvent().

◆ theMeasurementTrackerEventToken_

edm::EDGetTokenT<MeasurementTrackerEvent> FastTSGFromPropagation::theMeasurementTrackerEventToken_

private

Definition at line 171 of file FastTSGFromPropagation.h.

Referenced by setEvent().

◆ theNavigation

std::unique_ptr<const DirectTrackerNavigation> FastTSGFromPropagation::theNavigation

private

Definition at line 135 of file FastTSGFromPropagation.h.

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

◆ thePropagatorName

std::string FastTSGFromPropagation::thePropagatorName

private

Definition at line 159 of file FastTSGFromPropagation.h.

Referenced by init(), and propagator().

◆ theResetMethod

std::string FastTSGFromPropagation::theResetMethod

private

Definition at line 155 of file FastTSGFromPropagation.h.

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

◆ theSelectStateFlag

bool FastTSGFromPropagation::theSelectStateFlag

private

Definition at line 157 of file FastTSGFromPropagation.h.

Referenced by init(), and passSelection().

◆ theService

const MuonServiceProxy* FastTSGFromPropagation::theService

private

Definition at line 139 of file FastTSGFromPropagation.h.

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

◆ theSigmaZ

double FastTSGFromPropagation::theSigmaZ

private

Definition at line 165 of file FastTSGFromPropagation.h.

Referenced by init(), and passSelection().

◆ theSimTrackCollectionToken_

edm::EDGetTokenT<edm::SimTrackContainer> FastTSGFromPropagation::theSimTrackCollectionToken_

private

Definition at line 168 of file FastTSGFromPropagation.h.

Referenced by setEvent().

◆ theSimTracks

edm::Handle<edm::SimTrackContainer> FastTSGFromPropagation::theSimTracks

private

Definition at line 174 of file FastTSGFromPropagation.h.

Referenced by setEvent().

◆ theTracker

edm::ESHandle<GeometricSearchTracker> FastTSGFromPropagation::theTracker

private

Definition at line 131 of file FastTSGFromPropagation.h.

Referenced by init(), and setEvent().

◆ theTTRHBuilder

edm::ESHandle<TransientTrackingRecHitBuilder> FastTSGFromPropagation::theTTRHBuilder

private

Definition at line 177 of file FastTSGFromPropagation.h.

Referenced by init(), and trackerSeeds().

◆ theUpdateStateFlag

bool FastTSGFromPropagation::theUpdateStateFlag

private

Definition at line 153 of file FastTSGFromPropagation.h.

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

◆ theUpdator

std::unique_ptr<const TrajectoryStateUpdator> FastTSGFromPropagation::theUpdator

private

Definition at line 141 of file FastTSGFromPropagation.h.

Referenced by init(), and updator().

◆ theUseVertexStateFlag

bool FastTSGFromPropagation::theUseVertexStateFlag

private

Definition at line 151 of file FastTSGFromPropagation.h.

Referenced by init(), and outerTkState().

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ FastTSGFromPropagation() [1/2]

◆ FastTSGFromPropagation() [2/2]

◆ ~FastTSGFromPropagation()

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()

◆ stateOnDet()

◆ trackerSeeds()

◆ updator()

◆ validMeasurements()

◆ zDis()

Member Data Documentation

◆ beamSpot_

◆ recHitCombinations

◆ recHitCombinationsToken_

◆ theAdjustAtIp

◆ theBeamSpot

◆ theCacheId_MT

◆ theCacheId_TG

◆ theCategory

◆ theConfig

◆ theErrorMatrixAdjuster

◆ theEstimator

◆ theFixedErrorRescaling

◆ theFlexErrorRescaling

◆ theGeometry

◆ theMaxChi2

◆ theMeasTracker

◆ theMeasTrackerEvent

◆ theMeasurementTrackerEventToken_

◆ theNavigation

◆ thePropagatorName

◆ theResetMethod

◆ theSelectStateFlag

◆ theService

◆ theSigmaZ

◆ theSimTrackCollectionToken_

◆ theSimTracks

◆ theTracker

◆ theTTRHBuilder

◆ theUpdateStateFlag

◆ theUpdator

◆ theUseVertexStateFlag