#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 void	setEvent (const edm::Event &event) const
	set Event for the internal MeasurementTracker data member More...

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	~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

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, ExpressReco_HICollisions_FallBack::id, TrackingGeometry::idToDet(), insideOut, prof2calltree::l, TempTrajectory::lastLayer(), TempTrajectory::lastMeasurement(), LogDebug, Propagator::magneticField(), DetLayer::nextLayers(), outsideIn, TrajectoryMeasurement::recHit(), TrajectoryStateOnSurface::rescaleError(), query::result, mathSSE::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();
       TrajectoryStateTransform tsTransform;
       TrajectoryStateOnSurface currentState(tsTransform.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);
 
 }