FFTJetAdjusters.h

Go to the documentation of this file.

#include <cassert>
 
#include "JetMETCorrections/FFTJetObjects/interface/AbsFFTJetAdjuster.h"
#include "FWCore/Utilities/interface/Exception.h"
 
template <class MyJet, class Adjustable>
struct FFTSimpleScalingAdjuster : public AbsFFTJetAdjuster<MyJet, Adjustable> {
  inline ~FFTSimpleScalingAdjuster() override {}
 
  void adjust(const MyJet& /* jet */,
              const Adjustable& in,
              const double* factors,
              const unsigned lenFactors,
              Adjustable* out) const override {
    if (lenFactors != 1U)
      throw cms::Exception("FFTJetBadConfig") << "In FFTSimpleScalingAdjuster::adjust: wrong number of "
                                              << "scales (expected 1, got " << lenFactors << ")\n";
    assert(factors);
    assert(out);
    *out = in;
    *out *= factors[0];
  }
};
 
template <class MyJet, class Adjustable>
struct FFTUncertaintyAdjuster : public AbsFFTJetAdjuster<MyJet, Adjustable> {
  inline ~FFTUncertaintyAdjuster() override {}
 
  void adjust(const MyJet& /* jet */,
              const Adjustable& in,
              const double* factors,
              const unsigned lenFactors,
              Adjustable* out) const override {
    if (lenFactors != 1U)
      throw cms::Exception("FFTJetBadConfig") << "In FFTUncertaintyAdjuster::adjust: wrong number of "
                                              << "scales (expected 1, got " << lenFactors << ")\n";
    assert(factors);
    assert(out);
    *out = in;
    const double s = factors[0];
    out->setVariance(in.variance() + s * s);
  }
};
 
template <class MyJet, class Adjustable>
struct FFTScalingAdjusterWithUncertainty : public AbsFFTJetAdjuster<MyJet, Adjustable> {
  inline ~FFTScalingAdjusterWithUncertainty() override {}
 
  void adjust(const MyJet& /* jet */,
              const Adjustable& in,
              const double* factors,
              const unsigned lenFactors,
              Adjustable* out) const override {
    if (lenFactors != 2U)
      throw cms::Exception("FFTJetBadConfig") << "In FFTScalingAdjusterWithUncertainty::adjust: wrong "
                                              << "number of scales (expected 2, got " << lenFactors << ")\n";
    assert(factors);
    assert(out);
    *out = in;
    *out *= factors[0];
    const double s = factors[1];
    out->setVariance(in.variance() + s * s);
  }
};