GsfTrajectorySmoother.h

Go to the documentation of this file.

#ifndef GsfTrajectorySmoother_H_
#define GsfTrajectorySmoother_H_

#include "TrackingTools/PatternTools/interface/TrajectorySmoother.h"
#include "TrackingTools/PatternTools/interface/TrajectoryStateUpdator.h"
#include "TrackingTools/GeomPropagators/interface/Propagator.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "TrackingTools/PatternTools/interface/MeasurementEstimator.h"
#include "TrackingTools/PatternTools/interface/TrajectoryMeasurement.h"
#include "TrackingTools/GsfTracking/interface/GsfPropagatorWithMaterial.h"
#include "TrackingTools/GsfTools/interface/GsfPropagatorAdapter.h"
#include "TrackingTools/GsfTracking/interface/FullConvolutionWithMaterial.h"
#include "TrackingTools/DetLayers/interface/DetLayerGeometry.h"

class MultiTrajectoryStateMerger;

class GsfTrajectorySmoother : public TrajectorySmoother {

private:

  typedef TrajectoryStateOnSurface TSOS;
  typedef TrajectoryMeasurement TM;

public:

  GsfTrajectorySmoother(const GsfPropagatorWithMaterial& aPropagator,
            const TrajectoryStateUpdator& aUpdator,
            const MeasurementEstimator& aEstimator,
            const MultiTrajectoryStateMerger& merger,
            float errorRescaling,
            const bool materialBeforeUpdate = true,
            const DetLayerGeometry* detLayerGeometry=0);

  virtual ~GsfTrajectorySmoother();

  virtual std::vector<Trajectory> trajectories(const Trajectory& aTraj) const;
  const Propagator* propagator() const {
    if ( thePropagator) return thePropagator;
    else  return theGeomPropagator;
  }
  const TrajectoryStateUpdator* updator() const {return theUpdator;}
  const MeasurementEstimator* estimator() const {return theEstimator;}

  virtual GsfTrajectorySmoother* clone() const
  {
    return new GsfTrajectorySmoother(*thePropagator,*theUpdator,*theEstimator,
                     *theMerger,theErrorRescaling,theMatBeforeUpdate,theGeometry);
  }

private:
  GsfPropagatorWithMaterial* thePropagator;
  const GsfPropagatorAdapter* theGeomPropagator;
  const FullConvolutionWithMaterial* theConvolutor;
  const TrajectoryStateUpdator* theUpdator;
  const MeasurementEstimator* theEstimator;
  const MultiTrajectoryStateMerger* theMerger;
  
  bool theTiming;
  bool theMatBeforeUpdate;
  float theErrorRescaling;

  const DetLayerGeometry dummyGeometry;
  const DetLayerGeometry* theGeometry;

};

#endif //TR_GsfTrajectorySmoother_H_