#include <MuonCkfTrajectoryBuilder.h>

Inheritance diagram for MuonCkfTrajectoryBuilder:

Public Member Functions
	MuonCkfTrajectoryBuilder (const edm::ParameterSet &conf, const TrajectoryStateUpdator updator, const Propagator propagatorAlong, const Propagator propagatorOpposite, const Propagator propagatorProximity, const Chi2MeasurementEstimatorBase estimator, const TransientTrackingRecHitBuilder RecHitBuilder, const MeasurementTracker measurementTracker, const TrajectoryFilter filter)

virtual	~MuonCkfTrajectoryBuilder ()

Public Member Functions inherited from CkfTrajectoryBuilder
	CkfTrajectoryBuilder (const edm::ParameterSet &conf, const TrajectoryStateUpdator updator, const Propagator propagatorAlong, const Propagator propagatorOpposite, const Chi2MeasurementEstimatorBase estimator, const TransientTrackingRecHitBuilder recHitBuilder, const MeasurementTracker measurementTracker, const TrajectoryFilter *filter)

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, const TrajectoryStateUpdator updator, const Propagator propagatorAlong, const Propagator propagatorOpposite, const Chi2MeasurementEstimatorBase estimator, const TransientTrackingRecHitBuilder RecHitBuilder, const MeasurementTracker measurementTracker, const TrajectoryFilter filter, const TrajectoryFilter inOutFilter=0)

virtual void	setDebugger (CkfDebugger *dbg) const

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

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 TempTrajectory &traj, std::vector< TrajectoryMeasurement > &result) const

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

void	limitedCandidates (TempTrajectoryContainer &candidates, TrajectoryContainer &result) const

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

Protected Member Functions inherited from BaseCkfTrajectoryBuilder
void	addToResult (TempTrajectory &traj, TrajectoryContainer &result, bool inOut=false) const

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

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

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

TempTrajectory	createStartingTrajectory (const TrajectorySeed &seed) const

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

StateAndLayers	findStateAndLayers (const TempTrajectory &traj) const

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

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

Protected Attributes
Chi2MeasurementEstimatorBase *	theEtaPhiEstimator

const Propagator *	theProximityPropagator

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 Propagator *	theBackwardPropagator

const Chi2MeasurementEstimatorBase *	theEstimator

const Propagator *	theForwardPropagator

const LayerMeasurements *	theLayerMeasurements

const MeasurementTracker *	theMeasurementTracker

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 FreeTrajectoryState	FTS

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

typedef TrajectoryMeasurement	TM

typedef TrajectoryStateOnSurface	TSOS

Detailed Description

Definition at line 6 of file MuonCkfTrajectoryBuilder.h.

Constructor & Destructor Documentation

MuonCkfTrajectoryBuilder::MuonCkfTrajectoryBuilder	(	const edm::ParameterSet &	conf,
		const TrajectoryStateUpdator *	updator,
		const Propagator *	propagatorAlong,
		const Propagator *	propagatorOpposite,
		const Propagator *	propagatorProximity,
		const Chi2MeasurementEstimatorBase *	estimator,
		const TransientTrackingRecHitBuilder *	RecHitBuilder,
		const MeasurementTracker *	measurementTracker,
		const TrajectoryFilter *	filter
	)

Definition at line 16 of file MuonCkfTrajectoryBuilder.cc.

References dPhi(), edm::ParameterSet::getParameter(), BaseCkfTrajectoryBuilder::theEstimator, theEtaPhiEstimator, theRescaleErrorIfFail, and theUseSeedLayer.

                                                                         : 
   CkfTrajectoryBuilder(conf,updator,propagatorAlong,propagatorOpposite,estimator,RecHitBuilder,measurementTracker,filter),
   theProximityPropagator(propagatorProximity)
 {
   //and something specific to me ?
   theUseSeedLayer = conf.getParameter<bool>("useSeedLayer");
   theRescaleErrorIfFail = conf.getParameter<double>("rescaleErrorIfFail");
   double dEta=conf.getParameter<double>("deltaEta");
   double dPhi=conf.getParameter<double>("deltaPhi");
   if (dEta>0 && dPhi>0)
     theEtaPhiEstimator = new EtaPhiEstimator(dEta,dPhi,theEstimator);
   else theEtaPhiEstimator = (Chi2MeasurementEstimatorBase*)0;
 
 
 }

MuonCkfTrajectoryBuilder::~MuonCkfTrajectoryBuilder ( )

virtual

Definition at line 40 of file MuonCkfTrajectoryBuilder.cc.

References theEtaPhiEstimator.

 {
   if (theEtaPhiEstimator) delete theEtaPhiEstimator;
 }

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

References PrintoutHelper::dumpMeasurements(), TransverseImpactPointExtrapolator::extrapolate(), TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), LogDebug, LayerMeasurements::measurements(), BaseCkfTrajectoryBuilder::theEstimator, theEtaPhiEstimator, BaseCkfTrajectoryBuilder::theLayerMeasurements, 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;
     }
 
     std::vector<TM> tmp =
       theLayerMeasurements->measurements((**il),stateToUse, *prop, *theEstimator);
     
     if (tmp.size()==1 && theEtaPhiEstimator){
       LogDebug("CkfPattern")<<"only an invalid hit is found. trying differently";
       tmp = theLayerMeasurements->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 TempTrajectory &	traj,
		std::vector< TrajectoryMeasurement > &	result
	)		const

protectedvirtual

Reimplemented from CkfTrajectoryBuilder.

Definition at line 129 of file MuonCkfTrajectoryBuilder.cc.

References alongMomentum, collectMeasurement(), PTrajectoryStateOnDet::detId(), GeometricSearchTracker::detLayer(), TempTrajectory::direction(), TempTrajectory::empty(), TrajectoryStateOnSurface::freeState(), g, MeasurementTracker::geometricSearchTracker(), MeasurementTracker::geomTracker(), i, errorMatrix2Lands_multiChannel::id, TrackingGeometry::idToDet(), insideOut, prof2calltree::l, TempTrajectory::lastLayer(), TempTrajectory::lastMeasurement(), LogDebug, Propagator::magneticField(), DetLayer::nextLayers(), outsideIn, TrajectoryMeasurement::recHit(), TrajectoryStateOnSurface::rescaleError(), query::result, hitfit::return, TempTrajectory::seed(), python.multivaluedict::sort(), TrajectorySeed::startingState(), BaseCkfTrajectoryBuilder::theForwardPropagator, BaseCkfTrajectoryBuilder::theMeasurementTracker, 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 =traj.seed().startingState();
       DetId id(ptod.detId());
       const GeomDet * g = theMeasurementTracker->geomTracker()->idToDet(id);
       const Surface * surface=&g->surface();
       
       TrajectoryStateOnSurface currentState(trajectoryStateTransform::transientState(ptod,surface,theForwardPropagator->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 = l->nextLayers(*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 = l->nextLayers(((traj.direction()==alongMomentum)?insideOut:outsideIn));
           }
           invalidHits=0;
           collectMeasurement(l,nl,currentState,result,invalidHits,theForwardPropagator);
         }
 
       //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,theForwardPropagator);
         }
 
     }
   else //regular case
     {
 
       TSOS currentState( traj.lastMeasurement().updatedState());
 
       nl = traj.lastLayer()->nextLayers( *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,theForwardPropagator);
     }
 
 
   // 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);
 
 }