SimpleClusterSmearer.h

#ifndef SimpleClusterSmearer_h
#define SimpleClusterSmearer_h 1

#include "marlin/MarlinConfig.h"

#ifdef MARLIN_CLHEP  // only if CLHEP is available !

#include "marlin/IFourVectorSmearer.h"
#include <vector>

namespace marlin{

  struct ClusterResolution {
    
    ClusterResolution() :
      A(0.) ,
      B(0.) ,
      ThMin(0.) ,
      ThMax(0.) {}
    
    ClusterResolution(float a, float b, float thMin, float thMax) :
      A(a) ,
      B(b) ,
      ThMin(thMin) ,
      ThMax(thMax) {}
    
    float A ;
    float B ;
    float ThMin ;
    float ThMax ;
    
  } ;

  class SimpleClusterSmearer : public IFourVectorSmearer {
    
  
    typedef std::vector<ClusterResolution> ResVec ;

  public:
    
    SimpleClusterSmearer(const std::vector<float>& resVec ) ;
    
    
    virtual ~SimpleClusterSmearer() {} 
    
    virtual HepLorentzVector smearedFourVector( const HepLorentzVector& v, int pdgCode ) ;
    
  protected:

    ResVec _resVec ;

  } ;
  


} // end namespace 

#endif // MARLIN_CLHEP  
#endif // SimpleClusterSmearer_h