QCDandCmsLikeSmear.h

#ifndef FastJetStability_QCDandCmsLikeSmear_h
#define FastJetStability_QCDandCmsLikeSmear_h

#include "JetStabilityDefs.h"
#include "StabilityUserInfo.h"
#include "AbsSmearingAgent.h"
#include "CmsLikeJetSmear.h"
#include "QCDSmearSequence.h"
#include "QCDSmear.h"

namespace stab {

class QCDandCmsLikeSmear: public AbsSmearingAgent {

public:
  inline QCDandCmsLikeSmear()
    :AbsSmearingAgent() 
    ,m_CmsSmear1(false,true)
    ,m_CmsSmear2(true,false)
    ,m_QCDSmear(m_alphaS)
    { }

  virtual inline ~QCDandCmsLikeSmear() { }

  inline virtual void initialize(const Event &ev, const StabilityEventInfo &info) 
    { 
      m_event.clear();
      m_CmsSmear1.initialize(ev,info);
      m_CmsSmear1.smear(m_event);
      assert(m_event.size() <= ev.size() );
      m_CmsSmear2.initialize(m_event,info);
      //m_QCDSmear.initialize(m_event);
    }

  virtual void smear(Event &);

private:

  Event m_event;

  AlphaStrong m_alphaS;
  CmsLikeJetSmear m_CmsSmear1;
  CmsLikeJetSmear m_CmsSmear2;
  QCDSmearSequence m_QCDSmear; 

};

} // end stab namespace

#endif