File indexing completed on 2023-03-17 11:10:35
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& ,
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