#include <GsfTrajectoryFitter.h>

Inheritance diagram for GsfTrajectoryFitter:

Public Member Functions
virtual GsfTrajectoryFitter *	clone () const

const MeasurementEstimator *	estimator () const

virtual std::vector< Trajectory >	fit (const Trajectory &aTraj) const

virtual std::vector< Trajectory >	fit (const TrajectorySeed &aSeed, const RecHitContainer &hits) const

virtual std::vector< Trajectory >	fit (const TrajectorySeed &aSeed, const RecHitContainer &hits, const TSOS &firstPredTsos) const

	GsfTrajectoryFitter (const Propagator &aPropagator, const TrajectoryStateUpdator &aUpdator, const MeasurementEstimator &aEstimator, const MultiTrajectoryStateMerger &aMerger, const DetLayerGeometry *detLayerGeometry=0)

const MultiTrajectoryStateMerger *	merger () const

const Propagator *	propagator () const

const TrajectoryStateUpdator *	updator () const

virtual	~GsfTrajectoryFitter ()

Public Member Functions inherited from TrajectoryFitter
virtual std::vector< Trajectory >	fit (const Trajectory &traj, fitType type) const

virtual std::vector< Trajectory >	fit (const TrajectorySeed &seed, const RecHitContainer &hits, fitType type) const

virtual std::vector< Trajectory >	fit (const TrajectorySeed &seed, const RecHitContainer &hits, const TrajectoryStateOnSurface &tsos, fitType type) const

virtual	~TrajectoryFitter ()

Private Types
typedef TrajectoryMeasurement	TM

typedef TrajectoryStateOnSurface	TSOS

Private Attributes
const DetLayerGeometry	dummyGeometry

const MeasurementEstimator *	theEstimator

const DetLayerGeometry *	theGeometry

const MultiTrajectoryStateMerger *	theMerger

const Propagator *	thePropagator

bool	theTiming

const TrajectoryStateUpdator *	theUpdator

Additional Inherited Members
Public Types inherited from TrajectoryFitter
enum	fitType { standard =0, looper =1 }

typedef Trajectory::RecHitContainer	RecHitContainer

typedef TrajectoryFitterRecord	Record

Detailed Description

A GSF fitter, similar to KFTrajectoryFitter.

Definition at line 20 of file GsfTrajectoryFitter.h.

Member Typedef Documentation

typedef TrajectoryMeasurement GsfTrajectoryFitter::TM

private

Definition at line 24 of file GsfTrajectoryFitter.h.

typedef TrajectoryStateOnSurface GsfTrajectoryFitter::TSOS

private

Definition at line 23 of file GsfTrajectoryFitter.h.

Constructor & Destructor Documentation

GsfTrajectoryFitter::GsfTrajectoryFitter	(	const Propagator &	aPropagator,
		const TrajectoryStateUpdator &	aUpdator,
		const MeasurementEstimator &	aEstimator,
		const MultiTrajectoryStateMerger &	aMerger,
		const DetLayerGeometry *	detLayerGeometry = `0`
	)

Constructor with explicit components for propagation, update, chi2 calculation, merging and flag for merging before / after the update (i.e. fully configured)

Definition at line 12 of file GsfTrajectoryFitter.cc.

References dummyGeometry, and theGeometry.

Referenced by clone().

                                                                :
   thePropagator(aPropagator.clone()),
   theUpdator(aUpdator.clone()),
   theEstimator(aEstimator.clone()),
   theMerger(aMerger.clone()),
   theGeometry(detLayerGeometry)
 {
   if(!theGeometry) theGeometry = &dummyGeometry;
   //   static SimpleConfigurable<bool> timeConf(false,"GsfTrajectoryFitter:activateTiming");
   //   theTiming = timeConf.value();
 }

GsfTrajectoryFitter::~GsfTrajectoryFitter ( )

virtual

Definition at line 28 of file GsfTrajectoryFitter.cc.

References theEstimator, theMerger, thePropagator, and theUpdator.

                                           {
   delete thePropagator;
   delete theUpdator;
   delete theEstimator;
   delete theMerger;
 }

Member Function Documentation

virtual GsfTrajectoryFitter* GsfTrajectoryFitter::clone ( void ) const

inlinevirtual

Implements TrajectoryFitter.

Definition at line 51 of file GsfTrajectoryFitter.h.

References GsfTrajectoryFitter(), theEstimator, theGeometry, theMerger, thePropagator, and theUpdator.

   {
     return new GsfTrajectoryFitter(*thePropagator,*theUpdator,*theEstimator,*theMerger,theGeometry);
   }

const MeasurementEstimator* GsfTrajectoryFitter::estimator ( void ) const

inline

Definition at line 48 of file GsfTrajectoryFitter.h.

References theEstimator.

Referenced by fit().

48 {return theEstimator;}

GsfTrajectoryFitter::theEstimator

const MeasurementEstimator * theEstimator

Definition: GsfTrajectoryFitter.h:59

std::vector< Trajectory > GsfTrajectoryFitter::fit ( const Trajectory & aTraj ) const

virtual

Implements TrajectoryFitter.

Definition at line 35 of file GsfTrajectoryFitter.cc.

References Trajectory::empty(), Trajectory::firstMeasurement(), Trajectory::recHits(), Trajectory::seed(), and TrajectoryMeasurement::updatedState().

 {  
   if(aTraj.empty()) return std::vector<Trajectory>();
  
   TM firstTM = aTraj.firstMeasurement();
   TSOS firstTsos = TrajectoryStateWithArbitraryError()(firstTM.updatedState());
   
   return fit(aTraj.seed(), aTraj.recHits(), firstTsos);
 }

std::vector< Trajectory > GsfTrajectoryFitter::fit	(	const TrajectorySeed &	aSeed,
		const RecHitContainer &	hits
	)		const

virtual

Implements TrajectoryFitter.

Definition at line 45 of file GsfTrajectoryFitter.cc.

                                                             {
 
   edm::LogError("GsfTrajectoryFitter") 
     << "GsfTrajectoryFitter::fit(TrajectorySeed, vector<RecHit>) not implemented";
 
   return std::vector<Trajectory>();
 }

std::vector< Trajectory > GsfTrajectoryFitter::fit	(	const TrajectorySeed &	aSeed,
		const RecHitContainer &	hits,
		const TSOS &	firstPredTsos
	)		const

virtual

!! no invalid hits on cylinders anymore??

Implements TrajectoryFitter.

Definition at line 54 of file GsfTrajectoryFitter.cc.

References MeasurementEstimator::estimate(), estimator(), Trajectory::foundHits(), DetLayerGeometry::idToLayer(), TrajectoryStateOnSurface::isValid(), LogDebug, MultiTrajectoryStateMerger::merge(), merger(), Propagator::propagate(), propagator(), Trajectory::push(), edm::second(), TrajectoryStateOnSurface::surface(), theGeometry, TrajectoryStateUpdator::update(), and updator().

                                                           {
 
   //   static TimingReport::Item* propTimer =
   //     &(*TimingReport::current())[string("GsfTrajectoryFitter:propagation")];
   //   propTimer->switchCPU(false);
   //   if ( !theTiming )  propTimer->switchOn(false);
   //   static TimingReport::Item* updateTimer =
   //     &(*TimingReport::current())[string("GsfTrajectoryFitter:update")];
   //   updateTimer->switchCPU(false);
   //   if ( !theTiming )  updateTimer->switchOn(false);
 
   if(hits.empty()) return std::vector<Trajectory>();
 
   Trajectory myTraj(aSeed, propagator()->propagationDirection());
 
   TSOS predTsos(firstPredTsos);
   if(!predTsos.isValid()) {
     edm::LogInfo("GsfTrajectoryFitter") 
       << "GsfTrajectoryFitter: predicted tsos of first measurement not valid!";
     return std::vector<Trajectory>();
   } 
 
   TSOS currTsos;
   if(hits.front()->isValid()) {
     //update
     TransientTrackingRecHit::RecHitPointer preciseHit = hits.front()->clone(predTsos);
     {
       //       TimeMe t(*updateTimer,false);
       currTsos = updator()->update(predTsos, *preciseHit);
     }
     if (!predTsos.isValid() || !currTsos.isValid()){
       edm::LogError("InvalidState")<<"first hit";
       return std::vector<Trajectory>();
     }
     myTraj.push(TM(predTsos, currTsos, preciseHit, 0., theGeometry->idToLayer(preciseHit->geographicalId() )),
         estimator()->estimate(predTsos, *preciseHit).second);
   } else {
     currTsos = predTsos;
     if (!predTsos.isValid()){
       edm::LogError("InvalidState")<<"first invalid hit";
       return std::vector<Trajectory>();
     }
     myTraj.push(TM(predTsos, *hits.begin(),0., theGeometry->idToLayer((*hits.begin())->geographicalId()) ));
   }
   
   for(RecHitContainer::const_iterator ihit = hits.begin() + 1; 
       ihit != hits.end(); ihit++) {    
     //
     // temporary protection copied from KFTrajectoryFitter.
     //
     if ((**ihit).isValid() == false && (**ihit).det() == 0) {
       LogDebug("GsfTrajectoryFitter") << " Error: invalid hit with no GeomDet attached .... skipping";
       continue;
     }
 
     //     //
     //     // check type of surface in case of invalid hit
     //     // (in this version only propagations to planes are
     //     // supported for multi trajectory states)
     //     //
     //     if ( !(**ihit).isValid() ) {
     //       const BoundPlane* plane = 
     //  dynamic_cast<const BoundPlane*>(&(**ihit).det().surface());
     //       //
     //       // no plane: insert invalid measurement
     //       //
     //       if ( plane==0 ) {
     //  myTraj.push(TM(TrajectoryStateOnSurface(),&(**ihit)));
     //  continue;
     //       }
     //     }
     {
       //       TimeMe t(*propTimer,false);
       predTsos = propagator()->propagate(currTsos,
                      (**ihit).det()->surface());
     }
     if(!predTsos.isValid()) {
       if ( myTraj.foundHits()>=3 ) {
     edm::LogInfo("GsfTrajectoryFitter") 
       << "GsfTrajectoryFitter: predicted tsos not valid! \n"
       << "Returning trajectory with " << myTraj.foundHits() << " found hits.";
     return std::vector<Trajectory>(1,myTraj);
       }
       else {
       edm::LogInfo("GsfTrajectoryFitter") 
     << "GsfTrajectoryFitter: predicted tsos not valid after " << myTraj.foundHits()
     << " hits, discarding candidate!";
     return std::vector<Trajectory>();
       }
     }
     if ( merger() ) predTsos = merger()->merge(predTsos);
     
     if((**ihit).isValid()) {
       //update
       TransientTrackingRecHit::RecHitPointer preciseHit = (**ihit).clone(predTsos);
       {
     //  TimeMe t(*updateTimer,false);
     currTsos = updator()->update(predTsos, *preciseHit);
       }
       if (!predTsos.isValid() || !currTsos.isValid()){
     edm::LogError("InvalidState")<<"inside hit";
     return std::vector<Trajectory>();
       }
       myTraj.push(TM(predTsos, currTsos, preciseHit,
              estimator()->estimate(predTsos, *preciseHit).second,
              theGeometry->idToLayer(preciseHit->geographicalId() )));
     } else {
       currTsos = predTsos;
       if (!predTsos.isValid()){
       edm::LogError("InvalidState")<<"inside invalid hit";
       return std::vector<Trajectory>();
       }
       myTraj.push(TM(predTsos, *ihit,0., theGeometry->idToLayer( (*ihit)->geographicalId()) ));
     }
   }
   return std::vector<Trajectory>(1, myTraj);
 }