WilsonMatchingPatch.h

Go to the documentation of this file.

#ifndef HYPERISO_WILSON_MATCHING_PATCH_H
#define HYPERISO_WILSON_MATCHING_PATCH_H
 
#include <functional>
#include <string>
#include <unordered_set>
 
#include "Include.h"
#include "Math.h"
#include "SourcesView.h"
 
using WilsonMatchingPatchFunction = std::function<scalar_t(const ParamSrc&)>;
 
struct WilsonMatchingPatch {
    WGroupId group{};
    WCoefId coefficient{};
    QCDOrder order{QCDOrder::LO};
    ContributionType contribution{ContributionType::SM};
    std::unordered_set<ParamId> sources{};
    WilsonMatchingPatchFunction compute{};
    std::string label{};
    bool enabled{true};
 
    bool marty_only{false};
};
 
inline WilsonMatchingPatch make_hyperiso_c9_sm_photon_patch() {
    WilsonMatchingPatch patch;
    patch.group = GroupMapper::to_id(WGroup::B);
    patch.coefficient = WCoefMapper::to_id(WCoef::C9);
    patch.order = QCDOrder::LO;
    patch.contribution = ContributionType::SM;
    patch.sources = {
        ParamId{ParameterType::WILSON, "WPARAM_MATCH_SM", LhaID(2, 1)},
        ParamId{ParameterType::WILSON, "WPARAM_MATCH_SM", LhaID(3)}
    };
    patch.compute = [](const ParamSrc& src) -> scalar_t {
        const double xt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", LhaID(2, 1));
        const scalar_t L = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", LhaID(3));
        return scalar_t(-D0t(xt) + 38. / 27.) - scalar_t(4. / 9.) * L;
    };
    patch.label = "Hyperiso:C9:SM-photon";
    patch.marty_only = true;
    return patch;
}
 
// Backward-compatible alias for older user code.  The old name only described
// the light u/c part; the helper now returns the full SM photon matching piece.
inline WilsonMatchingPatch make_hyperiso_c9_light_photon_sm_patch() {
    return make_hyperiso_c9_sm_photon_patch();
}
 
inline WilsonMatchingPatch make_hyperiso_c9_thdm_photon_patch() {
    WilsonMatchingPatch patch;
    patch.group = GroupMapper::to_id(WGroup::B);
    patch.coefficient = WCoefMapper::to_id(WCoef::C9);
    patch.order = QCDOrder::LO;
    patch.contribution = ContributionType::BSM;
    patch.sources = {
        ParamId{ParameterType::WILSON, "WPARAM_SI_BSM", LhaID(7)},
        ParamId{ParameterType::WILSON, "WPARAM_MATCH_BSM", LhaID(1)}
    };
    patch.compute = [](const ParamSrc& src) -> scalar_t {
        const double lu = src.get_val(ParameterType::WILSON, "WPARAM_SI_BSM", LhaID(7));
        const double yt = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_BSM", LhaID(1));
        return scalar_t(-D9H0(yt, lu));
    };
    patch.label = "Hyperiso:C9:THDM-photon";
    patch.marty_only = true;
    return patch;
}
 
inline WilsonMatchingPatch make_hyperiso_cp9_thdm_photon_patch() {
    WilsonMatchingPatch patch;
    patch.group = GroupMapper::to_id(WGroup::BPrime);
    patch.coefficient = WCoefMapper::to_id(WCoef::CP9);
    patch.order = QCDOrder::LO;
    patch.contribution = ContributionType::BSM;
    patch.sources = {
        ParamId{ParameterType::WILSON, "WPARAM_SI_BSM", LhaID(8)},
        ParamId{ParameterType::WILSON, "WPARAM_MATCH_BSM", LhaID(1)},
        ParamId{ParameterType::WILSON, "WPARAM_MATCH_SM", LhaID(6)},
        ParamId{ParameterType::WILSON, "WPARAM_MATCH_SM", LhaID(5, 1)},
        ParamId{ParameterType::SM, "MASS", LhaID(3)}
    };
    patch.compute = [](const ParamSrc& src) -> scalar_t {
        const double ld     = src.get_val(ParameterType::WILSON, "WPARAM_SI_BSM", LhaID(8));
        const double yt     = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_BSM", LhaID(1));
        const double mt_muW = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", LhaID(6));
        const double mb_muW = src.get_val(ParameterType::WILSON, "WPARAM_MATCH_SM", LhaID(5, 1));
        const double ms     = src.get_val(ParameterType::SM, "MASS", LhaID(3));
        return scalar_t(-ms * mb_muW / (mt_muW * mt_muW) * D9H0(yt, ld));
    };
    patch.label = "Hyperiso:CP9:THDM-photon";
    patch.marty_only = true;
    return patch;
}
 
#endif // HYPERISO_WILSON_MATCHING_PATCH_H