---

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 ( float mass ) :
    theMass(mass),
    theDeltaP(0.),
    theDeltaCov() {}
  //  MaterialEffectsUpdator ();

  virtual ~MaterialEffectsUpdator () {}

  virtual TrajectoryStateOnSurface updateState (const TrajectoryStateOnSurface& TSoS, 
                                                const PropagationDirection propDir) const;
  virtual double deltaP (const TrajectoryStateOnSurface& TSoS, const PropagationDirection propDir) const {
    // check for material
    if ( !TSoS.surface().mediumProperties() )  return 0.;
    // check for change (avoid using compute method if possible)
    if ( newArguments(TSoS,propDir) )  compute(TSoS,propDir);
    return theDeltaP;
  }


  virtual const AlgebraicSymMatrix55 &deltaLocalError (const TrajectoryStateOnSurface& TSoS, 
                                              const PropagationDirection propDir) const {
    // check for material
    if ( !TSoS.surface().mediumProperties() )  return theNullMatrix;
    // check for change (avoid using compute method if possible)
    if ( newArguments(TSoS,propDir) )  compute(TSoS,propDir);
    return theDeltaCov;
  }  
  inline float 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;

 protected:
  // check of arguments for use with cached values
  virtual bool newArguments (const TrajectoryStateOnSurface&, const PropagationDirection) const {
    return true;
  }
  // storage of arguments for later use
  virtual void storeArguments (const TrajectoryStateOnSurface&, const PropagationDirection) const {
  }
  
 private:
  float theMass;


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

#endif

---