SMParamSetter.h

Go to the documentation of this file.

#ifndef SMPARAMSETTER_H
#define SMPARAMSETTER_H
 
#include <cmath>
#include <set>
#include <string>
#include <iostream>
#include <vector>
#include <optional>
 
#include "config.hpp"
#include "LhaID.h"
#include "Interpreter.h"
#include "IMartyParameterProxy.h"
#include "CinematicExtractor.h"
 
class SMParamSetter {
public:
    SMParamSetter(const std::string& model,
                  std::set<std::string> special_blocks,
                  std::shared_ptr<IMartyParameterProxy<std::string, LhaID>> sm_proxy,
                  std::shared_ptr<IMartyParameterProxy<std::string, LhaID>> bsm_proxy = nullptr,
                  std::string cinematic_template = "");
 
    std::unordered_map<std::string, double> setParam(const std::string& name, const InterpretedParam& interpretedParam);
 
private:
    scalar_t calculateValue(const InterpretedParam& interpretedParam);
 
    scalar_t calculateKinematicInvariant(const LhaID& code) const;
    scalar_t calculateOneToThreeInvariant(const LhaID& code, const std::vector<scalar_t>& masses) const;
    scalar_t calculateOneToTwoInvariant(const LhaID& code, const std::vector<scalar_t>& masses) const;
    scalar_t legacyKinematicInvariant(const LhaID& code) const;
 
    std::vector<scalar_t> extractMassesForCurrentProcess() const;
    scalar_t massValueForParticle(const std::string& particle_name) const;
 
    Model model_type;
    std::set<std::string> special_blocks;
 
    std::shared_ptr<IMartyParameterProxy<std::string, LhaID>> sm_proxy;
    std::shared_ptr<IMartyParameterProxy<std::string, LhaID>> bsm_proxy;
 
    std::string cinematic_template;
    std::optional<CinematicProcess> cinematic_process;
};
 
#endif // SMPARAMSETTER_H