FFTJetObjects/interface/L2ResScaleCalculator.h

0001 #ifndef JetMETCorrections_FFTJetObjects_L2ResScaleCalculator_h
0002 #define JetMETCorrections_FFTJetObjects_L2ResScaleCalculator_h
0003
0004 #include <cassert>
0005
0006 #include "JetMETCorrections/FFTJetObjects/interface/AbsFFTSpecificScaleCalculator.h"
0007 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0008 #include "FWCore/Utilities/interface/Exception.h"
0009
0010 class L2ResScaleCalculator : public AbsFFTSpecificScaleCalculator {
0011 public:
0012   inline explicit L2ResScaleCalculator(const edm::ParameterSet& ps)
0013       : m_radiusFactor(ps.getParameter<double>("radiusFactor")) {}
0014
0015   inline ~L2ResScaleCalculator() override {}
0016
0017   inline void mapFFTJet(const reco::Jet& /* jet */,
0018                         const reco::FFTJet<float>& fftJet,
0019                         const math::XYZTLorentzVector& current,
0020                         double* buf,
0021                         const unsigned dim) const override {
0022     if (dim != 2)
0023       throw cms::Exception("FFTJetBadConfig") << "In L2ResScaleCalculator::mapFFTJet: "
0024                                               << "invalid table dimensionality: " << dim << std::endl;
0025     assert(buf);
0026     const double radius = fftJet.f_recoScale();
0027     buf[0] = radius * m_radiusFactor;
0028     buf[1] = current.eta();
0029   }
0030
0031 private:
0032   double m_radiusFactor;
0033 };
0034
0035 #endif  // JetMETCorrections_FFTJetObjects_L2ResScaleCalculator_h