CmsLikeJetSmear.cc

#include <cmath>
#include <cassert>

#include "generator/pythia6.h"

#include "CmsLikeJetSmear.h"

namespace stab {

  CmsLikeJetSmear::CmsLikeJetSmear(const bool makePtPositive,const bool runEfficiency) 
        : m_ecalStochasticTerm(0.028),
          m_ecalNoiseTerm(0.125),
          m_ecalConstantTerm(0.003),
          m_hcalStochasticTerm(0.85),
          m_hcalNoiseTerm(0.0),
          m_hcalConstantTerm(0.075),
          m_trackPtResolution(0.02),
          m_angularResolutionPhoton(0.005),
          m_angularResolutionNeutralHadron(0.025),
          m_angularResolutionTrack(0.0015),
          m_onlyPositivePt(makePtPositive),
      m_runEfficiency(runEfficiency),
      m_eff(0.25)
    {
    }

    void CmsLikeJetSmear::processPhoton(const PseudoJet& in, const StabilityUserInfo &info,
                                        ParticleCollection* out
                    ) const
    {
        const double sigmaPt = ecalEtUncertainty(in.pt());
        PseudoJet newP4 = smear4Mom(in, sigmaPt, m_angularResolutionPhoton,
                                        m_onlyPositivePt);

    //StabilityUserInfo *info = new StabilityUserInfo( in.user_info<StabilityUserInfo>() );
    //newP4.set_user_info(info);
        out->push_back( newP4 );
    }

    void CmsLikeJetSmear::processTrack(const PseudoJet& in, const StabilityUserInfo &info,
                                       ParticleCollection* out
                    ) const
    {
        const double sigmaPt = m_trackPtResolution*in.pt();
        PseudoJet newP4 = smear4Mom(in, sigmaPt, m_angularResolutionTrack,
                                        m_onlyPositivePt);

    //StabilityUserInfo *info = new StabilityUserInfo( in.user_info<StabilityUserInfo>() );
    //newP4.set_user_info(info);
        out->push_back( newP4 );
    }     
          
    void CmsLikeJetSmear::processNeutralHadron(const PseudoJet& in, const StabilityUserInfo &info,
                                               ParticleCollection* out
                           ) const
    {     
        const double sigmaPt = hcalEtUncertainty(in.pt());
        PseudoJet newP4 = smear4Mom(in, sigmaPt,
                                        m_angularResolutionNeutralHadron,
                                        m_onlyPositivePt);

    //StabilityUserInfo *info = new StabilityUserInfo( in.user_info<StabilityUserInfo>() );
    //newP4.set_user_info( info );
        out->push_back( newP4 );
    }


} // end stab namespace