#include <MuonCkfTrajectoryBuilder.h>

Inheritance diagram for MuonCkfTrajectoryBuilder:

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

virtual	~MuonCkfTrajectoryBuilder ()

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

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

	CkfTrajectoryBuilder (const edm::ParameterSet &conf, TrajectoryFilter *filter)

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

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

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

	~CkfTrajectoryBuilder ()

Public Member Functions inherited from BaseCkfTrajectoryBuilder
	BaseCkfTrajectoryBuilder (const edm::ParameterSet &conf, TrajectoryFilter filter, TrajectoryFilter inOutFilter=0)

	BaseCkfTrajectoryBuilder (const BaseCkfTrajectoryBuilder &)=delete

const TransientTrackingRecHitBuilder *	hitBuilder () const

BaseCkfTrajectoryBuilder &	operator= (const BaseCkfTrajectoryBuilder &)=delete

virtual void	setDebugger (CkfDebugger *dbg) const

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

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

void	setNavigationSchool (NavigationSchool const *nv)

virtual void	unset () const

virtual	~BaseCkfTrajectoryBuilder ()

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

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

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

Protected Member Functions inherited from CkfTrajectoryBuilder
void	limitedCandidates (const TrajectorySeed &seed, TempTrajectory &startingTraj, TrajectoryContainer &result) const

void	limitedCandidates (const boost::shared_ptr< const TrajectorySeed > &sharedSeed, TempTrajectoryContainer &candidates, 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 (boost::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 (Trajectory &traj, const std::vector< TrajectoryMeasurement > &meas, const MeasurementTrackerEvent theMeasurementTracker, const Propagator theForwardPropagator, const Chi2MeasurementEstimatorBase theEstimator, const TransientTrackingRecHitBuilder theTTRHBuilder) const

virtual bool	analyzeMeasurementsDebugger (TempTrajectory &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

StateAndLayers	findStateAndLayers (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

const MeasurementTrackerEvent *	theMeasurementTracker

const NavigationSchool *	theNavigationSchool = 0

const Propagator *	thePropagatorAlong

const Propagator *	thePropagatorOpposite

const TransientTrackingRecHitBuilder *	theTTRHBuilder

const TrajectoryStateUpdator *	theUpdator

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 TrajectoryFilter *	createTrajectoryFilter (const edm::ParameterSet &pset, edm::ConsumesCollector &iC)

Detailed Description

Definition at line 7 of file MuonCkfTrajectoryBuilder.h.

Constructor & Destructor Documentation

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

Definition at line 18 of file MuonCkfTrajectoryBuilder.cc.

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

                                                                                                        :
   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");
 }

MuonCkfTrajectoryBuilder::~MuonCkfTrajectoryBuilder ( )

virtual

Definition at line 31 of file MuonCkfTrajectoryBuilder.cc.

32 {

33 }

Member Function Documentation

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 78 of file MuonCkfTrajectoryBuilder.cc.

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

Referenced by findCompatibleMeasurements().

                                                         {
   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);
 }

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

protectedvirtual

Reimplemented from CkfTrajectoryBuilder.

Definition at line 133 of file MuonCkfTrajectoryBuilder.cc.

References alongMomentum, collectMeasurement(), PTrajectoryStateOnDet::detId(), GeometricSearchTracker::detLayer(), TempTrajectory::direction(), TempTrajectory::empty(), BaseCkfTrajectoryBuilder::forwardPropagator(), TrajectoryStateOnSurface::freeState(), g, MeasurementTrackerEvent::geometricSearchTracker(), MeasurementTrackerEvent::geomTracker(), i, TrackingGeometry::idToDet(), insideOut, prof2calltree::l, TempTrajectory::lastLayer(), TempTrajectory::lastMeasurement(), LogDebug, NavigationSchool::nextLayers(), outsideIn, TrajectoryMeasurement::recHit(), TrajectoryStateOnSurface::rescaleError(), query::result, mathSSE::return(), python.multivaluedict::sort(), TrajectorySeed::startingState(), BaseCkfTrajectoryBuilder::theMeasurementTracker, BaseCkfTrajectoryBuilder::theNavigationSchool, theProximityPropagator, theRescaleErrorIfFail, theUseSeedLayer, trajectoryStateTransform::transientState(), and TrajectoryMeasurement::updatedState().

 {
   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.size()!=0)
           {
         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.size()==0 || (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.size()==0){
             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.size()!=0 && (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);
 
 }

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

overrideprotectedvirtual

Implements BaseCkfTrajectoryBuilder.

Definition at line 35 of file MuonCkfTrajectoryBuilder.cc.

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

                                                                                             {
   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.reset(new EtaPhiEstimator(theDeltaEta, theDeltaPhi, theEstimator));
 }

Member Data Documentation

const double MuonCkfTrajectoryBuilder::theDeltaEta

protected

Definition at line 22 of file MuonCkfTrajectoryBuilder.h.

Referenced by setEvent_().

const double MuonCkfTrajectoryBuilder::theDeltaPhi

protected

Definition at line 23 of file MuonCkfTrajectoryBuilder.h.

Referenced by setEvent_().

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

protected

Definition at line 26 of file MuonCkfTrajectoryBuilder.h.

Referenced by setEvent_().

std::unique_ptr<Chi2MeasurementEstimatorBase> MuonCkfTrajectoryBuilder::theEtaPhiEstimator

protected

Definition at line 27 of file MuonCkfTrajectoryBuilder.h.

Referenced by collectMeasurement(), and setEvent_().

const Propagator* MuonCkfTrajectoryBuilder::theProximityPropagator

protected

Definition at line 25 of file MuonCkfTrajectoryBuilder.h.

Referenced by findCompatibleMeasurements(), and setEvent_().

const std::string MuonCkfTrajectoryBuilder::theProximityPropagatorName

protected

Definition at line 24 of file MuonCkfTrajectoryBuilder.h.

Referenced by setEvent_().

double MuonCkfTrajectoryBuilder::theRescaleErrorIfFail

protected

Definition at line 21 of file MuonCkfTrajectoryBuilder.h.

Referenced by findCompatibleMeasurements(), and MuonCkfTrajectoryBuilder().

bool MuonCkfTrajectoryBuilder::theUseSeedLayer

protected

Definition at line 20 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

Member Function Documentation

Member Data Documentation