/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_7/src/JetMETCorrections/FFTJetObjects/interface/FFTJetScaleCalculators.h

Go to the documentation of this file.

#ifndef JetMETCorrections_FFTJetObjects_FFTJetScaleCalculators_h
#define JetMETCorrections_FFTJetObjects_FFTJetScaleCalculators_h

#include <cassert>
#include <cmath>

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

template<class MyJet, class Adjustable>
class FFTEtaLogPtConeRadiusMapper : 
    public AbsFFTJetScaleCalculator<MyJet, Adjustable>
{
public:
    inline explicit FFTEtaLogPtConeRadiusMapper(
        CPP11_shared_ptr<npstat::AbsMultivariateFunctor> f)
        : AbsFFTJetScaleCalculator<MyJet, Adjustable>(f) {}

private:
    inline void map(const MyJet& jet,
                    const Adjustable& current,
                    double* buf, const unsigned dim) const
    {
        assert(buf);
        if (dim != 3)
            throw cms::Exception("FFTJetBadConfig")
                << "In FFTEtaLogPtConeRadiusMapper::map: "
                << "invalid table dimensionality: "
                << dim << std::endl;
        buf[0] = current.vec().eta();
        buf[1] = log(current.vec().pt());
        buf[2] = jet.getFFTSpecific().f_recoScale();
    }
};

template<class MyJet, class Adjustable>
class FFTSpecificScaleCalculator : 
    public AbsFFTJetScaleCalculator<MyJet, Adjustable>
{
public:
    //
    // This class will assume the ownership of the
    // AbsFFTSpecificScaleCalculator object provided
    // in the constructor
    //
    inline FFTSpecificScaleCalculator(
        CPP11_shared_ptr<npstat::AbsMultivariateFunctor> f,
        const AbsFFTSpecificScaleCalculator* p)
        : AbsFFTJetScaleCalculator<MyJet, Adjustable>(f), calc_(p) {assert(p);}

    inline virtual ~FFTSpecificScaleCalculator() {delete calc_;}

private:
    inline void map(const MyJet& jet,
                    const Adjustable& current,
                    double* buf, const unsigned dim) const
    {
        return calc_->mapFFTJet(jet, jet.getFFTSpecific(),
                                current.vec(), buf, dim);
    }

    const AbsFFTSpecificScaleCalculator* calc_;
};

#endif // JetMETCorrections_FFTJetObjects_FFTJetScaleCalculators_h