L2AbsScaleCalculator.h

Go to the documentation of this file.

#ifndef JetMETCorrections_FFTJetObjects_L2AbsScaleCalculator_h
#define JetMETCorrections_FFTJetObjects_L2AbsScaleCalculator_h

#include <cmath>
#include <cassert>

#include "JetMETCorrections/FFTJetObjects/interface/AbsFFTSpecificScaleCalculator.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/Exception.h"

class L2AbsScaleCalculator : public AbsFFTSpecificScaleCalculator {
public:
  inline explicit L2AbsScaleCalculator(const edm::ParameterSet& ps)
      : m_radiusFactor(ps.getParameter<double>("radiusFactor")),
        m_zeroPtLog(ps.getParameter<double>("zeroPtLog")),
        m_takePtLog(ps.getParameter<bool>("takePtLog")) {}

  inline ~L2AbsScaleCalculator() override {}

  inline void mapFFTJet(const reco::Jet& /* jet */,
                        const reco::FFTJet<float>& fftJet,
                        const math::XYZTLorentzVector& current,
                        double* buf,
                        const unsigned dim) const override {
    if (dim != 2)
      throw cms::Exception("FFTJetBadConfig") << "In L2AbsScaleCalculator::mapFFTJet: "
                                              << "invalid table dimensionality: " << dim << std::endl;
    assert(buf);
    const double radius = fftJet.f_recoScale();
    const double pt = current.pt();
    buf[0] = radius * m_radiusFactor;
    if (m_takePtLog) {
      if (pt > 0.0)
        buf[1] = log(pt);
      else
        buf[1] = m_zeroPtLog;
    } else
      buf[1] = pt;
  }

private:
  double m_radiusFactor;
  double m_zeroPtLog;
  bool m_takePtLog;
};

#endif  // JetMETCorrections_FFTJetObjects_L2AbsScaleCalculator_h