#include <MuonCkfTrajectoryBuilder.h>

Inheritance diagram for MuonCkfTrajectoryBuilder:

Public Member Functions
	MuonCkfTrajectoryBuilder (const edm::ParameterSet &conf, edm::ConsumesCollector &iC)

	~MuonCkfTrajectoryBuilder () override

Public Member Functions inherited from CkfTrajectoryBuilder
TempTrajectory	buildTrajectories (const TrajectorySeed &, TrajectoryContainer &ret, unsigned int &nCandPerSeed, const TrajectoryFilter *) const override

	CkfTrajectoryBuilder (const edm::ParameterSet &conf, edm::ConsumesCollector &iC)

	CkfTrajectoryBuilder (const edm::ParameterSet &conf, std::unique_ptr< TrajectoryFilter > filter)

void	rebuildTrajectories (TempTrajectory const &startingTraj, const TrajectorySeed &, TrajectoryContainer &result) const override

TrajectoryContainer	trajectories (const TrajectorySeed &seed) const override
	trajectories building starting from a seed More...

void	trajectories (const TrajectorySeed &seed, TrajectoryContainer &ret) const override
	trajectories building starting from a seed More...

	~CkfTrajectoryBuilder () override

Public Member Functions inherited from BaseCkfTrajectoryBuilder
	BaseCkfTrajectoryBuilder (const BaseCkfTrajectoryBuilder &)=delete

	BaseCkfTrajectoryBuilder (const edm::ParameterSet &conf, std::unique_ptr< TrajectoryFilter > filter, std::unique_ptr< TrajectoryFilter > inOutFilter=nullptr)

const TransientTrackingRecHitBuilder *	hitBuilder () const

BaseCkfTrajectoryBuilder &	operator= (const BaseCkfTrajectoryBuilder &)=delete

virtual void	setDebugger (CkfDebugger *dbg) const

void	setEvent (const edm::Event &event) const override

void	setEvent (const edm::Event &iEvent, const edm::EventSetup &iSetup, const MeasurementTrackerEvent *data)

void	setNavigationSchool (NavigationSchool const *nv)

void	unset () const override

	~BaseCkfTrajectoryBuilder () override

Public Member Functions inherited from TrajectoryBuilder
virtual void	rebuildSeedingRegion (const TrajectorySeed &, TrajectoryContainer &result) const

virtual	~TrajectoryBuilder ()

Protected Member Functions
void	collectMeasurement (const DetLayer layer, const std::vector< const DetLayer > &nl, const TrajectoryStateOnSurface &currentState, std::vector< TM > &result, int &invalidHits, const Propagator *) const

void	findCompatibleMeasurements (const TrajectorySeed &seed, const TempTrajectory &traj, std::vector< TrajectoryMeasurement > &result) const override

void	setEvent_ (const edm::Event &iEvent, const edm::EventSetup &iSetup) override

Protected Member Functions inherited from CkfTrajectoryBuilder
unsigned int	limitedCandidates (const std::shared_ptr< const TrajectorySeed > &sharedSeed, TempTrajectoryContainer &candidates, TrajectoryContainer &result) const

unsigned int	limitedCandidates (const TrajectorySeed &seed, TempTrajectory &startingTraj, TrajectoryContainer &result) const

void	setEvent_ (const edm::Event &iEvent, const edm::EventSetup &iSetup) override

void	updateTrajectory (TempTrajectory &traj, TM &&tm) const

Protected Member Functions inherited from BaseCkfTrajectoryBuilder
void	addToResult (std::shared_ptr< const TrajectorySeed > const &seed, TempTrajectory &traj, TrajectoryContainer &result, bool inOut=false) const

void	addToResult (TempTrajectory const &traj, TempTrajectoryContainer &result, bool inOut=false) const

virtual bool	analyzeMeasurementsDebugger (TempTrajectory &traj, const std::vector< TrajectoryMeasurement > &meas, const MeasurementTrackerEvent theMeasurementTracker, const Propagator theForwardPropagator, const Chi2MeasurementEstimatorBase theEstimator, const TransientTrackingRecHitBuilder theTTRHBuilder) const

virtual bool	analyzeMeasurementsDebugger (Trajectory &traj, const std::vector< TrajectoryMeasurement > &meas, const MeasurementTrackerEvent theMeasurementTracker, const Propagator theForwardPropagator, const Chi2MeasurementEstimatorBase theEstimator, const TransientTrackingRecHitBuilder theTTRHBuilder) const

const Propagator *	backwardPropagator (const TrajectorySeed &seed) const

TempTrajectory	createStartingTrajectory (const TrajectorySeed &seed) const

virtual void	fillSeedHistoDebugger (std::vector< TrajectoryMeasurement >::iterator begin, std::vector< TrajectoryMeasurement >::iterator end) const

StateAndLayers	findStateAndLayers (const TrajectorySeed &seed, const TempTrajectory &traj) const

const Propagator *	forwardPropagator (const TrajectorySeed &seed) const

void	moveToResult (TempTrajectory &&traj, TempTrajectoryContainer &result, bool inOut=false) const

bool	qualityFilter (const TempTrajectory &traj, bool inOut=false) const

void	setData (const MeasurementTrackerEvent *data)

bool	toBeContinued (TempTrajectory &traj, bool inOut=false) const

Protected Attributes
const double	theDeltaEta

const double	theDeltaPhi

edm::ESWatcher< BaseCkfTrajectoryBuilder::Chi2MeasurementEstimatorRecord >	theEstimatorWatcher

std::unique_ptr< Chi2MeasurementEstimatorBase >	theEtaPhiEstimator

const Propagator *	theProximityPropagator

const std::string	theProximityPropagatorName

double	theRescaleErrorIfFail

bool	theUseSeedLayer

Protected Attributes inherited from CkfTrajectoryBuilder
bool	theAlwaysUseInvalidHits

bool	theIntermediateCleaning

float	theLostHitPenalty

int	theMaxCand
	set Event for the internal MeasurementTracker data member More...

Protected Attributes inherited from BaseCkfTrajectoryBuilder
const Chi2MeasurementEstimatorBase *	theEstimator = nullptr

const MeasurementTrackerEvent *	theMeasurementTracker = nullptr

const NavigationSchool *	theNavigationSchool = nullptr

const Propagator *	thePropagatorAlong = nullptr

const Propagator *	thePropagatorOpposite = nullptr

const TransientTrackingRecHitBuilder *	theTTRHBuilder = nullptr

const TrajectoryStateUpdator *	theUpdator = nullptr

Additional Inherited Members
Public Types inherited from CkfTrajectoryBuilder
typedef std::vector< TempTrajectory >	TempTrajectoryContainer

typedef std::vector< Trajectory >	TrajectoryContainer

Public Types inherited from BaseCkfTrajectoryBuilder
typedef std::vector< TempTrajectory >	TempTrajectoryContainer

typedef std::vector< Trajectory >	TrajectoryContainer

typedef TrajectoryContainer::iterator	TrajectoryIterator

Public Types inherited from TrajectoryBuilder
typedef std::vector< Trajectory >	TrajectoryContainer

typedef TrajectoryContainer::iterator	TrajectoryIterator

Protected Types inherited from BaseCkfTrajectoryBuilder
typedef TrackingComponentsRecord	Chi2MeasurementEstimatorRecord

typedef FreeTrajectoryState	FTS

typedef std::pair< TSOS, std::vector< const DetLayer * > >	StateAndLayers

typedef TrajectoryMeasurement	TM

typedef TrajectoryStateOnSurface	TSOS

Static Protected Member Functions inherited from BaseCkfTrajectoryBuilder
static std::unique_ptr< TrajectoryFilter >	createTrajectoryFilter (const edm::ParameterSet &pset, edm::ConsumesCollector &iC)

Detailed Description

Definition at line 7 of file MuonCkfTrajectoryBuilder.h.

Constructor & Destructor Documentation

◆ MuonCkfTrajectoryBuilder()

MuonCkfTrajectoryBuilder::MuonCkfTrajectoryBuilder	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &	iC
	)

Definition at line 18 of file MuonCkfTrajectoryBuilder.cc.

     : CkfTrajectoryBuilder(conf, iC),
       theDeltaEta(conf.getParameter<double>("deltaEta")),
       theDeltaPhi(conf.getParameter<double>("deltaPhi")),
       theProximityPropagatorName(conf.getParameter<std::string>("propagatorProximity")),
       theProximityPropagator(nullptr),
       theEtaPhiEstimator(nullptr) {
   //and something specific to me ?
   theUseSeedLayer = conf.getParameter<bool>("useSeedLayer");
   theRescaleErrorIfFail = conf.getParameter<double>("rescaleErrorIfFail");
 }

References edm::ParameterSet::getParameter(), theRescaleErrorIfFail, and theUseSeedLayer.

◆ ~MuonCkfTrajectoryBuilder()

MuonCkfTrajectoryBuilder::~MuonCkfTrajectoryBuilder ( )

override

Definition at line 30 of file MuonCkfTrajectoryBuilder.cc.

30 {}

Member Function Documentation

◆ collectMeasurement()

void MuonCkfTrajectoryBuilder::collectMeasurement	(	const DetLayer *	layer,
		const std::vector< const DetLayer * > &	nl,
		const TrajectoryStateOnSurface &	currentState,
		std::vector< TM > &	result,
		int &	invalidHits,
		const Propagator *	prop
	)		const

protected

Definition at line 74 of file MuonCkfTrajectoryBuilder.cc.

                                                                                 {
   for (std::vector<const DetLayer*>::const_iterator il = nl.begin(); il != nl.end(); il++) {
     TSOS stateToUse = currentState;
  
     if (layer == (*il)) {
       LogDebug("CkfPattern") << " self propagating in findCompatibleMeasurements.\n from: \n" << stateToUse;
       //self navigation case
       // go to a middle point first
       TransverseImpactPointExtrapolator middle;
       GlobalPoint center(0, 0, 0);
       stateToUse = middle.extrapolate(stateToUse, center, *prop);
  
       if (!stateToUse.isValid())
         continue;
       LogDebug("CkfPattern") << "to: " << stateToUse;
     }
  
     LayerMeasurements layerMeasurements(theMeasurementTracker->measurementTracker(), *theMeasurementTracker);
     std::vector<TM> tmp = layerMeasurements.measurements((**il), stateToUse, *prop, *theEstimator);
  
     if (tmp.size() == 1 && theEtaPhiEstimator) {
       LogDebug("CkfPattern") << "only an invalid hit is found. trying differently";
       tmp = layerMeasurements.measurements((**il), stateToUse, *prop, *theEtaPhiEstimator);
     }
     LogDebug("CkfPattern") << tmp.size() << " measurements returned by LayerMeasurements";
  
     if (!tmp.empty()) {
       // FIXME durty-durty-durty cleaning: never do that please !
       /*      for (vector<TM>::iterator it = tmp.begin(); it!=tmp.end(); ++it)
               {if (it->recHit()->det()==0) it=tmp.erase(it)--;}*/
  
       if (result.empty())
         result = tmp;
       else {
         // keep one dummy TM at the end, skip the others
         result.insert(result.end() - invalidHits, tmp.begin(), tmp.end());
       }
       invalidHits++;
     }
   }
  
   LogDebug("CkfPattern") << "starting from:\n"
                          << "x: " << currentState.globalPosition() << "\n"
                          << "p: " << currentState.globalMomentum() << "\n"
                          << PrintoutHelper::dumpMeasurements(result);
 }

References PrintoutHelper::dumpMeasurements(), TransverseImpactPointExtrapolator::extrapolate(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), LogDebug, MeasurementTrackerEvent::measurementTracker(), mps_fire::result, BaseCkfTrajectoryBuilder::theEstimator, theEtaPhiEstimator, BaseCkfTrajectoryBuilder::theMeasurementTracker, and createJobs::tmp.

Referenced by findCompatibleMeasurements().

◆ findCompatibleMeasurements()

void MuonCkfTrajectoryBuilder::findCompatibleMeasurements	(	const TrajectorySeed &	seed,
		const TempTrajectory &	traj,
		std::vector< TrajectoryMeasurement > &	result
	)		const

overrideprotectedvirtual

Reimplemented from CkfTrajectoryBuilder.

Definition at line 126 of file MuonCkfTrajectoryBuilder.cc.

                                                                                                           {
   int invalidHits = 0;
  
   std::vector<const DetLayer*> nl;
  
   if (traj.empty()) {
     LogDebug("CkfPattern") << "using JR patch for no measurement case";
     //what if there are no measurement on the Trajectory
  
     //set the currentState to be the one from the trajectory seed starting point
     PTrajectoryStateOnDet ptod = seed.startingState();
     DetId id(ptod.detId());
     const GeomDet* g = theMeasurementTracker->geomTracker()->idToDet(id);
     const Surface* surface = &g->surface();
  
     TrajectoryStateOnSurface currentState(
         trajectoryStateTransform::transientState(ptod, surface, forwardPropagator(seed)->magneticField()));
  
     //set the next layers to be that one the state is on
     const DetLayer* l = theMeasurementTracker->geometricSearchTracker()->detLayer(id);
  
     if (theUseSeedLayer) {
       {
         //get the measurements on the layer first
         LogDebug("CkfPattern") << "using the layer of the seed first.";
         nl.push_back(l);
         collectMeasurement(l, nl, currentState, result, invalidHits, theProximityPropagator);
       }
  
       //if fails: try to rescale locally the state to find measurements
       if ((unsigned int)invalidHits == result.size() && theRescaleErrorIfFail != 1.0 && !result.empty()) {
         result.clear();
         LogDebug("CkfPattern") << "using a rescale by " << theRescaleErrorIfFail << " to find measurements.";
         TrajectoryStateOnSurface rescaledCurrentState = currentState;
         rescaledCurrentState.rescaleError(theRescaleErrorIfFail);
         invalidHits = 0;
         collectMeasurement(l, nl, rescaledCurrentState, result, invalidHits, theProximityPropagator);
       }
     }
  
     //if fails: go to next layers
     if (result.empty() || (unsigned int)invalidHits == result.size()) {
       result.clear();
       LogDebug("CkfPattern") << "Need to go to next layer to get measurements";
       //the following will "JUMP" the first layer measurements
       nl = theNavigationSchool->nextLayers(*l, *currentState.freeState(), traj.direction());
       if (nl.empty()) {
         LogDebug("CkfPattern") << " there was no next layer with wellInside. Use the next with no check.";
         //means you did not get any compatible layer on the next 1/2 tracker layer.
         // use the next layers with no checking
         nl = theNavigationSchool->nextLayers(*l, ((traj.direction() == alongMomentum) ? insideOut : outsideIn));
       }
       invalidHits = 0;
       collectMeasurement(l, nl, currentState, result, invalidHits, forwardPropagator(seed));
     }
  
     //if fails: this is on the next layers already, try rescaling locally the state
     if (!result.empty() && (unsigned int)invalidHits == result.size() && theRescaleErrorIfFail != 1.0) {
       result.clear();
       LogDebug("CkfPattern") << "using a rescale by " << theRescaleErrorIfFail
                              << " to find measurements on next layers.";
       TrajectoryStateOnSurface rescaledCurrentState = currentState;
       rescaledCurrentState.rescaleError(theRescaleErrorIfFail);
       invalidHits = 0;
       collectMeasurement(l, nl, rescaledCurrentState, result, invalidHits, forwardPropagator(seed));
     }
  
   } else  //regular case
   {
     TSOS currentState(traj.lastMeasurement().updatedState());
  
     nl = theNavigationSchool->nextLayers(*traj.lastLayer(), *currentState.freeState(), traj.direction());
     if (nl.empty()) {
       LogDebug("CkfPattern") << " no next layers... going " << traj.direction() << "\n from: \n"
                              << currentState
                              << "\n from detId: " << traj.lastMeasurement().recHit()->geographicalId().rawId();
       return;
     }
  
     collectMeasurement(traj.lastLayer(), nl, currentState, result, invalidHits, forwardPropagator(seed));
   }
  
   // sort the final result, keep dummy measurements at the end
   if (result.size() > 1) {
     sort(result.begin(), result.end() - invalidHits, TrajMeasLessEstim());
   }
  
 #ifdef DEBUG_INVALID
   bool afterInvalid = false;
   for (std::vector<TM>::const_iterator i = result.begin(); i != result.end(); i++) {
     if (!i->recHit().isValid())
       afterInvalid = true;
     if (afterInvalid && i->recHit().isValid()) {
       edm::LogError("CkfPattern") << "CkfTrajectoryBuilder error: valid hit after invalid!";
     }
   }
 #endif
  
   //analyseMeasurements( result, traj);
 }

References alongMomentum, collectMeasurement(), PTrajectoryStateOnDet::detId(), GeometricSearchTracker::detLayer(), TempTrajectory::direction(), TempTrajectory::empty(), BaseCkfTrajectoryBuilder::forwardPropagator(), TrajectoryStateOnSurface::freeState(), g, MeasurementTrackerEvent::geometricSearchTracker(), MeasurementTrackerEvent::geomTracker(), mps_fire::i, triggerObjects_cff::id, TrackingGeometry::idToDet(), insideOut, createfilelist::int, cmsLHEtoEOSManager::l, TempTrajectory::lastLayer(), TempTrajectory::lastMeasurement(), LogDebug, HLT_2018_cff::magneticField, NavigationSchool::nextLayers(), outsideIn, TrajectoryMeasurement::recHit(), TrajectoryStateOnSurface::rescaleError(), mps_fire::result, SurveyInfoScenario_cff::seed, BaseCkfTrajectoryBuilder::theMeasurementTracker, BaseCkfTrajectoryBuilder::theNavigationSchool, theProximityPropagator, theRescaleErrorIfFail, theUseSeedLayer, trajectoryStateTransform::transientState(), and TrajectoryMeasurement::updatedState().

◆ setEvent_()

void MuonCkfTrajectoryBuilder::setEvent_	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprotectedvirtual

Implements BaseCkfTrajectoryBuilder.

Definition at line 32 of file MuonCkfTrajectoryBuilder.cc.

                                                                                             {
   CkfTrajectoryBuilder::setEvent_(iEvent, iSetup);
  
   edm::ESHandle<Propagator> propagatorProximityHandle;
   iSetup.get<TrackingComponentsRecord>().get(theProximityPropagatorName, propagatorProximityHandle);
   theProximityPropagator = propagatorProximityHandle.product();
  
   // theEstimator is set for this event in the base class
   if (theEstimatorWatcher.check(iSetup) && theDeltaEta > 0 && theDeltaPhi > 0)
     theEtaPhiEstimator = std::make_unique<EtaPhiEstimator>(theDeltaEta, theDeltaPhi, theEstimator);
 }

References edm::EventSetup::get(), get, iEvent, edm::ESHandle< T >::product(), CkfTrajectoryBuilder::setEvent_(), theDeltaEta, theDeltaPhi, BaseCkfTrajectoryBuilder::theEstimator, theEstimatorWatcher, theEtaPhiEstimator, theProximityPropagator, and theProximityPropagatorName.

Member Data Documentation

◆ theDeltaEta

const double MuonCkfTrajectoryBuilder::theDeltaEta

protected

Definition at line 29 of file MuonCkfTrajectoryBuilder.h.

Referenced by setEvent_().

◆ theDeltaPhi

const double MuonCkfTrajectoryBuilder::theDeltaPhi

protected

Definition at line 30 of file MuonCkfTrajectoryBuilder.h.

Referenced by setEvent_().

◆ theEstimatorWatcher

edm::ESWatcher<BaseCkfTrajectoryBuilder::Chi2MeasurementEstimatorRecord> MuonCkfTrajectoryBuilder::theEstimatorWatcher

protected

Definition at line 33 of file MuonCkfTrajectoryBuilder.h.

Referenced by setEvent_().

◆ theEtaPhiEstimator

std::unique_ptr<Chi2MeasurementEstimatorBase> MuonCkfTrajectoryBuilder::theEtaPhiEstimator

protected

Definition at line 34 of file MuonCkfTrajectoryBuilder.h.

Referenced by collectMeasurement(), and setEvent_().

◆ theProximityPropagator

const Propagator* MuonCkfTrajectoryBuilder::theProximityPropagator

protected

Definition at line 32 of file MuonCkfTrajectoryBuilder.h.

Referenced by findCompatibleMeasurements(), and setEvent_().

◆ theProximityPropagatorName

const std::string MuonCkfTrajectoryBuilder::theProximityPropagatorName

protected

Definition at line 31 of file MuonCkfTrajectoryBuilder.h.

Referenced by setEvent_().

◆ theRescaleErrorIfFail

double MuonCkfTrajectoryBuilder::theRescaleErrorIfFail

protected

Definition at line 28 of file MuonCkfTrajectoryBuilder.h.

Referenced by findCompatibleMeasurements(), and MuonCkfTrajectoryBuilder().

◆ theUseSeedLayer

bool MuonCkfTrajectoryBuilder::theUseSeedLayer

protected

Definition at line 27 of file MuonCkfTrajectoryBuilder.h.

Referenced by findCompatibleMeasurements(), and MuonCkfTrajectoryBuilder().

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ MuonCkfTrajectoryBuilder()

◆ ~MuonCkfTrajectoryBuilder()

Member Function Documentation

◆ collectMeasurement()

◆ findCompatibleMeasurements()

◆ setEvent_()

Member Data Documentation

◆ theDeltaEta

◆ theDeltaPhi

◆ theEstimatorWatcher

◆ theEtaPhiEstimator

◆ theProximityPropagator

◆ theProximityPropagatorName

◆ theRescaleErrorIfFail

◆ theUseSeedLayer