#include <GsfTrajectoryFitter.h>

Inheritance diagram for GsfTrajectoryFitter:

Public Member Functions
virtual std::unique_ptr < TrajectoryFitter >	clone () const override

const MeasurementEstimator *	estimator () const

Trajectory	fitOne (const Trajectory &t, fitType type) const

Trajectory	fitOne (const TrajectorySeed &aSeed, const RecHitContainer &hits, const TrajectoryStateOnSurface &firstPredTsos, fitType type) const

Trajectory	fitOne (const TrajectorySeed &aSeed, const RecHitContainer &hits, fitType type) 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

virtual void	setHitCloner (TkCloner const *hc)

const TrajectoryStateUpdator *	updator () const

virtual	~GsfTrajectoryFitter ()

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

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

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

virtual	~TrajectoryFitter ()

Private Types
typedef TrajectoryMeasurement	TM

typedef TrajectoryStateOnSurface	TSOS

Private Attributes
const DetLayerGeometry	dummyGeometry

const MeasurementEstimator *	theEstimator

const DetLayerGeometry *	theGeometry

TkCloner const *	theHitCloner = 0

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 17 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 33 of file GsfTrajectoryFitter.cc.

References theEstimator, theMerger, thePropagator, and theUpdator.

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

Member Function Documentation

virtual std::unique_ptr<TrajectoryFitter> GsfTrajectoryFitter::clone ( void ) const

inlineoverridevirtual

Implements TrajectoryFitter.

Definition at line 54 of file GsfTrajectoryFitter.h.

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

   {
     return std::unique_ptr<TrajectoryFitter>(
         new GsfTrajectoryFitter(*thePropagator,
                                 *theUpdator,
                                 *theEstimator,
                                 *theMerger,
                                 theGeometry));
   }

const MeasurementEstimator* GsfTrajectoryFitter::estimator ( void ) const

inline

Definition at line 51 of file GsfTrajectoryFitter.h.

References theEstimator.

Referenced by fitOne().

51 {return theEstimator;}

GsfTrajectoryFitter::theEstimator

const MeasurementEstimator * theEstimator

Definition: GsfTrajectoryFitter.h:72

Trajectory GsfTrajectoryFitter::fitOne	(	const Trajectory &	t,
		fitType	type
	)		const

virtual

Implements TrajectoryFitter.

Definition at line 40 of file GsfTrajectoryFitter.cc.

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

                                                                                   {  
   if(aTraj.empty()) return Trajectory();
  
   TM const & firstTM = aTraj.firstMeasurement();
   TSOS firstTsos = TrajectoryStateWithArbitraryError()(firstTM.updatedState());
   
   return fitOne(aTraj.seed(), aTraj.recHits(), firstTsos,type);
 }

Trajectory GsfTrajectoryFitter::fitOne	(	const TrajectorySeed &	aSeed,
		const RecHitContainer &	hits,
		const TrajectoryStateOnSurface &	firstPredTsos,
		fitType	type
	)		const

virtual

!! no invalid hits on cylinders anymore??

Implements TrajectoryFitter.

Definition at line 58 of file GsfTrajectoryFitter.cc.

References assert(), MeasurementEstimator::estimate(), estimator(), Trajectory::foundHits(), DetLayerGeometry::idToLayer(), TrajectoryStateOnSurface::isValid(), LogDebug, TkCloner::makeShared(), MultiTrajectoryStateMerger::merge(), merger(), Propagator::propagate(), propagator(), Trajectory::push(), edm::second(), TrajectoryStateOnSurface::surface(), theGeometry, theHitCloner, 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 Trajectory();
 
   Trajectory myTraj(aSeed, propagator()->propagationDirection());
 
   TSOS predTsos(firstPredTsos);
   if(!predTsos.isValid()) {
     edm::LogInfo("GsfTrajectoryFitter") 
       << "GsfTrajectoryFitter: predicted tsos of first measurement not valid!";
     return Trajectory();
   } 
 
   TSOS currTsos;
   if(hits.front()->isValid()) {
      auto const & ihit = hits.front();
     //update
      assert( (!(ihit)->canImproveWithTrack()) | (nullptr!=theHitCloner));
      assert( (!(ihit)->canImproveWithTrack()) | (nullptr!=dynamic_cast<BaseTrackerRecHit const*>(ihit.get())));
      auto preciseHit = theHitCloner->makeShared(ihit,predTsos);
     {
       //       TimeMe t(*updateTimer,false);
       currTsos = updator()->update(predTsos, *preciseHit);
     }
     if (!predTsos.isValid() || !currTsos.isValid()){
       edm::LogError("InvalidState")<<"first hit";
       return 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 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 myTraj;
       }
       else {
       edm::LogInfo("GsfTrajectoryFitter") 
     << "GsfTrajectoryFitter: predicted tsos not valid after " << myTraj.foundHits()
     << " hits, discarding candidate!";
     return Trajectory();
       }
     }
     if ( merger() ) predTsos = merger()->merge(predTsos);
     
     if((**ihit).isValid()) {
       //update
        assert( (!(*ihit)->canImproveWithTrack()) | (nullptr!=theHitCloner));
        assert( (!(*ihit)->canImproveWithTrack()) | (nullptr!=dynamic_cast<BaseTrackerRecHit const*>((*ihit).get())));
        auto preciseHit = theHitCloner->makeShared(*ihit,predTsos);
       {
     //  TimeMe t(*updateTimer,false);
     currTsos = updator()->update(predTsos, *preciseHit);
       }
       if (!predTsos.isValid() || !currTsos.isValid()){
     edm::LogError("InvalidState")<<"inside hit";
     return 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 Trajectory();
       }
       myTraj.push(TM(predTsos, *ihit,0., theGeometry->idToLayer( (*ihit)->geographicalId()) ));
     }
   }
   return myTraj;
 }

Trajectory GsfTrajectoryFitter::fitOne	(	const TrajectorySeed &	aSeed,
		const RecHitContainer &	hits,
		fitType	type
	)		const

virtual

Implements TrajectoryFitter.

Definition at line 49 of file GsfTrajectoryFitter.cc.

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

const MultiTrajectoryStateMerger* GsfTrajectoryFitter::merger ( ) const

inline

Definition at line 52 of file GsfTrajectoryFitter.h.

References theMerger.

Referenced by fitOne().

52 {return theMerger;}

GsfTrajectoryFitter::theMerger

const MultiTrajectoryStateMerger * theMerger

Definition: GsfTrajectoryFitter.h:73

const Propagator* GsfTrajectoryFitter::propagator ( void ) const

inline

Definition at line 49 of file GsfTrajectoryFitter.h.

References thePropagator.

Referenced by fitOne().

49 {return thePropagator;}

GsfTrajectoryFitter::thePropagator

const Propagator * thePropagator

Definition: GsfTrajectoryFitter.h:70

virtual void GsfTrajectoryFitter::setHitCloner ( TkCloner const * hc )

inlinevirtual

Implements TrajectoryFitter.

Definition at line 64 of file GsfTrajectoryFitter.h.

References AnalysisDataFormats_SUSYBSMObjects::hc, and theHitCloner.

                                                  {
      theHitCloner = hc;
   }

const TrajectoryStateUpdator* GsfTrajectoryFitter::updator ( ) const

inline

Definition at line 50 of file GsfTrajectoryFitter.h.

References theUpdator.

Referenced by fitOne().

50 {return theUpdator;}

GsfTrajectoryFitter::theUpdator

const TrajectoryStateUpdator * theUpdator

Definition: GsfTrajectoryFitter.h:71

Member Data Documentation

const DetLayerGeometry GsfTrajectoryFitter::dummyGeometry

private

Definition at line 75 of file GsfTrajectoryFitter.h.

Referenced by GsfTrajectoryFitter().

const MeasurementEstimator* GsfTrajectoryFitter::theEstimator

private

Definition at line 72 of file GsfTrajectoryFitter.h.

Referenced by clone(), estimator(), and ~GsfTrajectoryFitter().

const DetLayerGeometry* GsfTrajectoryFitter::theGeometry

private

Definition at line 76 of file GsfTrajectoryFitter.h.

Referenced by clone(), fitOne(), and GsfTrajectoryFitter().

TkCloner const* GsfTrajectoryFitter::theHitCloner = 0

private

Definition at line 74 of file GsfTrajectoryFitter.h.

Referenced by fitOne(), and setHitCloner().

const MultiTrajectoryStateMerger* GsfTrajectoryFitter::theMerger

private

Definition at line 73 of file GsfTrajectoryFitter.h.

Referenced by clone(), merger(), and ~GsfTrajectoryFitter().

const Propagator* GsfTrajectoryFitter::thePropagator

private

Definition at line 70 of file GsfTrajectoryFitter.h.

Referenced by clone(), propagator(), and ~GsfTrajectoryFitter().

bool GsfTrajectoryFitter::theTiming

private

Definition at line 78 of file GsfTrajectoryFitter.h.

const TrajectoryStateUpdator* GsfTrajectoryFitter::theUpdator

private

Definition at line 71 of file GsfTrajectoryFitter.h.

Referenced by clone(), updator(), and ~GsfTrajectoryFitter().

Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation