MaterialEffectsUpdator.h

Go to the documentation of this file.

#ifndef _CR_MATERIALEFFECTSUPDATOR_H_
#define _CR_MATERIALEFFECTSUPDATOR_H_

#include "DataFormats/GeometrySurface/interface/Surface.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
#include "DataFormats/TrajectorySeed/interface/PropagationDirection.h"
class MaterialEffectsUpdator
{  
public:
  MaterialEffectsUpdator ( double mass );
  virtual ~MaterialEffectsUpdator ();

  virtual TrajectoryStateOnSurface updateState (const TrajectoryStateOnSurface& TSoS, 
                        const PropagationDirection propDir) const;

  virtual bool updateStateInPlace (TrajectoryStateOnSurface& TSoS, 
                   const PropagationDirection propDir) const;

 
  virtual double deltaP (const TrajectoryStateOnSurface& TSoS, const PropagationDirection propDir) const;


  virtual const AlgebraicSymMatrix55 &deltaLocalError (const TrajectoryStateOnSurface& TSoS, 
                               const PropagationDirection propDir) const;

  inline double mass () const {
    return theMass;
  }

  virtual MaterialEffectsUpdator* clone()  const = 0;

 private:
  // here comes the actual computation of the values
  virtual void compute (const TrajectoryStateOnSurface&, const PropagationDirection) const = 0;

  // check of arguments for use with cached values
  bool newArguments (const TrajectoryStateOnSurface & TSoS, PropagationDirection  propDir) const;
  
 private:
  double theMass;

  // chache previous call state
  mutable double theLastOverP;
  mutable double theLastDxdz;
  mutable float  theLastRL;
  mutable PropagationDirection theLastPropDir;


protected:  
  mutable double theDeltaP;
  mutable AlgebraicSymMatrix55 theDeltaCov;
  static  AlgebraicSymMatrix55  theNullMatrix;
};

#endif