lha_elements.h

Go to the documentation of this file.

#ifndef HYPERISO_LHA_ELEMENTS_H
#define HYPERISO_LHA_ELEMENTS_H
 
#include <string>
#include <optional>
#include <vector>
#include <memory>
#include <map>
#include <sstream>
#include <iostream>
#include <limits>
 
#include "Logger.h"
#include "LhaID.h"
#include "LhaBlockPrototype.h"
#include "DBNode.h"
 
enum class RenormalizationScheme {
    POLE, MSBAR, DRBAR, ONE_S, KIN, INVARIANT, MOM, SMOM, NONE
};
 
class AbstractElement {
protected:
    const LhaID id;  
public:
    inline explicit AbstractElement(const LhaID& id_) : id(id_) {}
 
    inline LhaID getId() const { return this->id; }
 
    virtual std::pair<double, double> getBinning() const = 0;
 
    virtual double getScale() const = 0;
 
    virtual std::string toString() const = 0;
 
    virtual std::shared_ptr<DBNode> toDBNode() const = 0;
 
    virtual ~AbstractElement() = default;
 
};
 
template <typename T>
class LhaElement : public AbstractElement {
private:
    T value;                                        
    std::optional<RenormalizationScheme> rScheme;   
    std::optional<double> Q;                        
    std::optional<std::pair<double, double>> bin;   
    LhaID encodeId(const Prototype& prototype, const std::vector<std::string>& line);
 
public:
    LhaElement(const Prototype& prototype, const std::vector<std::string>& line);
 
    inline T getValue() const { return this->value; }
 
    inline RenormalizationScheme getScheme() const { return this->rScheme.has_value() ? this->rScheme.value() : RenormalizationScheme::NONE; }
    
    inline double getScale() const override { return this->Q.has_value() ? this->Q.value() : 0; }
 
    inline std::pair<double, double> getBinning() const override { return this->bin.value_or(std::pair(-1.0, -1.0)); }
    
    std::string toString() const override;
 
    std::shared_ptr<DBNode> toDBNode() const override;
    
};
 
class LhaElementFactory {
    public:
        static std::shared_ptr<AbstractElement> createElement(const Prototype& prototype, const std::vector<std::string>& line);
    };
 
#endif // HYPERISO_LHA_ELEMENTS_H