RCModule.h

/*
 * RCModule.h
 *
 *  Created on: Feb 12, 2021
 *      Author: Pawel Sznajder (NCBJ)
 */
 
#ifndef MODULES_RADIATIVECORRECTIONS_RCMODULE_H_
#define MODULES_RADIATIVECORRECTIONS_RCMODULE_H_
 
#include <ElementaryUtils/logger/CustomException.h>
#include <ElementaryUtils/string_utils/Formatter.h>
#include <partons/beans/channel/ChannelType.h>
#include <stddef.h>
#include <string>
#include <tuple>
#include <vector>
 
#include "../../beans/containers/ExperimentalConditions.h"
#include "../../beans/containers/KinematicRange.h"
#include "../EpicModuleObject.h"
 
namespace EPIC {
class Event;
} /* namespace EPIC */
 
namespace EPIC {
class ExperimentalConditions;
} /* namespace EPIC */
 
namespace EPIC {
 
template<typename ObsKinType>
class RCModule: public EpicModuleObject {
 
public:
 
    virtual ~RCModule() {
    }
 
    size_t getNVariables() const;
 
    const std::vector<KinematicRange>& getVariableRanges() const {
        return m_variableRanges;
    }
 
    virtual std::tuple<double, ExperimentalConditions, ObsKinType> evaluate(
            const ExperimentalConditions& experimentalConditions,
            const ObsKinType& obsKin, const std::vector<double>& par) const = 0;
 
    virtual void updateEvent(Event& event, const std::vector<double>& par) const = 0;
 
protected:
 
    RCModule(const std::string &className,
            PARTONS::ChannelType::Type channelType) :
            EpicModuleObject(className, channelType), m_nVariables(0) {
    }
 
    RCModule(const RCModule &other) :
            EpicModuleObject(other), m_nVariables(other.m_nVariables), m_variableRanges(
                    other.m_variableRanges) {
    }
 
protected:
 
    virtual void checkParameters(const std::vector<double>& par) const {
 
        if (par.size() != m_nVariables) {
            throw ElemUtils::CustomException(getClassName(), __func__,
                    ElemUtils::Formatter() << "Wrong number of variables, is: "
                            << par.size() << ", expected: " << m_nVariables);
        }
    }
 
    void setVariableRanges(const std::vector<KinematicRange>& variableRanges) {
 
        m_nVariables = variableRanges.size();
        m_variableRanges = variableRanges;
    }
 
private:
 
    size_t m_nVariables; 
    std::vector<KinematicRange> m_variableRanges; 
};
}
 
#endif /* MODULES_RADIATIVECORRECTIONS_RCMODULE_H_ */