CmsLikeJetSmear.h

#ifndef FastJetStability_CmsLikeJetSmear_h
#define FastJetStability_CmsLikeJetSmear_h

#include "AbsSmearJet.hh"
#include "CmsLikeJetEfficiency.h"

namespace stab {

class CmsLikeJetSmear: public AbsSmearJet {

public:
  CmsLikeJetSmear(bool makePtPositive=false,bool runEfficiency=true);
  inline virtual ~CmsLikeJetSmear() {} 

  inline virtual void initialize(const Event &ev,const StabilityEventInfo &info) 
    {
      this->getEvent().clear();
      this->getInfo().clear();
      if ( m_runEfficiency ) {
        m_eff.initialize(ev,info);
        m_eff.setFriendlyInfo( this->getInfo() );
        m_eff.smear(this->getEvent());
      } else {
    this->getEvent() = ev;  
    this->getInfo() = info;
      }
    }


  // be friends with QCDandCmsLikeSmear
  friend class QCDandCmsLikeSmear;

private:

  inline double ecalEtUncertainty(const double Et) const
        {
            return calorimeterEtUncertainty(m_ecalStochasticTerm, m_ecalNoiseTerm,
                                            m_ecalConstantTerm, Et);
        }
     
  inline double hcalEtUncertainty(const double Et) const
        {
            return calorimeterEtUncertainty(m_hcalStochasticTerm, m_hcalNoiseTerm,
                                            m_hcalConstantTerm, Et);
        }

  // Pure virtual methods from the base class
  virtual void processPhoton(const PseudoJet& in, const StabilityUserInfo &info,
                             ParticleCollection* out
                 ) const;
  virtual void processTrack(const PseudoJet& in, const StabilityUserInfo &info,
                            ParticleCollection* out
                ) const;
  virtual void processNeutralHadron(const PseudoJet& in, const StabilityUserInfo &info,
                                    ParticleCollection* out
                    ) const;


  // ----------------------- member data --------------------------
  // ECAL resolution parameters (photons)
  double m_ecalStochasticTerm;
  double m_ecalNoiseTerm;
  double m_ecalConstantTerm;

  // HCAL resolution parameters (neutral hadrons, missed tracks)
  double m_hcalStochasticTerm;
  double m_hcalNoiseTerm;
  double m_hcalConstantTerm;

  // Track Pt resolution
  double m_trackPtResolution;

  // Angular resolutions
  double m_angularResolutionPhoton;
  double m_angularResolutionNeutralHadron;
  double m_angularResolutionTrack;

  bool m_onlyPositivePt;
  bool m_runEfficiency;

  CmsLikeJetEfficiency m_eff;


};  // end CmsLikeJetSmear class declaration

} // end stab namespace

#endif